Gene duplication and stress genomics in Brassicas: Current understanding and future prospects

Lineage-specific gene duplication and expansion of DUF1216 gene family in Brassicaceae

Proteins containing domain of unknown function (DUF) are prevalent in eukaryotic genome. The DUF1216 proteins possess a conserved DUF1216 domain resembling to the mediator protein of Arabidopsis RNA polymerase II transcriptional subunit-like protein. The DUF1216 family are specifically existed in Brassicaceae, however, no comprehensive evolutionary analysis of DUF1216 genes have been performed. We performed a first comprehensive genome-wide analysis of DUF1216 proteins in Brassicaceae. Totally 284 DUF1216 genes were identified in 27 Brassicaceae species and classified into four subfamilies on the basis of phylogenetic analysis. The analysis of gene structure and conserved motifs revealed that DUF1216 genes within the same subfamily exhibited similar intron/exon patterns and motif composition. The majority members of DUF1216 genes contain a signal peptide in the N-terminal, and the ninth position of the signal peptide in most DUF1216 is cysteine. Synteny analysis revealed that segmental duplication is a major mechanism for expanding of DUF1216 genes in Brassica oleracea, Brassica juncea, Brassica napus, Lepidium meyneii, and Brassica carinata, while in Arabidopsis thaliana and Capsella rubella, tandem duplication plays a major role in the expansion of the DUF1216 gene family. The analysis of Ka/Ks (non-synonymous substitution rate/synonymous substitution rate) ratios for DUF1216 paralogous indicated that most of gene pairs underwent purifying selection. DUF1216 genes displayed a specifically high expression in reproductive tissues in most Brassicaceae species, while its expression in Brassica juncea was specifically high in root. Our studies offered new insights into the phylogenetic relationships, gene structures and expressional patterns of DUF1216 members in Brassicaceae, which provides a foundation for future functional analysis.

Combined genomic and transcriptomic analysis reveals the contribution of tandem duplication genes to low-temperature adaptation in perennial ryegrass

Perennial ryegrass (Lolium perenne L.) is an agronomically important cool-season grass species that is widely used as forage for ruminant animal production and cultivated in temperate regions for the establishment of lawns. However, the underlying genetic mechanism of the response of L. perenne to low temperature is still unclear. In the present study, we performed a comprehensive study and identified 3,770 tandem duplication genes (TDGs) in L. perenne, and evolutionary analysis revealed that L. perenne might have undergone a duplication event approximately 7.69 Mya. GO and KEGG pathway functional analyses revealed that these TDGs were mainly enriched in photosynthesis, hormone-mediated signaling pathways and responses to various stresses, suggesting that TDGs contribute to the environmental adaptability of L. perenne. In addition, the expression profile analysis revealed that the expression levels of TDGs were highly conserved and significantly lower than those of all genes in different tissues, while the frequency of differentially expressed genes (DEGs) from TDGs was much higher than that of DEGs from all genes in response to low-temperature stress. Finally, in-depth analysis of the important and expanded gene family indicated that the members of the ELIP subfamily could rapidly respond to low temperature and persistently maintain higher expression levels during all low temperature stress time points, suggesting that ELIPs most likely mediate low temperature responses and help to facilitate adaptation to low temperature in L. perenne. Our results provide evidence for the genetic underpinning of low-temperature adaptation and valuable resources for practical application and genetic improvement for stress resistance in L. perenne.

Comprehensive analysis of cucumber RAV family genes and functional characterization of CsRAV1 in salt and ABA tolerance in cucumber

The RAV (related to ABI3 and VP1) transcription factors are specific and exist in plants, which contain a B3 DNA binding domain and/or an APETALA2 (AP2) DNA binding domain. RAVs have been extensively studied in plants, and more and more evidences show that RAVs are involved in various aspects of plant growth and development, stress resistance and hormone signal transduction. However, the systematic analysis of RAV family in cucumber is rarely reported. In this study, eight CsRAV genes were identified in cucumber genome and we further comprehensively analyzed their protein physicochemical properties, conserved domains, gene structure and phylogenetic relationships. The synteny analysis and gene duplications of CsRAV genes were also analysed. Cis-element analysis revealed that the CsRAVs promoter contained several elements related to plant hormones and abiotic stress. Expression analysis showed that NaCl and ABA could significantly induce CsRAV genes expression. Subcellular localization revealed that all CsRAVs were localized in the nucleus. In addition, 35S:CsRAV1 transgenic Arabidopsis and cucumber seedlings enhanced NaCl and ABA tolerance, revealing CsRAV1 may be an important regulator of abiotic stress response. In conclusion, comprehensive analysis of CsRAVs would provide certain reference for understanding the evolution and function of the CsRAV genes.

