Genetic Signature of Rice Domestication Shown by a Variety of Genes

Cultivated rice was domesticated from common wild rice. However, little is known about genetic adaptation under domestication. We investigated the nucleotide variation of both cultivated rice and its wild progenitors at 22 R-gene and 10 non–R-gene loci. A significant regression was observed between wild rice and rice cultivars in their polymorphic levels, particularly in their nonsynonymous substitutions (θ _a ). Our data also showed that a similar proportion (approximately 60%) of nucleotide variation in wild rice was retained in cultivated rice in both R-genes and non–R-genes. Interestingly, the slope always was >1 and the intercept always >0 in linear regressions when a cultivar’s polymorphism was x-axis. The slope and intercept values can provide a basis by which to estimate the founder effect and the strength of artificial direct selection. A larger founder effect than previously reported and a strong direct-selection effect were shown in rice genes. In addition, two-directional selection was commonly found in differentiated genes between indica and japonica rice subspecies. This kind of selection may explain the mosaic origins of indica and japonica rice subspecies. Furthermore, in most R-genes, no significant differentiation between cultivated and wild rice was detected. We found evidence for genetic introgression from wild rice, which may have played an important role during the domestication of rice R-genes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. Yuanli Zhang and Jiao Wang contributed equally to this work.     

Genome-wide investigation on the genetic variations of rice disease resistance genes

Exploitation of plant disease resistance (R) gene in breeding programs has been proven to be the most efficient strategy for coping with the threat of pathogens. An understanding of R-gene variation is the basis for this strategy. Here we report a genome-wide investigation on the variation of NBS-LRR-encoding genes, the common type of R genes, between two sequenced rice genomes, Oryza sativa L. var. Nipponbare and 93–11. We show that the allelic nucleotide diversity in 65.0% of 397 least-divergent pairs is not high (0.344% on average), while the remaining 35% display a greater diversity (5.4% on average). The majority of conserved R genes is single-copy and/or located as a singleton. The clustered, particularly the complex-clustered, R-genes contribute greatly to the rich genetic variation. Surprisingly only 11.2% of R-genes have remarkably high ratios of non-synonymous to synonymous rates, which is much less than the 17.4% observed between Arabidopsis genomes. Noticeable “artificially selective sweeping” could be detected in a large proportion of the conserved R-genes, a scenario described in the “arms race” co-evolutionary model. Based on our study, a variation pattern of R-genes is proposed and confirmed by the analysis of R-genes from other rice lines, indicating that the observed variation pattern may be common in all rice lines.Electronic Supplementary Material Supplementary material is available for this article at     

Genetic variation of NBS-LRR class resistance genes in rice lines

The use of plant disease resistance (R) genes in breeding programs needs an understanding of their variation patterns. In our current study, we investigated the polymorphisms of 44 NBS-LRR class R-genes among 21 rice cultivars and 14 wild rice populations. Our data suggested that there were four basic types of variations: conserved, diversified, intermediate-diversified, and present/absent patterns. Common characteristics at a locus of conserved R-genes were: copy-number uniformity, clear divergence (long branches) with other paralogs, and highly identical alleles. On the other hand, copy-number variability, a nearly equal and non-zero branch lengths, and high levels of nucleotide diversity were observed at the loci of highly diversified R-genes. Research suggests that the ratio of diverse alleles to the total number of genes at a locus is one of the best criteria to characterize the variation pattern of an R-gene. Our data suggested that a significant genetic reduction was detected only in four present/absent R-genes, compared with the variation observed in wild rice. In general, no difference was detected between wild rice and cultivars, japonica and indica rice, or between lines from different geographic regions. Our results also suggested that R-genes were under strong selection, which shaped R-gene variation patterns. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Sihai Yang and Tingting Gu contributed equally to this work.     

Patterns of DNA-Sequence Divergence Between <Emphasis Type="Italic">Drosophila miranda</Emphasis> and <Emphasis Type="Italic">D. pseudoobscura</Emphasis>

Contrary to the classical view, a large amount of non-coding DNA seems to be selectively constrained in Drosophila and other species. Here, using Drosophila miranda BAC sequences and the Drosophila pseudoobscura genome sequence, we aligned coding and non-coding sequences between D. pseudoobscura and D. miranda, and investigated their patterns of evolution. We found two patterns that have previously been observed in comparisons between Drosophila melanogaster and its relatives. First, there is a negative correlation between intron divergence and intron length, suggesting that longer non-coding sequences may contain more regulatory elements than shorter sequences. Our other main finding is a negative correlation between the rate of non-synonymous substitutions (d _N) and codon usage bias (F _op), showing that fast-evolving genes have a lower codon usage bias, consistent with strong positive selection interfering with weak selection for codon usage.     

