Positive and negative selection in murine ultraconserved noncoding elements

There are many more selectively constrained noncoding than coding nucleotides in the mammalian genome, but most mammalian noncoding DNA is subject to weak selection, on average. One of the most striking discoveries to have emerged from comparisons among mammalian genomes is the hundreds of noncoding elements of more than 200 bp in length that show absolute conservation among mammalian orders. These elements represent the tip of the iceberg of a much larger class of conserved noncoding elements (CNEs). Much evidence suggests that CNEs are selectively constrained and not mutational cold-spots, and there is evidence that some CNEs play a role in the regulation of development. Here, we quantify negative and positive selection acting in murine CNEs by analyzing within-species nucleotide variation and between-species divergence of CNEs that we identified using a phylogenetically independent comparison. The distribution of fitness effects of new mutations in CNEs, inferred from within-species polymorphism, suggests that CNEs receive a higher number of strongly selected deleterious mutations and many fewer nearly neutral mutations than amino acid sites of protein-coding genes or regulatory elements close to genes. However, we also show that CNEs experience a far higher proportion of adaptive substitutions than any known category of genomic sites in murids. The absolute rate of adaptation of CNEs is similar to that of amino acid sites of proteins. This result suggests that there is widespread adaptation in mammalian conserved noncoding DNA elements, some of which have been implicated in the regulation of crucially important processes, including development.     

Positive and negative selection on noncoding DNA close to protein-coding genes in wild house mice

During the past two decades, evidence has accumulated of adaptive evolution within protein-coding genes in a variety of species. However, with the exception of Drosophila and humans, little is known about the extent of adaptive evolution in noncoding DNA. Here, we study regions upstream and downstream of protein-coding genes in the house mouse Mus musculus castaneus, a species that has a much larger effective population size (N(e)) than humans. We analyze polymorphism data for 78 genes from 15 wild-caught M. m. castaneus individuals and divergence to a closely related species, Mus famulus. We find high levels of nucleotide diversity and moderate levels of selective constraint in upstream and downstream regions compared with nonsynonymous sites of protein-coding genes. From the polymorphism data, we estimate the distribution of fitness effects (DFE) of new mutations and infer that most new mutations in upstream and downstream regions behave as effectively neutral and that only a small fraction is strongly negatively selected. We also estimate the fraction of substitutions that have been driven to fixation by positive selection (α) and the ratio of adaptive to neutral divergence (ω(α)). We find that α for upstream and downstream regions (～ 10%) is much lower than α for nonsynonymous sites (～ 50%). However, ω(α) estimates are very similar for nonsynonymous sites (～ 10%) and upstream and downstream regions (～ 5%). We conclude that negative selection operating in upstream and downstream regions of M. m. castaneus is weak and that the low values of α for upstream and downstream regions relative to nonsynonymous sites are most likely due to the presence of a higher proportion of neutrally evolving sites and not due to lower absolute rates of adaptive substitution.     

Positive selection versus demography: evolutionary inferences based on an unusual haplotype structure in Drosophila simulans

Coalescent simulations were used to investigate the possible role of population subdivision and history in shaping nucleotide variation in a recombining 88-kb genomic fragment of Drosophila simulans displaying an unusual large-scale haplotype structure. The multilocus analysis, based on summary statistics using specific demographic null models under recombination, indicates that the observed levels of linkage disequilibrium differed significantly from the values expected under different bottleneck and population admixture scenarios. These results indicate that demography alone may not account for the observed pattern of variation and support the previous claim that the data are better described by a model in which an adaptive mutation has not yet gone to fixation.     

Testing models of selection and demography in Drosophila simulans

We analyze patterns of nucleotide variability at 15 X-linked loci and 14 autosomal loci from a North American population of Drosophila simulans. We show that there is significantly more linkage disequilibrium on the X chromosome than on chromosome arm 3R and much more linkage disequilibrium on both chromosomes than expected from estimates of recombination rates, mutation rates, and levels of diversity. To explore what types of evolutionary models might explain this observation, we examine a model of recurrent, nonoverlapping selective sweeps and a model of a recent drastic bottleneck (e.g., founder event) in the demographic history of North American populations of D. simulans. The simple sweep model is not consistent with the observed patterns of linkage disequilibrium nor with the observed frequencies of segregating mutations. Under a restricted range of parameter values, a simple bottleneck model is consistent with multiple facets of the data. While our results do not exclude some influence of selection on X vs. autosome variability levels, they suggest that demography alone may account for patterns of linkage disequilibrium and the frequency spectrum of segregating mutations in this population of D. simulans.     

