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

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

Plant self-incompatibility in natural populations: a critical assessment of recent theoretical and empirical advances

Self-incompatibility systems in plants are genetic systems that prevent self-fertilization in hermaphrodites through recognition and rejection of pollen expressing the same allelic specificity as that expressed in the pistils. The evolutionary properties of these self-recognition systems have been revealed through a fascinating interplay between empirical advances and theoretical developments. In 1939, Wright suggested that the main evolutionary force driving the genetic and molecular properties of these systems was strong negative frequency-dependent selection acting on pollination success. The empirical observation of high allelic diversity at the self-incompatibility locus in several species, followed by the discovery of very high molecular divergence among alleles in all plant families where the locus has been identified, supported Wright's initial theoretical predictions as well as many of its later developments. In the last decade, however, advances in the molecular characterization of the incompatibility reaction and in the analysis of allelic frequencies and allelic divergence from natural populations have stimulated new theoretical investigations that challenged some important assumptions of Wright's model of gametophytic self-incompatibility. We here review some of these recent empirical and theoretical advances that investigated: (i) the hypothesis that S-alleles are selectively equivalent, and the evolutionary consequences of genetic interactions between alleles; (ii) the occurrence of frequency-dependent selection in female fertility; (iii) the evolutionary genetics of self-incompatibility systems in subdivided populations; (iv) the evolutionary implications of the self-incompatibility locus's genetic architecture; and (v) of its interactions with the genomic environment.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Can we continue to neglect genomic variation in introgression rates when inferring the history of speciation? A case study in a Mytilus hybrid zone

The use of molecular data to reconstruct the history of divergence and gene flow between populations of closely related taxa represents a challenging problem. It has been proposed that the long‐standing debate about the geography of speciation can be resolved by comparing the likelihoods of a model of isolation with migration and a model of secondary contact. However, data are commonly only fit to a model of isolation with migration and rarely tested against the secondary contact alternative. Furthermore, most demographic inference methods have neglected variation in introgression rates and assume that the gene flow parameter (Nm) is similar among loci. Here, we show that neglecting this source of variation can give misleading results. We analysed DNA sequences sampled from populations of the marine mussels, Mytilus edulis and M. galloprovincialis, across a well‐studied mosaic hybrid zone in Europe and evaluated various scenarios of speciation, with or without variation in introgression rates, using an Approximate Bayesian Computation (ABC) approach. Models with heterogeneous gene flow across loci always outperformed models assuming equal migration rates irrespective of the history of gene flow being considered. By incorporating this heterogeneity, the best‐supported scenario was a long period of allopatric isolation during the first three‐quarters of the time since divergence followed by secondary contact and introgression during the last quarter. By contrast, constraining migration to be homogeneous failed to discriminate among any of the different models of gene flow tested. Our simulations thus provide statistical support for the secondary contact scenario in the European Mytilus hybrid zone that the standard coalescent approach failed to confirm. Our results demonstrate that genomic variation in introgression rates can have profound impacts on the biological conclusions drawn from inference methods and needs to be incorporated in future studies.     

Spatial covariation of mutation and nonsynonymous substitution rates in vertebrate mitochondrial genomes

Mitochondrial genomes encode fundamental subunits of the basic energy producing machinery of eukaryotic cells that are under strong functional constraint. Paradoxically, these genes evolve rapidly in general, and there is substantial variation in evolutionary rates among genes within genomes. In order to investigate spatial variation in selection intensity, we conducted tests of neutrality using ratios of synonymous to nonsynonymous substitutions (dN/dS = omega) on numerous protein gene segments from fishes and mammals. Values of omega were very low for nearly all genomic regions. However, values of both omega and dN varied in a clinal pattern with increasing distance from the light-strand origin of replication. Spatial heterogeneity of nonsynonymous substitution rates exhibits a significantly positive correlation with variation in mutation rates that are related to the mode of mitochondrial DNA replication. The finding that nonsynonymous substitution rates are proportional to mutation rates is expected if a majority of substitutions are selectively neutral or slightly deleterious. Spatial patterns of among-gene variation in nonsynonymous rates were highly similar between fishes and mammals, suggesting that forces governing mitochondrial gene evolution have remained relatively constant over 450 Myr of vertebrate evolution. Conservation of substitution patterns despite major shifts in thermal habit and metabolic demands among taxa implicates a conserved replication mechanism controlling relative mutation rates as a major determinant of mitochondrial protein evolution.     

Modelling contemporary evolution in stickleback