Role of seed germination in adaptation and reproductive isolation in Arabidopsis lyrata

Seed germination is an important developmental and life history stage. Yet, the evolutionary impact of germination has mainly been studied in the context of dormancy, or for its role in reproductive isolation between species. Here, we aim to examine multiple consequences of genetic divergence on germination traits between two Arabidopsis lyrata subspecies: ssp. petraea (Eurasia) and ssp. lyrata (North America). Postdormancy germination time, a potentially adaptive trait, showed differentiation between the populations, and quantitative trait loci (QTL) mapping revealed that the trait variation is mainly controlled by two antagonistic loci. These QTL areas contain several candidate genes with known function in postdormancy germination in A. thaliana. The sequence variation of three genes was consistent with differential selection, and they also included fixed nonsynonymous substitutions with potential to account for the phenotypic differentiation. We further show that the divergence between the subspecies has led to a slight but significant reduction in hybrid germination proportions, indicating incipient reproductive isolation. Comparison of reciprocal F₁ and F₂ progenies suggests that Bateson–Dobzhansky–Muller incompatibilities likely act through uniparentally inherited factors. Examination of genomewide transmission ratio distortion further revealed that cytonuclear interactions cause substantial pregermination inviability in the hybrids. These results confirm that seed germination has adaptive potential beyond the dormancy stage and that hybrid seed inviability can be one of the first reproductive barriers to arise during divergence.     

Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.     

The complete mitochondrial genome of dhole <Emphasis Type="Italic">Cuon alpinus</Emphasis>: phylogenetic analysis and dating evolutionary divergence within canidae

The dhole (Cuon alpinus) is the only existent species in the genus Cuon (Carnivora: Canidae). In the present study, the complete mitochondrial genome of the dhole was sequenced. The total length is 16672 base pairs which is the shortest in Canidae. Sequence analysis revealed that most mitochondrial genomic functional regions were highly consistent among canid animals except the CSB domain of the control region. The difference in length among the Canidae mitochondrial genome sequences is mainly due to the number of short segments of tandem repeated in the CSB domain. Phylogenetic analysis was progressed based on the concatenated data set of 14 mitochondrial genes of 8 canid animals by using maximum parsimony (MP), maximum likelihood (ML) and Bayesian (BI) inference methods. The genera Vulpes and Nyctereutes formed a sister group and split first within Canidae, followed by that in the Cuon. The divergence in the genus Canis was the latest. The divarication of domestic dogs after that of the Canis lupus laniger is completely supported by all the three topologies. Pairwise sequence divergence data of different mitochondrial genes among canid animals were also determined. Except for the synonymous substitutions in protein-coding genes, the control region exhibits the highest sequence divergences. The synonymous rates are approximately two to six times higher than those of the non-synonymous sites except for a slightly higher rate in the non-synonymous substitution between Cuon alpinus and Vulpes vulpes. 16S rRNA genes have a slightly faster sequence divergence than 12S rRNA and tRNA genes. Based on nucleotide substitutions of tRNA genes and rRNA genes, the times since divergence between dhole and other canid animals, and between domestic dogs and three subspecies of wolves were evaluated. The result indicates that Vulpes and Nyctereutes have a close phylogenetic relationship and the divergence of Nyctereutes is a little earlier. The Tibetan wolf may be an archaic pedigree within wolf subspecies. The genetic distance between wolves and domestic dogs is less than that among different subspecies of wolves. The domestication of dogs was about 1.56–1.92 million years ago or even earlier.     

