Testing the neutral theory of molecular evolution with genomic data from Drosophila

Although positive selection has been detected in many genes, its overall contribution to protein evolution is debatable. If the bulk of molecular evolution is neutral, then the ratio of amino-acid (A) to synonymous (S) polymorphism should, on average, equal that of divergence. A comparison of the A/S ratio of polymorphism in Drosophila melanogaster with that of divergence from Drosophila simulans shows that the A/S ratio of divergence is twice as high---a difference that is often attributed to positive selection. But an increase in selective constraint owing to an increase in effective population size could also explain this observation, and, if so, all genes should be affected similarly. Here we show that the difference between polymorphism and divergence is limited to only a fraction of the genes, which are also evolving more rapidly, and this implies that positive selection is responsible. A higher A/S ratio of divergence than of polymorphism is also observed in other species, which suggests a rate of adaptive evolution that is far higher than permitted by the neutral theory of molecular evolution.     

Insights into hominid evolution from the gorilla genome sequence

Gorillas are humans' closest living relatives after chimpanzees, and are of comparable importance for the study of human origins and evolution. Here we present the assembly and analysis of a genome sequence for the western lowland gorilla, and compare the whole genomes of all extant great ape genera. We propose a synthesis of genetic and fossil evidence consistent with placing the human-chimpanzee and human-chimpanzee-gorilla speciation events at approximately 6 and 10 million years ago. In 30% of the genome, gorilla is closer to human or chimpanzee than the latter are to each other; this is rarer around coding genes, indicating pervasive selection throughout great ape evolution, and has functional consequences in gene expression. A comparison of protein coding genes reveals approximately 500 genes showing accelerated evolution on each of the gorilla, human and chimpanzee lineages, and evidence for parallel acceleration, particularly of genes involved in hearing. We also compare the western and eastern gorilla species, estimating an average sequence divergence time 1.75 million years ago, but with evidence for more recent genetic exchange and a population bottleneck in the eastern species. The use of the genome sequence in these and future analyses will promote a deeper understanding of great ape biology and evolution.     

Adaptive evolution drives divergence of a hybrid inviability gene between two species of Drosophila

Speciation--the splitting of one species into two--occurs by the evolution of any of several forms of reproductive isolation between taxa, including the intrinsic sterility and inviability of hybrids. Abundant evidence shows that these hybrid fitness problems are caused by incompatible interactions between loci: new alleles that become established in one species are sometimes functionally incompatible with alleles at interacting loci from another species. However, almost nothing is known about the genes involved in such hybrid incompatibilities or the evolutionary forces that drive their divergence. Here we identify a gene that causes epistatic inviability in hybrids between two fruitfly species, Drosophila melanogaster and D. simulans. Our population genetic analysis reveals that this gene--which encodes a nuclear pore protein--evolved by positive natural selection in both species' lineages. These results show that a lethal hybrid incompatibility has evolved as a by-product of adaptive protein evolution.     

MHC II DRB variation and trans-species polymorphism in the golden snub-nosed monkey (Rhinopithecus roxellana)

Genetic variation is generally believed to be important in studying endangered species’ adaptive potential.Early studies assessed genetic diversity using nearly neutral markers,such as microsatellite loci and mitochondrial DNA(mtDNA),which are very informative for phylogenetic and phylogeographic reconstructions.However,the variation at these loci cannot provide direct information on selective processes involving the interaction of individuals with their environment,or on the capability to resist continuously evolving pathogens and parasites.The importance of genetic diversity at informative adaptive markers,such as major histocompatibility complex(MHC) genes,is increasingly being realized,especially in endangered,isolated species.Small population size and isolation make the golden snub-nosed monkey(Rhinopithecus roxellana) particularly susceptible to genetic variation losses through inbreeding and restricted gene flow.In this study,we compared the genetic variation and population structure of microsatellites,mtDNA,and the most relevant adaptive region of the MHC II-DRB genes in the golden snub-nosed monkey.We examined three Chinese R.roxellana populations and found the same variation patterns in all gene regions,with the population from Shennongjia population,Hubei Province,showing the lowest polymorphism among three populations.Genetic drift that outweighed balancing selection and the founder effect in these populations may explain the similar genetic variation pattern found in these neutral and adaptive genes.     

Human-specific loss of regulatory DNA and the evolution of human-specific traits

