Evolution of stress-regulated gene expression in duplicate genes of Arabidopsis thaliana

Unique trajectory of gene family evolution from genomic analysis of nearly all known species in an ancient yeast lineage

Gene gains and losses are a major driver of genome evolution; their precise characterization can provide insights into the origin and diversification of major lineages. Here, we examined gene family evolution of 1154 genomes from nearly all known species in the medically and technologically important yeast subphylum Saccharomycotina. We found that yeast gene family evolution differs from that of plants, animals, and filamentous ascomycetes, and is characterized by smaller overall gene numbers yet larger gene family sizes for a given gene number. Faster-evolving lineages (FELs) in yeasts experienced significantly higher rates of gene losses-commensurate with a narrowing of metabolic niche breadth-but higher speciation rates than their slower-evolving sister lineages (SELs). Gene families most often lost are those involved in mRNA splicing, carbohydrate metabolism, and cell division and are likely associated with intron loss, metabolic breadth, and non-canonical cell cycle processes. Our results highlight the significant role of gene family contractions in the evolution of yeast metabolism, genome function, and speciation, and suggest that gene family evolutionary trajectories have differed markedly across major eukaryotic lineages.

Functional innovation through new genes as a general evolutionary process

In the past decade, our understanding of how new genes originate in diverse organisms has advanced substantially, and more than a dozen molecular mechanisms for generating initial gene structures were identified, in addition to gene duplication. These new genes have been found to integrate into and modify pre-existing gene networks primarily through mutation and selection, revealing new patterns and rules with stable origination rates across various organisms. This progress has challenged the prevailing belief that new proteins evolve from pre-existing genes, as new genes may arise de novo from noncoding DNA sequences in many organisms, with high rates observed in flowering plants. New genes have important roles in phenotypic and functional evolution across diverse biological processes and structures, with detectable fitness effects of sexual conflict genes that can shape species divergence. Such knowledge of new genes can be of translational value in agriculture and medicine.

The mixed auto-/allooctoploid genome of Japanese knotweed (Reynoutria japonica) provides insights into its polyploid origin and invasiveness

Reynoutria japonica Houtt. (Polygonaceae), a traditional Chinese medicine, is one of the top 100 most destructive invasive species worldwide due to its aggressive growth and strong adaptability. Here, we report an 8.04 Gb chromosome-scale assembly of R. japonica with 88 chromosomes across eight homologous sets. Through a combined phylogenetic and genomic analysis, we demonstrate that R. japonica is a mixed auto-/allooctoploid (AAAABBBB). Subgenome A (SubA) exhibited a close phylogenetic relationship with the related species Fallopia multiflora. We also unveiled the origin and evolutionary history of octoploid R. japonica based on resequencing data from Reynoutria species with different ploidy. Comparative genomics analysis revealed the genetic basis of R. japonica's invasivity and adaptability. The auxin response factor (ARF) gene family was significantly expanded in R. japonica, and these genes were highly expressed in rhizomes. We also investigated the collaboration and differentiation of the duplicated genes resulting from auto- and allo-polyploidization at the genomic variation, gene family evolution, and gene expression levels in R. japonica. Transcriptomic analysis of stem internodes and apices at different developmental stages revealed that the octuplication and significant expansion of the SAUR19 and SAUR63 subfamilies due to tandem replication in SubA, and the high expression of these genes in stems, likely contribute to the rapid growth of R. japonica. Our study provides important clues into adaptive evolution and polyploidy dominant traits in invasive plants, and will also provide important guidance for the breeding of polyploid crops.

Overlooked and misunderstood: can glutathione conjugates be clues to understanding plant glutathione transferases?

Plant glutathione transferases (GSTs) constitute a large and diverse family of enzymes that are involved in plant stress response, metabolism and defence, yet their physiological functions remain largely elusive. Consistent with the traditional view on GSTs across organisms as detoxification enzymes, in vitro most plant GSTs catalyse glutathionylation, conjugation of the tripeptide glutathione (GSH; γ-Glu-Cys-Gly) onto reactive molecules. However, when it comes to elucidating GST functions, it remains a key challenge that the endogenous plant glutathione conjugates (GS-conjugates) that would result from such glutathionylation reactions are rarely reported. Furthermore, GSTs often display high substrate promiscuity, and their proposed substrates are prone to spontaneous chemical reactions with GSH; hence, single-gene knockouts rarely provide clear chemotypes or phenotypes. In a few cases, GS-conjugates are demonstrated to be biosynthetic intermediates that are rapidly further metabolized towards a pathway end product, explaining their low abundance and rare detection. In this review, we summarize the current knowledge of plant GST functions and how and possibly why evolution has resulted in a broad and extensive expansion of the plant GST family. Finally, we demonstrate that endogenous GS-conjugates are more prevalent in plants than assumed and suggest they are overlooked as clues towards the identification of plant GST functions. This article is part of the theme issue 'The evolution of plant metabolism'.

Decrease in purifying selection pressures on wheat homoeologous genes: tetraploidization versus hexaploidization

A series of polyploidizations in higher-order polyploids is the main event affecting gene content in a genome. Each polyploidization event can lead to massive functional divergence because of the subsequent decrease in selection pressure on duplicated genes; however, the causal relationship between multiple rounds of polyploidization and the functional divergence of duplicated genes is poorly understood. We focused on the Triticum-Aegilops complex lineage and compared selection pressure before and after tetraploidization and hexaploidization events. Although both events led to decreased selection pressure on homoeologous gene pairs (compared with diploids and tetraploids), the initial tetraploidization had a greater impact on selection pressure on homoeologous gene pairs than did subsequent hexaploidization. Consistent with this, selection pressure on expression patterns for the initial event relaxed more than those for the subsequent event. Surprisingly, the decreased selection pressure on these homoeologous genes was independent of the existence of in-paralogs within the same subgenome. Wheat homoeologous pairs had different evolutionary consequences compared with orthologs related to other mechanisms (ancient allopolyploidization, ancient autopolyploidization, and small-scale duplication). Furthermore, tetraploidization and hexaploidization also seemed to have different evolutionary consequences. This suggests that homoeologous genes retain unique functions, including functions that are unlikely to be preserved in genes generated by the other duplication mechanisms. We found that their unique functions differed between tetraploidization and hexaploidization (e.g., reproductive and chromosome segregation processes). These findings imply that the substantial number of gene pairs resulting from multiple allopolyploidization events, especially initial tetraploidization, may have been a unique source of functional divergence.

Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms

Following whole-genome duplication (WGD), duplicate gene pairs (homoeologs) can evolve varying degrees of expression divergence. However, the determinants influencing these relative expression level differences (RFPKM) between homoeologs remain elusive. In this study, we analyzed the RFPKM between homoeologs in 3 angiosperms, Nymphaea colorata, Nelumbo nucifera, and Acorus tatarinowii, all having undergone a single WGD since the origin of angiosperms. Our results show significant positive correlations in RFPKM of homoeologs among tissues within the same species, and among orthologs across these 3 species, indicating convergent expression balance/bias between homoeologous gene copies following independent WGDs. We linked RFPKM between homoeologs to gene attributes associated with dosage-balance constraints, such as protein-protein interactions, lethal-phenotype scores in Arabidopsis (Arabidopsis thaliana) orthologs, domain numbers, and expression breadth. Notably, homoeologs with lower RFPKM often had more interactions and higher lethal-phenotype scores, indicating selective pressures favoring balanced expression. Also, homoeologs with lower RFPKM were more likely to be retained after WGDs in angiosperms. Within Nelumbo, greater RFPKM between homoeologs correlated with increased cis- and trans-regulatory differentiation between species, highlighting the ongoing escalation of gene expression divergence. We further found that expression degeneration in 1 copy of homoeologs is inclined toward nonfunctionalization. Our research highlights the importance of balanced expression, shaped by dosage-balance constraints, in the evolutionary retention of homoeologs in plants.

Dynamics of ER stress-induced gene regulation in plants

