Evolution of Cis-Regulatory Elements and Regulatory Networks in Duplicated Genes of Arabidopsis

Evolutionary Consequences of Functional and Regulatory Divergence of HD-Zip I Transcription Factors as a Source of Diversity in Protein Interaction Networks in Plants

The HD superfamily has been studied in detail for several decades. The plant-specific HD-Zip I subfamily attracts the most attention because of its involvement in plant development and stress responses. In this review, we provide a comprehensive insight into the evolutionary events responsible for the functional redundancy and diversification of the HD-Zip I genes in regulating various biological processes. We summarized the evolutionary history of the HD-Zip family, highlighting the important role of WGDs in its expansion and divergence of retained duplicates in the genome. To determine the relationship between the evolutionary origin and functional conservation of HD-Zip I in different species, we performed a phylogenetic analysis, compared their expression profiles in different tissues and under stress and traced the role of orthologs and paralogs in regulating developmental processes. We found that HD-Zip I from different species have similar gene structures with a highly conserved HD and Zip, bind to the same DNA sequences and are involved in similar biological processes. However, they exhibit a functional diversity, which is manifested in altered expression patterns. Some of them are involved in the regulation of species-specific leaf morphology and phenotypes. Here, we discuss the role of changes in functional domains involved in DNA binding and protein interaction of HD-Zip I and in cis-regulated regions of its target genes in promoting adaptive innovations through the formation of de novo regulatory systems. Understanding the role of the HD-Zip I subfamily in organism-environment interactions remains a challenge for evolutionary developmental biology (evo-devo).

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

Enrichment of intrinsically disordered residues in ohnologs facilitates abiotic stress resilience in Brassica rapa

Arabidopsis thaliana and Brassica rapa are in the same evolutionary lineage, although the latter experienced an additional whole genome triplication event. Therefore, it would be intriguing to investigate the traits that gene duplication imposes to mediate plant stress tolerance. Here, we noticed that B. rapa abiotic stress resistance (ASR) genes which code at least one stress responsive domain have a significantly higher number of paralogs than A. thaliana. Analysing the disordered content of the ASR genes in both species, we found that intrinsically disordered residues (IDR) are specifically enriched in whole genome duplication (WGD) derived paralogs. Subsequently, domain similarity analysis between WGD pairs of both species has revealed that majority of WGD pairs in B. rapa did not share domains with each other. Furthermore, domain enrichment analysis has shown that B. rapa paralogs contain 36 distinct stress responsive enriched domains, significantly higher than A. thaliana paralogs. Next, we performed MSA to investigate the domain conservation between orthologs and ohnologs pairs, we found that 80.13% of B. rapa ohnologs acquire new domains, depicting the fact that ohnologs play a significant role in stress-related behaviours. The average IDR content of the ohnologs enriching new domains after gene duplication in B. rapa (0.19), is also significantly higher than A. thaliana (0.04). Interestingly, we also found that all of these attributes i.e., exhibiting higher number of WGD paralogs and enhancement of IDR in ASR genes of B. rapa compared to A. thaliana is exclusive for ASR genes only. No such significant differences were observed in randomly selected non-ASR genes between the two species. Together these results provide strong support for the hypothesis that augmentation of IDR content followed by a whole genome duplication event imposes the stress resistance potentiality in B. rapa. This research will shed light on the mechanism of how B. rapa is able to successfully adapt to stress over the evolutionary timescale.

Genome-wide characterization of the xyloglucan endotransglucosylase/hydrolase gene family in Solanum lycopersicum L. and gene expression analysis in response to arbuscular mycorrhizal symbiosis

