Divergence of conserved non-coding sequences: rate estimates and relative rate tests

In many eukaryotic genomes only a small fraction of the DNA codes for proteins, but the non-protein coding DNA harbors important genetic elements directing the development and the physiology of the organisms, like promoters, enhancers, insulators, and micro-RNA genes. The molecular evolution of these genetic elements is difficult to study because their functional significance is hard to deduce from sequence information alone. Here we propose an approach to the study of the rate of evolution of functional non-coding sequences at a macro-evolutionary scale. We identify functionally important non-coding sequences as Conserved Non-Coding Nucleotide (CNCN) sequences from the comparison of two outgroup species. The CNCN sequences so identified are then compared to their homologous sequences in a pair of ingroup species, and we monitor the degree of modification these sequences suffered in the two ingroup lineages. We propose a method to test for rate differences in the modification of CNCN sequences among the two ingroup lineages, as well as a method to estimate their rate of modification. We apply this method to the full sequences of the HoxA clusters from six gnathostome species: a shark, Heterodontus francisci; a basal ray finned fish, Polypterus senegalus; the amphibian, Xenopus tropicalis; as well as three mammalian species, human, rat and mouse. The results show that the evolutionary rate of CNCN sequences is not distinguishable among the three mammalian lineages, while the Xenopus lineage has a significantly increased rate of evolution. Furthermore the estimates of the rate parameters suggest that in the stem lineage of mammals the rate of CNCN sequence evolution was more than twice the rate observed within the placental amniotes clade, suggesting a high rate of evolution of cis-regulatory elements during the origin of amniotes and mammals. We conclude that the proposed methods can be used for testing hypotheses about the rate and pattern of evolution of putative cis-regulatory elements.     

Divergence at neutral and non-neutral loci in Drosophila buzzatii populations and their hybrids

The impact of intraspecific hybridisation on fitness and morphological traits depends on the history of natural selection and genetic drift, which may have led to differently coadapted gene-complexes in the parental populations. The divergence at neutral and non-neutral loci between populations can be evaluated by estimating F_ST and Q_ST respectively, and hence give an estimate of drift and selection in the populations. Here we investigate (1) whether divergence between populations in quantitative traits (wing size and shape) can be attributed to selection or drift alone, (2) The impact of intraspecific hybridisation on estimators for divergence at neutral (F_ST) and non-neutral loci (Q_ST) in hybrids, (3) If measurement of shape is more informative than size in order to detect divergence in quantitative traits between populations. The aims were addressed by performing two hybridisations between three populations of Drosophila buzzatii, one between populations from Argentina and the Canary Islands (separated for 200 years), and the other between populations from Argentina and Australia (separated for 80 years). We observed the highest divergence at neutral loci between the Argentinean and Canary Island populations, but highest morphological divergence between the Argentinean and Australian populations, indicating that natural selection is acting on the wings. Divergence based on Q_ST measures in the hybrids was sensitive towards increased phenotypic variance (σ²p) within groups and should be used with care when σ²p of populations differ. Our results indicate that measures of shape give a better estimate of divergence at the underlying quantitative traits loci than measures of size.     

Estimation of the neutral mutation rate in a finite population from DNA sequence data

The sampling theory for the infinite site model taking into account the phylogenetic relationship between the alleles is developed for those cases in which two or three alleles are observed in the sample. From this theory a maximum likelihood estimate of θ = 4Nμ can be obtained. Unlike the maximum likelihood estimate of θ based on the infinite allele model or the number of segregating sites, this estimate of θ is a function of the frequencies of the alleles. This method is used to estimate θ for mitochondrial DNA in Drosophila melanogaster and D. virilis from data obtained by Shah and Langley (1979. Nature (London)281, 696–699) using restriction endonucleases.     

The spontaneous mutation rate in Drosophila simulans