Humans differ from other animals in many aspects of anatomy, physiology, and behaviour; however, the genotypic basis of most human-specific traits remains unknown. Recent whole-genome comparisons have made it possible to identify genes with elevated rates of amino acid change or divergent expression in humans, and non-coding sequences with accelerated base pair changes. Regulatory alterations may be particularly likely to produce phenotypic effects while preserving viability, and are known to underlie interesting evolutionary differences in other species. Here we identify molecular events particularly likely to produce significant regulatory changes in humans: complete deletion of sequences otherwise highly conserved between chimpanzees and other mammals. We confirm 510 such deletions in humans, which fall almost exclusively in non-coding regions and are enriched near genes involved in steroid hormone signalling and neural function. One deletion removes a sensory vibrissae and penile spine enhancer from the human androgen receptor (AR) gene, a molecular change correlated with anatomical loss of androgen-dependent sensory vibrissae and penile spines in the human lineage. Another deletion removes a forebrain subventricular zone enhancer near the tumour suppressor gene growth arrest and DNA-damage-inducible, gamma (GADD45G), a loss correlated with expansion of specific brain regions in humans. Deletions of tissue-specific enhancers may thus accompany both loss and gain traits in the human lineage, and provide specific examples of the kinds of regulatory alterations and inactivation events long proposed to have an important role in human evolutionary divergence.     

Genetic evidence for complex speciation of humans and chimpanzees

The genetic divergence time between two species varies substantially across the genome, conveying important information about the timing and process of speciation. Here we develop a framework for studying this variation and apply it to about 20 million base pairs of aligned sequence from humans, chimpanzees, gorillas and more distantly related primates. Human-chimpanzee genetic divergence varies from less than 84% to more than 147% of the average, a range of more than 4 million years. Our analysis also shows that human-chimpanzee speciation occurred less than 6.3 million years ago and probably more recently, conflicting with some interpretations of ancient fossils. Most strikingly, chromosome X shows an extremely young genetic divergence time, close to the genome minimum along nearly its entire length. These unexpected features would be explained if the human and chimpanzee lineages initially diverged, then later exchanged genes before separating permanently.     

Phylogenetic origins and adaptive evolution of avian and mammalian haemoglobin genes

Recent years have seen rapid growth in amino acid sequence data on globins and nucleotide sequence data on haemoglobin genes and pseudogenes, and cladistic analysis of these data continues to reveal new facets of globin evolution. Our present findings demonstrate: (1) avian and mammalian embryonic alpha genes (pi and xi, respectively) had a monophyletic origin involving an alpha locus duplication about 400 Myr ago soon after the duplication which separated alpha and beta genes; (2) much later in phylogeny, independent beta-gene duplications produced the embryonic rho locus of birds and embryonic epsilon and fetal gamma loci of mammals. This parallels the earlier finding that myoglobins evolved more than once from generalized globin ancestors. Here we support the view that such globin evolution resulted from natural selection acting on mutations in duplicated genes. Thus, our evidence contradicts the neutralist view in which almost all amino acid substitutions in descent to extant globins evaded positive selection.     

Complex speciation of humans and chimpanzees

Genetic data from two or more species provide information about the process of speciation. In their analysis of DNA from humans, chimpanzees, gorillas, orangutans and macaques (HCGOM), Patterson et al. suggest that the apparently short divergence time between humans and chimpanzees on the X chromosome is explained by a massive interspecific hybridization event in the ancestry of these two species. However, Patterson et al. do not statistically test their own null model of simple speciation before concluding that speciation was complex, and--even if the null model could be rejected--they do not consider other explanations of a short divergence time on the X chromosome. These include natural selection on the X chromosome in the common ancestor of humans and chimpanzees, changes in the ratio of male-to-female mutation rates over time, and less extreme versions of divergence with gene flow (see ref. 2, for example). I therefore believe that their claim of hybridization is unwarranted.     

禾本科植物β-淀粉酶基因家族分子进化及响应非生物胁迫的表达模式分析

 β-淀粉酶(beta-amylase,BAM)是一类关键的淀粉水解酶,在禾谷类作物生长发育过程中起着重要作用,与植物多种非生物胁迫响应相关.本研究通过系统发育分析,将水稻、玉米、高粱、谷子、二穗短柄草5 种禾本科植物中共54 个BAM 基因分为10 个同源基因簇,每个同源基因簇都涵盖了这5 个物种,因此推测在禾本科祖先物种中至少含有10 个BAM 基因,并且在禾本科植物分化后没有发生明显的基因丢失事件.基于对编码蛋白质序列的功能分化分析,表明同源基因簇间存在明显的进化速率的差异.对10 个同源基因簇进行了适应性进化检测,发现有3 个同源簇在禾本科植物的进化过程中经历了适应性进化.此外,对水稻β-淀粉酶的表达分析发现,一些β-淀粉酶具有组织特异性表达特征,并且至少有5 个水稻的β-淀粉酶基因具有受到非生物逆境的胁迫而表现出不同的表达模式.本研究结果为进一步探讨禾本科BAM 基因的生物学功能提供了一定的理论基础.     