Use of Yucca (Yucca schidigera) Extracts as Biostimulants to Promote Germination and Early Vigor and as Natural Fungicides

Climate change is increasing drought and salinity in many cultivated areas, therefore threatening food production. There is a great demand for novel agricultural inputs able to maintain yield under the conditions imposed by the anthropogenic global warming. Biostimulants have been proposed as a useful tool to achieve this objective. We have investigated the biostimulant effect of different yucca (Yucca schidigera) extracts on plant growth at different stages of development under different abiotic stress conditions. The extracts were tested in the model plant Arabidopsis thaliana, and in three different crops; tomato (Solanum lycopersicum var microtom), broccoli (Brassica oleracea var. italica) and lettuce (Lactuca sativa var romana). We have found that the investigated extracts are able to promote germination and early vigor under drought/osmotic and salt stress induced either by sodium chloride or lithium chloride. This effect is particularly strong in Arabidopsis thaliana and in the Brassicaceae broccoli. We have also determined using antibiograms against the model yeast Saccharomyces cerevisiae that the evaluated extracts may be used also as a natural fungicide. The results in this report show that yucca extracts may be used to enhance early vigor in some crops and as a natural fungicide, providing a new and useful tool for farmers.

Evolutionary analysis and expression profiling of the HSP70 gene family in response to abiotic stresses in tomato (Solanum lycopersicum)

Heat shock protein 70 (HSP70) genes play essential roles in guarding plants against abiotic stresses, including heat, drought, and salt. In this study, the SlHSP70 gene family in tomatoes has been characterized using bioinformatic tools. 25 putative SlHSP70 genes in the tomato genome were found and classified into five subfamilies, with multi-subcellular localizations. Twelve pairs of gene duplications were identified, and segmental events were determined as the main factor for the gene family expansion. Based on public RNA-seq data, gene expression analysis identified the majority of genes expressed in the examined organelles. Further RNA-seq analysis and then quantitative RT-PCR validation showed that many SlHSP70 members are responsible for cellular feedback to heat, drought, and salt treatments, in which, at least five genes might be potential key players in the stress response. Our results provided a thorough overview of the SlHSP70 gene family in the tomato, which may be useful for the evolutionary and functional analysis of SlHSP70 under abiotic stress conditions.

Genome-wide identification, comprehensive characterization of transcription factors, cis-regulatory elements, protein homology, and protein interaction network of DREB gene family in Solanum lycopersicum

Tomato is a drought-sensitive crop which has high susceptibility to adverse climatic changes. Dehydration-responsive element-binding (DREB) are significant plant transcription factors that have a vital role in regulating plant abiotic stress tolerance by networking with DRE/CRT cis-regulatory elements in response to stresses. In this study, bioinformatics analysis was performed to conduct the genome-wide identification and characterization of DREB genes and promoter elements in Solanum lycopersicum. In genome-wide coverage, 58 SlDREB genes were discovered on 12 chromosomes that justified the criteria of the presence of AP2 domain as conserved motifs. Intron-exon organization and motif analysis showed consistency with phylogenetic analysis and confirmed the absence of the A3 class, thus dividing the SlDREB genes into five categories. Gene expansion was observed through tandem duplication and segmental duplication gene events in SlDREB genes. Ka/Ks values were calculated in ortholog pairs that indicated divergence time and occurrence of purification selection during the evolutionary period. Synteny analysis demonstrated that 32 out of 58 and 47 out of 58 SlDREB genes were orthologs to Arabidopsis and Solanum tuberosum, respectively. Subcellular localization predicted that SlDREB genes were present in the nucleus and performed primary functions in DNA binding to regulate the transcriptional processes according to gene ontology. Cis-acting regulatory element analysis revealed the presence of 103 motifs in 2.5-kbp upstream promoter sequences of 58 SlDREB genes. Five representative SlDREB proteins were selected from the resultant DREB subgroups for 3D protein modeling through the Phyre2 server. All models confirmed about 90% residues in the favorable region through Ramachandran plot analysis. Moreover, active catalytic sites and occurrence in disorder regions indicated the structural and functional flexibility of SlDREB proteins. Protein association networks through STRING software suggested the potential interactors that belong to different gene families and are involved in regulating similar functional and biological processes. Transcriptome data analysis has revealed that the SlDREB gene family is engaged in defense response against drought and heat stress conditions in tomato. Overall, this comprehensive research reveals the identification and characterization of SlDREB genes that provide potential knowledge for improving abiotic stress tolerance in tomato.

Genome-Scale Investigation of GARP Family Genes Reveals Their Pivotal Roles in Nutrient Stress Resistance in Allotetraploid Rapeseed