Endoplasmic reticulum (ER) stress is a potentially lethal condition that is induced by the abnormal accumulation of unfolded or misfolded secretory proteins in the ER. In eukaryotes, ER stress is managed by the unfolded protein response (UPR) through a tightly regulated, yet highly dynamic, reprogramming of gene transcription. Although the core principles of the UPR are similar across eukaryotes, unique features of the plant UPR reflect the adaptability of plants to their ever-changing environments and the need to balance the demands of growth and development with the response to environmental stressors. The past decades have seen notable progress in understanding the mechanisms underlying ER stress sensing and signalling transduction pathways, implicating the UPR in the effects of physiological and induced ER stress on plant growth and crop yield. Facilitated by sequencing technologies and advances in genetic and genomic resources, recent efforts have driven the discovery of transcriptional regulators and elucidated the mechanisms that mediate the dynamic and precise gene regulation in response to ER stress at the systems level.

BnaABF3 and BnaMYB44 regulate the transcription of zeaxanthin epoxidase genes in carotenoid and abscisic acid biosynthesis

Zeaxanthin epoxidase (ZEP) is a key enzyme that catalyzes the conversion of zeaxanthin to violaxanthin in the carotenoid and abscisic acid (ABA) biosynthesis pathways. The rapeseed (Brassica napus) genome has 4 ZEP (BnaZEP) copies that are suspected to have undergone subfunctionalization, yet the 4 genes' underlying regulatory mechanisms remain unknown. Here, we genetically confirmed the functional divergence of the gene pairs BnaA09.ZEP/BnaC09.ZEP and BnaA07.ZEP/BnaC07.ZEP, which encode enzymes with tissue-specific roles in carotenoid and ABA biosynthesis in flowers and leaves, respectively. Molecular and transgenic experiments demonstrated that each BnaZEP pair is transcriptionally regulated via ABA-responsive element-binding factor 3 s (BnaABF3s) and BnaMYB44s as common and specific regulators, respectively. BnaABF3s directly bound to the promoters of all 4 BnaZEPs and activated their transcription, with overexpression of individual BnaABF3s inducing BnaZEP expression and ABA accumulation under drought stress. Conversely, loss of BnaABF3s function resulted in lower expression of several genes functioning in carotenoid and ABA metabolism and compromised drought tolerance. BnaMYB44s specifically targeted and repressed the expression of BnaA09.ZEP/BnaC09.ZEP but not BnaA07.ZEP/BnaC07.ZEP. Overexpression of BnaA07.MYB44 resulted in increased carotenoid content and an altered carotenoid profile in petals. Additionally, RNA-seq analysis indicated that BnaMYB44s functions as a repressor in phenylpropanoid and flavonoid biosynthesis. These findings provide clear evidence for the subfunctionalization of duplicated genes and contribute to our understanding of the complex regulatory network involved in carotenoid and ABA biosynthesis in B. napus.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

Dynamics of cis-regulatory sequences and transcriptional divergence of duplicated genes in soybean

Transcriptional divergence of duplicated genes after whole genome duplication (WGD) has been described in many plant lineages and is often associated with subgenome dominance, a genome-wide mechanism. However, it is unknown what underlies the transcriptional divergence of duplicated genes in polyploid species that lack subgenome dominance. Soybean is a paleotetraploid with a WGD that occurred 5 to 13 Mya. Approximately 50% of the duplicated genes retained from this WGD exhibit transcriptional divergence. We developed accessible chromatin region (ACR) datasets from leaf, flower, and seed tissues using MNase-hypersensitivity sequencing. We validated enhancer function of several ACRs associated with known genes using CRISPR/Cas9-mediated genome editing. The ACR datasets were used to examine and correlate the transcriptional patterns of 17,111 pairs of duplicated genes in different tissues. We demonstrate that ACR dynamics are correlated with divergence of both expression level and tissue specificity of individual gene pairs. Gain or loss of flanking ACRs and mutation of cis-regulatory elements (CREs) within the ACRs can change the balance of the expression level and/or tissue specificity of the duplicated genes. Analysis of DNA sequences associated with ACRs revealed that the extensive sequence rearrangement after the WGD reshaped the CRE landscape, which appears to play a key role in the transcriptional divergence of duplicated genes in soybean. This may represent a general mechanism for transcriptional divergence of duplicated genes in polyploids that lack subgenome dominance.

Review: WRKY transcription factors: Understanding the functional divergence

WRKY transcription factors (TFs) play crucial roles in the growth and development of plants and their response to environmental changes. WRKY TFs have been detected in sequenced plant genomes. The functions and regulatory networks of many WRKY TFs, especially from Arabidopsis thaliana (AtWRKY TFs), have been revealed, and the origin of WRKY TFs in plants is clear. Nonetheless, the relationship between WRKY TFs function and classification is unclear. Furthermore, the functional divergence of homologous WRKY TFs in plants is unclear. In this review, WRKY TFs were explored based on WRKY-related literature published from 1994 to 2022. WRKY TFs were identified in 234 species at the genome and transcriptome levels. The biological functions of ∼ 71 % of AtWRKY TFs were uncovered. Although functional divergence occurred in homologous WRKY TFs, different WRKY TF groups had no preferential function.

Identification of CDK gene family and functional analysis of CqCDK15 under drought and salt stress in quinoa

as one of the oldest cultivated crops in the world, quinoa has been widely valued for its rich nutritional value and green health. In this study, 22 CDK genes (CqCDK01-CqCDK22) were identified from quinoa genome using bioinformatics method. The number of amino acids was 173-811, the molecular weight was 19,554.89 Da-91,375.70 Da, and the isoelectric point was 4.57-9.77. The phylogenetic tree divided 21 CqCDK genes into six subfamilies, the gene structure showed that 12 (54.5%) CqCDK genes (CqCDK03, CqCDK04, CqCDK05, CqCDK06, CqCDK07, CqCDK11, CqCDK14, CqCDK16, CqCDK18, CqCDK19, CqCDK20 and CqCDK21) had UTR regions at 5' and 3' ends. Each CDK protein had different motifs (3-9 motifs), but the genes with the same motifs were located in the same branch. Promoter analysis revealed 41 cis-regulatory elements related to plant hormones, abiotic stresses, tissue-specific expression and photoresponse. The results of real-time fluorescence quantitative analysis showed that the expression level of some CDK genes was higher under drought and salt stress, which indicated that CDK genes could help plants to resist adverse environmental effects. Subcellular localization showed that CqCDK15 gene was localized to the nucleus and cytoplasm, and transgenic plants overexpressing CqCDK15 gene showed higher drought and salt tolerance compared to the controls. Therefore, CDK genes are closely related to quinoa stress resistance. In this study, the main functions of quinoa CDK gene family and its expression level in different tissues and organs were analyzed in detail, which provided some theoretical support for quinoa stress-resistant breeding. Meanwhile, this study has important implications for further understanding the function of the CDK gene family in quinoa and our understanding of the CDK family in vascular plant.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Unravelling structural, functional, evolutionary and genetic basis of SWEET transporters regulating abiotic stress tolerance in maize

Sugars Will Eventually be Exported Transporters (SWEETs) are the novel sugar transporters widely distributed among living systems. SWEETs play a crucial role in various bio-physiological processes, viz., plant developmental, nectar secretion, pollen development, and regulation of biotic and abiotic stresses, in addition to their prime sugar-transporting activity. Thus, in-depth structural, evolutionary, and functional characterization of maize SWEET transporters was performed for their utility in maize improvement. The mining of SWEET genes in the latest maize genome release (v.5) showed an uneven distribution of 20 ZmSWEETs. The comprehensive structural analyses and docking of ZmSWEETs with four sugars, viz., fructose, galactose, glucose, and sucrose, revealed frequent amino acid residues forming hydrogen (asparagine, valine, serine) and hydrophobic (tryptophan, glycine, and phenylalanine) interactions. Evolutionary analyses of SWEETs showed a mixed lineage with 50-100 % commonality of ortho-groups and -sequences evolved under strong purifying selection (Ka/Ks < 0.5). The duplication analysis showed non-functionalization (ZmSWEET18 in B73) and neo- and sub-functionalization (ZmSWEET3, ZmSWEET6, ZmSWEET9, ZmSWEET19, and ZmSWEET20) events in maize. Functional analyses of ZmSWEET genes through co-expression, in silico expression and qRT-PCR assays showed the relevance of ZmSWEETs expression in regulating drought, heat, and waterlogging stress tolerances in maize. The first ever ZmSWEET-regulatory network revealed 286 direct (ZmSWEET-TF: 140 ZmSWEET-miRNA: 146) and 1226 indirect (TF-TF: 597; TF-miRNA: 629) edges. The present investigation has given new insights into the complex transcriptional and post-transcriptional regulation and the regulatory and functional relevance of ZmSWEETs in assigning stress tolerance in maize.

