Rapid evolution in a pair of recent duplicate segments of rice

Gene duplication has been considered the most important way of generating genetic novelties. The subsequent evolution right after gene duplication is critical for new function to occur. Here we analyzed the evolutionary pattern for a recently duplicated segment between rice chromosomes 11 and 12. This duplication event was estimated to occur about 6 million years ago, during the divergence of the B- and C-genome rice species. The duplicate segment in chromosome 12 has significantly higher frequency of sequence rearrangement rate than non-duplicated regions. The rearrangement rate is approximately 6.5 breakages/Mb per million years, about six times higher than the fastest rate ever reported in eukaryotes. The genes within both segments experienced accelerated nucleotide substitution rates revealed by synonymous (Ks) and non-synonymous divergence (Ka) between Oryza sativa indica and O. sativa japonica. Analysis using EST data also implicates rapid divergence in expression between these segmental duplicate genes. These overall rapid changes from different perspective for the first time provide evidence that relaxation of selection also occurs in large-scale duplications.     

The IQD Gene Family in Soybean: Structure,Phylogeny, Evolution and Expression

Members of the plant-specific IQ67-domain (IQD) protein family are involved in plant development and the basal defense response. Although systematic characterization of this family has been carried out in Arabidopsis, tomato (Solanum lycopersicum), Brachypodium distachyon and rice (Oryza sativa), systematic analysis and expression profiling of this gene family in soybean (Glycine max) have not previously been reported. In this study, we identified and structurally characterized IQD genes in the soybean genome. A complete set of 67 soybean IQD genes (GmIQD1–67) was identified using Blast search tools, and the genes were clustered into four subfamilies (IQD I–IV) based on phylogeny. These soybean IQD genes are distributed unevenly across all 20 chromosomes, with 30 segmental duplication events, suggesting that segmental duplication has played a major role in the expansion of the soybean IQD gene family. Analysis of the Ka/Ks ratios showed that the duplicated genes of the GmIQD family primarily underwent purifying selection. Microsynteny was detected in most pairs: genes in clade 1–3 might be present in genome regions that were inverted, expanded or contracted after the divergence; most gene pairs in clade 4 showed high conservation with little rearrangement among these gene-residing regions. Of the soybean IQD genes examined, six were most highly expressed in young leaves, six in flowers, one in roots and two in nodules. Our qRT-PCR analysis of 24 soybean IQD III genes confirmed that these genes are regulated by MeJA stress. Our findings present a comprehensive overview of the soybean IQD gene family and provide insights into the evolution of this family. In addition, this work lays a solid foundation for further experiments aimed at determining the biological functions of soybean IQD genes in growth and development.     

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.     

Expansion Mechanisms and Functional Divergence of the Glutathione S-Transferase Family in Sorghum and Other Higher Plants

Glutathione S-transferases (GSTs) exist in various eukaryotes and function in detoxification of xenobiotics and in response to abiotic and biotic stresses. We have carried out a genome-wide survey of this gene family in 10 plant genomes. Our data show that tandem duplication has been regarded as the major expansion mechanism and both monocot and dicot plants may have practiced different expansion and evolutionary history. Non-synonymous substitutions per site (Ka) and synonymous substitutions per site (Ks) analyses showed that N- and C-terminal functional domains of GSTs (GST_N and GST_C) seem to have evolved under a strong purifying selection (Ka/Ks < 1) under different selective pressures. Differential evolutionary rates between GST_N and GST_C and high degree of expression divergence have been regarded as the major drivers for the retention of duplicated genes and the adaptability to various stresses. Expression profiling also indicated that the gene family plays a role not only in stress-related biological processes but also in the sugar-signalling pathway. Our survey provides additional annotation of the plant GST gene family and advance the understanding of plant GSTs in lineage-specific expansion and species diversification.     

Species-Specific Expansion and Molecular Evolution of the 3-hydroxy-3-methylglutaryl Coenzyme A Reductase (HMGR) Gene Family in Plants

