A new statistical method for estimating divergence dates of species from DNA sequence data by a molecular clock approach is developed. This method takes into account effectively the information contained in a set of DNA sequence data. The molecular clock of mitochondrial DNA (mtDNA) was calibrated by setting the date of divergence between primates and ungulates at the Cretaceous-Tertiary boundary (65 million years ago), when the extinction of dinosaurs occurred. A generalized least-squares method was applied in fitting a model to mtDNA sequence data, and the clock gave dates of 92.3 +/- 11.7, 13.3 +/- 1.5, 10.9 +/- 1.2, 3.7 +/- 0.6, and 2.7 +/- 0.6 million years ago (where the second of each pair of numbers is the standard deviation) for the separation of mouse, gibbon, orangutan, gorilla, and chimpanzee, respectively, from the line leading to humans. Although there is some uncertainty in the clock, this dating may pose a problem for the widely believed hypothesis that the pipedal creature Australopithecus afarensis, which lived some 3.7 million years ago at Laetoli in Tanzania and at Hadar in Ethiopia, was ancestral to man and evolved after the human-ape splitting. Another likelier possibility is that mtDNA was transferred through hybridization between a proto-human and a proto-chimpanzee after the former had developed bipedalism.

Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.

In a previous paper, we introduced MUSCLE, a new program for creating multiple alignments of protein sequences, giving a brief summary of the algorithm and showing MUSCLE to achieve the highest scores reported to date on four alignment accuracy benchmarks. Here we present a more complete discussion of the algorithm, describing several previously unpublished techniques that improve biological accuracy and / or computational complexity. We introduce a new option, MUSCLE-fast, designed for high-throughput applications. We also describe a new protocol for evaluating objective functions that align two profiles. We compare the speed and accuracy of MUSCLE with CLUSTALW, Progressive POA and the MAFFT script FFTNS1, the fastest previously published program known to the author. Accuracy is measured using four benchmarks: BAliBASE, PREFAB, SABmark and SMART. We test three variants that offer highest accuracy (MUSCLE with default settings), highest speed (MUSCLE-fast), and a carefully chosen compromise between the two (MUSCLE-prog). We find MUSCLE-fast to be the fastest algorithm on all test sets, achieving average alignment accuracy similar to CLUSTALW in times that are typically two to three orders of magnitude less. MUSCLE-fast is able to align 1,000 sequences of average length 282 in 21 seconds on a current desktop computer. MUSCLE offers a range of options that provide improved speed and / or alignment accuracy compared with currently available programs. MUSCLE is freely available at http://www.drive5.com/muscle
 .

/pdf/muscle-a-multiple-sequence-alignment-method-with-reduced-9d0g1qzdlk.pdf

MUSCLE: a multiple sequence alignment method with reduced time and space complexity

Our capacity to sequence human genomes has exceeded our ability to interpret genetic variation. Current genomic annotations tend to exploit a single information type (e.g. conservation) and/or are restricted in scope (e.g. to missense changes). Here, we describe Combined Annotation Dependent Depletion (CADD), a framework that objectively integrates many diverse annotations into a single, quantitative score. We implement CADD as a support vector machine trained to differentiate 14.7 million high-frequency human derived alleles from 14.7 million simulated variants. We pre-compute “C-scores” for all 8.6 billion possible human single nucleotide variants and enable scoring of short insertions/deletions. C-scores correlate with allelic diversity, annotations of functionality, pathogenicity, disease severity, experimentally measured regulatory effects, and complex trait associations, and highly rank known pathogenic variants within individual genomes. The ability of CADD to prioritize functional, deleterious, and pathogenic variants across many functional categories, effect sizes and genetic architectures is unmatched by any current annotation.

/pdf/a-general-framework-for-estimating-the-relative-zgrfr3u3u8.pdf

A general framework for estimating the relative pathogenicity of human genetic variants

Preservation of duplicate genes by complementary, degenerative mutations.

Summary
1The role of phenotypic plasticity in evolution has historically been a contentious issue because of debate over whether plasticity shields genotypes from selection or generates novel opportunities for selection to act. Because plasticity encompasses diverse adaptive and non-adaptive responses to environmental variation, no single conceptual framework adequately predicts the diverse roles of plasticity in evolutionary change.
2Different types of phenotypic plasticity can uniquely contribute to adaptive evolution when populations are faced with new or altered environments. Adaptive plasticity should promote establishment and persistence in a new environment, but depending on how close the plastic response is to the new favoured phenotypic optimum dictates whether directional selection will cause adaptive divergence between populations. Further, non-adaptive plasticity in response to stressful environments can result in a mean phenotypic response being further away from the favoured optimum or alternatively increase the variance around the mean due to the expression of cryptic genetic variation. The expression of cryptic genetic variation can facilitate adaptive evolution if by chance it results in a fitter phenotype.
3We conclude that adaptive plasticity that places populations close enough to a new phenotypic optimum for directional selection to act is the only plasticity that predictably enhances fitness and is most likely to facilitate adaptive evolution on ecological time-scales in new environments. However, this type of plasticity is likely to be the product of past selection on variation that may have been initially non-adaptive.
4We end with suggestions on how future empirical studies can be designed to better test the importance of different kinds of plasticity to adaptive evolution.

/pdf/adaptive-versus-non-adaptive-phenotypic-plasticity-and-the-55680hih4y.pdf

Adaptive versus non‐adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments

IN 1966, I became interested in the amino acid sequences of cytochrome c molecules ([Jukes 1966][1]). I noted that these sequences differed in the cytochromes c of various species to an extent that seemed unnecessary from the standpoint of their function. I stated, “The changes produced in

https://www.genetics.org/content/genetics/154/3/956.full.pdf

The neutral theory of molecular evolution.

Estimation of evolutionary changes in certain homologous polypeptide chains.

Deleterious mutations and neutral substitutions, discussing molecular evolution model for DNA and proteins

Deleterious Mutations and Neutral Substitutions

We have investigated the relationship between the G + C content of silent (synonymous) sites in codons and the amino acid composition of encoded proteins for approximately 1,600 human genes. There are positive correlations between silent site G + C and the proportions of codons for Arg, Pro, Ala, Trp, His, Gln, and Leu and negative ones for Tyr, Phe, Asn, Ile, Lys, Asp, Thr, and Glu. The median proteins coded by groups of genes that differ in silent-site G + C content also differ in amino acid composition, as do some proteins coded by homologous genes. The pattern of compositional change can be largely explained by directional mutation pressure, the genetic code, and differences in the frequencies of accepted amino acid substitutions; the shifts in protein composition are likely to be selectively neutral.

Relationship between G + C in silent sites of codons and amino acid composition of human proteins.

Schultz and Yarus (J. Mol. Biol. 235:1377–1380, 1994) have proposed that reassignment of codons in the genetic code passes through a stage in which the codons are ambiguously translated. In contrast, we state that such ambiguity would be deleterious, and that, to be reassigned, a codon, together with the tRNA that translates the codon, must first disappear from coding sequences, after which a tRNA appears with a mutated anticodon, and this enables the codon to reappear with a changed meaning. In the case of a stop codon, the relevant release factor must change so as not to recognize it.

Thomas H. Jukes

Papers

The neutral theory of molecular evolution.

Estimation of evolutionary changes in certain homologous polypeptide chains.

Deleterious Mutations and Neutral Substitutions

Relationship between G + C in silent sites of codons and amino acid composition of human proteins.

On codon reassignment.