The GARP genes are plant-specific transcription factors (TFs) and play key roles in regulating plant development and abiotic stress resistance. However, few systematic analyses of GARPs have been reported in allotetraploid rapeseed (Brassica napus L.) yet. In the present study, a total of 146 BnaGARP members were identified from the rapeseed genome based on the sequence signature. The BnaGARP TFs were divided into five subfamilies: ARR, GLK, NIGT1/HRS1/HHO, KAN, and PHL subfamilies, and the members within the same subfamilies shared similar exon-intron structures and conserved motif configuration. Analyses of the Ka/Ks ratios indicated that the GARP family principally underwent purifying selection. Several cis-acting regulatory elements, essential for plant growth and diverse biotic and abiotic stresses, were identified in the promoter regions of BnaGARPs. Further, 29 putative miRNAs were identified to be targeting BnaGARPs. Differential expression of BnaGARPs under low nitrate, ammonium toxicity, limited phosphate, deficient boron, salt stress, and cadmium toxicity conditions indicated their potential involvement in diverse nutrient stress responses. Notably, BnaA9.HHO1 and BnaA1.HHO5 were simultaneously transcriptionally responsive to these nutrient stresses in both hoots and roots, which indicated that BnaA9.HHO1 and BnaA1.HHO5 might play a core role in regulating rapeseed resistance to nutrient stresses. Therefore, this study would enrich our understanding of molecular characteristics of the rapeseed GARPs and will provide valuable candidate genes for further in-depth study of the GARP-mediated nutrient stress resistance in rapeseed.

Genome-Wide Identification of Brassicaceae Hormone-Related Transcription Factors and Their Roles in Stress Adaptation and Plant Height Regulation in Allotetraploid Rapeseed

Phytohormone-related transcription factors (TFs) are involved in regulating stress responses and plant growth. However, systematic analysis of these TFs in Brassicaceae is limited, and their functions in stress adaptation and plant height (PH) regulation remain unclear. In this study, 2115 hormone-related TFs were identified in nine Brassicaceae species. Specific domains were found in several Brassicaceae hormone-related TFs, which may be associated with diverse functions. Syntenic analysis indicated that expansion of these genes was mainly caused by segmental duplication, with whole-genome duplication occurring in some species. Differential expression analysis and gene co-expression network analysis identified seven phytohormone-related TFs (BnaWRKY7, 21, 32, 38, 52, BnaGL3-4, and BnaAREB2-5) as possible key genes for cadmium (Cd) toxicity, salinity stress, and potassium (K) and nitrogen (N) deficiencies. Furthermore, BnaWRKY42 and BnaARR21 may play essential roles in plant height. Weighted gene co-expression network analysis (WGCNA) identified 15 phytohormone-related TFs and their potential target genes regulating stress adaptation and plant height. Among the above genes, BnaWRKY56 and BnaWRKY60 responded to four different stresses simultaneously, and BnaWRKY42 was identified in two dwarf rapeseeds. In summary, several candidate genes for stress resistance (BnaWRKY56 and BnaWRKY60) and plant height (BnaWRKY42) were identified. These findings should help elucidate the biological roles of Brassicaceae hormone-related TFs, and the identified candidate genes should provide a genetic resource for the potential development of stress-tolerant and dwarf oilseed plants.

Genome-Wide Identification and Expressional Profiling of the Metal Tolerance Protein Gene Family in Brassica napus

The Cation Diffusion Facilitator (CDF) family, also named Metal Tolerance Protein (MTP), is one of the gene families involved in heavy metal transport in plants. However, a comprehensive study of MTPs in Brassica napus has not been reported yet. In the present study, we identified 33 BnMTP genes from the rapeseed genome using bioinformatic analyses. Subsequently, we analyzed the phylogenetic relationship, gene structure, chromosome distribution, conserved domains, and motifs of the BnMTP gene family. The 33 BnMTPs were phylogenetically divided into three major clusters (Zn-CDFs, Fe/Zn-CDFs, and Mn-CDFs) and seven groups (group 1, 5, 6, 7, 8, 9, and 12). The structural characteristics of the BnMTP members were similar in the same group, but different among groups. Evolutionary analysis indicated that the BnMTP gene family mainly expanded through whole-genome duplication (WGD) and segmental duplication events. Moreover, the prediction of cis-acting elements and microRNA target sites suggested that BnMTPs might be involved in plant growth, development, and stress responses. In addition, we found the expression of 24 BnMTPs in rapeseed leaves or roots could respond to heavy metal ion treatments. These results provided an important basis for clarifying the biological functions of BnMTPs, especially in heavy metal detoxification, and will be helpful in the phytoremediation of heavy metal pollution in soil.

Identification and characterization of G protein-coupled receptors in Spodoptera frugiperda (Insecta: Lepidoptera)