The terpene compounds represent the largest and most diverse class of plant secondary metabolites which are important in plant growth and development. The 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR; EC 1.1.1.34) is one of the key enzymes contributed to terpene biosynthesis. To better understand the basic characteristics and evolutionary history of the HMGR gene family in plants, a genome-wide analysis of HMGR genes from 20 representative species was carried out. A total of 56 HMGR genes in the 14 land plant genomes were identified, but no genes were found in all 6 algal genomes. The gene structure and protein architecture of all plant HMGR genes were highly conserved. The phylogenetic analysis revealed that the plant HMGRs were derived from one ancestor gene and finally developed into four distinct groups, two in the monocot plants and two in dicot plants. Species-specific gene duplications, caused mainly by segmental duplication, led to the limited expansion of HMGR genes in Zea mays, Gossypium raimondii, Populus trichocarpa and Glycine max after the species diverged. The analysis of Ka/Ks ratios and expression profiles indicated that functional divergence after the gene duplications was restricted. The results suggested that the function and evolution of HMGR gene family were dramatically conserved throughout the plant kingdom.     

A statistical analysis of nucleotide substitutions in the Drosophila Adh region reflects irregularities in molecular clocks

Substitutions rates are expected to be rather constant when a gene is compared between species. To analyze this feature, Ka/Ks ratios have been studied for Alcohol dehydrogenase (Adh) and Alcohol dehydrogenase duplication (Adh-dup) genes in Drosophila species. Adh Ka/Ks values are lower in intrasubgenus comparisons involving species of the Sophophora group than when these species are compared to the D. immigrans and S. lebanonensis, and this difference does not occur in the Adh-dup comparisons.     

Reductive genome evolution from the mother of Rickettsia

The Rickettsia genus is a group of obligate intracellular α-proteobacteria representing a paradigm of reductive evolution. Here, we investigate the evolutionary processes that shaped the genomes of the genus. The reconstruction of ancestral genomes indicates that their last common ancestor contained more genes, but already possessed most traits associated with cellular parasitism. The differences in gene repertoires across modern Rickettsia are mainly the result of differential gene losses from the ancestor. We demonstrate using computer simulation that the propensity of loss was variable across genes during this process. We also analyzed the ratio of nonsynonymous to synonymous changes (Ka/Ks) calculated as an average over large sets of genes to assay the strength of selection acting on the genomes of Rickettsia, Anaplasmataceae, and free-living γ-proteobacteria. As a general trend, Ka/Ks were found to decrease with increasing divergence between genomes. The high Ka/Ks for closely related genomes are probably due to a lag in the removal of slightly deleterious nonsynonymous mutations by natural selection. Interestingly, we also observed a decrease of the rate of gene loss with increasing divergence, suggesting a similar lag in the removal of slightly deleterious pseudogene alleles. For larger divergence (Ks > 0.2), Ka/Ks converge toward similar values indicating that the levels of selection are roughly equivalent between intracellular α-proteobacteria and their free-living relatives. This contrasts with the view that obligate endocellular microorganisms tend to evolve faster as a consequence of reduced effectiveness of selection, and suggests a major role of enhanced background mutation rates on the fast protein divergence in the obligate intracellular α-proteobacteria.     

Genomic,molecular evolution,and expression analysis of NOX genes in soybean (Glycine max)

Reactive oxygen species (ROS) are versatile signaling molecules in sensing stresses and play critical roles in signaling and development. Plasma membrane NADPH oxidases (NOXs) are key producers of ROS, and play important roles in the regulation of plant-pathogen interactions. Here, we performed a comprehensive analysis of the NOX gene family in the soybean genome (Glycine max) and 17 NOX (GmNOX) genes were identified. Structural analysis revealed that the GmNOX proteins in soybean were as conserved as those in other plants. 8 duplicated gene pairs were formed by a Glycine-specific whole-genome duplication (WGD) event approximately 13 million years ago (Mya). The Ka/Ks ratios of GmNOX genes ranged from 0.04 to 0.28, suggesting that the GmNOX family had undergone purifying selection in soybean. Gene expression patterns showed different expression of these duplicate genes, suggesting that the GmNOXs were retained by substantial subfunctionalization during the soybean evolutionary processes. Subsequently, the expression of GmNOXs in response to drought and phytohormones were characterized via qPCR. Importantly, four GmNOXs showed strong expression in nodules, pointing to their probable involvement in nodulation. Thus, our results shed light on the evolutionary history of this family in soybean and contribute to the functional characterization of GmNOX genes in soybean.     

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.     

Characterization of duplicated Dunaliella viridis SPT1 genes provides insights into early gene divergence after duplication