The cost of inbreeding in Arabidopsis

Population geneticists have long sought to estimate the distribution of selection intensities among genes of diverse function across the genome. Only recently have DNA sequencing and analytical techniques converged to make this possible. Important advances have come from comparing genetic variation within species (polymorphism) with fixed differences between species (divergence). These approaches have been used to examine individual genes for evidence of selection. Here we use the fact that the time since species divergence allows combination of data across genes. In a comparison of amino-acid replacements among species of the mustard weed Arabidopsis with those among species of the fruitfly Drosophila, we find evidence for predominantly beneficial gene substitutions in Drosophila but predominantly detrimental substitutions in Arabidopsis. We attribute this difference to the Arabidopsis mating system of partial self-fertilization, which corroborates a prediction of population genetics theory that species with a high frequency of inbreeding are less efficient in eliminating deleterious mutations owing to their reduced effective population size.     

HLA-A and B polymorphisms predate the divergence of humans and chimpanzees

Major histocompatibility complex (MHC) glycoproteins bind processed fragments of proteins and present them to the receptors of T lymphocytes. The extraordinary polymorphism of class I MHC molecules in man (HLA-A, B and C) and mouse (H-2 K, D and L) poses many questions concerning their diversification and evolution. Comparison of allelic sequences within a species suggests diversity is generated by the assortment of point mutations into varied combinations by mechanisms of recombination and gene conversion. We have now compared class I MHC alleles in two closely related species: humans (Homo sapiens) and chimpanzees (Pan troglodytes). Chimpanzee homologues of HLA-A, HLA-B and a non-classical gene have been identified. No features distinguishing human and chimpanzee alleles could be found. Individual HLA-A or B alleles are more closely related to individual chimpanzee alleles than to other HLA-A or B alleles. These results show that a considerable proportion of contemporary HLA-A and B polymorphism existed before divergence of the chimpanzee and human lines. The stability of the polymorphism indicates that hyper-mutational mechanisms are not necessary to account for HLA-A, B and C diversity.     

Prediction of sites under adaptive evolution in flavin-containing monooxygenases: Selection pattern revisited

Flavin-containing monooxygenase (FMO), like cytochrome P450 (CYP), is a monooxygenase that uses the reducing equivalents of NADPH to reduce one atom of molecular oxygen to water, while the other atom is used to oxidize the substrate. Recently, it was shown that some CYP isoforms have been subject to positive selection. However, it is unknown whether the highly conserved phase I detoxification enzyme, FMO, has undergone similar positive Darwinian selection. We used maximum-likelihood models of codon substitution, evolutionary fingerprinting, and cross species comparison to investigate the occurrence of adaptive evolution in FMO sequences. We used recent genomic data from a range of species, including vertebrates and invertebrates. We present the evidence for the occurrence of adaptive evolution in mammalian FMO 3, 4, 5, and fugu FMOs but not in mammalian FMO 1, FMO 2, frog FMOs, other fish FMOs and invertebrate FMOs. The sites under adaptive evolution were significantly associated with the insertion domain in mammalian FMO 5. We identified specific amino acid sites in FMOs 3–5 that are likely targets for selection based on the patterns of parallel amino acid change. The most likely role of adaptive evolution is the repair of mutations that permitted optimal NADP+ binding and improved catalytic efficiency. The occurrence of positive selection during the evolution of phase I detoxification enzymes such as FMOs 3–5 and fugu FMO suggests the occurrence of both high selection pressure acting on species within their unique habitats and significant changes in intensity and direction (forms of xenobiotics and drugs) resulting from changes in microhabitat and food.     

Speciation through sensory drive in cichlid fish

Theoretically, divergent selection on sensory systems can cause speciation through sensory drive. However, empirical evidence is rare and incomplete. Here we demonstrate sensory drive speciation within island populations of cichlid fish. We identify the ecological and molecular basis of divergent evolution in the cichlid visual system, demonstrate associated divergence in male colouration and female preferences, and show subsequent differentiation at neutral loci, indicating reproductive isolation. Evidence is replicated in several pairs of sympatric populations and species. Variation in the slope of the environmental gradients explains variation in the progress towards speciation: speciation occurs on all but the steepest gradients. This is the most complete demonstration so far of speciation through sensory drive without geographical isolation. Our results also provide a mechanistic explanation for the collapse of cichlid fish species diversity during the anthropogenic eutrophication of Lake Victoria.     

Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection

