Molecular population genetics and the search for adaptive evolution in plants

The first papers on plant molecular population genetics were published approximately 10 years ago. Since that time, well over 50 additional studies of plant nucleotide polymorphism have been published, and many of these studies focused on detecting the signature of balancing or positive selection at a locus. In this review, we discuss some of the theoretical and statistical issues surrounding the detection of selection, with focus on plant populations, and we also summarize the empirical plant molecular population genetics literature. At face value, the literature suggests that a history of balancing or positive selection in plant genes is rampant. In two well-studied taxa (maize and Arabidopsis) over 20% of studied genes have been interpreted as containing the signature of selection. We argue that this is probably an overstatement of the prevalence of natural selection in plant genomes, for two reasons. First, demographic effects are difficult to incorporate and have generally not been well integrated into the plant population genetics literature. Second, the genes studied to date are not a random sample, so selected genes may be overrepresented. The next generation of studies in plant molecular population genetics requires additional sampling of local populations, explicit comparisons among loci, and improved theoretical methods to control for demography. Eventually, candidate loci should be confirmed by explicit consideration of phenotypic effects.     

Fertility selection,genetic selection and evolution

The relationship between fertility selection as measured by the correlation in progeny number between parents and offspring, and selection at individual loci is investigated in humans. Estimates for the magnitude of fertility selection (0.1) and the rate of gene substitution (0.5 gene substitutions per generation per genome) are used in various mathematical models for selection. It is found that the observed magnitude of fertility selection cannot be explained by non‐epistatic directional selection at individual loci. A symmetric quantitative directional selection model is consistent with the observed data. But it is possible that fertility selection does not have a genetic basis.     

Adaptive evolution of metabolic pathways in Drosophila

The adaptive significance of enzyme variation has been of central interest in population genetics. Yet, how natural selection operates on enzymes in the larger context of biochemical pathways has not been broadly explored. A basic expectation is that natural selection on metabolic phenotypes will target enzymes that control metabolic flux, but how adaptive variation is distributed among enzymes in metabolic networks is poorly understood. Here, we use population genetic methods to identify enzymes responding to adaptive selection in the pathways of central metabolism in Drosophila melanogaster and Drosophila simulans. We report polymorphism and divergence data for 17 genes that encode enzymes of 5 metabolic pathways that converge at glucose-6-phosphate (G6P). Deviations from neutral expectations were observed at five loci. Of the 10 genes that encode the enzymes of glycolysis, only aldolase (Ald) deviated from neutrality. The other 4 genes that were inconsistent with neutral evolution (glucose-6-phosphate dehydrogenase [G6pd]), phosphoglucomutase [Pgm], trehalose-6-phosphate synthetase [Tps1], and glucose-6phosphatase [G6pase] encode G6P branch point enzymes that catalyze reactions at the entry point to the pentose-phosphate, glycogenic, trehalose synthesis, and gluconeogenic pathways. We reconcile these results with population genetics theory and existing arguments on metabolic regulation and propose that the incidence of adaptive selection in this system is related to the distribution of flux control. The data suggest that adaptive evolution of G6P branch point enzymes may have special significance in metabolic adaptation.     

Adaptive molecular evolution of HINTW,a female-specific gene in birds

It is well established that many genes on the male-specific Y chromosome of organisms such as mammals are involved in male reproduction and may evolve rapidly because of positive selection on male reproductive traits. In contrast, very little is known about the function and evolution of W-linked genes restricted to the female genome of organisms with female heterogamety. For birds (males ZZ, females ZW), only one W-linked gene (HINTW) is sufficiently different from its Z-linked homolog to indicate a female-specific function. Here, we report that HINTW shows evidence of adaptive molecular evolution, implying strong positive selection for new functional properties in female birds. Moreover, because HINTW is expressed in the gonads of female birds just before sexual differentiation and is thus a candidate for sex determination, it suggests adaptive evolution related to female development. This provides the first example of Darwinian evolution of a gene restricted to the female genome of any organism. Given that HINTW exists in multiple copies on W, similar to some testis-specific genes amplified on mammalian Y, avian HINTW may thus potentially represent a female parallel to the organization and evolution of Y chromosome genes involved in male reproduction and development.     