Xyloglucan endotransglucosylase/hydrolases (XTHs) are a glycoside hydrolase protein family involved in the biosynthesis of xyloglucans, with essential roles in the regulation of plant cell wall extensibility. By taking advantage of the whole genome sequence in Solanum lycopersicum, 37 SlXTHs were identified in the present work. SlXTHs were classified into four subfamilies (ancestral, I/II, III-A, III-B) when aligned to XTHs of other plant species. Gene structure and conserved motifs showed similar compositions in each subfamily. Segmental duplication was the primary mechanism accounting for the expansion of SlXTH genes. In silico expression analysis showed that SlXTH genes exhibited differential expression in several tissues. GO analysis and 3D protein structure indicated that all 37 SlXTHs participate in cell wall biogenesis and xyloglucan metabolism. Promoter analysis revealed that some SlXTHs have MeJA- and stress-responsive elements. qRT-PCR expression analysis of nine SlXTHs in leaves and roots of mycorrhizal colonized vs. non-colonized plants showed that eight of these genes were differentially expressed in leaves and four in roots, suggesting that SlXTHs might play roles in plant defense induced by arbuscular mycorrhiza. Our results provide valuable insight into the function of XTHs in S. lycopersicum, in addition to the response of plants to mycorrhizal colonization.

Genome-wide identification and comparative analysis of CLE family in rapeseed and its diploid progenitors

Crop genomics and breeding CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) proteins belong to a small peptide family in plants. During plant development, CLE gene family members play a pivotal role in regulating cell-to-cell communication and stem cell maintenance. However, the evolutionary process and functional importance of CLEs are unclear in Brassicaceae. In this study, a total of 70 BnCLEs were identified in Brassica napus (2n = 4x = 38, A_nC_n): 32 from the A_n subgenome, 36 from the C_n subgenome, and 2 from the unanchored subgenome. Meanwhile, 29 BrCLE and 32 BoCLE genes were explored in Brassica rapa (2n = 2x = 20, A_r) and Brassica oleracea (2n = 2x = 18, C_o). Phylogenetic analysis revealed that 163 CLEs derived from three Brassica species and Arabidopsis thaliana can be divided into seven subfamilies. Homology and synteny analyses indicated whole-genome triplication (WGT) and segmental duplication may be the major contributors to the expansion of CLE family. In addition, RNA-seq and qPCR analysis indicated that 19 and 16 BnCLEs were more highly expressed in immature seeds and roots than in other tissues. Some CLE gene pairs exhibited different expression patterns in the same tissue, which indicated possible functional divergence. Furthermore, genetic variations and regional association mapping analysis indicated that 12 BnCLEs were potential genes for regulating important agronomic traits. This study provided valuable information to understand the molecular evolution and biological function of CLEs in B. napus and its diploid progenitors, which will be helpful for genetic improvement of high-yield breeding in B. napus.

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.

Increased ACS Enzyme Dosage Causes Initiation of Climacteric Ethylene Production in Tomato

Fruits of wild tomato species show different ethylene-dependent ripening characteristics, such as variations in fruit color and whether they exhibit a climacteric or nonclimacteric ripening transition. 1-Aminocyclopropane-1-carboxylic acid (ACC) synthase (ACS) and ACC oxidase (ACO) are key enzymes in the ethylene biosynthetic pathway encoded by multigene families. Gene duplication is a primary driver of plant diversification and angiosperm evolution. Here, interspecific variations in the molecular regulation of ethylene biosynthesis and perception during fruit ripening in domesticated and wild tomatoes were investigated. Results showed that the activated ACS genes were increased in number in red-ripe tomato fruits than in green-ripe tomato fruits; therefore, elevated dosage of ACS enzyme promoted ripening ethylene production. Results showed that the expression of three ACS isogenes ACS1A, ACS2, and ACS4, which are involved in autocatalytic ethylene production, was higher in red-ripe tomato fruits than in green-ripe tomato fruits. Elevated ACS enzyme dosage promoted ethylene production, which corresponded to the climacteric response of red-ripe tomato fruits. The data suggest that autoinhibitory ethylene production is common to all tomato species, while autocatalytic ethylene production is specific to red-ripe species. The essential regulators Non-ripening (NOR) and Ripening-Inhibitor (RIN) have experienced gene activation and overlapped with increasing ACS enzyme dosage. These complex levels of transcript regulation link higher ethylene production with spatiotemporal modulation of gene expression in red-ripe tomato species. Taken together, this study shows that bursts in ethylene production that accompany fruit color changes in red-ripe tomatoes are likely to be an evolutionary adaptation for seed dispersal.