The purpose of this investigation was to determine whether there exists a correlation in Drosophila between the spontaneous mutation rate and the amount of dispersed middle repetitive (mobile) DNA sequences. The amount of these sequences is 7 times less in Drosophila simulans as compared to Drosophila melanogaster. Therefore, if a correlation exists, the spontaneous mutation rate in Drosophila simulans should be 7 times lower than that in Drosophila melanogaster. We isolated an X-chromosome inversion after X-irradiation of wild-type Drosophila simulans males, that reduced crossing-over between white and forked, two X-linked visible markers, to less than 1%. This inversion was subsequently used to determine the sex-linked recessive lethal mutation rate in Drosophila simulans males of a laboratory strain marked with white. The frequency of these lethal mutations found is not different from that observed in Drosophila melanogaster males of laboratory strains.     

Sojourn times and substitution rate at overdominant and linked neutral loci

The sojourn times until fixation of an overdominant allele were investigated based on the diffusion equation. Furthermore, the rate of accumulation of mutations, or the substitution rate, was predicted from the mean extinction time of a common overdominant allele. The substitution rate calculated theoretically agreed well with that determined by computer simulation. Overdominant selection enhances the polymorphism at linked loci, while its effect on the sojourn times and the substitution rate at linked loci has not been studied yet. To solve these problems, a model that assumed two linked loci, each with infinite alleles, was examined by computer simulation. A decrease in the recombination rate between two loci markedly changed the distribution of sojourn times of a neutral allele. Although overdominant selection obviously increased the sojourn times and the polymorphism at a linked locus, the rate of nucleotide substitution at the neutral locus was not influenced significantly even if complete linkage was assumed. These results suggest that, in regions containing overdominant genes, linked neutral loci will exhibit elevated levels of polymorphism, but their rate of molecular evolution remains that predicted by neutral theory.     

The rate and spectrum of microsatellite mutation in Caenorhabditis elegans and Daphnia pulex

The effective use of microsatellite loci as tools for microevolutionary analysis requires knowledge of the factors influencing the rate and pattern of mutation, much of which is derived from indirect inference from population samples. Interspecific variation in microsatellite stability also provides a glimpse into aspects of phylogenetic constancy of mutational processes. Using long-term series of mutation-accumulation lines, we have obtained direct estimates of the spectrum of microsatellite mutations in two model systems: the nematode Caenorhabditis elegans and the microcrustacean Daphnia pulex. Although the scaling of the mutation rate with the number of tandem repeats is highly consistent across distantly related species, including yeast and human, the per-cell-division mutation rate appears to be elevated in multicellular species. Contrary to the expectations under the stepwise mutation model, most microsatellite mutations in C. elegans and D. pulex involve changes of multiple repeat units, with expansions being much more common than contractions.     

Characterizing the time dependency of human mitochondrial DNA mutation rate estimates

Previous research has established a discrepancy of nearly anorder of magnitude between pedigree-based and phylogeny-based(human vs. chimpanzee) estimates of the mitochondrial DNA (mtDNA)control region mutation rate. We characterize the time dependencyof the human mitochondrial hypervariable region one mutationrate by generating 14 new phylogeny-based mutation rate estimatesusing within-human comparisons and archaeological dates. Rateestimates based on population events between 15,000 and 50,000years ago are at least 2-fold lower than pedigree-based estimates.These within-human estimates are also higher than estimatesgenerated from phylogeny-based human–chimpanzee comparisons.Our new estimates establish a rapid decay in evolutionary mutationrate between approximately 2,500 and 50,000 years ago and aslow decay from 50,000 to 6 Ma. We then extend this analysisto the mtDNA-coding region. Our within-human coding region mutationrate estimates display a similar, though less rapid, time-dependentdecay. We explore the possibility that multiple hits explainthe discrepancy between pedigree-based and phylogeny-based mutationrates. We conclude that whereas nucleotide substitution modelsincorporating multiple hits do provide a possible explanationfor the discrepancy between pedigree-based and human–chimpanzeemutation rate estimates, they do not explain the rapid declineof within-human rate estimates. We propose that demographicprocesses such as serial bottlenecks prior to the Holocene couldexplain the difference between rates estimated before and after15,000 years ago. Our findings suggest that human mtDNA estimatesof dates of population and phylogenetic events should be adjustedin light of this time dependency of the mutation rate estimates.     