Genome-Wide Comprehensive Analysis of the Nitrogen Metabolism Toolbox Reveals Its Evolution and Abiotic Stress Responsiveness in Rice (Oryza sativa L.)

Nitrogen metabolism (NM) plays an essential role in response to abiotic stresses for plants. Enzyme activities have been extensively studied for nitrogen metabolism-associated pathways, but the knowledge of nitrogen metabolism-associated genes involved in stress response is still limited, especially for rice. In this study, we performed the genome-wide characterization of the genes putatively involved in nitrogen metabolism. A total of 1110 potential genes were obtained to be involved in nitrogen metabolism from eight species (Arabidopsis thaliana (L.) Heynh., Glycine max (L.) Merr., Brassica napus L., Triticum aestivum L., Sorghum bicolor L., Zea mays L., Oryza sativa L. and Amborella trichopoda Baill.), especially 104 genes in rice. The comparative phylogenetic analysis of the superfamily revealed the complicated divergence of different NM genes. The expression analysis among different tissues in rice indicates the NM genes showed diverse functions in the pathway of nitrogen absorption and assimilation. Distinct expression patterns of NM genes were observed in rice under drought stress, heat stress, and salt stress, indicating that the NM genes play a curial role in response to abiotic stress. Most NM genes showed a down-regulated pattern under heat stress, while complicated expression patterns were observed for different genes under salt stress and drought stress. The function of four representative NM genes (OsGS2, OsGLU, OsGDH2, and OsAMT1;1) was further validated by using qRT-PCR analysis to confirm their responses to these abiotic stresses. Based on the predicted transcription factor binding sites (TFBSs), we built a co-expression regulatory network containing transcription factors (TFs) and NM genes, of which the constructed ERF and Dof genes may act as the core genes to respond to abiotic stresses. This study provides novel sights to the interaction between nitrogen metabolism and the response to abiotic stresses.

Evolutionary analysis of the LORELEI gene family in plants reveals regulatory subfunctionalization

A signaling complex comprising members of the LORELEI (LRE)-LIKE GPI-anchored protein (LLG) and Catharanthus roseus RECEPTOR-LIKE KINASE 1-LIKE (CrRLK1L) families perceive RAPID ALKALINIZATION FACTOR (RALF) peptides and regulate growth, reproduction, immunity, and stress responses in Arabidopsis (Arabidopsis thaliana). Genes encoding these proteins are members of multigene families in most angiosperms and could generate thousands of signaling complex variants. However, the links between expansion of these gene families and the functional diversification of this critical signaling complex as well as the evolutionary factors underlying the maintenance of gene duplicates remain unknown. Here, we investigated LLG gene family evolution by sampling land plant genomes and explored the function and expression of angiosperm LLGs. We found that LLG diversity within major land plant lineages is primarily due to lineage-specific duplication events, and that these duplications occurred both early in the history of these lineages and more recently. Our complementation and expression analyses showed that expression divergence (i.e. regulatory subfunctionalization), rather than functional divergence, explains the retention of LLG paralogs. Interestingly, all but one monocot and all eudicot species examined had an LLG copy with preferential expression in male reproductive tissues, while the other duplicate copies showed highest levels of expression in female or vegetative tissues. The single LLG copy in Amborella trichopoda is expressed vastly higher in male compared to in female reproductive or vegetative tissues. We propose that expression divergence plays an important role in retention of LLG duplicates in angiosperms.

Genome-wide analysis of tandem duplicated genes and their expression under salt stress in seashore paspalum

Seashore paspalum (Paspalum vaginatum) is a halophytic, warm-season grass which is closely related to various grain crops. Gene duplication plays an important role in plant evolution, conferring significant plant adaptation at the genomic level. Here, we identified 2,542 tandem duplicated genes (TDGs) in the P. vaginatum genome and estimated the divergence time of pairs of TDGs based on synonymous substitution rates (Ks). Expression of P. vaginatum TDGs resulted in enrichment in many GO terms and KEGG pathways when compared to four other closely-related species. The GO terms included: "ion transmembrane transporter activity," "anion transmembrane transporter activity" and "cation transmembrane transport," and KEGG pathways included "ABC transport." RNA-seq analysis of TDGs showed tissue-specific expression under salt stress, and we speculated that P. vaginatum leaves became adapted to salt stress in the earlier whole-genome duplication (WGD; ~83.3 million years ago; Ma), whereas the entire P. vaginatum plant acquired a large number of TDGs related to salt stress in the second WGD (~23.3 Ma). These results can be used as a reference resource to accelerate salt-resistance research in other grasses and crops.

De novo genes in Arachis hypogaea cv. Tifrunner: systematic identification, molecular evolution, and potential contributions to cultivated peanut

De novo genes are derived from non-coding sequences, and they can play essential roles in organisms. Cultivated peanut (Arachis hypogaea) is a major oil and protein crop derived from a cross between Arachis duranensis and Arachis ipaensis. However, few de novo genes have been documented in Arachis. Here, we identified 381 de novo genes in A. hypogaea cv. Tifrunner based on comparison with five closely related Arachis species. There are distinct differences in gene expression patterns and gene structures between conserved and de novo genes. The identified de novo genes originated from ancestral sequence regions associated with metabolic and biosynthetic processes, and they were subsequently integrated into existing regulatory networks. De novo paralogs and homoeologs were identified in A. hypogaea cv. Tifrunner. De novo paralogs and homoeologs with conserved expression have mismatching cis-acting elements under normal growth conditions. De novo genes potentially have pluripotent functions in responses to biotic stresses as well as in growth and development based on quantitative trait locus data. This work provides a foundation for future research examining gene birth processes and gene function in Arachis and related taxa.

Duplication of a manganese/cadmium transporter gene reduces cadmium accumulation in rice grain

Global contamination of soils with toxic cadmium (Cd) is a serious health threat. Here we found that a tandem duplication of a gene encoding a manganese/Cd transporter, OsNramp5, was responsible for low-Cd accumulation in Pokkali, an old rice cultivar. This duplication doubled the expression of OsNramp5 gene but did not alter its spatial expression pattern and cellular localization. Higher expression of OsNramp5 increased uptake of Cd and Mn into the root cells but decreased Cd release to the xylem. Introgression of this allele into Koshihikari, an elite rice cultivar, through backcrossing significantly reduced Cd accumulation in the grain when cultivated in soil heavily contaminated with Cd but did not affect both grain yield and eating quality. This study not only reveals the molecular mechanism underlying low-Cd accumulation but also provides a useful target for breeding rice cultivars with low-Cd accumulation.

Genome-Wide Identification of GDSL-Type Esterase/Lipase Gene Family in Dasypyrum villosum L. Reveals That DvGELP53 Is Related to BSMV Infection

GDSL-type esterase/lipase proteins (GELPs) characterized by a conserved GDSL motif at their N-terminus belong to the lipid hydrolysis enzyme superfamily. In plants, GELPs play an important role in plant growth, development and stress response. The studies of the identification and characterization of the GELP gene family in Triticeae have not been reported. In this study, 193 DvGELPs were identified in Dasypyrum villosum and classified into 11 groups (clade A–K) by means of phylogenetic analysis. Most DvGELPs contain only one GDSL domain, only four DvGELPs contain other domains besides the GDSL domain. Gene structure analysis indicated 35.2% DvGELP genes have four introns and five exons. In the promoter regions of the identified DvGELPs, we detected 4502 putative cis-elements, which were associated with plant hormones, plant growth, environmental stress and light responsiveness. Expression profiling revealed 36, 44 and 17 DvGELPs were highly expressed in the spike, the root and the grain, respectively. Further investigation of a root-specific expressing GELP, DvGELP53, indicated it was induced by a variety of biotic and abiotic stresses. The knockdown of DvGELP53 inhibited long-distance movement of BSMV in the tissue of D. villosum. This research provides a genome-wide glimpse of the D. villosum GELP genes and hints at the participation of DvGELP53 in the interaction between virus and plants.