Retrieving a disrupted gene encoding phospholipase A for fibre enhancement in allotetraploid cultivated cotton

After polyploidization originated from one interspecific hybridization event in Gossypium, Gossypium barbadense evolved to produce extra-long staple fibres than Gossypium hirsutum (Upland cotton), which produces a higher fibre yield. The genomic diversity between G. barbadense and G. hirsutum thus provides a genetic basis for fibre trait variation. Recently, rapid accumulation of gene disruption or deleterious mutation was reported in allotetraploid cotton genomes, with unknown impacts on fibre traits. Here, we identified gene disruptions in allotetraploid G. hirsutum (18.14%) and G. barbadense (17.38%) through comparison with their presumed diploid progenitors. Relative to conserved genes, these disrupted genes exhibited faster evolution rate, lower expression level and altered gene co-expression networks. Within a module regulating fibre elongation, a hub gene experienced gene disruption in G. hirsutum after polyploidization, with a 2-bp deletion in the coding region of GhNPLA1D introducing early termination of translation. This deletion was observed in all of the 34 G. hirsutum landraces and 36 G. hirsutum cultivars, but not in 96% of 57 G. barbadense accessions. Retrieving the disrupted gene GhNPLA1D using its homoeolog GhNPLA1A achieved longer fibre length in G. hirsutum. Further enzyme activity and lipids analysis confirmed that GhNPLA1A encodes a typical phospholipase A and promotes cotton fibre elongation via elevating intracellular levels of linolenic acid and 34:3 phosphatidylinositol. Our work opens a strategy for identifying disrupted genes and retrieving their functions in ways that can provide valuable resources for accelerating fibre trait enhancement in cotton breeding.

Cis-regulatory sequences in plants: Their importance, discovery, and future challenges

The identification and characterization of cis-regulatory DNA sequences and how they function to coordinate responses to developmental and environmental cues is of paramount importance to plant biology. Key to these regulatory processes are cis-regulatory modules (CRMs), which include enhancers and silencers. Despite the extraordinary advances in high-quality sequence assemblies and genome annotations, the identification and understanding of CRMs, and how they regulate gene expression, lag significantly behind. This is especially true for their distinguishing characteristics and activity states. Here, we review the current knowledge on CRMs and breakthrough technologies enabling identification, characterization, and validation of CRMs; we compare the genomic distributions of CRMs with respect to their target genes between different plant species, and discuss the role of transposable elements harboring CRMs in the evolution of gene expression. This is an exciting time to study cis-regulomes in plants; however, significant existing challenges need to be overcome to fully understand and appreciate the role of CRMs in plant biology and in crop improvement.

Genome-wide identification of alcohol dehydrogenase (ADH) gene family under waterlogging stress in wheat (Triticum aestivum)

Background

Alcohol dehydrogenase (ADH) plays an important role in plant survival under anaerobic conditions. Although some research about ADH in many plants have been carried out, the bioinformatics analysis of the ADH gene family from Triticum aestivum and their response to abiotic stress is unclear.

Methods

A total of 22 ADH genes were identified from the wheat genome, and these genes could be divided into two subfamilies (subfamily I and subfamily II). All TaADH genes belonged to the Medium-chain ADH subfamily. Sequence alignment analysis showed that all TaADH proteins contained a conservative GroES-like domain and Zinc-binding domain. A total of 64 duplicated gene pairs were found, and the Ka/Ks value of these gene pairs was less than 1, which indicated that these genes were relatively conservative and did not change greatly in the process of duplication.

Results

The organizational analysis showed that nine TaADH genes were highly expressed in all organs, and the rest of TaADH genes had tissue specificity. Cis-acting element analysis showed that almost all of the TaADH genes contained an anaerobic response element. The expression levels of ADH gene in waterlogging tolerant and waterlogging sensitive wheat seeds were analyzed by quantitative real-time PCR (qRT-PCR). This showed that some key ADH genes were significantly responsive to waterlogging stress at the seed germination stage, and the response of waterlogging tolerant and waterlogging sensitive wheat seeds to waterlogging stress was regulated by different ADH genes. The results may be helpful to further study the function of TaADH genes and to determine the candidate gene for wheat stress resistance breeding.