The major histocompatibility complex (MHC) loci are known to be highly polymorphic in humans, mice and certain other mammals, with heterozygosity as high as 80-90% (ref. 1). Four different hypotheses have been proposed to explain this high degree of polymorphism: (1) a high mutation rate, (2) gene conversion or interlocus genetic exchange, (3) over dominant selection and (4) frequency-dependent selection. In an attempt to establish which of these hypotheses is correct, we examined the pattern of nucleotide substitution between polymorphic alleles in the region of the antigen recognition site (ARS) and other regions of human and mouse class I MHC genes. The results indicate that in ARS the rate of nonsynonymous (amino acid altering) substitution is significantly higher than that of synonymous substitution in both humans and mice, whereas in other regions the reverse is true. This observation, together with a theoretical study and other considerations, supports the hypothesis of overdominant selection (heterozygote advantage).     

Rapid evolution of male reproductive genes in the descent of man

A diverse body of morphological and genetic evidence has suggested that traits pertaining to male reproduction may have evolved much more rapidly than other types of character. Recently, DNA sequence comparisons have also shown a very high level of divergence in male reproductive proteins between closely related Drosophila species, among marine invertebrates and between mouse and rat. Here we show that rapid evolution of male reproductive genes is observable in primates and is quite notable in the lineages to human and chimpanzee. Nevertheless, rapid evolution by itself is not necessarily an indication of positive darwinian selection; relaxation of negative selection is often equally compatible with the DNA sequence data. By taking three statistical approaches, we show that positive darwinian selection is often the driving force behind this rapid evolution. These results open up opportunities to test the hypothesis that sexual selection plays some role in the molecular evolution of higher primates.     

High mutation rate and predominance of insertions in the Caenorhabditis elegans nuclear genome

Mutations have pivotal functions in the onset of genetic diseases and are the fundamental substrate for evolution. However, present estimates of the spontaneous mutation rate and spectrum are derived from indirect and biased measurements. For instance, mutation rate estimates for Caenorhabditis elegans are extrapolated from observations on a few genetic loci with visible phenotypes and vary over an order of magnitude. Alternative approaches in mammals, relying on phylogenetic comparisons of pseudogene loci and fourfold degenerate codon positions, suffer from uncertainties in the actual number of generations separating the compared species and the inability to exclude biases associated with natural selection. Here we provide a direct and unbiased estimate of the nuclear mutation rate and its molecular spectrum with a set of C. elegans mutation-accumulation lines that reveal a mutation rate about tenfold higher than previous indirect estimates and an excess of insertions over deletions. Because deletions dominate patterns of C. elegans pseudogene variation, our observations indicate that natural selection might be significant in promoting small genome size, and challenge the prevalent assumption that pseudogene divergence accurately reflects the spontaneous mutation spectrum.     

Independent evolution of bitter-taste sensitivity in humans and chimpanzees

It was reported over 65 years ago that chimpanzees, like humans, vary in taste sensitivity to the bitter compound phenylthiocarbamide (PTC). This was suggested to be the result of a shared balanced polymorphism, defining the first, and now classic, example of the effects of balancing selection in great apes. In humans, variable PTC sensitivity is largely controlled by the segregation of two common alleles at the TAS2R38 locus, which encode receptor variants with different ligand affinities. Here we show that PTC taste sensitivity in chimpanzees is also controlled by two common alleles of TAS2R38; however, neither of these alleles is shared with humans. Instead, a mutation of the initiation codon results in the use of an alternative downstream start codon and production of a truncated receptor variant that fails to respond to PTC in vitro. Association testing of PTC sensitivity in a cohort of captive chimpanzees confirmed that chimpanzee TAS2R38 genotype accurately predicts taster status in vivo. Therefore, although Fisher et al.'s observations were accurate, their explanation was wrong. Humans and chimpanzees share variable taste sensitivity to bitter compounds mediated by PTC receptor variants, but the molecular basis of this variation has arisen twice, independently, in the two species.     

Detecting positive darwinian selection in brain-expressed genes during human evolution