Comparative Genomics of Mitochondrial DNA in Drosophila simulans

The current study compares the nucleotide variation among 22 complete mitochondrial genomes of the three distinct Drosophila simulans haplotypes with intron 1 of the alcohol dehydrogenase-related locus. This is the first study to investigate the sequence variation of multiple complete mitochondrial genomes within distinct mitochondrial haplotypes of a single species. Patterns of variation suggest distinct forces are influencing the evolution of mitochondrial DNA (mtDNA) and autosomal DNA in D. simulans. First, there is little variation within each mtDNA haplotype but strong differentiation among them. In contrast, there is no support for differentiation of the mitochondrial haplotypes at the autosomal locus. Second, there is a significant deficiency of mitochondrial variation in each haplotype relative to the autosomal locus. Third, the ratio of nonsynonymous to synonymous substitutions is not equal in all branches of the well-resolved phylogeny. There is an excess of nonsynonymous substitutions relative to synonymous substitutions within each D. simulans haplotype. This result is similar to that previously observed within the mtDNA of distinct species. A single evolutionary force may be causally linked to the observed patterns of mtDNA variation—a rickettsia-like microorganism, Wolbachia pipientis, which is known to directly influence mitochondrial evolution but have a less direct influence on autosomal loci. Received: 16 September 1999 / Accepted: 14 March 2000     

Positive and negative selection on the mitochondrial genome

Several recent studies have confirmed that mitochondrial DNA variation and evolution are not consistent with the neutral theory of molecular evolution and might be inappropriate for estimating effective population sizes. Evidence for the action of both positive and negative selection on mitochondrial genes has been put forward, and the complex genetics of mitochondrial DNA adds to the challenge of resolving this debate. The solution could lie in distinguishing genetic drift from 'genetic draft' and in dissecting the physiology of mitochondrial fitness.     

Positive and negative selection on mammalian Y chromosomes

Y chromosomes are genetically degenerate in most organisms studied. The loss of genes from Y chromosomes is thought to be due to the inefficiency of purifying selection in nonrecombining regions, which leads to the accumulation of deleterious mutations via the processes of hitchhiking, background selection, and Muller's ratchet. As the severity of these processes depends on the number of functional genes linked together on the nonrecombining Y, it is not clear whether these processes are still at work on the old, gene-poor mammalian Y chromosomes. If purifying selection is indeed less efficient in the Y-linked, compared to the X-linked genes, deleterious nonsynonymous substitutions are expected to accumulate faster on the Y chromosome. However, positive selection on Y-linked genes could also increase the rate of amino acid-changing substitutions. Thus, the previous reports of an elevated nonsynonymous substitution rate in Y-linked genes are still open to interpretation. Here, we report evidence for positive selection in two out of three studied mammalian Y-linked genes, suggesting that adaptive Darwinian evolution may be common on mammalian Y chromosomes. Taking positive selection into account, we demonstrate that purifying selection is less efficient in mammalian Y-linked genes compared to their X-linked homologues, suggesting that these genes continue to degenerate.     

Purifying selection maintains highly conserved noncoding sequences in Drosophila