Comparative Molecular Evolution of Primary (<Emphasis Type="BoldItalic">Buchnera</Emphasis>) and Secondary Symbionts of Aphids Based on Two Protein-Coding Genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary (Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other γ-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.     

Abundance of plastid DNA insertions in nuclear genomes of rice and Arabidopsis

Pairwise comparison of whole plastid and draft nuclear genomic sequences of Arabidopsis thaliana and Oryza sativa L. ssp. indica shows that rice nuclear genomic sequences contain homologs of plastid DNA covering about 94 kb (83%) of plastid genome and including one or more full-length intact (without mutations resulting in premature stop codons) homologues of 26 known protein-coding (KPC) plastid genes. By contrast, only about 20 kb (16%) of chloroplast DNA, including a single intact plastid-derived KPC gene, is presented in the nucleus of A. thaliana. Sixteen rice plastid genes have at least one nuclear copy without any mutation or with only synonymous substitutions. Nuclear copies for other ten plastid genes contain both synonymous and non-synonymous substitutions. Multiple ESTs for 25 out of 26 KPC genes were also found, as well as putative promoters for some of them. The study of substitutions pattern shows that some of nuclear homologues of plastid genes may be functional and/or are under the pressure of the positive natural selection. The similar comparative analysis performed on rice chromosome 1 revealed 27 contigs containing plastid-derived sequences, totalling about 84 kb and covering two thirds of chloroplast DNA, with the intact nuclear copies of 26 different KPC genes. One of these contigs, AP003280, includes almost 57 kb (45%) of chloroplast genome with the intact copies of 22 KPC genes. At the same time, we observed that relative locations of homologues in plastid DNA and the nuclear genome are significantly different.     

Genetic signatures of plant resistance genes with known function within and between species

Plant disease resistance (R) genes have undergone significant evolutionary divergence to cope with rapid changes in pathogens. These highly variable evolutionary patterns may have contributed to diversity in R gene protein families or structures. Here, the evolutionary patterns of 76 identified R genes and their homologs were investigated within and between plant species. Results demonstrated that nucleotide binding sites and leucine-rich-repeat genes located in loci with complex evolutionary histories tended to evolve rapidly, have high variation in copy numbers, exhibit high levels of nucleotide variation and frequent gene conversion events, and also exhibit high non-synonymous to synonymous substitution ratios in LRR regions. However, non-NBS-LRR R genes are relatively well conserved with constrained variation and are more likely to participate in the basic defense system of hosts. In addition, both conserved and highly divergent evolutionary patterns were observed for the same R genes and were consistent with inter- and intra-specific distributions of some R genes. These results thus indicate either continuous or altered evolutionary patterns between and within species. The present investigation is the first attempt to investigate evolutionary patterns among all clearly functional R genes. The results reported here thus provide a foundation for future plant disease studies.     

Evolution of the number of LRRs in plant disease resistance genes

The largest group of plant resistance (R) genes contain the regions that encode the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains (NBS-LRR genes). To gain new resistance, amino acid substitutions and changes in number of the LRRs that recognize the presence of pathogens are considered important. In this study, we focus on the evolution of the number of LRRs and analyze the genome data of five plant species, Arabidopsis thaliana, Oryza sativa, Medicago truncatula, Lotus japonicus and Populus trichocarpa. We first categorized the NBS-LRR genes in each species into groups and subgroups based on the phylogenetic relationships of their NBS domain sequences. Then we estimated the evolutionary rate of the number of LRRs relative to the synonymous divergence in the NBS domain sequences by a maximum likelihood method assuming the single stepwise mutation model. The estimates ranged from 4.5 to 600 and differed between groups in the same species or between species. This indicated different roles played by different groups of the NBS-LRR genes within a species or the effects of various life history characteristics, such as generation time, of the species. We also tested the fit of the model to the data using the variance of number of LRRs in each subgroup. In some subgroups in some plants (16 out of 174 subgroups), the results of simulation using the estimated rates significantly deviated from the observed data. Those subgroups may have undergone different modes of selection from the other subgroups.     

Sexually dimorphic gene expression and transcriptome evolution provide mixed evidence for a fast‐Z effect in Heliconius

Sex chromosomes have different evolutionary properties compared to autosomes due to their hemizygous nature. In particular, recessive mutations are more readily exposed to selection, which can lead to faster rates of molecular evolution. Here, we report patterns of gene expression and molecular evolution for a group of butterflies. First, we improve the completeness of the Heliconius melpomene reference annotation, a neotropical butterfly with a ZW sex determination system. Then, we analyse RNA from male and female whole abdomens and sequence female ovary and gut tissue to identify sex‐ and tissue‐specific gene expression profiles in H. melpomene. Using these expression profiles, we compare (a) sequence divergence and polymorphism; (b) the strength of positive and negative selection; and (c) rates of adaptive evolution, for Z and autosomal genes between two species of Heliconius butterflies, H. melpomene and H. erato. We show that the rate of adaptive substitutions is higher for Z than autosomal genes, but contrary to expectation, it is also higher for male‐biased than female‐biased genes. Additionally, we find no significant increase in the rate of adaptive evolution or purifying selection on genes expressed in ovary tissue, a heterogametic‐specific tissue. Our results contribute to a growing body of literature from other ZW systems that also provide mixed evidence for a fast‐Z effect where hemizygosity influences the rate of adaptive substitutions.     

Increased constraints on MC4R during primate and human evolution

The melanocortin 4 receptor (MC4R) is routinely investigated for the role it plays in human obesity, as mutations in MC4R are the most common dominantly inherited form of the disease. As little is known about the evolutionary history of this locus, we investigated patterns of variation at MC4R in a worldwide sample of 1,015 humans from 51 populations, and in 8 central chimpanzees. There is a significant paucity of diversity at MC4R in humans, but not in chimpanzees. The spectrum of mutations in humans, combined with the overall low level of diversity, suggests that most (if not all) of the observed non-synonymous polymorphisms are likely to be transient deleterious mutations. The MC4R coding region was resequenced in 12 primate species and sequences from an additional 29 vertebrates were included in molecular evolutionary analyses. MC4R is highly conserved throughout vertebrate evolution, and has apparently been subject to high levels of continuous purifying selection that increased approximately threefold during primate evolution. Furthermore, the strong selection extends to codon usage bias, where most silent mutations are expected to be either quickly fixed or removed from the population, which may help explain the unusually low levels of silent polymorphisms in humans. Finally, there is a significant tendency for non-synonymous mutations that impact MC4R function to occur preferentially at sites that are identified by evolutionary analyses as being subject to very strong purifying selection. The information from this study should help inform future epidemiological investigations of MC4R. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The effects of recombination,mutation and selection on the evolution of the Rp1 resistance genes in grasses

Plant immune genes, or resistance genes, are involved in a co‐evolutionary arms race with a diverse range of pathogens. In agronomically important grasses, such R genes have been extensively studied because of their role in pathogen resistance and in the breeding of resistant cultivars. In this study, we evaluate the importance of recombination, mutation and selection on the evolution of the R gene complex Rp1 of Sorghum, Triticum, Brachypodium, Oryza and Zea. Analyses show that recombination is widespread, and we detected 73 independent instances of sequence exchange, involving on average 1567 of 4692 nucleotides analysed (33.4%). We were able to date 24 interspecific recombination events and found that four occurred postspeciation, which suggests that genetic introgression took place between different grass species. Other interspecific events seemed to have been maintained over long evolutionary time, suggesting the presence of balancing selection. Significant positive selection (i.e. a relative excess of nonsynonymous substitutions (d_N/d_S>1)) was detected in 17–95 codons (0.42–2.02%). Recombination was significantly associated with areas with high levels of polymorphism but not with an elevated d_N/d_S ratio. Finally, phylogenetic analyses show that recombination results in a general overestimation of the divergence time (mean = 14.3%) and an alteration of the gene tree topology if the tree is not calibrated. Given that the statistical power to detect recombination is determined by the level of polymorphism of the amplicon as well as the number of sequences analysed, it is likely that many studies have underestimated the importance of recombination relative to the mutation rate.     

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide‐toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H₂S) and compared two population pairs of sulphide‐adapted and ancestral fish by sequencing population pools of >200 individuals (Pool‐Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection‐mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide‐adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of J_corr = 0.0032 at the level of SNPs, divergent genome regions (J_corr = 0.0061) and genes therein (J_corr = 0.0091). At the level of metabolic pathways, the overlap was J_corr = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.     

Testing the neutral theory of molecular evolution using genomic data: a comparison of the human and bovine transcriptome

Despite growing evidence of rapid evolution in protein coding genes, the contribution of positive selection to intra- and interspecific differences in protein coding regions of the genome is unclear. We attempted to see if genes coding for secreted proteins and genes with narrow expression, specifically those preferentially expressed in the mammary gland, have diverged at a faster rate between domestic cattle (Bos taurus) and humans (Homo sapiens) than other genes and whether positive selection is responsible. Using a large data set, we identified groups of genes based on secretion and expression patterns and compared them for the rate of nonsynonymous (dN) and synonymous (dS) substitutions per site and the number of radical (Dr) and conservative (Dc) amino acid substitutions. We found evidence of rapid evolution in genes with narrow expression, especially for those expressed in the liver and mammary gland and for genes coding for secreted proteins. We compared common human polymorphism data with human-cattle divergence and found that genes with high evolutionary rates in human-cattle divergence also had a large number of common human polymorphisms. This argues against positive selection causing rapid divergence in these groups of genes. In most cases dN/dS ratios were lower in human-cattle divergence than in common human polymorphism presumably due to differences in the effectiveness of purifying selection between long-term divergence and short-term polymorphism.     

Evolutionary aspects of intraspecific polymorphism of the <Emphasis Type="Italic">Arabidopsis thaliana</Emphasis> genes coding for subunit I of the Mg-chelatase complex

Intraspecific polymorphism of CHLI 2, which codes for subunit I of the Mg-chelatase complex involved in Mg-protoporphyrin IX synthesis, was investigated in 19 ecotypes of Arabidopsis thaliana. Sequence divergence by 35 nucleotides was found. Nucleotide substitutions in 12 out of the 35 positions result in amino acid changes in the functionally significant C-terminal domain of the protein. Of the two divergent sequence haplogroups identified, haplogroup Col displayed an excess of low-frequency polymorphisms, suggesting the effect of purifying selection and a functional significance of CHLI 2. The results suggest a difference in evolutionary dynamics for CHLI 1 and CHLI 2, which determine different forms of subunit I of the Mg-chelatase complex in A. thaliana.