During the past decade, two lines of research have advanced our understanding of micro‐evolution. On the one hand, a number of studies have generated evidence for strong selection on phenotypes ( Kingsolver et al. 2001 ) and the contemporary (sometimes deemed ‘rapid’) evolution of phenotypic traits ( Hendry & Kinnison 1999 ). On the other hand, other studies have sought to identify the genes that underlie ecologically important traits ( Ungerer et al. 2008 ). Over the next decade, micro‐evolutionists might expect considerable progress from the study of contemporary evolution at both the phenotypic and genetic level simultaneously. In this issue of Molecular Ecology, Le Rouzic et al. (2011) present a teaser for this approach. They examined contemporary evolution of an adaptive trait with a well‐studied genetic basis, the number of lateral plates, in threespine stickleback (Gasterosteus aculeatus L.). A time series of 20 years of change for this trait after introduction into a pond in Norway was compared with a similar time series of 12 years following the invasion of a lake in Alaska. Using a modelling approach, the authors then teased apart selection acting upon the phenotype and selection acting on a major effect gene. In both time series, selection was strong and consistent. The models suggested that selection could act directly on the phenotype, or through the gene’s pleiotropic effects.     

Effects of metabolic rate on protein evolution

Since the modern evolutionary synthesis was first proposed early in the twentieth century, attention has focused on assessing the relative contribution of mutation versus natural selection on protein evolution. Here we test a model that yields general quantitative predictions on rates of protein evolution by combining principles of individual energetics with Kimura's neutral theory. The model successfully predicts much of the heterogeneity in rates of protein evolution for diverse eukaryotes (i.e. fishes, amphibians, reptiles, birds, mammals) from different thermal environments. Data also show that the ratio of non-synonymous to synonymous nucleotide substitution is independent of body size, and thus presumably of effective population size. These findings indicate that rates of protein evolution are largely controlled by mutation rates, which in turn are strongly influenced by individual metabolic rate.     

Phylogenetic network and physicochemical properties of nonsynonymous mutations in the protein-coding genes of human mitochondrial DNA

Theories on molecular evolution predict that phylogenetically recent nonsynonymous mutations should contain more non-neutral amino acid replacements than ancient mutations. We analyzed 840 complete coding-region human mitochondrial DNA (mtDNA) sequences for nonsynonymous mutations and evaluated the mutations in terms of the physicochemical properties of the amino acids involved. We identified 465 distinct missense and 6 nonsense mutations. 48% of the amino acid replacements changed polarity, 26% size, 8% charge, 32% aliphaticity, 13% aromaticity, and 44% hydropathy. The reduced-median networks of the amino acid changes revealed relatively few differences between the major continent-specific haplogroups, but a high variation and highly starlike phylogenies within the haplogroups. Some 56% of the mutations were private, and 25% were homoplasic. Nonconservative changes were more common than expected among the private mutations but less common among the homoplasic mutations. The asymptotic maximum of the number of nonsynonymous mutations in European mtDNA was estimated to be 1,081. The results suggested that amino acid replacements in the periphery of phylogenetic networks are more deleterious than those in the central parts, indicating that purifying selection prevents the fixation of some alleles.     

Comparative genomics and the study of evolution by natural selection

Genomics profoundly affects most areas of biology, including ecology and evolutionary biology. By examining genome sequences from multiple species, comparative genomics offers new insight into genome evolution and the way natural selection moulds DNA sequence evolution. Functional divergence, as manifested in the accumulation of nonsynonymous substitutions in protein-coding genes, differs among lineages in a manner seemingly related to population size. For example, the ratio of nonsynonymous to synonymous substitution (d_N/d_S) is higher in apes than in rodents, compatible with Ohta's nearly neutral theory of molecular evolution, which suggests that the fixation of slightly deleterious mutations contributes to protein evolution at an extent negatively correlated with effective population size. While this supports the idea that functional evolution is not necessarily adaptive, comparative genomics is uncovering a role for positive Darwinian selection in 10–40% of all genes in different lineages, estimates that are likely to increase when the addition of more genomes gives increased power. Again, population size seems to matter also in this context, with a higher proportion of fixed amino acid changes representing advantageous mutations in large populations. Genes that are particularly prone to be driven by positive selection include those involved with reproduction, immune response, sensory perception and apoptosis. Genetic innovations are also frequently obtained by the gain or loss of complete gene sequences. Moreover, it is increasingly realized, from comparative genomics, that purifying selection conserves much more than just the protein-coding part of the genome, and this points at an important role for regulatory elements in trait evolution. Finally, genome sequencing using outbred or multiple individuals has provided a wealth of polymorphism data that gives information on population history, demography and marker evolution.