To understand the genetic basis that underlies the phenotypic divergence between human and non- human primates, we screened a total of 7176 protein-coding genes expressed in the human brain and compared them with the chimpanzee orthologs to identify genes that show evidence of rapid evolution in the human lineage. Our results showed that the nonsynonymous/synonymous substitution (Ka/Ks) ratio for genes expressed in the brain of human and chimpanzee is 0.3854, suggesting that the brain-expressed genes are under functional constraint. The X-linked human brain-expressed genes evolved more rapidly than autosomal ones. We further dissected the molecular evolutionary patterns of 34 candidate genes by sequencing representative primate species to identify lineage-specific adaptive evolution. Fifteen out of the 34 candidate genes showed evidence of positive Darwinian selection in human and/or chimpanzee lineages. These genes are predicted to play diverse functional roles in em- bryonic development, spermatogenesis and male fertility, signal transduction, sensory nociception, and neural function. This study together with others demonstrated the usefulness and power of phy- logenetic comparison of multiple closely related species in detecting lineage-specific adaptive evolu- tion, and the identification of the positively selected brain-expressed genes may add new knowledge to the understanding of molecular mechanism of human origin.     

Gene expression divergence recapitulates the developmental hourglass model

The observation that animal morphology tends to be conserved during the embryonic phylotypic period (a period of maximal similarity between the species within each animal phylum) led to the proposition that embryogenesis diverges more extensively early and late than in the middle, known as the hourglass model. This pattern of conservation is thought to reflect a major constraint on the evolution of animal body plans. Despite a wealth of morphological data confirming that there is often remarkable divergence in the early and late embryos of species from the same phylum, it is not yet known to what extent gene expression evolution, which has a central role in the elaboration of different animal forms, underpins the morphological hourglass pattern. Here we address this question using species-specific microarrays designed from six sequenced Drosophila species separated by up to 40 million years. We quantify divergence at different times during embryogenesis, and show that expression is maximally conserved during the arthropod phylotypic period. By fitting different evolutionary models to each gene, we show that at each time point more than 80% of genes fit best to models incorporating stabilizing selection, and that for genes whose evolutionarily optimal expression level is the same across all species, selective constraint is maximized during the phylotypic period. The genes that conform most to the hourglass pattern are involved in key developmental processes. These results indicate that natural selection acts to conserve patterns of gene expression during mid-embryogenesis, and provide a genome-wide insight into the molecular basis of the hourglass pattern of developmental evolution.     

SPL家族基因复制及功能分化分析

【目的】基因复制及随后的功能分化是基因组和物种演化的重要驱动力。植物特有的转录因子家族SPL(SQUAMOSA-promoter binding protein like)广泛参与调控植物生长发育及响应逆境胁迫,为研究重复基因的起源方式和进化命运提供了良好的研究系统。本研究对葡萄(Vitis vinifera)、番木瓜(Carica papaya)、毛果杨(Populus trichocarpa)和拟南芥(Arabidopsis thaliana)4种模式植物的SPL基因家族开展基因复制及功能分化分析,为进一步研究SPL基因功能、预测种属特异性的功能基因提供系统进化角度的参考。【方法】利用SBP特征结构域,鉴定葡萄、番木瓜、毛果杨和拟南芥4种模式植物中SPL基因家族成员,并利用最大似然法构建系统进化树。基于物种内、物种间基因组共线性,分析SPL基因家族发生基因复制的方式及差异保留情况,并计算保留的SPL直系和旁系同源基因的同义、非同义替换率,分析功能分化情况。【结果】在4种模式植物中共鉴定出SPL基因73个,其中42个是miR156的靶基因。系统进化分析显示:73个SPL基因聚类为9个主要分支,miR156靶向SPL基因成簇聚集在6个主要分支;Clade I中SPL基因编码的2个锌指结构基序为C4和C₂HC,而其余8个分支中SPL基因的锌指结构基序由C3H和C₂HC组成。大规模基因组复制事件(片段复制或全基因组复制)是SPL基因家族发生基因重复的主要方式。根据基因组复制事件推算,15个古基因位点理论上应复制出的360个位点中,83.6%的重复位点发生丢失或演化成非SPL基因。本研究鉴定出旁系同源基因17对,直系同源基因27对,且所有旁系和直系同源基因的K_a/K_s(非同义替换率和同义替换率之比)值均小于1。【结论】在不同物种中保留下来的SPL直系同源基因受到较强的纯化选择,在功能上具有保守性;同一物种中保留下来的SPL旁系同源基因在进化过程中维持部分功能冗余,但在组织表达偏好性和蛋白功能上已呈现出不同形式的分化。