The sodium-dependent phosphate transporter gene from unicellular green algae Dunaliella viridis, DvSPT1, shares similarity with members of Pi transporter family. Sequencing analysis of D. viridis BAC clone containing the DvSPT1 gene revealed two inverted duplicated copies of this gene (DvSPT1 and DvSPT1-2 respectively). The duplication covered most of both genes except for their 3' downstream region. The duplicated genomic sequences exhibited 97.9% identity with a synonymous divergence of Ks=0.0126 in the coding region. This data indicated very recent gene duplication in D. viridis genome, providing an excellent opportunity to investigate sequence and expression divergence of duplicated genes at an early stage. Scatted point mutations and length polymorphism of simple sequence repeats (SSRs) were predominant among the sequence divergence soon after gene duplication. Due to sequence divergence in the 5' regulatory regions and a swap of the entire 3' downstream regions (3'-UTR), DvSPT1 and DvSPT1-2 showed expression divergence in response to extra-cellular NaCl concentration changes. According to their expression patterns, the two diverged gene copies would provide better adaptation to a broader range of extra-cellular NaCl concentration. Furthermore, Southern blot analysis indicated that there might be a large phosphate transporter gene family in D. viridis.     

重复基因的进化一回顾与进展

基因重复是普遍存在的生物学现象, 是基因组和遗传系统多样化的重要推动力量, 在生物进化过程中发挥着极其重要的作用。基因重复有何利弊, 基因发生重复后, 2个重复子拷贝的保留在基因功能方面是否存在偏好性, 子拷贝在表达和进化速率上如何分化, 以及重复基因为什么会被保留下来一直是进化生物学领域研究的热点问题之一。该文对以上重复基因研究的热点问题进行了介绍, 并对重复基因的进化机制和理论模型及其近年来的一些主要研究进展进行了综述。     

重复基因的进化－－回顾与进展

基因重复是普遍存在的生物学现象, 是基因组和遗传系统多样化的重要推动力量, 在生物进化过程中发挥着极其重要的作用。基因重复有何利弊, 基因发生重复后, 2个重复子拷贝的保留在基因功能方面是否存在偏好性, 子拷贝在表达和进化速率上如何分化, 以及重复基因为什么会被保留下来一直是进化生物学领域研究的热点问题之一。该文对以上重复基因研究的热点问题进行了介绍, 并对重复基因的进化机制和理论模型及其近年来的一些主要研究进展进行了综述。     

Asymmetric evolution of duplicate genes encoding the CCAAT-binding factor NF-Y in plant genomes

NF-Y is a ubiquitous CCAAT-binding factor composed of NF-YA, NF-YB and NF-YC. Multiple genes encoding NF-Y subunits have been identified in plant genomes. It remains unclear whether the duplicate genes underwent different evolutionary patterns. Likelihood-ratio tests were used to examine whether the amino acid substitution rates are the same between duplicate genes. The influences of selection on evolution were evaluated by comparing the conservative and radical amino acid substitution rates, as well as maximum-likelihood analysis. Some NF-YB and NF-YC duplicates showed significant evidence of asymmetric evolution but not the NF-YA duplicates. Most amino acid replacements in the NF-YB and NF-YC duplicates result in changes in hydropathy, polar requirement and polarity. The physicochemical changes in the sequences of NF-YB seem to be coupled to asymmetric divergence in gene function. Plant NF-Y genes have evolved in different patterns. Relaxed selective constraints following gene duplication are most likely responsible for the unequal evolutionary rates and distinct divergence patterns of duplicate NF-Y genes. Positive selection may have promoted amino acid hydropathy changes in the NF-YC duplicates.     

Positive selection signatures in the TLR7 family

While adaptive immunity genes evolve rapidly under the influence of positive selection, innate immune system genes are known to evolve slowly due to strong purifying selection. Among the sensors of the innate immune system, Toll-like receptors (TLRs) are particularly important due to their ability to recognize and respond to pathogen-associated molecular patterns (PAMP), such as lipopolysaccharides, peptidoglycans, and nucleic acids from bacteria or viruses. In the present study, we examine the evolutionary process that has operated on the TLR7 family genes TLR7, TLR8, and TLR9. The results demonstrate that the average Ka/Ks (the ratio between nonsynonymous and synonymous substitution rates) of each TLR family gene is far lower than one regardless of estimating methods, supporting previous observations of strong purifying selection in this gene family. Interestingly, however, analysis of Ka/Ks ratios along the coding regions of TLR7 family genes by sliding-window analysis reveals a few narrow high peaks (Ka/Ks > 1). The most prominent peak corresponds to a specific region in the ectodomain, which exists only in the TLR7 family, suggesting that this unique structure of the TLR7 family might have been a target of positive selection in a variety of lineages. Furthermore, maximum likelihood model tests suggest that positive selection is the best explanation for a certain fraction of the amino acid substitutions in the TLR9.     

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.