The majority of metazoan genomes consist of nonprotein-coding regions, although the functional significance of most noncoding DNA sequences remains unknown. Highly conserved noncoding sequences (CNSs) have proven to be reliable indicators of functionally constrained sequences such as cis-regulatory elements and noncoding RNA genes. However, CNSs may arise from nonselective evolutionary processes such as genomic regions with extremely low mutation rates known as mutation "cold spots." Here we combine comparative genomic data from recently completed insect genome projects with population genetic data in Drosophila melanogaster to test predictions of the mutational cold spot model of CNS evolution in the genus Drosophila. We find that point mutations in intronic and intergenic CNSs exhibit a significant reduction in levels of divergence relative to levels of polymorphism, as well as a significant excess of rare derived alleles, compared with either the nonconserved spacer regions between CNSs or with 4-fold silent sites in coding regions. Controlling for the effects of purifying selection, we find no evidence of positive selection acting on Drosophila CNSs, although we do find evidence for the action of recurrent positive selection in the spacer regions between CNSs. We estimate that approximately 85% of sites in Drosophila CNSs are under constraint with selection coefficients (N(e)s) on the order of 10-100, and thus, the estimated strength and number of sites under purifying selection is greater for Drosophila CNSs relative to those in the human genome. These patterns of nonneutral molecular evolution are incompatible with the mutational cold spot hypothesis to explain the existence of CNSs in Drosophila and, coupled with similar findings in mammals, argue against the general likelihood that CNSs are generated by mutational cold spots in any metazoan genome.     

Direct and correlated responses to selection on age at physiological maturity in Drosophila simulans

Biologists who study the timing of development in insects have focused on variation in duration of pre‐adult stages almost without exception. However, development is not complete until adults are not only morphologically mature, but also reproductively mature. Here we describe an experiment in the fruit fly, Drosophila simulans, in which we used artificial selection to create lines with shortened and lengthened duration from eclosion to the age when the first egg was laid. We found significant genetic variation for this trait. The response to selection on age when the first egg was laid was due to variation among females. Delayed adult development was correlated with rapid pre‐adult development and longer life span in females. The approach we use here resolves some difficulties with previous approaches used to study the genetics of senescence, and provides an opportunity to study the hitherto unexamined predictions derived from classic models for the evolution of senescence.     

Positive and negative DNA elements of the Drosophila grimshawi s18 chorion gene assayed in Drosophila melanogaster.

Germ line transformation has been used to map the cis regulatory DNA elements responsible for the precise and evolutionarily stable developmental expression of the s18 chorion gene. Constructs containing chimeric combinations of Drosophila melanogaster and D. grimshawi DNA regions, as well as D. grimshawi sequences alone, can direct expression in the follicular epithelium, in an s18-specific temporal and spatial pattern. The results indicate that both positive and negative regulatory elements can function when transferred from D. grimshawi to D. melanogaster. The first ca. 100 bp of the 5'-flanking DNA region constitute a minimal, developmentally regulated promoter, expression of which is inhibited by the next 100-bp DNA segment and activated by positive elements located further upstream. Expression of the minimal promoter can also be enhanced by more distant chorion regulatory elements, provided the inhibitory DNA segment is absent.     

Positive and negative selection on the human genome.

The distinction between deleterious, neutral, and adaptive mutations is a fundamental problem in the study of molecular evolution. Two significant quantities are the fraction of DNA variation in natural populations that is deleterious and destined to be eliminated and the fraction of fixed differences between species driven by positive Darwinian selection. We estimate these quantities using the large number of human genes for which there are polymorphism and divergence data. The fraction of amino acid mutations that is neutral is estimated to be 0.20 from the ratio of common amino acid (A) to synonymous (S) single nucleotide polymorphisms (SNPs) at frequencies of > or =15%. Among the 80% of amino acid mutations that are deleterious at least 20% of them are only slightly deleterious and often attain frequencies of 1-10%. We estimate that these slightly deleterious mutations comprise at least 3% of amino acid SNPs in the average individual or at least 300 per diploid genome. This estimate is not sensitive to human population history. The A/S ratio of fixed differences is greater than that of common SNPs and suggests that a large fraction of protein divergence is adaptive and driven by positive Darwinian selection.     

Phylogenetic models and model selection for noncoding DNA