Gene expression and molecular evolution

The combination of complete genome sequence information and estimates of mRNA abundances have begun to reveal causes of both silent and protein sequence evolution. Translational selection appears to explain patterns of synonymous codon usage in many prokaryotes as well as a number of eukaryotic model organisms (with the notable exception of vertebrates). Relationships between gene length and codon usage bias, however, remain unexplained. Intriguing correlations between expression patterns and protein divergence suggest some general mechanisms underlying protein evolution.     

Insights into chromosomal evolution of Cicadomorpha using fluorochrome staining and mapping 18S rRNA and H3 histone genes

The infraorder Cicadomorpha (Hemiptera) is a cosmopolitan species‐rich lineage of phytophagous insects. They have holocentric chromosomes and vary greatly in diploid number across families, with X0 as the predominant sex male mechanism. Here, we advance the understanding of chromosome mapping of repetitive elements of four families of cicadomorphan insects, the spittlebugs (Cercopidae), leafhoppers (Cicadellidae), and treehoppers (Aetalionidae and Membracidae). Sampled individuals from 19 species show considerable variation in diploid number, which may have originated from fusions between autosomes or between autosomes and the ancient X. The distribution of CMA₃⁺ blocks, primarily observed in low numbers in autosomal regions, was a conserved trait. Likewise, fluorescence in situ hybridization (FISH) mapping revealed mainly one locus per haploid genome for the 18S rRNA gene and for H3, each of which is located on distinct chromosomes. Despite the extensive variation in the number of autosomes and sex systems, the number of loci of ribosomal and H3 genes remained stable and may reflect the ancestral genome organizations in these groups. These results shed light on the chromosomal‐level organization in Cicadomorpha and provide new insights into the evolutionary history of karyotypes and repetitive elements in this diverse insect lineage.     

Adaptive protein evolution and regulatory divergence in Drosophila

Two recent studies demonstrated a positive correlation between divergence in gene expression and protein sequence in Drosophila. This correlation could be driven by positive selection or variation in functional constraint. To distinguish between these alternatives, we compared patterns of molecular evolution for 1,862 genes with two previously reported estimates of expression divergence in Drosophila. We found a slight negative trend (nonsignificant) between positive selection on protein sequence and divergence in expression levels between Drosophila melanogaster and Drosophila simulans. Conversely, shifts in expression patterns during Drosophila development showed a positive association with adaptive protein evolution, though as before the relationship was weak and not significant. Overall, we found no strong evidence for an increase in the incidence of positive selection on protein-coding regions in genes with divergent expression in Drosophila, suggesting that the previously reported positive association between protein and regulatory divergence primarily reflects variation in functional constraint.     

Selection and demographic history shape the molecular evolution of the gamete compatibility protein bindin in Pisaster sea stars

Reproductive compatibility proteins have been shown to evolve rapidly under positive selection leading to reproductive isolation, despite the potential homogenizing effects of gene flow. This process has been implicated in both primary divergence among conspecific populations and reinforcement during secondary contact; however, these two selective regimes can be difficult to discriminate from each other. Here, we describe the gene that encodes the gamete compatibility protein bindin for three sea star species in the genus Pisaster. First, we compare the full‐length bindin‐coding sequence among all three species and analyze the evolutionary relationships between the repetitive domains of the variable second bindin exon. The comparison suggests that concerted evolution of repetitive domains has an effect on bindin divergence among species and bindin variation within species. Second, we characterize population variation in the second bindin exon of two species: We show that positive selection acts on bindin variation in Pisaster ochraceus but not in Pisaster brevispinus, which is consistent with higher polyspermy risk in P. ochraceus. Third, we show that there is no significant genetic differentiation among populations and no apparent effect of sympatry with congeners that would suggest selection based on reinforcement. Fourth, we combine bindin and cytochrome c oxidase 1 data in isolation‐with‐migration models to estimate gene flow parameter values and explore the historical demographic context of our positive selection results. Our findings suggest that positive selection on bindin divergence among P. ochraceus alleles can be accounted for in part by relatively recent northward population expansions that may be coupled with the potential homogenizing effects of concerted evolution.     