The Transposable Element Environment of Human Genes Differs According to Their Duplication Status and Essentiality

Transposable elements (TEs) are major components of eukaryotic genomes and represent approximately 45% of the human genome. TEs can be important sources of novelty in genomes and there is increasing evidence that TEs contribute to the evolution of gene regulation in mammals. Gene duplication is an evolutionary mechanism that also provides new genetic material and opportunities to acquire new functions. To investigate how duplicated genes are maintained in genomes, here, we explored the TE environment of duplicated and singleton genes. We found that singleton genes have more short-interspersed nuclear elements and DNA transposons in their vicinity than duplicated genes, whereas long-interspersed nuclear elements and long-terminal repeat retrotransposons have accumulated more near duplicated genes. We also discovered that this result is highly associated with the degree of essentiality of the genes with an unexpected accumulation of short-interspersed nuclear elements and DNA transposons around the more-essential genes. Our results underline the importance of taking into account the TE environment of genes to better understand how duplicated genes are maintained in genomes.

Primary Metabolism co-Opted for Defensive Chemical Production in the Carabid Beetle, Harpalus pensylvanicus

Of the approximately one million described insect species, ground beetles (Coleoptera: Carabidae) have long captivated the attention of evolutionary biologists due to the diversity of defensive compounds they synthesize. Produced using defensive glands in the abdomen, ground beetle chemicals represent over 250 compounds including predator-deterring formic acid, which has evolved as a defensive strategy at least three times across Insecta. Despite being a widespread method of defense, formic acid biosynthesis is poorly understood in insects. Previous studies have suggested that the folate cycle of one-carbon (C1) metabolism, a pathway involved in nucleotide biosynthesis, may play a key role in defensive-grade formic acid production in ants. Here, we report on the defensive gland transcriptome of the formic acid-producing ground beetle Harpalus pensylvanicus. The full suite of genes involved in the folate cycle of C1 metabolism are significantly differentially expressed in the defensive glands of H. pensylvanicus when compared to gene expression profiles in the rest of the body. We also find support for two additional pathways potentially involved in the biosynthesis of defensive-grade formic acid, the kynurenine pathway and the methionine salvage cycle. Additionally, we have found an array of differentially expressed genes in the secretory lobes involved in the biosynthesis and transport of cofactors necessary for formic acid biosynthesis, as well as genes presumably involved in the detoxification of secondary metabolites including formic acid. We also provide insight into the evolution of the predominant gene family involved in the folate cycle (MTHFD) and suggest that high expression of folate cycle genes rather than gene duplication and/or neofunctionalization may be more important for defensive-grade formic acid biosynthesis in H. pensylvanicus. This provides the first evidence in Coleoptera and one of a few examples in Insecta of a primary metabolic process being co-opted for defensive chemical biosynthesis. Our results shed light on potential mechanisms of formic acid biosynthesis in the defensive glands of a ground beetle and provide a foundation for further studies into the evolution of formic acid-based chemical defense strategies in insects.

An Overview of Duplicated Gene Detection Methods: Why the Duplication Mechanism Has to Be Accounted for in Their Choice

Gene duplication is an important evolutionary mechanism allowing to provide new genetic material and thus opportunities to acquire new gene functions for an organism, with major implications such as speciation events. Various processes are known to allow a gene to be duplicated and different models explain how duplicated genes can be maintained in genomes. Due to their particular importance, the identification of duplicated genes is essential when studying genome evolution but it can still be a challenge due to the various fates duplicated genes can encounter. In this review, we first describe the evolutionary processes allowing the formation of duplicated genes but also describe the various bioinformatic approaches that can be used to identify them in genome sequences. Indeed, these bioinformatic approaches differ according to the underlying duplication mechanism. Hence, understanding the specificity of the duplicated genes of interest is a great asset for tool selection and should be taken into account when exploring a biological question.

Identification, genetic diversity and recombination analysis of Watermelon Mosaic Virus isolates