Differential Regulation of Maize and Sorghum Orthologs in Response to the Fungal Pathogen Exserohilum turcicum

Pathogens that infect more than one host offer an opportunity to study how resistance mechanisms have evolved across different species. Exserohilum turcicum infects both maize and sorghum and the isolates are host-specific, offering a unique system to examine both compatible and incompatible interactions. We conducted transcriptional analysis of maize and sorghum in response to maize-specific and sorghum-specific E. turcicum isolates and identified functionally related co-expressed modules. Maize had a more robust transcriptional response than sorghum. E. turcicum responsive genes were enriched in core orthologs in both crops, but only up to 16% of core orthologs showed conserved expression patterns. Most changes in gene expression for the core orthologs, including hub genes, were lineage specific, suggesting a role for regulatory divergent evolution. We identified several defense-related shared differentially expressed (DE) orthologs with conserved expression patterns between the two crops, suggesting a role for parallel evolution of those genes in both crops. Many of the differentially expressed genes (DEGs) during the incompatible interaction were related to quantitative disease resistance (QDR). This work offers insights into how different hosts with relatively recent divergence interact with a common pathogen. Our results are important for developing resistance to this critical pathogen and understanding the evolution of host-pathogen interactions.

Phylogenetic and ion-response analyses reveal a relationship between gene expansion and functional divergence in the Ca2+/cation antiporter family in Angiosperms

Plant CaCA superfamily genes with higher tendency to retain after WGD are more gene expression and function differentiated in ion-response. Plants and animals face different environmental stresses but share conserved Ca²⁺ signaling pathways, such as Ca²⁺/Cation transport. The Ca²⁺/cation antiporters superfamily (CaCAs) is an ancient and widespread family of ion-coupled cation transporters found in all kingdoms of life. We analyzed the molecular evolution progress of the family through comparative genomics and phylogenetics of CaCAs genes from plants and animals, grouping these genes into several families and clades, and identified multiple gene duplication retention events, particularly in the CAX (H⁺/cation exchanger), CCX (cation/Ca²⁺ exchanger), and NCL (Na⁺/Ca²⁺ exchanger-like) families. The tendency of duplication retention differs between families and gene clades. The gene duplication events were probably the result of whole-genome duplication (WGD) in plants and might have led to functional divergence. Tissue and ion-response expression analyses revealed that CaCAs genes with more highly differentiated expression patterns are more likely to be retained as duplicates than those with more conserved expression profiles. Phenotype of Arabidopsis thaliana mutants showed that loss of genes with a greater tendency to be retained after duplication resulted in more severe growth deficiency. CaCAs genes in salt-tolerant species tended to inherit the expression characteristics of their most recent common ancestral genes, with conservative ion-response expression. This study indicates a possible evolutionary scheme for cation transport and illustrates distinct fates and a mechanism for the evolution of gene duplicates. The increased copy numbers of genes and divergences in expression might have contributed to the divergent functions of CaCAs protein, allowing plants to cope with environmental stresses and adapt to a larger number of ecological niches.

Genes, pathways and networks responding to drought stress in oil palm roots

Oil palm is the most productive oilseed crop and its oil yield is seriously affected by frequent drought stress. However, little is known about the molecular responses of oil palm to drought stress. We studied the root transcriptomic responses of oil palm seedlings under 14-day drought stress. We identified 1293 differentially expressed genes (DEGs), involved in several molecular processes, including cell wall biogenesis and functions, phenylpropanoid biosynthesis and metabolisms, ion transport and homeostasis and cellular ketone metabolic process, as well as small molecule biosynthetic process. DEGs were significantly enriched into two categories: hormone regulation and metabolism, as well as ABC transporters. In addition, three protein–protein interaction networks: ion transport, reactive nitrogen species metabolic process and nitrate assimilation, were identified to be involved in drought stress responses. Finally, 96 differentially expressed transcription factors were detected to be associated with drought stress responses, which were classified into 28 families. These results provide not only novel insights into drought stress responses, but also valuable genomic resources to improve drought tolerance of oil palm by both genetic modification and selective breeding.

Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress

BACKGROUND

Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea).

RESULTS

In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content.

CONCLUSIONS

Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

Molecular basis for neofunctionalization of duplicated E3 ubiquitin ligases underlying adaptation to drought tolerance in Arabidopsis thaliana

Multigene families in plants expanded from ancestral genes via gene duplication mechanisms constitute a significant fraction of the coding genome. Although most duplicated genes are lost over time, many are retained in the genome. Clusters of tandemly arrayed genes are commonly found in the plant genome where they can promote expansion of gene families. In the present study, promoter fusion to the GUS reporter gene was used to examine the promoter architecture of duplicated E3 ligase genes that are part of group C in the Arabidopsis thaliana ATL family. Acquisition of gene expression by AtATL78, possibly generated from defective AtATL81 expression, is described. AtATL78 expression was purportedly enhanced by insertion of a TATA box within the core promoter region after a short tandem duplication that occurred during evolution of Brassicaceae lineages. This gene is associated with an adaptation to drought tolerance of A. thaliana. These findings also suggest duplicated genes could serve as a reservoir of tacit genetic information, and expression of these duplicated genes is activated upon acquisition of core promoter sequences. Remarkably, drought transcriptome profiling in response to rehydration suggests that ATL78-dependent gene expression predominantly affects genes with root-specific activities.

Transcriptional Rewiring, Adaptation, and the Role of Gene Duplication in the Metabolism of Ethanol of Saccharomyces cerevisiae

Ethanol is the main by-product of yeast sugar fermentation that affects microbial growth parameters, being considered a dual molecule, a nutrient and a stressor. Previous works demonstrated that the budding yeast arose after an ancient hybridization process resulted in a tier of duplicated genes within its genome, many of them with implications in this ethanol "produce-accumulate-consume" strategy. The evolutionary link between ethanol production, consumption, and tolerance versus ploidy and stability of the hybrids is an ongoing debatable issue. The implication of ancestral duplicates in this metabolic rewiring, and how these duplicates differ transcriptionally, remains unsolved. Here, we study the transcriptomic adaptive signatures to ethanol as a nonfermentative carbon source to sustain clonal yeast growth by experimental evolution, emphasizing the role of duplicated genes in the adaptive process. As expected, ethanol was able to sustain growth but at a lower rate than glucose. Our results demonstrate that in asexual populations a complete transcriptomic rewiring was produced, strikingly by downregulation of duplicated genes, mainly whole-genome duplicates, whereas small-scale duplicates exhibited significant transcriptional divergence between copies. Overall, this study contributes to the understanding of evolution after gene duplication, linking transcriptional divergence with duplicates' fate in a multigene trait as ethanol tolerance.IMPORTANCE Gene duplication events have been related with increasing biological complexity through the tree of life, but also with illnesses, including cancer. Early evolutionary theories indicated that duplicated genes could explore alternative functions due to relaxation of selective constraints in one of the copies, as the other remains as ancestral-function backup. In unicellular eukaryotes like yeasts, it has been demonstrated that the fate and persistence of duplicates depend on duplication mechanism (whole-genome or small-scale events), shaping their actual genomes. Although it has been shown that small-scale duplicates tend to innovate and whole-genome duplicates specialize in ancestral functions, the implication of duplicates' transcriptional plasticity and transcriptional divergence on environmental and metabolic responses remains largely obscure. Here, by experimental adaptive evolution, we show that Saccharomyces cerevisiae is able to respond to metabolic stress (ethanol as nonfermentative carbon source) due to the persistence of duplicated genes. These duplicates respond by transcriptional rewiring, depending on their transcriptional background. Our results shed light on the mechanisms that determine the role of duplicates, and on their evolvability.

A genome-wide survey of copy number variations reveals an asymmetric evolution of duplicated genes in rice

Background

Copy number variations (CNVs) are an important type of structural variations in the genome that usually affect gene expression levels by gene dosage effect. Understanding CNVs as part of genome evolution may provide insights into the genetic basis of important agricultural traits and contribute to the crop breeding in the future. While available methods to detect CNVs utilizing next-generation sequencing technology have helped shed light on prevalence and effects of CNVs, the complexity of crop genomes poses a major challenge and requires development of additional tools.