Awareness of the complex structure and evolutionary dynamics of noncoding DNA has improved both noncoding sequence alignment and the use of microstructural changes as characters in phylogenetic analysis. The next step is to consider improvements in the use and selection of phylogenetic models for noncoding sequence data. Models of character evolution are central to phylogeny estimation, but the use of an inadequate model can mislead topology selection and branch length estimations. This is particularly likely when sequence divergence is either limited (nearly invariable, as in population-level or species-level studies) or extreme (nearly saturated, as in deep-level studies that focus on conserved secondary structures). Noncoding data sets are often at these extremes, and they can be particularly awkward for model definition and model selection. This paper introduces the goals of model use in phylogenetics and identifies ten issues that arise from the application of models to noncoding sequence data. It is concluded that most of these issues derive from small data set sizes, very low or very high sequence variability, limitations of current phylogenetic models, and possibly character definition and nonindependence. Recommendations are made that should help to improve alignment, character quality, model selection, and phylogeny estimation based on noncoding sequence data.     

Linkage disequilibrium and recent selection at three immunity receptor loci in Drosophila simulans

Immune system genes in a California population sample of Drosophila simulans were shown to bear several hallmarks of the effects of past directional selection. One potential effect of directional selection is an increase in linkage disequilibrium among the polymorphic sites that are linked to the site under selection. In this study, we focus on three D. simulans immunity loci, Hmu, Sr-CI/Sr-CIII, and Tehao, for which the polymorphic sites are in nearly perfect linkage disequilibrium, an unusual finding even with respect to other immunity genes sampled from the same lines. The most likely explanation for this finding is that, at each locus, two divergent alleles have been selected to intermediate frequencies in the recent past. The extent to which the linkage disequilibrium extends to the flanks of each of the immunity genes is minimal, suggesting that the favored mutations actually occurred within the immunity genes themselves. Furthermore, the excess linkage disequilibrium found in the California population is not found in an African D. simulans population sample and may be a result of novel pathogen-mediated selection pressures encountered during establishment of non-African populations.     

Environmental heterogeneity,multivariate sexual selection and genetic constraints on cuticular hydrocarbons in Drosophila simulans

Sexual selection is responsible for the evolution of many elaborate traits, but sexual trait evolution could be influenced by opposing natural selection as well as genetic constraints. As such, the evolution of sexual traits could depend heavily on the environment if trait expression and attractiveness vary between environments. Here, male Drosophila simulans were reared across a range of diets and temperatures, and we examined differences between these environments in terms of (i) the expression of male cuticular hydrocarbons (CHCs) and (ii) which male CHC profiles were most attractive to females. Temperature had a strong effect on male CHC expression, whereas the effect of diet was weaker. Male CHCs were subject to complex patterns of directional, quadratic and correlational sexual selection, and we found differences between environments in the combination of male CHCs that were most attractive to females, with clearer differences between diets than between temperatures. We also show that genetic covariance between environments is likely to cause a constraint on independent CHC evolution between environments. Our results demonstrate that even across the narrow range of environmental variation studied here, predicting the outcome of sexual selection can be extremely complicated, suggesting that studies ignoring multiple traits or environments may provide an over‐simplified view of the evolution of sexual traits.     

Insertional DNA and spontaneous mutation at the white locus in Drosophila simulans

A large body of data on molecular analyses of several multiallelic loci in Drosophila melanogaster has demonstrated a high incidence of mobile DNA element insertions among spontaneous mutations. In the sibling species D. simulans, the dispersed, middle repetitive, nomadic sequences are reduced to about one-seventh that of its sibling species (Dowsett and Young 1982). Does this reduced amount of middle repetitive DNA (or mobile DNA sequences) mean that in D. simulans the occurrence of insertion mutants will be rare compared with that of D. melanogaster? To test this possibility, we collected seven different spontaneous white mutants of D. simulans and studied their molecular gene structures. Five out of seven mutants had insertion sequences which varied in length from 0.4 kb to 16 kb. One bore a deletion spanning the w region and another showed no gross structural alteration. Thus the proportion of insertional mutations at the white locus in D. simulans is equivalent to that observed in D. melanogaster. Among the five insertional mutants, one, wmky, showed genetic instability; the other four were stable. wmky was found to mutate at a frequency of 2.1 x 10(-5) in meiotic cells and may also be unstable in somatic cells.     

Positive and negative selection in the DAZ gene family