Watermelon mosaic virus (WMV) is an important virus causing adverse effects on cucurbits throughout the world. In this study, we recorded WMV infection in the watermelon (Citrullus lanatus)-growing area of Alwar and Sikar in districts of Rajasthan, India. The RT-PCR-based detection was performed to confirm the presence of WMV, by using potyvirus-degenerated coat protein primers. Further, the complete genome sequences of two WMV isolates were compared with previously reported genome sequences. The complete genome of each isolate was 10,030 nt long, excluding the poly-A tails. Phylogeny relationships of the WMV isolates in the present study revealed the presence of uneven evolutionary pressure among the different WMV viral genomic segments. The analysis revealed that all the WMV isolates were divided into three clusters and the Indian WMV isolates cluster together with the French isolate. Recombination analysis of WMV exhibited significant recombination hotspots in the P1, NIa-Pro and Nib-CP regions. Our finding highlights the importance of genetic variability and recombination analysis to provide a better understanding of WMV molecular diversity.

Comprehensive phylogenomic analysis of ERF genes in sorghum provides clues to the evolution of gene functions and redundancy among gene family members

APETALA2/Ethylene-Responsive transcription factors (AP2/ERF), with their multifunctional roles in plant development, hormone signaling and stress tolerance, are important candidates for engineering crop plants. Here, we report identification and analysis of gene structure, phylogenetic distribution, expression, chromosomal localization and cis-acting promoter analysis of AP2/ERF genes in the C4 crop plant sorghum. We identified 158 ERF genes in sorghum with 52 of them encoding dehydration-responsive binding elements (DREB) while 106 code for ERF subfamily proteins. Phylogenetic analysis organized sorghum ERF proteins into 11 distinct groups exhibiting clade-specific expansion. About 68% ERF genes have paralogs indicating gene duplications as major cause of expansion of ERF family in sorghum. Analysis of spatiotemporal expression patterns using publicly available data revealed their tissue/genotype-preferential accumulation. In addition, 40 ERF genes exhibited differential accumulation in response to heat and/or drought stress. About 25% of the segmental gene pairs and eleven tandem duplicated genes exhibited high correlation (> 0.7) in their expression patterns indicating genetic redundancy. Comparative phylogenomic analysis of sorghum ERFs with 74 genetically characterized ERF genes from other plant species provided significant clues to sorghum ERF functions. Overall data generated here provides an overview of evolutionary relationship among ERF gene family members in sorghum and with respect to previously characterized ERF genes from other plant species. This information will be instrumental in initiating functional genomic studies of ERF candidates in sorghum.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein–protein, genetic and homologous interactions, and directed protein–DNA, regulatory and miRNA–mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

Asymmetric Redundancy of ZERZAUST and ZERZAUST HOMOLOG in Different Accessions of Arabidopsis thaliana

Divergence among duplicate genes is one of the important sources of evolutionary innovation. But, the contribution of duplicate divergence to variation in Arabidopsis accessions is sparsely known. Recently, we studied the role of a cell wall localized protein, ZERZAUST (ZET), in Landsberg erecta (Ler) accession, lack of which results in aberrant plant morphology. Here, we present the study of ZET in Columbia (Col) accession, which not only showed differential expression patterns in comparison to Ler, but also revealed its close homolog, ZERZAUST HOMOLOG (ZETH). Although, genetic analysis implied redundancy, expression analysis revealed divergence, with ZETH showing minimal expression in both Col and Ler. In addition, ZETH shows relatively higher expression levels in Col compared to Ler. Our data also reveal compensatory up-regulation of ZETH in Col, but not in Ler, implying it is perhaps dispensable in Ler. However, a novel CRISPR/Cas9-induced zeth allele confirmed that ZETH has residual activity in Ler. Finally, the synergistic interaction of the receptor-like kinase gene, ERECTA with ZET in ameliorating morphological defects suggests crucial role of modifiers on plant phenotype. The results provide genetic evidence for accession-specific differences in compensation mechanism and asymmetric gene contribution. Thus, our work reveals a novel example for how weakly expressed homologs contribute to diversity among accessions.

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.

