Differential contributions to the transcriptome of duplicated genes in response to abiotic stresses in natural and synthetic polyploids

Polyploidy has occurred throughout plant evolution and can result in considerable changes to gene expression when it takes place and over evolutionary time. Little is known about the effects of abiotic stress conditions on duplicate gene expression patterns in polyploid plants. We examined the expression patterns of 60 duplicated genes in leaves, roots and cotyledons of allotetraploid Gossypium hirsutum in response to five abiotic stress treatments (heat, cold, drought, high salt and water submersion) using single-strand conformation polymorphism assays, and 20 genes in a synthetic allotetraploid. Over 70% of the genes showed stress-induced changes in the relative expression levels of the duplicates under one or more stress treatments with frequent variability among treatments. Twelve pairs showed opposite changes in expression levels in response to different abiotic stress treatments. Stress-induced expression changes occurred in the synthetic allopolyploid, but there was little correspondence in patterns between the natural and synthetic polyploids. Our results indicate that abiotic stress conditions can have considerable effects on duplicate gene expression in a polyploid, with the effects varying by gene, stress and organ type. Differential expression in response to environmental stresses may be a factor in the preservation of some duplicated genes in polyploids.     

Divergence in expression between duplicated genes in Arabidopsis

New genes may arise through tandem duplication, dispersed small-scale duplication, and polyploidy, and patterns of divergence between duplicated genes may vary among these classes. We have examined patterns of gene expression and coding sequence divergence between duplicated genes in Arabidopsis thaliana. Due to the simultaneous origin of polyploidy-derived gene pairs, we can compare covariation in the rates of expression divergence and sequence divergence within this group. Among tandem and dispersed duplicates, much of the divergence in expression profile appears to occur at or shortly after duplication. Contrary to findings from other eukaryotic systems, there is little relationship between expression divergence and synonymous substitutions, whereas there is a strong positive relationship between expression divergence and nonsynonymous substitutions. Because this pattern is pronounced among the polyploidy-derived pairs, we infer that the strength of purifying selection acting on protein sequence and expression pattern is correlated. The polyploidy-derived pairs are somewhat atypical in that they have broader expression patterns and are expressed at higher levels, suggesting differences among polyploidy- and nonpolyploidy-derived duplicates in the types of genes that revert to single copy. Finally, within many of the duplicated pairs, 1 gene is expressed at a higher level across all assayed conditions, which suggests that the subfunctionalization model for duplicate gene preservation provides, at best, only a partial explanation for the patterns of expression divergence between duplicated genes.     

Homeologous gene expression in response to growing temperature in a recent Allopolyploid (Coffea arabica L.)

Allopolyploidy is considered as a major factor contributing to speciation, diversification, and plant ecological adaptation. In particular, the expression of duplicate genes (homeologs) can be altered leading to functional plasticity and to phenotypic novelty. This study investigated the influence of growing temperatures on homeologous gene expression in Coffea arabica L., a recent allopolyploid involving 2 closely related diploid parental species. The relative expression of homeologs of 13 genes all located in the same genomic region was analyzed using an SNP ratio quantification method based on dideoxy-terminated sequences of cDNA amplicons. The relative expression of homeologous genes varied depending on the gene, the organ, and the growing condition. Nevertheless, expression of both homeologs was always detected (i.e., no silencing). Although the growing conditions were suitable for one or other of the parental species, neither subgenome appeared preferentially expressed. Furthermore, relative homeologous expression showed moderate variations across organs and conditions and appeared uncorrelated between adjacent genes. These results indicate the absence of signs of subfunctionalization suggesting C. arabica has not undergone noticeable diploidization. Furthermore, these results suggest that the expression of homeologous genes in C. arabica is regulated by a shared trans-regulation mechanism acting similarly on the 2 subgenomes and that the observed biases in the relative homeolog expression may result from cis fine-scale factors.     

Polyploidy and genome evolution in plants

Genome doubling (polyploidy) has been and continues to be a pervasive force in plant evolution. Modern plant genomes harbor evidence of multiple rounds of past polyploidization events, often followed by massive silencing and elimination of duplicated genes. Recent studies have refined our inferences of the number and timing of polyploidy events and the impact of these events on genome structure. Many polyploids experience extensive and rapid genomic alterations, some arising with the onset of polyploidy. Survivorship of duplicated genes are differential across gene classes, with some duplicate genes more prone to retention than others. Recent theory is now supported by evidence showing that genes that are retained in duplicate typically diversify in function or undergo subfunctionalization. Polyploidy has extensive effects on gene expression, with gene silencing accompanying polyploid formation and continuing over evolutionary time.     