Because a microdeletion containing the DAZ gene is the most frequently observed deletion in infertile men, the DAZ gene was considered a strong candidate for the azoospermia factor. A recent evolutionary analysis, however, suggested that DAZ was free from functional constraints and consequently played little or no role in human spermatogenesis. The major evidence for this surprising conclusion is that the nonsynonymous substitution rate is similar to the synonymous rate and to the rate in introns. In this study, we reexamined the evolution of the DAZ gene family by using maximum-likelihood methods, which accommodate variable selective pressures among sites or among branches. The results suggest that DAZ is not free from functional constraints. Most amino acids in DAZ are under strong selective constraint, while a few sites are under diversifying selection with nonsynonymous/ synonymous rate ratios (d(N)/d(S)) well above 1. As a result, the average d(N)/d(S) ratio over sites is not a sensible measure of selective pressure on the protein. Lineage-specific analysis indicated that human members of this gene family were evolving by positive Darwinian selection, although the evidence was not strong.     

Factors affecting mitochondrial DNA quality from museum preserved Drosophila simulans

In this study we investigate how traditional killing methods and storage regimes affected mitochondrial DNA quality in Drosophila simulans. Here we define quality with three criteria: (1) size of extracted DNA, (2) extraction yield, and (3) ability to amplify from four target regions. Killing methods had a significant effect on extraction yield, but not on PCR success. Highest DNA yields were extracted from specimens exposed to cyanide, while the lowest were from specimens killed in 70% ethanol. Specimens stored for two years contained badly sheared DNA, which translated into a significant decrease in PCR success compared to freshly assayed specimens. The most dramatic decrease in PCR success occurred in the 1822 bp and 1332 bp amplicons, compared to the 959 bp and 291 bp fragments. Naphthalene did not affect any aspect of DNA quality; time of storage affected PCR success regardless of naphthalene environment. This study serves to further refine our understanding of DNA degradation.     

Patterns of DNA sequence variation suggest the recent action of positive selection in the janus-ocnus region of Drosophila simulans

Levels of nucleotide polymorphism in three paralogous Drosophila simulans genes, janusA (janA), janusB (janB), and ocnus (ocn), were surveyed by DNA sequencing. The three genes lie in tandem within a 2.5-kb region of chromosome arm 3R. In a sample of eight alleles from a worldwide distribution we found a significant departure from neutrality by several statistical tests. The most striking feature of this sample was that in a 1.7-kb region containing the janA and janB genes, 30 out of 31 segregating sites contained variants present only once in the sample, and 29 of these unique variants were found in the same allele. A restriction survey of an additional 28 lines of D. simulans revealed strong linkage disequilibrium over the janA-janB region and identified six more alleles matching the rare haplotype. Among the rare alleles, the level of DNA sequence variation was typical for D. simulans autosomal genes and showed no departure from neutrality. In addition, the rare haplotype was more similar to the D. melanogaster sequence, indicating that it was the ancestral form. These results suggest that the derived haplotype has risen to high worldwide frequency relatively recently, most likely as a result of natural selection.     

Natural selection and ribosomal DNA in Drosophila

Natural populations of Drosophila mercatorum are variable for the number of X-linked 28S ribosomal genes bearing a 5-kilobase insert. A separate polymorphic X-linked gene controls whether 28S repeats bearing the insert are preferentially underreplicated during the formation of polytene tissue. Female flies having at least a third of their 28S genes bearing the insert and lacking the ability to preferentially underreplicate inserted repeats display the abnormal abdomen syndrome. The syndrome is characterized by retention of juvenile abdominal cuticle into the adult, a slowdown in larval developmental time, and an increase in early female fecundity. The life history traits are expressed in nature and provide a basis for strong natural selection. The abnormal abdomen syndrome should be favored whenever the adult age structure is skewed towards young individuals, and field studies confirm this prediction. The closely related species, Drosophila hydei, also bears these inserts and appears to be subject to similar selection. However, D. mercatorum responds to this selection primarily through the allelic variation that controls preferential underreplication, whereas D. hydei responds primarily through adjustment of the proportion of inserted 28S genes. This is interpreted to mean that the evolution of a multigene family arises from the interaction of population-level and DNA-level processes.