Results

Here, we generated genomic and transcriptomic data of 93 rice (Oryza sativa L.) accessions and developed a comprehensive pipeline to call CNVs in this large-scale dataset. We analyzed the correlation between CNVs and gene expression levels and found that approximately 13% of the identified genes showed a significant correlation between their expression levels and copy numbers. Further analysis showed that about 36% of duplicate pairs were involved in pseudogenetic events while only 5% of them showed functional differentiation. Moreover, the offspring copy mainly contributed to the expression levels and seemed more likely to become a pseudogene, whereas the parent copy tended to maintain the function of ancestral gene.

Conclusion

We provide a high-accuracy CNV dataset that will contribute to functional genomics studies and molecular breeding in rice. We also showed that gene dosage effect of CNVs in rice is not exponential or linear. Our work demonstrates that the evolution of duplicated genes is asymmetric in both expression levels and gene fates, shedding a new insight into the evolution of duplicated genes.

The Evolution of Gene Duplicates in Angiosperms and the Impact of Protein-Protein Interactions and the Mechanism of Duplication

Gene duplicates, generated through either whole genome duplication (WGD) or small-scale duplication (SSD), are prominent in angiosperms and are believed to play an important role in adaptation and in generating evolutionary novelty. Previous studies reported contrasting evolutionary and functional dynamics of duplicate genes depending on the mechanism of origin, a behavior that is hypothesized to stem from constraints to maintain the relative dosage balance between the genes concerned and their interaction context. However, the mechanisms ultimately influencing loss and retention of gene duplicates over evolutionary time are not yet fully elucidated. Here, by using a robust classification of gene duplicates in Arabidopsis thaliana, Solanum lycopersicum, and Zea mays, large RNAseq expression compendia and an extensive protein-protein interaction (PPI) network from Arabidopsis, we investigated the impact of PPIs on the differential evolutionary and functional fate of WGD and SSD duplicates. In all three species, retained WGD duplicates show stronger constraints to diverge at the sequence and expression level than SSD ones, a pattern that is also observed for shared PPI partners between Arabidopsis duplicates. PPIs are preferentially distributed among WGD duplicates and specific functional categories. Furthermore, duplicates with PPIs tend to be under stronger constraints to evolve than their counterparts without PPIs regardless of their mechanism of origin. Our results support dosage balance constraint as a specific property of genes involved in biological interactions, including physical PPIs, and suggest that additional factors may be differently influencing the evolution of genes following duplication, depending on the species, time, and mechanism of origin.

Genome-Wide Analysis of ROS Antioxidant Genes in Resurrection Species Suggest an Involvement of Distinct ROS Detoxification Systems during Desiccation

Abiotic stress is one of the major threats to plant crop yield and productivity. When plants are exposed to stress, production of reactive oxygen species (ROS) increases, which could lead to extensive cellular damage and hence crop loss. During evolution, plants have acquired antioxidant defense systems which can not only detoxify ROS but also adjust ROS levels required for proper cell signaling. Ascorbate peroxidase (APX), glutathione peroxidase (GPX), catalase (CAT) and superoxide dismutase (SOD) are crucial enzymes involved in ROS detoxification. In this study, 40 putative APX, 28 GPX, 16 CAT, and 41 SOD genes were identified from genomes of the resurrection species Boea hygrometrica, Selaginella lepidophylla, Xerophyta viscosa, and Oropetium thomaeum, and the mesophile Selaginellamoellendorffii. Phylogenetic analyses classified the APX, GPX, and SOD proteins into five clades each, and CAT proteins into three clades. Using co-expression network analysis, various regulatory modules were discovered, mainly involving glutathione, that likely work together to maintain ROS homeostasis upon desiccation stress in resurrection species. These regulatory modules also support the existence of species-specific ROS detoxification systems. The results suggest molecular pathways that regulate ROS in resurrection species and the role of APX, GPX, CAT and SOD genes in resurrection species during stress.

Genome wide identification and comparative analysis of glutathione transferases (GST) family genes in Brassica napus

Glutathione transferases (GSTs) are multifunctional enzymes that play important roles in plant development and responses to biotic and abiotic stress. However, a systematic analysis of GST family members in Brassica napus has not yet been reported. In this study, we identified 179 full-length GST genes in B. napus, 44.2% of which are clustered on various chromosomes. In addition, we identified 141 duplicated GST gene pairs in B. napus. Molecular evolutionary analysis showed that speciation and whole-genome triplication played important roles in the divergence of the B. napus GST duplicated genes. Transcriptome analysis of 21 tissues at different developmental stages showed that 47.6% of duplicated GST gene pairs have divergent expression patterns, perhaps due to structural divergence. We constructed a GST gene coexpression network with genes encoding various transcription factors (NAC, MYB, WRKY and bZIP) and identified six modules, including genes expressed during late seed development (after 40 days; BnGSTU19, BnGSTU20 and BnGSTZ1) and in the seed coat (BnGSTF6 and BnGSTF12), stamen and anther (BnGSTF8), root and stem (BnGSTU21), leaves and funiculus, as well as during the late stage of pericarp development (after 40 days; BnGSTU12 and BnGSTF2) and in the radicle during seed germination (BnGSTF14, BnGSTU1, BnGSTU28, and BnGSTZ1). These findings lay the foundation for elucidating the roles of GSTs in B. napus.

Old and young duplicate genes reveal different responses to environmental changes in Arachis duranensis

Old and young duplicate genes have been reported in some organisms. However, little is known about the properties of old and young duplicate genes in Arachis. Here, we have identified old and young duplicate genes in Arachis duranensis, and analyzed the evolution, gene complexity, gene expression pattern, and functional divergence between old and young duplicate genes. Our results showed different evolutionary, gene complexity and gene expression patterns, as well as differing correlations between old and young duplicate genes. Gene ontology results showed that old duplicate genes play a crucial role in lipid and amino acid biosynthesis and the oxidation-reduction process and that young duplicate genes are preferentially involved in photosynthesis and response to biotic stimulus. Transcriptome data sets revealed that most old and young duplicate genes had asymmetric function, and only a few duplicate genes exhibited symmetric function under drought and nematode stress. We found that old duplicate genes are preferentially involved in lipid and amino acid metabolism and response to abiotic stress, while young duplicate genes are likely to participate in photosynthesis and response to biotic stress. This work provides a better understanding of the evolution and functional divergence of old and young duplicate genes in A. duranensis.

Fine-Tuning the Expression of Duplicate Genes by Translational Regulation in Arabidopsis and Maize

Plant genomes are extensively shaped by various types of gene duplication. However, in this active area of investigation, the vast majority of studies focus on the sequence and transcription of duplicate genes, leaving open the question of how translational regulation impacts the expression and evolution of duplicate genes. We explored this issue by analyzing the ribo- and mRNA-seq data sets across six tissue types and stress conditions in Arabidopsis thaliana and maize (Zea mays). We dissected the relative contributions of transcriptional and translational regulation to the divergence in the abundance of ribosome footprint (RF) for different types of duplicate genes. We found that the divergence in RF abundance was largely programmed at the transcription level and that translational regulation plays more of a modulatory role. Intriguingly, translational regulation is characterized by its strong directionality, with the divergence in translational efficiency (TE) globally counteracting the divergence in mRNA abundance, indicating partial buffering of the transcriptional divergence between paralogs by translational regulation. Divergence in TE was associated with several sequence features. The faster-evolving copy in a duplicate pair was more likely to show lower RF abundance, which possibly results from relaxed purifying selection compared with its paralog. A considerable proportion of duplicates displayed differential TE across tissue types and stress conditions, most of which were enriched in photosynthesis, energy production, and translation-related processes. Additionally, we constructed a database TDPDG-DB (http://www.plantdupribo.tk), providing an online platform for data exploration. Overall, our study illustrates the roles of translational regulation in fine-tuning duplicate gene expression in plants.

Expression and regulatory asymmetry of retained Arabidopsis thaliana transcription factor genes derived from whole genome duplication

Background

Transcription factors (TFs) play a key role in regulating plant development and response to environmental stimuli. While most genes revert to single copy after whole genome duplication (WGD) event, transcription factors are retained at a significantly higher rate. Little is known about how TF duplicates have diverged in their expression and regulation, the answer to which may contribute to a better understanding of the elevated retention rate among TFs.

Results

Here we assessed what features may explain differences in the retention of TF duplicates and other genes using Arabidopsis thaliana as a model. We integrated 34 expression, sequence, and conservation features to build a linear model for predicting the extent of duplicate retention following WGD events among TFs and 19 groups of genes with other functions. We found that TFs was the least well predicted, demonstrating the features of TFs are substantially deviated from duplicate genes in other function groups. Consistent with this, the evolution of TF expression patterns and cis-regulatory cites favors the partitioning of ancestral states among the resulting duplicates: one “ancestral” TF duplicate retains most ancestral expression and cis-regulatory sites, while the “non-ancestral” duplicate is enriched for novel regulatory sites. By modeling the retention of ancestral expression and cis-regulatory states in duplicate pairs using a system of differential equations, we found that TF duplicate pairs in a partitioned state are preferentially maintained.

Conclusions

These TF duplicates with asymmetrically partitioned ancestral states are likely maintained because one copy retains ancestral functions while the other, at least in some cases, acquires novel cis-regulatory sites that may be important for novel, adaptive traits.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1398-z) contains supplementary material, which is available to authorized users.