Differential expression of duplicated peroxidase genes in the allotetraploid Brassica napus

Gene redundancy due to polyploidization provides a selective advantage for plant adaptation. We examined the expression patterns of two peroxidase genes (BnPOX1 and BnPOX2) in the natural allotetraploid Brassica napus and the model diploid progenitors Brassica rapa (Br) and Brassica oleracea (Bo) in response to the fungal pathogen Sclerotinia sclerotiorum. We demonstrated that the Bo homeolog of BnPOX1 was up-regulated after infection, while both BnPOX2 homeologs were down-regulated. A bias toward reciprocal expression of the homeologs of BnPOX1 in different organs in the natural allotetraploid of B. napus was also observed. These results suggest that subfunctionalization of the duplicated BnPOX genes after B. napus polyploidization as well as subneofunctionalization of the homeologs in response to this specific biotic stress has occurred. Retention of expression patterns in the diploid progenitors and the natural allotetraploid in some organs indicates that the function of peroxidase genes has been conserved during evolution.     

Gene duplication within the Green Lineage: the case of TEL genes

Recent years have witnessed a breathtaking increase in the availability of genome sequence data, providing evidence of the highly duplicate nature of eukaryotic genomes. Plants are exceptional among eukaryotic organisms in that duplicate loci compose a large fraction of their genomes, partly because of the frequent occurrence of polyploidy (or whole-genome duplication) events. Tandem gene duplication and transposition have also contributed to the large number of duplicated genes in plant genomes. Evolutionary analyses allowed the dynamics of duplicate gene evolution to be studied and several models were proposed. It seems that, over time, many duplicated genes were lost and some of those that were retained gained new functions and/or expression patterns (neofunctionalization) or subdivided their functions and/or expression patterns between them (subfunctionalization). Recent studies have provided examples of genes that originated by duplication with successive diversification within plants. In this review, we focused on the TEL (TERMINAL EAR1-like) genes to illustrate such mechanisms. Emerged from the mei2 gene family, these TEL genes are likely to be land plant-specific. Phylogenetic analyses revealed one or two TEL copies per diploid genome. TEL gene degeneration and loss in several Angiosperm species such as in poplar and maize seem to have occurred. In Arabidopsis thaliana, whose genome experienced at least three polyploidy events followed by massive gene loss and genomic reorganization, two TEL genes were retained and two new shorter TEL-like (MCT) genes emerged. Molecular and expression analyses suggest for these genes sub- and neofunctionalization events, but confirmation will come from their functional characterization.     

Expression pattern shifts following duplication indicative of subfunctionalization and neofunctionalization in regulatory genes of Arabidopsis

Gene duplication plays an important role in the evolution of diversity and novel function and is especially prevalent in the nuclear genomes of flowering plants. Duplicate genes may be maintained through subfunctionalization and neofunctionalization at the level of expression or coding sequence. In order to test the hypothesis that duplicated regulatory genes will be differentially expressed in a specific manner indicative of regulatory subfunctionalization and/or neofunctionalization, we examined expression pattern shifts in duplicated regulatory genes in Arabidopsis. A two-way analysis of variance was performed on expression data for 280 phylogenetically identified paralogous pairs. Expression data were extracted from global expression profiles for wild-type root, stem, leaf, developing inflorescence, nearly mature flower buds, and seedpod. Gene, organ, and gene by organ interaction (G x O) effects were examined. Results indicate that 85% of the paralogous pairs exhibited a significant G x O effect indicative of regulatory subfunctionalization and/or neofunctionalization. A significant G x O effect was associated with complementary expression patterns in 45% of pairwise comparisons. No association was detected between a G x O effect and a relaxed evolutionary constraint as detected by the ratio of nonsynonymous to synonymous substitutions. Ancestral gene expression patterns inferred across a Type II MADS-box gene phylogeny suggest several cases of regulatory neofunctionalization and organ-specific nonfunctionalization. Complete linkage clustering of gene expression levels across organs suggests that regulatory modules for each organ are independent or ancestral genes had limited expression. We propose a new classification, regulatory hypofunctionalization, for an overall decrease in expression level in one member of a paralogous pair while still having a significant G x O effect. We conclude that expression divergence specifically indicative of subfunctionalization and/or neofunctionalization contributes to the maintenance of most if not all duplicated regulatory genes in Arabidopsis and hypothesize that this results in increasing expression diversity or specificity of regulatory genes after each round of duplication.     