Direct estimation of the mitochondrial DNA mutation rate in Drosophila melanogaster

Mitochondrial DNA (mtDNA) variants are widely used in evolutionary genetics as markers for population history and to estimate divergence times among taxa. Inferences of species history are generally based on phylogenetic comparisons, which assume that molecular evolution is clock-like. Between-species comparisons have also been used to estimate the mutation rate, using sites that are thought to evolve neutrally. We directly estimated the mtDNA mutation rate by scanning the mitochondrial genome of Drosophila melanogaster lines that had undergone approximately 200 generations of spontaneous mutation accumulation (MA). We detected a total of 28 point mutations and eight insertion-deletion (indel) mutations, yielding an estimate for the single-nucleotide mutation rate of 6.2 × 10⁻⁸ per site per fly generation. Most mutations were heteroplasmic within a line, and their frequency distribution suggests that the effective number of mitochondrial genomes transmitted per female per generation is about 30. We observed repeated occurrences of some indel mutations, suggesting that indel mutational hotspots are common. Among the point mutations, there is a large excess of G→A mutations on the major strand (the sense strand for the majority of mitochondrial genes). These mutations tend to occur at nonsynonymous sites of protein-coding genes, and they are expected to be deleterious, so do not become fixed between species. The overall mtDNA mutation rate per base pair per fly generation in Drosophila is estimated to be about 10× higher than the nuclear mutation rate, but the mitochondrial major strand G→A mutation rate is about 70× higher than the nuclear rate. Silent sites are substantially more strongly biased towards A and T than nonsynonymous sites, consistent with the extreme mutation bias towards A+T. Strand-asymmetric mutation bias, coupled with selection to maintain specific nonsynonymous bases, therefore provides an explanation for the extreme base composition of the mitochondrial genome of Drosophila.     

The rate of spontaneous mutation for life-history traits in Caenorhabditis elegans.

Spontaneous mutations were accumulated in 100 replicate lines of Caenorhabditis elegans over a period of approximately 50 generations. Periodic assays of these lines and comparison to a frozen control suggest that the deleterious mutation rate for typical life-history characters in this species is at least 0.05 per diploid genome per generation, with the average mutational effect on the order of 14% or less in the homozygous state and the average mutational heritability approximately 0.0034. While the average mutation rate per character and the average mutational heritability for this species are somewhat lower than previous estimates for Drosophila, these differences can be reconciled to a large extent when the biological differences between these species are taken into consideration.     

Time dependency of molecular rate estimates and systematic overestimation of recent divergence times

Studies of molecular evolutionary rates have yielded a wide range of rate estimates for various genes and taxa. Recent studies based on population-level and pedigree data have produced remarkably high estimates of mutation rate, which strongly contrast with substitution rates inferred in phylogenetic (species-level) studies. Using Bayesian analysis with a relaxed-clock model, we estimated rates for three groups of mitochondrial data: avian protein-coding genes, primate protein-coding genes, and primate d-loop sequences. In all three cases, we found a measurable transition between the high, short-term (< 1-2 Myr) mutation rate and the low, long-term substitution rate. The relationship between the age of the calibration and the rate of change can be described by a vertically translated exponential decay curve, which may be used for correcting molecular date estimates. The phylogenetic substitution rates in mitochondria are approximately 0.5% per million years for avian protein-coding sequences and 1.5% per million years for primate protein-coding and d-loop sequences. Further analyses showed that purifying selection offers the most convincing explanation for the observed relationship between the estimated rate and the depth of the calibration. We rule out the possibility that it is a spurious result arising from sequence errors, and find it unlikely that the apparent decline in rates over time is caused by mutational saturation. Using a rate curve estimated from the d-loop data, several dates for last common ancestors were calculated: modern humans and Neandertals (354 ka; 222-705 ka), Neandertals (108 ka; 70-156 ka), and modern humans (76 ka; 47-110 ka). If the rate curve for a particular taxonomic group can be accurately estimated, it can be a useful tool for correcting divergence date estimates by taking the rate decay into account. Our results show that it is invalid to extrapolate molecular rates of change across different evolutionary timescales, which has important consequences for studies of populations, domestication, conservation genetics, and human evolution.     