Angiosperm-Wide and Family-Level Analyses of AP2/ERF Genes Reveal Differential Retention and Sequence Divergence After Whole-Genome Duplication

Plants are immobile and often face stressful environmental conditions, prompting the evolution of genes regulating environmental responses. Such evolution is achieved largely through gene duplication and subsequent divergence. One of the most important gene families involved in regulating plant environmental responses and development is the AP2/ERF superfamily; however, the evolutionary history of these genes is unclear across angiosperms and in major angiosperm families adapted to various ecological niches. Specifically, the impact on gene copy number of whole-genome duplication events occurring around the time of the origins of several plant families is unknown. Here, we present the first angiosperm-wide comparative study of AP2/ERF genes, identifying 75 Angiosperm OrthoGroups (AOGs), each derived from an ancestral angiosperm gene copy. Among these AOGs, 21 retain duplicates with increased copy number in many angiosperm lineages, while the remaining 54 AOGs tend to maintain low copy number. Further analyses of multiple species in the Brassicaceae family indicated that family-specific duplicates experienced differential selective pressures in coding regions, with some paralogs showing signs of positive selection. Further, cis regulatory elements also exhibit extensive divergence between duplicates in Arabidopsis. Moreover, comparison of expression levels suggested that AP2/ERF genes with frequently retained duplicates are enriched for broad expression patterns, offering increased opportunities for functional diversification via changes in expression patterns, and providing a mechanism for repeated duplicate retention in some AOGs. Our results represent the most comprehensive evolutionary history of the AP2/ERF gene family, and support the hypothesis that AP2/ERF genes with broader expression patterns are more likely to be retained as duplicates than those with narrower expression profiles, which could lead to a higher chance of duplicate gene subfunctionalization. The greater tendency of some AOGs to retain duplicates, allowing expression and functional divergence, may facilitate the evolution of complex signaling networks in response to new environmental conditions.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Dissecting the Genomic Diversification of Late Embryogenesis Abundant (LEA) Protein Gene Families in Plants

Late embryogenesis abundant (LEA) proteins include eight multigene families that are expressed in response to water loss during seed maturation and in vegetative tissues of desiccation tolerant species. To elucidate LEA proteins evolution and diversification, we performed a comprehensive synteny and phylogenetic analyses of the eight gene families across 60 complete plant genomes. Our integrated comparative genomic approach revealed that synteny conservation and diversification contributed to LEA family expansion and functional diversification in plants. We provide examples that: 1) the genomic diversification of the Dehydrin family contributed to differential evolution of amino acid sequences, protein biochemical properties, and gene expression patterns, and led to the appearance of a novel functional motif in angiosperms; 2) ancient genomic diversification contributed to the evolution of distinct intrinsically disordered regions of LEA_1 proteins; 3) recurrent tandem-duplications contributed to the large expansion of LEA_2; and 4) dynamic synteny diversification played a role on the evolution of LEA_4 and its function on plant desiccation tolerance. Taken together, these results show that multiple evolutionary mechanisms have not only led to genomic diversification but also to structural and functional plasticity among LEA proteins which have jointly contributed to the adaptation of plants to water-limiting environments.

Comparative analysis of selection mode reveals different evolutionary rate and expression pattern in Arachis duranensis and Arachis ipaënsis duplicated genes

Our results reveal that K_s is a determining factor affecting selective pressure and different evolution and expression patterns are detected between PSGs and NSGs in wild Arachis duplicates. Selective pressure, including purifying (negative) and positive selection, can be detected in organisms. However, studies on comparative evolutionary rates, gene expression patterns and gene features between negatively selected genes (NSGs) and positively selected genes (PSGs) are lagging in paralogs of plants. Arachis duranensis and Arachis ipaënsis are ancestors of the cultivated peanut, an important oil and protein crop. Here, we carried out a series of systematic analyses, comparing NSG and PSG in paralogs, using genome sequences and transcriptome datasets in A. duranensis and A. ipaënsis. We found that synonymous substitution rate (K_s) is a determining factor affecting selective pressure in A. duranensis and A. ipaënsis duplicated genes. Lower expression level, lower gene expression breadth, higher codon bias and shorter polypeptide length were found in PSGs and not in NSGs. The correlation analyses showed that gene expression breadth was positively correlated with polypeptide length and GC content at the first codon site (GC1) in PSGs and NSGs, respectively. There was a negative correlation between expression level and polypeptide length in PSGs. In NSGs, the K_s was positively correlated with expression level, gene expression breadth, GC1, and GC content at the third codon site (GC3), but selective pressure was negatively correlated with expression level, gene expression breadth, polypeptide length, GC1, and GC3 content. The function of most duplicated gene pairs was divergent under drought and nematode stress. Taken together, our results show that different evolution and expression patterns occur between PSGs and NSGs in paralogs of two wild Arachis species.

Whole Genome Analysis of Cyclin Dependent Kinase (CDK) Gene Family in Cotton and Functional Evaluation of the Role of CDKF4 Gene in Drought and Salt Stress Tolerance in Plants

Cotton (Gossypium spp.) is the number one crop cultivated for fiber production and the cornerstone of the textile industry. Drought and salt stress are the major abiotic stresses, which can have a huge economic impact on cotton production; this has been aggravated with continued climate change, and compounded by pollution. Various survival strategies evolved by plants include the induction of various stress responsive genes, such as cyclin dependent kinases (CDKs). In this study, we performed a whole-genome identification and analysis of the CDK gene family in cotton. We identified 31, 12, and 15 CDK genes in G. hirsutum, G. arboreum, and G. raimondii respectively, and they were classified into 6 groups. CDK genes were distributed in 15, 10, and 9 linkage groups of AD, D, and A genomes, respectively. Evolutionary analysis revealed that segmental types of gene duplication were the primary force underlying CDK genes expansion. RNA sequence and RT-qPCR validation revealed that Gh_D12G2017 (CDKF4) was strongly induced by drought and salt stresses. The transient expression of Gh_D12G2017-GFP fusion protein in the protoplast showed that Gh_D12G2017 was localized in the nucleus. The transgenic Arabidopsis lines exhibited higher concentration levels of the antioxidant enzymes measured, including peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) concentrations under drought and salt stress conditions with very low levels of oxidants. Moreover, cell membrane stability (CMS), excised leaf water loss (ELWL), saturated leaf weight (SLW), and chlorophyll content measurements showed that the transgenic Arabidopsis lines were highly tolerant to either of the stress factors compared to their wild types. Moreover, the expression of the stress-related genes was also significantly up-regulated in Gh_D12G2017(CDKF4) transgenic Arabidopsis plants under drought and salt conditions. We infer that CDKF-4s and CDKG-2s might be the primary regulators of salt and drought responses in cotton.