Landscape of gene transposition-duplication within the Brassicaceae family

We developed the CLfinder-OrthNet pipeline that detects co-linearity among multiple closely related genomes, finds orthologous gene groups, and encodes the evolutionary history of each orthologue group into a representative network (OrthNet). Using a search based on network topology, we identified 1,394 OrthNets that included gene transposition-duplication (tr-d) events, out of 17,432 identified in six Brassicaceae genomes. Occurrences of tr-d shared by subsets of Brassicaceae genomes mirrored the divergence times between the genomes and their repeat contents. The majority of tr-d events resulted in truncated open reading frames (ORFs) in the duplicated loci. However, the duplicates with complete ORFs were significantly more frequent than expected from random events. These were derived from older tr-d events and had a higher chance of being expressed. We also found an enrichment of tr-d events with complete loss of intergenic sequence conservation between the original and duplicated loci. Finally, we identified tr-d events uniquely found in two extremophytes among the six Brassicaceae genomes, including tr-d of SALT TOLERANCE 32 and ZINC TRANSPORTER 3 that relate to their adaptive evolution. CLfinder-OrthNet provides a flexible toolkit to compare gene order, visualize evolutionary paths among orthologues as networks, and identify gene loci that share an evolutionary history.

Dual-targeting by CRISPR/Cas9 leads to efficient point mutagenesis but only rare targeted deletions in the rice genome

The present study investigated the efficiency of CRISPR/Cas9 in creating genomic deletions as the basis of its application in removing selection marker genes or the intergenic regions. Three loci, representing a transgene and two rice genes, were targeted at two sites each, in separate experiments, and the deletion of the defined fragments was investigated by PCR and sequencing. Genomic deletions were found at a low rate among the transformed callus lines that could be isolated, cultured, and regenerated into plants harboring the deletion. However, randomly regenerated plants showed mixed genomic effects, and generally did not harbor heritable genomic deletions. To determine whether point mutations occurred at each targeted site, a total of 114 plants consisting of primary transgenic lines and their progeny were analyzed. Ninety-three plants showed targeting, 60 of which were targeted at both sites. The presence of point mutations at both sites was correlated with the guide RNA efficiency. In summary, genomic deletions through dual-targeting by the paired-guide RNAs were generally observed in callus, while de novo point mutations at one or both sites occurred at high rates in transgenic plants and their progeny, generating a variety of insertion–deletions or single-nucleotide variations. In this study, point mutations were exceedingly favored over genomic deletions; therefore, for the recovery of plant lines harboring targeted deletions, identifying early transformed clones harboring the deletions, and isolating them for plant regeneration is recommended.

Evaluation of the diversity and phylogenetic implications of NAC transcription factor members of four reference species from the different embryophytic plant groups

NAC transcription factors (TFs) are one of the largest and important TF family that are involved in the regulation of plant growth and development. They are characterized by a highly conserved N-terminal domain and a variable C-terminal domain. In the present study, the amino acid sequences of NAC TFs from four embryophytic plant species viz. Arabidopsis thaliana (Angiosperm), Picea abies (Gymnosperm), Selaginella moellendorffii (Pteridophyte) and Physcomitrella patens (Bryophyte) as reference of the different plant groups were collected from the Plant Transcription Factor Database (PTFD) and the phylogenetic relationships were evaluated. The phylogenetic tree revealed that the majority of the NAC members were interspersed in the major subgroups that indicated the expansion of the NAC members predates the speciation events. Thirty one (31), five (05), one (1) and ten (10) paralog pairs were determined respectively for Arabidopsis, Picea, Selaginella and Physcomitrella. The structure-function relationship of paralog pairs were inferred from the phylogenetic tree of combined set of paralogous gene pairs by studying the prevalence of flanking regions and motif analysis of the NAC proteins. The motif analysis revealed the presence of an N-terminal conserved domain, a characteristic of the majority of NAC family proteins. Conserved motifs in the C-terminal region were absent in the majority of the protein sequences except few members in Arabidopsis and Physcomitrella. Also the time of gene duplication of the paralog pairs were calculated that revealed the duplication events occurred between 4.48 and 45.94 MYA Arabidopsis, 167.57-532.86 MYA in Picea, and 29.12-53.53 MYA in Physcomitrella.