The molecular evolution of acrosin in placental mammals

Acrosin is thought to fulfill several different roles in fertilization including that of a serine protease and in secondary zona pellucida (ZP) binding. However, acrosin's importance as a fertilization protein has been questioned. Especially since it was discovered that acrosin knockout mice are fertile. In this study, we explored the sites involved in serine protease activity and secondary binding. We also assessed conservation in functional sites across species and examined whether amino acid changes present in the human population have the potential to affect fertility. In addition, since many mammalian reproduction proteins have been found to evolve rapidly, we tested for positive selection. Sequences from 43 mammals from all 19 placental orders, which included a total of 828 nucleotides from acrosin exons 2, 3, 4, and a portion of exon 5, were obtained. We found that all sites of the serine catalytic triad as well as three other sites linked to catalytic activity were completely conserved. Five of six sites proposed to play a role in secondary binding were 100% conserved as basic residues. These results support an evolutionary conserved role for acrosin as a serine protease and secondary binding protein across placental mammals. We found statistically significant support for positive selection within acrosin, but no single amino acid site reached the significance level of P > 0.95 for inclusion within the category omega > 1. Based upon two amino acid mutation scoring systems, three out of seven human residue changing single nucleotide polymorphisms (SNPs) were found to be potentially protein-altering mutations.     

植物种群分子进化中对生境的适应——对荒漠植物遗传多样性研究结果的初步分析

长期以来,自然选择理论与中性理论对生物分子进化中的环境适应机理存在着激烈争论。目前,在植物种群分子进化中对生境适应的研究中正面临着一些难题：中性突变是分子水平进化的唯一原因,自然选择发挥主要作用的适应性进化是否存在于分子水平,选择与中性两种学说两种机制完全不同,如何才能将两者联系和统一起来,部分学者利用建立各种模型来描述自然选择对分子标记位点以及连锁序列的直接作用,如生态位宽度变异假设等。本研究小组以新疆阜康荒漠植物为研究对象,通过对两种重要荒漠植物遗传多样性的研究,分析两种植物各亚种群不同生境的生态因子与其遗传变异的关系,讨论生态位宽度变异假设,揭示遗传变异的产生与维持。中性论者与选择论者都试图解释生物环境适应与分子变异之间的关系。中性论和选择论是反映进化的两个侧面,它们不是绝对的,可以相互转化。     

Positive natural selection has driven the evolution of the Hsp70s in Diguetia spiders

Hsp70s are a ubiquitous family of highly conserved proteins. Hsp70s are chaperones and have important roles in both protein folding and thermotolerance. It has been widely assumed that Hsp70 sequence evolution is governed by the strong functional constraints imposed by its crucial cellular functions. In this study of cytosolic heat-inducible Hsp70s from three spider families, we have found clear evidence of positive natural selection altering Hsp70s in desert-dwelling and heat-loving Diguetidae spiders. These spiders are a small family restricted to deserts. They display heat-tolerant behaviours not seen in their closest relatives, the Pholcidae and Plectreuridae.     

Mitogenome evolution in ladybirds: Potential association with dietary adaptation

Dietary shifts can alter the relative availability of different nutrients and are therefore associated with metabolic adaptation in animals. The Coccinellidae (ladybirds) exhibits three major types of feeding habits and provides a useful model to study the effects of dietary changes on the evolution of mitogenomes, which encode proteins directly involved in energy metabolism. Here, mitogenomes of three coccinellid species were newly sequenced. These data were combined with other ten previously sequenced coccinellid mitogenomes to explore the relationship between mitogenome evolution and diets. Our results indicate that mitogenomic data can be effectively used to resolve phylogenetic relationships of Coccinellidae. Strong codon usage bias in coccinellid mitogenomes was predominantly determined by nucleotide composition. The 13 mitochondrial protein‐coding genes (PCGs) globally evolved under negative constraints, with some PCGs showing a stronger purifying selection. Six PCGs (nad3, nad4L, and nad5 from Complex I; cox1 and cox3 from Complex IV; and atp6 from Complex V) displayed signs of positive selection. Of these, adaptive changes in cox3 were potentially associated with metabolic differences resulting from dietary shifts in Coccinellidae. Our results provide insights into the adaptive evolution of coccinellid mitogenomes in response to both dietary shifts and other life history traits.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