Alterations of CorTFL1 and CorAP1 expression correlate with major evolutionary shifts of inflorescence architecture in Cornus (Cornaceae) - a proposed model for variation of closed inflorescence forms

TFL1-, AP1- and LFY-like genes are known to be key regulators of inflorescence development. However, it remains to be tested whether the evolutionary modifications of inflorescence morphology result from shifts in their expression patterns. We compared the spatiotemporal expression patterns of CorTFL1, CorAP1 and CorLFY in six closely related Cornus species that display four types of closed inflorescence morphology using quantitative real-time polymerase chain reaction (qRT-PCR) and RNA in situ hybridization. Character mapping on the phylogeny was conducted to identify evolutionary changes and to assess the correlation between changes in gene expression and inflorescence morphology. Results demonstrated variation of gene expression patterns among species and a strong correlation between CorTFL1 expression and the branch index of the inflorescence type. Evolutionary changes in CorTFL1 and CorAP1 expression co-occurred on the phylogeny with the morphological changes underpinning inflorescence divergence. The study found a clear correlation between the expression patterns of CorTFL1 and CorAP1 and the inflorescence architecture in a natural system displaying closed inflorescences. The results suggest a role for the alteration in CorTFL1 and CorAP1 expression during the evolutionary modification of inflorescences in Cornus. We propose that a TFL1-like and AP1-like gene-based model may explain variation of closed inflorescences in Cornus and other lineages.

The Tartary Buckwheat Genome Provides Insights into Rutin Biosynthesis and Abiotic Stress Tolerance

Tartary buckwheat (Fagopyrum tataricum) is an important pseudocereal crop that is strongly adapted to growth in adverse environments. Its gluten-free grain contains complete proteins with a well-balanced composition of essential amino acids and is a rich source of beneficial phytochemicals that provide significant health benefits. Here, we report a high-quality, chromosome-scale Tartary buckwheat genome sequence of 489.3 Mb that is assembled by combining whole-genome shotgun sequencing of both Illumina short reads and single-molecule real-time long reads, sequence tags of a large DNA insert fosmid library, Hi-C sequencing data, and BioNano genome maps. We annotated 33 366 high-confidence protein-coding genes based on expression evidence. Comparisons of the intra-genome with the sugar beet genome revealed an independent whole-genome duplication that occurred in the buckwheat lineage after they diverged from the common ancestor, which was not shared with rosids or asterids. The reference genome facilitated the identification of many new genes predicted to be involved in rutin biosynthesis and regulation, aluminum stress resistance, and in drought and cold stress responses. Our data suggest that Tartary buckwheat's ability to tolerate high levels of abiotic stress is attributed to the expansion of several gene families involved in signal transduction, gene regulation, and membrane transport. The availability of these genomic resources will facilitate the discovery of agronomically and nutritionally important genes and genetic improvement of Tartary buckwheat.

Co-expression network analysis of duplicate genes in maize (Zea mays L.) reveals no subgenome bias

Background

Gene duplication is prevalent in many species and can result in coding and regulatory divergence. Gene duplications can be classified as whole genome duplication (WGD), tandem and inserted (non-syntenic). In maize, WGD resulted in the subgenomes maize1 and maize2, of which maize1 is considered the dominant subgenome. However, the landscape of co-expression network divergence of duplicate genes in maize is still largely uncharacterized.

Results

To address the consequence of gene duplication on co-expression network divergence, we developed a gene co-expression network from RNA-seq data derived from 64 different tissues/stages of the maize reference inbred-B73. WGD, tandem and inserted gene duplications exhibited distinct regulatory divergence. Inserted duplicate genes were more likely to be singletons in the co-expression networks, while WGD duplicate genes were likely to be co-expressed with other genes. Tandem duplicate genes were enriched in the co-expression pattern where co-expressed genes were nearly identical for the duplicates in the network. Older gene duplications exhibit more extensive co-expression variation than younger duplications. Overall, non-syntenic genes primarily from inserted duplications show more co-expression divergence. Also, such enlarged co-expression divergence is significantly related to duplication age. Moreover, subgenome dominance was not observed in the co-expression networks – maize1 and maize2 exhibit similar levels of intra subgenome correlations. Intriguingly, the level of inter subgenome co-expression was similar to the level of intra subgenome correlations, and genes from specific subgenomes were not likely to be the enriched in co-expression network modules and the hub genes were not predominantly from any specific subgenomes in maize.

Conclusions

Our work provides a comprehensive analysis of maize co-expression network divergence for three different types of gene duplications and identifies potential relationships between duplication types, duplication ages and co-expression consequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3194-0) contains supplementary material, which is available to authorized users.

Evolutionary Fates and Dynamic Functionalization of Young Duplicate Genes in Arabidopsis Genomes

Gene duplication is a primary means to generate genomic novelties, playing an essential role in speciation and adaptation. Particularly in plants, a high abundance of duplicate genes has been maintained for significantly long periods of evolutionary time. To address the manner in which young duplicate genes were derived primarily from small-scale gene duplication and preserved in plant genomes and to determine the underlying driving mechanisms, we generated transcriptomes to produce the expression profiles of five tissues in Arabidopsis thaliana and the closely related species Arabidopsis lyrata and Capsella rubella Based on the quantitative analysis metrics, we investigated the evolutionary processes of young duplicate genes in Arabidopsis. We determined that conservation, neofunctionalization, and specialization are three main evolutionary processes for Arabidopsis young duplicate genes. We explicitly demonstrated the dynamic functionalization of duplicate genes along the evolutionary time scale. Upon origination, duplicates tend to maintain their ancestral functions; but as they survive longer, they might be likely to develop distinct and novel functions. The temporal evolutionary processes and functionalization of plant duplicate genes are associated with their ancestral functions, dynamic DNA methylation levels, and histone modification abundances. Furthermore, duplicate genes tend to be initially expressed in pollen and then to gain more interaction partners over time. Altogether, our study provides novel insights into the dynamic retention processes of young duplicate genes in plant genomes.

Evolutionary Fates and Dynamic Functionalization of Young Duplicate Genes in Arabidopsis Genomes

Gene duplication is a primary means to generate genomic novelties, playing an essential role in speciation and adaptation. Particularly in plants, a high abundance of duplicate genes has been maintained for significantly long periods of evolutionary time. To address the manner in which young duplicate genes were derived primarily from small-scale gene duplication and preserved in plant genomes and to determine the underlying driving mechanisms, we generated transcriptomes to produce the expression profiles of five tissues in Arabidopsis thaliana and the closely related species Arabidopsis lyrata and Capsella rubella Based on the quantitative analysis metrics, we investigated the evolutionary processes of young duplicate genes in Arabidopsis. We determined that conservation, neofunctionalization, and specialization are three main evolutionary processes for Arabidopsis young duplicate genes. We explicitly demonstrated the dynamic functionalization of duplicate genes along the evolutionary time scale. Upon origination, duplicates tend to maintain their ancestral functions; but as they survive longer, they might be likely to develop distinct and novel functions. The temporal evolutionary processes and functionalization of plant duplicate genes are associated with their ancestral functions, dynamic DNA methylation levels, and histone modification abundances. Furthermore, duplicate genes tend to be initially expressed in pollen and then to gain more interaction partners over time. Altogether, our study provides novel insights into the dynamic retention processes of young duplicate genes in plant genomes.

Evolution of Gene Duplication in Plants

Ancient duplication events and a high rate of retention of extant pairs of duplicate genes have contributed to an abundance of duplicate genes in plant genomes. These duplicates have contributed to the evolution of novel functions, such as the production of floral structures, induction of disease resistance, and adaptation to stress. Additionally, recent whole-genome duplications that have occurred in the lineages of several domesticated crop species, including wheat (Triticum aestivum), cotton (Gossypium hirsutum), and soybean (Glycine max), have contributed to important agronomic traits, such as grain quality, fruit shape, and flowering time. Therefore, understanding the mechanisms and impacts of gene duplication will be important to future studies of plants in general and of agronomically important crops in particular. In this review, we survey the current knowledge about gene duplication, including gene duplication mechanisms, the potential fates of duplicate genes, models explaining duplicate gene retention, the properties that distinguish duplicate from singleton genes, and the evolutionary impact of gene duplication.