Gene complexity and gene duplicability

Eukaryotic genes are on average more complex than prokaryotic genes in terms of expression regulation, protein length, and protein-domain structure [1-5]. Eukaryotes are also known to have a higher rate of gene duplication than prokaryotes do [6, 7]. Because gene duplication is the primary source of new genes [], the average gene complexity in a genome may have been increased by gene duplication if complex genes are preferentially duplicated. Here, we test this "gene complexity and gene duplicability" hypothesis with yeast genomic data. We show that, on average, duplicate genes from either whole-genome or individual-gene duplication have longer protein sequences, more functional domains, and more cis-regulatory motifs than singleton genes. This phenomenon is not a by-product of previously known mechanisms, such as protein function [10-13], evolutionary rate [14, 15], dosage [11], and dosage balance [16], that influence gene duplicability. Rather, it appears to have resulted from the sub-neo-functionalization process in duplicate-gene evolution [11]. Under this process, complex genes are more likely to be retained after duplication because they are prone to subfunctionalization, and gene complexity is regained via subsequent neofunctionalization. Thus, gene duplication increases both gene number and gene complexity, two important factors in the origin of genomic and organismal complexity.     

Novel patterns of gene expression in polyploid plants

Genome doubling, or polyploidy, is a major factor accounting for duplicate genes found in most eukaryotic genomes. Polyploidy has considerable effects on duplicate gene expression, including silencing and up- or downregulation of one of the duplicated genes. These changes can arise with the onset of polyploidization or within several generations after polyploid formation and they can have epigenetic causal factors. Many expression alterations are organ-specific. Specific genes can be independently and repeatedly silenced during polyploidization, whereas patterns for other genes appear to be more stochastic. Three recent reports have provided intriguing new insights into the patterns, timing and mechanisms of gene expression changes that accompany polyploidy in plants.     

番茄AT-hook基因家族的鉴定及胁迫条件下的表达分析

AT-hook蛋白家族在植物生长发育、器官构建及逆境胁迫和激素信号应答中发挥重要作用。本研究在番茄基因组范围内,利用生物信息学方法对番茄AT-hook基因家族的成员、分布、结构和功能进行分析。结果表明,番茄AT-hook家族包含32个成员,分为3种类型,其中类型Ⅰ含有13个成员;遗传进化分析表明番茄AT-hook基因成员与拟南芥家族基因具有相似分类。利用实时荧光定量PCR对番茄32个基因开展组织表达分析,结果表明AT-hook基因具有表达差异,主要在根和花中表达较高。氧化胁迫分析结果表明,32个基因受ABA、SA、盐、高温和低温诱导表达,其中部分基因显著上调或下调表达,很可能参与了番茄逆境胁迫条件下的防御应答反应。本研究结果将为番茄AT-hook家族基因的深入研究提供依据,为进一步解析番茄AT-hook基因的功能奠定基础。     

Analysing gene function after duplication

After gene duplication, mutations cause the gene copies to diverge. The classical model predicts that these mutations will generally lead to the loss of function of one gene copy; rarely, new functions will be created and both duplicate genes are conserved. In contrast, under the subfunctionalization model both duplicates are preserved due to the partition of different functions between the duplicates. A recent study provides support for the subfunctionalization model, identifying several expressed gene duplicates common to humans and mice that contain regions conserved in one duplicate but variable in the other (and vice versa). We discuss both the methodology used in this study and also how gene phylogeny may lead to additional evidence for the importance of subfunctionalization in the evolution of new genes.     

Comprehensive expression analysis of rice phospholipase D gene family during abiotic stresses and development