Adaptive evolution of genes underlying schizophrenia

Schizophrenia poses an evolutionary-genetic paradox because it exhibits strongly negative fitness effects and high heritability, yet it persists at a prevalence of approximately 1% across all human cultures. Recent theory has proposed a resolution: that genetic liability to schizophrenia has evolved as a secondary consequence of selection for human cognitive traits. This hypothesis predicts that genes increasing the risk of this disorder have been subject to positive selection in the evolutionary history of humans and other primates. We evaluated this prediction using tests for recent selective sweeps in human populations and maximum-likelihood tests for selection during primate evolution. Significant evidence for positive selection was evident using one or both methods for 28 of 76 genes demonstrated to mediate liability to schizophrenia, including DISC1, DTNBP1 and NRG1, which exhibit especially strong and well-replicated functional and genetic links to this disorder. Strong evidence of non-neutral, accelerated evolution was found for DISC1, particularly for exon 2, the only coding region within the schizophrenia-associated haplotype. Additionally, genes associated with schizophrenia exhibited a statistically significant enrichment in their signals of positive selection in HapMap and PAML analyses of evolution along the human lineage, when compared with a control set of genes involved in neuronal activities. The selective forces underlying adaptive evolution of these genes remain largely unknown, but these findings provide convergent evidence consistent with the hypothesis that schizophrenia represents, in part, a maladaptive by-product of adaptive changes during human evolution.     

Sexual selection and the adaptive evolution of mammalian ejaculate proteins

An elevated rate of substitution characterizes the molecular evolution of reproductive proteins from a wide range of taxa. Although the selective pressures explaining this rapid evolution are yet to be resolved, recent evidence implicates sexual selection as a potentially important explanatory factor. To investigate this hypothesis, we sought evidence of a high rate of adaptive gene evolution linked to postcopulatory sexual selection in muroid rodents, a model vertebrate group displaying a broad range of mating systems. Specifically, we sequenced 7 genes from diverse rodents that are expressed in the testes, prostate, or seminal vesicles, products of which have the potential to act in sperm competition. We inferred positive Darwinian selection in these genes by estimation of the ratio of nonsynonymous (d(N), amino acid changing) to synonymous (d(S), amino acid retaining) substitution rates (omega = d(N)/d(S)). Next, we tested whether variation in this ratio among lineages could be attributed to interspecific variation in mating systems, as inferred from the variation in these rodents' relative testis sizes (RTS). Four of the 7 genes examined (Prm1, Sva, Acrv1, and Svs2, but not Svp2, Msmb, or Spink3) exhibit unambiguous evidence of positive selection. One of these, the seminal vesicle-derived protein Svs2, also shows some evidence for a concentration of positive selection in those lineages in which sperm competition is common. However, this was not a general trend among all the rodent genes we examined. Using the same methods, we then reanalyzed previously published data on 2 primate genes, SEMG1 and SEMG2. Although SEMG2 also shows evidence of positive selection concentrated in lineages subject to high levels of sperm competition, no such trend was found for SEMG1. Overall, despite a high rate of positive selection being a feature of many ejaculate proteins, these results indicate that the action of sexual selection potentially responsible for elevated evolutionary rates may be difficult to detect on a gene-by-gene basis. Although the extreme diversity of reproductive phenotypes exhibited in nature attests to the power of sexual selection, the extent to which this force predominates in driving the rapid molecular evolution of reproductive genes therefore remains to be determined.