Deciphering the evolutionary interplay between subgenomes following polyploidy: A paleogenomics approach in grasses

How did plant species emerge from their most recent common ancestors (MRCAs) 250 million years ago? Modern plant genomes help to address such key questions in unveiling precise species genealogies. The field of paleogenomics is undergoing a paradigm shift for investigating species evolution from the study of ancestral genomes from extinct species to deciphering the evolutionary forces (in terms of duplication, fusion, fission, deletion, and translocation) that drove present-day plant diversity (in terms of chromosome/gene number and genome size). In this review, inferred ancestral karyotype genomes are shown to be powerful tools to (1) unravel the past history of extant species by recovering the variations of ancestral genomic compartments and (2) accelerate translational research by facilitating the transfer of genomic information from model systems to species of agronomic interest.

Purifying selection acts on coding and non-coding sequences of paralogous genes in Arabidopsis thaliana

Background

Whole-genome duplications in the ancestors of many diverse species provided the genetic material for evolutionary novelty. Several models explain the retention of paralogous genes. However, how these models are reflected in the evolution of coding and non-coding sequences of paralogous genes is unknown.

Results

Here, we analyzed the coding and non-coding sequences of paralogous genes in Arabidopsis thaliana and compared these sequences with those of orthologous genes in Arabidopsis lyrata. Paralogs with lower expression than their duplicate had more nonsynonymous substitutions, were more likely to fractionate, and exhibited less similar expression patterns with their orthologs in the other species. Also, lower-expressed genes had greater tissue specificity. Orthologous conserved non-coding sequences in the promoters, introns, and 3′ untranslated regions were less abundant at lower-expressed genes compared to their higher-expressed paralogs. A gene ontology (GO) term enrichment analysis showed that paralogs with similar expression levels were enriched in GO terms related to ribosomes, whereas paralogs with different expression levels were enriched in terms associated with stress responses.

Conclusions

Loss of conserved non-coding sequences in one gene of a paralogous gene pair correlates with reduced expression levels that are more tissue specific. Together with increased mutation rates in the coding sequences, this suggests that similar forces of purifying selection act on coding and non-coding sequences. We propose that coding and non-coding sequences evolve concurrently following gene duplication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2803-2) contains supplementary material, which is available to authorized users.

Transcriptional regulatory networks in Arabidopsis thaliana during single and combined stresses

Differentially evolved responses to various stress conditions in plants are controlled by complex regulatory circuits of transcriptional activators, and repressors, such as transcription factors (TFs). To understand the general and condition-specific activities of the TFs and their regulatory relationships with the target genes (TGs), we have used a homogeneous stress gene expression dataset generated on ten natural ecotypes of the model plant Arabidopsis thaliana, during five single and six combined stress conditions. Knowledge-based profiles of binding sites for 25 stress-responsive TF families (187 TFs) were generated and tested for their enrichment in the regulatory regions of the associated TGs. Condition-dependent regulatory sub-networks have shed light on the differential utilization of the underlying network topology, by stress-specific regulators and multifunctional regulators. The multifunctional regulators maintain the core stress response processes while the transient regulators confer the specificity to certain conditions. Clustering patterns of transcription factor binding sites (TFBS) have reflected the combinatorial nature of transcriptional regulation, and suggested the putative role of the homotypic clusters of TFBS towards maintaining transcriptional robustness against cis-regulatory mutations to facilitate the preservation of stress response processes. The Gene Ontology enrichment analysis of the TGs reflected sequential regulation of stress response mechanisms in plants.

The Impact of the Branched-Chain Ketoacid Dehydrogenase Complex on Amino Acid Homeostasis in Arabidopsis

The branched-chain amino acids (BCAAs) Leu, Ile, and Val are among nine essential amino acids that must be obtained from the diet of humans and other animals, and can be nutritionally limiting in plant foods. Despite genetic evidence of its importance in regulating seed amino acid levels, the full BCAA catabolic network is not completely understood in plants, and limited information is available regarding its regulation. In this study, transcript coexpression analyses revealed positive correlations among BCAA catabolism genes in stress, development, diurnal/circadian, and light data sets. A core subset of BCAA catabolism genes, including those encoding putative branched-chain ketoacid dehydrogenase subunits, is highly expressed during the night in plants on a diel cycle and in prolonged darkness. Mutants defective in these subunits accumulate higher levels of BCAAs in mature seeds, providing genetic evidence for their function in BCAA catabolism. In addition, prolonged dark treatment caused the mutants to undergo senescence early and overaccumulate leaf BCAAs. These results extend the previous evidence that BCAAs can be catabolized and serve as respiratory substrates at multiple steps. Moreover, comparison of amino acid profiles between mature seeds and dark-treated leaves revealed differences in amino acid accumulation when BCAA catabolism is perturbed. Together, these results demonstrate the consequences of blocking BCAA catabolism during both normal growth conditions and under energy-limited conditions.

An Ancestral Role for CONSTITUTIVE TRIPLE RESPONSE1 Proteins in Both Ethylene and Abscisic Acid Signaling

Land plants have evolved adaptive regulatory mechanisms enabling the survival of environmental stresses associated with terrestrial life. Here, we focus on the evolution of the regulatory CONSTITUTIVE TRIPLE RESPONSE1 (CTR1) component of the ethylene signaling pathway that modulates stress-related changes in plant growth and development. First, we compare CTR1-like proteins from a bryophyte, Physcomitrella patens (representative of early divergent land plants), with those of more recently diverged lycophyte and angiosperm species (including Arabidopsis [Arabidopsis thaliana]) and identify a monophyletic CTR1 family. The fully sequenced P. patens genome encodes only a single member of this family (PpCTR1L). Next, we compare the functions of PpCTR1L with that of related angiosperm proteins. We show that, like angiosperm CTR1 proteins (e.g. AtCTR1 of Arabidopsis), PpCTR1L modulates downstream ethylene signaling via direct interaction with ethylene receptors. These functions, therefore, likely predate the divergence of the bryophytes from the land-plant lineage. However, we also show that PpCTR1L unexpectedly has dual functions and additionally modulates abscisic acid (ABA) signaling. In contrast, while AtCTR1 lacks detectable ABA signaling functions, Arabidopsis has during evolution acquired another homolog that is functionally distinct from AtCTR1. In conclusion, the roles of CTR1-related proteins appear to have functionally diversified during land-plant evolution, and angiosperm CTR1-related proteins appear to have lost an ancestral ABA signaling function. Our study provides new insights into how molecular events such as gene duplication and functional differentiation may have contributed to the adaptive evolution of regulatory mechanisms in plants.

Molecular Evidence for Functional Divergence and Decay of a Transcription Factor Derived from Whole-Genome Duplication in Arabidopsis thaliana

Functional divergence between duplicate transcription factors (TFs) has been linked to critical events in the evolution of land plants and can result from changes in patterns of expression, binding site divergence, and/or interactions with other proteins. Although plant TFs tend to be retained post polyploidization, many are lost within tens to hundreds of million years. Thus, it can be hypothesized that some TFs in plant genomes are in the process of becoming pseudogenes. Here, we use a pair of salt tolerance-conferring transcription factors, DWARF AND DELAYED FLOWERING1 (DDF1) and DDF2, that duplicated through paleopolyploidy 50 to 65 million years ago, as examples to illustrate potential mechanisms leading to duplicate retention and loss. We found that the expression patterns of Arabidopsis thaliana (At)DDF1 and AtDDF2 have diverged in a highly asymmetric manner, and AtDDF2 has lost most inferred ancestral stress responses. Consistent with promoter disablement, the AtDDF2 promoter has fewer predicted cis-elements and a methylated repetitive element. Through comparisons of AtDDF1, AtDDF2, and their Arabidopsis lyrata orthologs, we identified significant differences in binding affinities and binding site preference. In particular, an AtDDF2-specific substitution within the DNA-binding domain significantly reduces binding affinity. Cross-species analyses indicate that both AtDDF1 and AtDDF2 are under selective constraint, but among A. thaliana accessions, AtDDF2 has a higher level of nonsynonymous nucleotide diversity compared with AtDDF1. This may be the result of selection in different environments or may point toward the possibility of ongoing functional decay despite retention for millions of years after gene duplication.