Estimates of the genomic mutation rate for detrimental alleles in Drosophila melanogaster

The net rate of mutation to deleterious but nonlethal alleles and the sizes of effects of these mutations are of great significance for many evolutionary questions. Here we describe three replicate experiments in which mutations have been accumulated on chromosome 3 of Drosophila melanogaster by means of single-male backcrosses of heterozygotes for a wild-type third chromosome. Egg-to-adult viability was assayed for nonlethal homozygous chromosomes. The rates of decline in mean and increase in variance (DM and DV, respectively) were estimated. Scaled up to the diploid whole genome, the mean DM for homozygous detrimental mutations over the three experiments was between 0.8 and 1.8%. The corresponding DV estimate was approximately 0.11%. Overall, the results suggest a lower bound estimate of at least 12% for the diploid per genome mutation rate for detrimentals. The upper bound estimates for the mean selection coefficient were between 2 and 10%, depending on the method used. Mutations with selection coefficients of at least a few percent must be the major contributors to the effects detected here and are likely to be caused mostly by transposable element insertions or indels.     

Monophyly, divergence times, and evolution of host plant use inferred from a revised phylogeny of the Drosophila repleta species group

We present a revised molecular phylogeny of the Drosophila repleta group including 62 repleta group taxa and nine outgroup species based on four mitochondrial and six nuclear DNA sequence fragments. With ca. 100 species endemic to the New World, the repleta species group represents one of the major species radiations in the genus Drosophila. Most repleta group species are associated with cacti in arid or semiarid regions. Contrary to previous results, maximum likelihood and Bayesian phylogenies of the 10-gene dataset strongly support the monophyly of the repleta group. Several previously described subdivisions in the group were also recovered, despite poorly resolved relationships between these clades. Divergence time estimates suggested that the repleta group split from its sister group about 21millionyears ago (Mya), although diversification of the crown group began ca. 16Mya. Character mapping of patterns of host plant use showed that flat leaf Opuntia use is common throughout the phylogeny and that shifts in host use from Opuntia to the more chemically complex columnar cacti occurred several times independently during the history of this group. Although some species retained the use of Opuntia after acquiring the use of columnar cacti, there were multiple, phylogenetically independent instances of columnar cactus specialization with loss of Opuntia as a host. Concordant with our proposed timing of host use shifts, these dates are consistent with the suggested times when the Opuntioideae originated in South America. We discuss the generally accepted South American origin of the repleta group.     

Low impact of germline transposition on the rate of mildly deleterious mutation in Caenorhabditis elegans

Little is known about the role of transposable element (TE) insertion in the production of mutations with mild effects on fitness, the class of mutations thought to be central to the evolution of many basic features of natural populations. We propagated mutation-accumulation (MA) lines of two RNAi-deficient strains of Caenorhabditis elegans that exhibit germline transposition. We show here that the impact of TE activity was to raise the level of mildly deleterious mutation by 2- to 8.5-fold, as estimated from fecundity, longevity, and body length measurements, compared to that observed in a parallel MA experiment with a control strain characterized by a lack of germline transposition. Despite this increase, the rate of mildly deleterious mutation was between one and two orders of magnitude lower than the rate of TE accumulation, which was approximately two new insertions per genome per generation. This study suggests that high rates of TE activity do not necessarily translate into high rates of detectable nonlethal mutation.     