Phospholipase D is one of the crucial enzymes involved in lipid mediated signaling, triggered during various developmental and physiological processes. Different members of PLD gene family have been known to be induced under different abiotic stresses and during developmental processes in various plant species. In this report, we are presenting a detailed microarray based expression analysis and expression profiles of entire set of PLD genes in rice genome, under three abiotic stresses (salt, cold and drought) and different developmental stages (3-vegetative stages and 11-reproductive stages). Seven and nine PLD genes were identified, which were expressed differentially under abiotic stresses and during reproductive developmental stages, respectively. PLD genes, which were expressed significantly under abiotic stresses exhibited an overlapping expression pattern and were also differentially expressed during developmental stages. Moreover, expression pattern for a set of stress induced genes was validated by real time PCR and it supported the microarray expression data. These findings emphasize the role of PLDs in abiotic stress signaling and development in rice. In addition, expression profiling for duplicated PLD genes revealed a functional divergence between the duplicated genes and signify the role of gene duplication in the evolution of this gene family in rice. This expressional study will provide an important platform in future for the functional characterization of PLDs in crop plants.     

Evolution of Stress-Regulated Gene Expression in Duplicate Genes of Arabidopsis thaliana

Due to the selection pressure imposed by highly variable environmental conditions, stress sensing and regulatory response mechanisms in plants are expected to evolve rapidly. One potential source of innovation in plant stress response mechanisms is gene duplication. In this study, we examined the evolution of stress-regulated gene expression among duplicated genes in the model plant Arabidopsis thaliana. Key to this analysis was reconstructing the putative ancestral stress regulation pattern. By comparing the expression patterns of duplicated genes with the patterns of their ancestors, duplicated genes likely lost and gained stress responses at a rapid rate initially, but the rate is close to zero when the synonymous substitution rate (a proxy for time) is >~0.8. When considering duplicated gene pairs, we found that partitioning of putative ancestral stress responses occurred more frequently compared to cases of parallel retention and loss. Furthermore, the pattern of stress response partitioning was extremely asymmetric. An analysis of putative cis-acting DNA regulatory elements in the promoters of the duplicated stress-regulated genes indicated that the asymmetric partitioning of ancestral stress responses are likely due, at least in part, to differential loss of DNA regulatory elements; the duplicated genes losing most of their stress responses were those that had lost more of the putative cis-acting elements. Finally, duplicate genes that lost most or all of the ancestral responses are more likely to have gained responses to other stresses. Therefore, the retention of duplicates that inherit few or no functions seems to be coupled to neofunctionalization. Taken together, our findings provide new insight into the patterns of evolutionary changes in gene stress responses after duplication and lay the foundation for testing the adaptive significance of stress regulatory changes under highly variable biotic and abiotic environments.     

Molecular evolution of two duplicated CDPK genes CPK7 and CPK12 in grass species: a case study in wheat (Triticum aestivum L.)

Gene duplication contributes to the expansion of gene families and subsequent functional diversification. Calcium-dependent protein kinases (CDPKs) are members of an important calcium sensor family involved in abiotic and biotic stress signaling in plants. We report here the molecular evolution and expression analysis of a pair of duplicated CDPK genes CPK7 and CPK12 that arose in the common ancestor of grass species. With higher nonsynonymous/synonymous ratios (dN/dS, or ω), CPK12 genes appear to diverge more rapidly than CPK7s, suggesting relaxed selection constraints on CPK12s. Sliding window analysis revealed increased dN and ω values at N-terminal regions and the calcium-binding EF hand loops. Likelihood analyses using various models in PAML 4.0 showed purifying selection on both CPK7 and CPK12 lineages. In addition to the divergence in cis-element combinations on their promoters, functional divergence of CPK7 and CPK12 genes was also observed in wheat where TaCPK7 was found to respond to drought (PEG), salt (NaCl), cold, and hydrogen peroxide (H(2)O(2)) while TaCPK12 responded only to the treatment of ABA, a feature that may complement or expand TaCPK7-mediated stress signaling networks of wheat. The contrasting expression patterns of CPK7 and CPK12 genes under stress conditions were also observed in rice, suggesting conservative functional evolution of these genes. Since no positive selection was detected between the two lineages, the divergence of CPK7 and CPK12 genes should be ascribed to subfunctionalization, rather than neofunctionalization. Thus, our work demonstrates another case of evolutionary employment of duplicated genes via subfunctionalization for better adaptation.