Expression Divergence as an Evolutionary Alternative Mechanism Adopted by Two Rice Subspecies Against Rice Blast Infection

BACKGROUND

Rice (Oryza sativa L.) is one of the most important crops that serves as staple food for ~ 50% of the human population worldwide. Some important agronomic traits that allow rice to cope with numerous abiotic and biotic stresses have been selected and fixed during domestication. Knowledge on how expression divergence of genes gradually contributes to phenotypic differentiation in response to biotic stress and their contribution to rice population speciation is still limited.

RESULTS

Here, we explored gene expression divergence between a japonica rice cultivar Nipponbare and an indica rice cultivar 93-11 in response to invasion by the filamentous ascomycete fungus Magnaporthe oryzae (Pyricularia oryzae), a plant pathogen that causes significant loss to rice production worldwide. We investigated differentially expressed genes in the two cultivars and observed that evolutionarily conserved orthologous genes showed highly variable expression patterns under rice blast infection. Analysis of promoter region of these differentially expressed orthologous genes revealed the existence of cis-regulatory elements associated with the differentiated expression pattern of these genes in the two rice cultivars. Further comparison of these regions in global rice population indicated their fixation and close relationship with rice population divergence.

CONCLUSION

We proposed that variation in the expression patterns of these orthologous genes mediated by cis-regulatory elements in the two rice cultivars, may constitute an alternative evolutionary mechanism that distinguishes these two genetically and ecologically divergent rice cultivars in response to M. oryzae infection.

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.

Bifunctional Chloroplastic DJ-1B from Arabidopsis thaliana is an Oxidation-Robust Holdase and a Glyoxalase Sensitive to H₂O₂

Members of the DJ-1 protein family are multifunctional enzymes whose loss increases the susceptibility of the cell to oxidative stress. However, little is known about the function of the plant DJ-1 homologs. Therefore, we analyzed the effect of oxidation on the structure and function of chloroplastic AtDJ-1B and studied the phenotype of T-DNA lines lacking the protein. In vitro oxidation of AtDJ-1B with H₂O₂ lowers its glyoxalase activity, but has no effect on its holdase chaperone function. Remarkably, upon oxidation, the thermostability of AtDJ-1B increases with no significant alteration of the overall secondary structure. Moreover, we found that AtDJ-1B transcript levels are invariable, and loss of AtDJ-1B does not affect plant viability, growth and stress response. All in all, two discrete functions of AtDJ-1B respond differently to H₂O₂, and AtDJ-1B is not essential for plant development under stress.

Unbiased subgenome evolution following a recent whole-genome duplication in pear (Pyrus bretschneideri Rehd.)

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear.

Cerium enhances germination and shoot growth, and alters mineral nutrient concentration in rice

Cerium (Ce) belongs to the rare earth elements (REEs), and although it is not essential for plants, it can stimulate growth and other physiological processes. The objective of this research was to evaluate the effect of Ce on seed germination, initial seedling growth, and vegetative growth in rice (Oryza sativa L.) cv. Morelos A-98. During the germination process, the seeds were treated with Ce concentrations of 0, 4, 8, and 12 μM; after 5 d, germination percentage was recorded and after 10 d seedling growth was measured. For vegetative growth, a hydroponic system was established where 14-d-old plants without previous Ce treatment were transferred into nutrient solution. After two weeks of acclimatizing, 0, 25, 50, and 100 μM Ce were added to the nutrient solution for 28 d. Ce significantly increased germination and the initial growth variables of the seedlings. During vegetative growth, Ce increased plant height, number of tillers, root volume, and shoot fresh and dry biomass, without affecting root biomass weight. With low Ce concentrations (25 and 50 μM), the concentrations of chlorophylls and amino acids in the shoots were similar to those in the control, like amino acid concentration in the roots at a concentration of 25 μM Ce. Conversely, the concentration of total sugars increased in the shoot with the application of 25, 50, and 100 μM Ce, and in the roots with the application of 50 μM Ce. Also, Ce did not affect the concentration of macro or micronutrients in the shoots. However, in the roots, the high Ce concentration decreased the concentrations of Ca, Fe, Mn, and Zn, while the Mg concentration increased. Our results indicate that Ce, at the right concentrations, can function as a biostimulant in rice germination and growth.