Chromosomal rearrangement inferred from comparisons of 12 Drosophila genomes

The availability of 12 complete genomes of various species of genus Drosophila provides a unique opportunity to analyze genome-scale chromosomal rearrangements among a group of closely related species. This article reports on the comparison of gene order between these 12 species and on the fixed rearrangement events that disrupt gene order. Three major themes are addressed: the conservation of syntenic blocks across species, the disruption of syntenic blocks (via chromosomal inversion events) and its relationship to the phylogenetic distribution of these species, and the rate of rearrangement events over evolutionary time. Comparison of syntenic blocks across this large genomic data set confirms that genetic elements are largely (95%) localized to the same Muller element across genus Drosophila species and paracentric inversions serve as the dominant mechanism for shuffling the order of genes along a chromosome. Gene-order scrambling between species is in accordance with the estimated evolutionary distances between them and we find it to approximate a linear process over time (linear to exponential with alternate divergence time estimates). We find the distribution of synteny segment sizes to be biased by a large number of small segments with comparatively fewer large segments. Our results provide estimated chromosomal evolution rates across this set of species on the basis of whole-genome synteny analysis, which are found to be higher than those previously reported. Identification of conserved syntenic blocks across these genomes suggests a large number of conserved blocks with varying levels of embryonic expression correlation in Drosophila melanogaster. On the other hand, an analysis of the disruption of syntenic blocks between species allowed the identification of fixed inversion breakpoints and estimates of breakpoint reuse and lineage-specific breakpoint event segregation.     

Pervasive Divergence of Transcriptional Gene Regulation in Caenorhabditis Nematodes

Because there is considerable variation in gene expression even between closely related species, it is clear that gene regulatory mechanisms evolve relatively rapidly. Because primary sequence conservation is an unreliable proxy for functional conservation of cis-regulatory elements, their assessment must be carried out in vivo. We conducted a survey of cis-regulatory conservation between C. elegans and closely related species C. briggsae, C. remanei, C. brenneri, and C. japonica. We tested enhancers of eight genes from these species by introducing them into C. elegans and analyzing the expression patterns they drove. Our results support several notable conclusions. Most exogenous cis elements direct expression in the same cells as their C. elegans orthologs, confirming gross conservation of regulatory mechanisms. However, the majority of exogenous elements, when placed in C. elegans, also directed expression in cells outside endogenous patterns, suggesting functional divergence. Recurrent ectopic expression of different promoters in the same C. elegans cells may reflect biases in the directions in which expression patterns can evolve due to shared regulatory logic of coexpressed genes. The fact that, despite differences between individual genes, several patterns repeatedly emerged from our survey, encourages us to think that general rules governing regulatory evolution may exist and be discoverable.     

Increase of the spontaneous mutation rate in a long-term experiment with Drosophila melanogaster

In a previous experiment, the effect of 255 generations of mutation accumulation (MA) on the second chromosome viability of Drosophila melanogaster was studied using 200 full-sib MA1 lines and a large C1 control, both derived from a genetically homogeneous base population. At generation 265, one of those MA1 lines was expanded to start 150 new full-sib MA2 lines and a new C2 large control. After 46 generations, the rate of decline in mean viability in MA2 was approximately 2.5 times that estimated in MA1, while the average degree of dominance of mutations was small and nonsignificant by generation 40 and moderate by generation 80. In parallel, the inbreeding depression rate for viability and the amount of additive variance for two bristle traits in C2 were 2-3 times larger than those in C1. The results are consistent with a mutation rate in the line from which MA2 and C2 were derived about 2.5 times larger than that in MA1. The mean viability of C2 remained roughly similar to that of C1, but the rate of MA2 line extinction increased progressively, leading to mutational collapse, which can be ascribed to accelerated mutation and/or synergy after important deleterious accumulation.