Spodoptera frugiperda (Insecta: Lepidoptera) is a destructive invasive pest feeding on various plants and causing serious damage to several economically-important crops. G protein-coupled receptors (GPCRs) are cellular receptors that coordinate diverse signaling processes, associated with many physiological processes and disease states. However, less information about GPCRs had been reported in S. frugiperda, limiting the recognition of signaling system and in-depth studies of this pest. Here, a total of 167 GPCRs were identified in S. frugiperda. Compared with other insects, the GPCRs of S. frugiperda were significantly expanded. A large of tandem duplication and segmental duplication events were observed, which may be the key factor to increase the size of GPCR family. In detail, these expansion events mainly concentrate on biogenic amine receptors, neuropeptide and protein hormone receptors, which may be involved in feeding, reproduction, life span, and tolerance of S. frugiperda. Additionally, 17 Mth/Mthl members were identified in S. frugiperda, which may be similar to the evolutionary pattern of 16 Mth/Mthl members in Drosophila. Moreover, the expression patterns across different developmental stages of all GPCR genes were also analyzed. Among these, most of the GPCR genes are poorly expressed in S. frugiperda and some highly expressed GPCR genes help S. frugiperda adapt to the environment better, such as Rh6 and AkhR. In this study, all GPCRs in S. frugiperda were identified for the first time, which provided a basis for further revealing the role of these receptors in the physiological and behavioral regulation of this pest.

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.

Genome-Wide Analysis of the Late Embryogenesis Abundant (LEA) and Abscisic Acid-, Stress-, and Ripening-Induced (ASR) Gene Superfamily from Canavalia rosea and Their Roles in Salinity/Alkaline and Drought Tolerance

Canavalia rosea (bay bean), distributing in coastal areas or islands in tropical and subtropical regions, is an extremophile halophyte with good adaptability to seawater and drought. Late embryogenesis abundant (LEA) proteins typically accumulate in response to various abiotic stresses, including dehydration, salinity, high temperature, and cold, or during the late stage of seed development. Abscisic acid-, stress-, and ripening-induced (ASR) genes are stress and developmentally regulated plant-specific genes. In this study, we reported the first comprehensive survey of the LEA and ASR gene superfamily in C. rosea. A total of 84 CrLEAs and three CrASRs were identified in C. rosea and classified into nine groups. All CrLEAs and CrASRs harbored the conserved motif for their family proteins. Our results revealed that the CrLEA genes were widely distributed in different chromosomes, and all of the CrLEA/CrASR genes showed wide expression features in different tissues in C. rosea plants. Additionally, we introduced 10 genes from different groups into yeast to assess the functions of the CrLEAs/CrASRs. These results contribute to our understanding of LEA/ASR genes from halophytes and provide robust candidate genes for functional investigations in plant species adapted to extreme environments.

Sequence and functional analysis of MIR319 promoter homologs from Brassica juncea reveals regulatory diversification and altered expression under stress

KEY MESSAGE

Extensive regulatory divergence during development, abiotic stress and ABA regime observed amongst promoter homologs and homeologs of MIR319 from Brassica juncea. Gene duplication followed by sub-functionalization, neo-functionalization, and pseudogenization are routes to functional and adaptive diversification. The influence of polyploidy on protein-coding genes is well investigated but little is known about their impact on transcriptional regulation of MIRNA gene family. The present study was therefore performed with an aim to uncover regulatory diversification of MIR319 homologs and homeologs in Brassica juncea. We employed comparative genomics to identify and isolate six promoter homologs of MIR319 from B. juncea. Regulatory diversification was studied using analysis of reporter activity driven by BjMIR319 promoters in a heterologous system employing promoter-reporter fusion constructs. MIR319 is known to play important roles in leaf and flower development, and multiple stress responses. Reporter activity was therefore monitored during development, hormonal and stress regimes. In-silico analyses revealed differential distribution of cis-regulatory motifs and functional analysis revealed distinct spatiotemporal expression patterns. The significance of presence of selected cis-regulatory motifs corresponding to heat, cold, salt and ABA stress were further functionally validated. It was observed that promoter of Bj -MIR319a-A01 was upregulated in response to cold and salt stress, while promoter of Bj -MIR319c-A04 (D1) and Bj -MIR319c-A05 (FL) were downregulated in response to high temperature. In summary, comparative analysis of homologous promoters from Brassica juncea, an allopolyploid revealed extensive sequence and functional diversity. Spatiotemporal activity of reporter gene driven by BjMIR319 promoter was distinct, and partially overlapping with from those reported previously for A. thaliana. The present study clearly demonstrates regulatory divergence amongst promoter homologs of MIR319 in Brassica juncea during development and stress response, and underlines the urgent need for dissection of promoter function and detailed characterization including identification of interacting trans-factors. Genbank accession numbers: MT379853-MT379858.