Different Modes of Gene Duplication Show Divergent Evolutionary Patterns and Contribute Differently to the Expansion of Gene Families Involved in Important Fruit Traits in Pear (Pyrus bretschneideri)

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30-45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

Gene-body CG methylation and divergent expression of duplicate genes in rice

Gene and genome duplication fosters genetic novelty, but redundant gene copies would undergo mutational decay unless preserved via selective or neutral forces. Molecular mechanisms mediating duplicate preservation remain incompletely understood. Several recent studies showed an association between DNA methylation and expression divergence of duplicated genes and suggested a role of epigenetic mechanism in duplicate retention. Here, we compare genome-wide gene-body CG methylation (BCGM) and duplicate gene expression between a rice mutant null for OsMet1-2(a major CG methytransferase in rice) and its isogenic wild-type. We demonstrate a causal link between BCGM divergence and expression difference of duplicate copies. Interestingly, the higher- and lower-expressing copies of duplicates as separate groups show broadly different responses with respect to direction of expression alteration upon loss of BCGM. A role for BCGM in conditioning expression divergence between copies of duplicates generally holds for duplicates generated by whole genome duplication (WGD) or by small-scale duplication processes. However, differences are evident among these categories, including a higher proportion of WGD duplicates manifesting expression alteration, and differential propensities to lose BCGM by the higher- and lower-expression copies in the mutant. Together, our results support the notion that differential epigenetic marking may facilitate long-term retention of duplicate genes.

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.

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.

MYB31/MYB42 Syntelogs Exhibit Divergent Regulation of Phenylpropanoid Genes in Maize, Sorghum and Rice

ZmMYB31 and ZmMYB42 are R2R3-MYB transcription factors implicated in the regulation of phenylpropanoid genes in maize. Here, we tested the hypothesis that the regulatory function of MYB31 and MYB42 is conserved in other monocots, specifically in sorghum and rice. We demonstrate that syntelogs of MYB31 and MYB42 do bind to phenylpropanoid genes that function in all stages of the pathway and in different tissues along the developmental gradient of seedling leaves. We found that caffeic acid O-methyltransferase (COMT1) is a common target of MYB31 and MYB42 in the mature leaf tissues of maize, sorghum and rice, as evidenced by Chromatin immunoprecipitation (ChIP) experiments. In contrast, 4-coumarate-CoA ligase (4CL2), ferulate-5-hydroxylase (F5H), and caffeoyl shikimate esterase (CSE), were targeted by MYB31 or MYB42, but in a more species-specific fashion. Our results revealed MYB31 and MYB42 participation in auto- and cross-regulation in all three species. Apart from a limited conservation of regulatory modules, MYB31 and MYB42 syntelogs appear to have undergone subfunctionalization following gene duplication and divergence of maize, sorghum, and rice. Elucidating the different regulatory roles of these syntelogs in the context of positive transcriptional activators may help guide attempts to alter the flux of intermediates towards lignin production in biofuel grasses.

Tools of the trade: studying molecular networks in plants

Driven by recent technological improvements, genes can be now studied in a larger biological context. Genes and their protein products rarely operate as a single entity and large-scale mapping by protein-protein interactions can unveil the molecular complexes that form in the cell to carry out various functions. Expression analysis under multiple conditions, supplemented with protein-DNA binding data can highlight when genes are active and how they are regulated. Representing these data in networks and finding strongly connected sub-graphs has proven to be a powerful tool to predict the function of unknown genes. As such networks are gradually becoming available for various plant species, it becomes possible to study how networks evolve. This review summarizes currently available network data and related tools for plants. Furthermore we aim to provide an outlook of future analyses that can be done in plants based on work done in other fields.