Expression partitioning between genes duplicated by polyploidy under abiotic stress and during organ development

Integrative analyses of long and short-read RNA sequencing reveal the spliced isoform regulatory network of seedling growth dynamics in upland cotton

The polyploid genome of cotton has significantly increased the transcript complexity. Recent advances in full-length transcript sequencing are now widely used to characterize the complete landscape of transcriptional events. Such studies in cotton can help us to explore the genetic mechanisms of the cotton seedling growth. Through long-read single-molecule RNA sequencing, this study compared the transcriptomes of three yield contrasting genotypes of upland cotton. Our analysis identified different numbers of spliced isoforms from 31,166, 28,716, and 28,713 genes in SJ48, Z98, and DT8 cotton genotypes, respectively, most of which were novel compared to previous cotton reference transcriptomes, and showed significant differences in the number of exon structures and coding sequence length due to intron retention. Quantification of isoform expression revealed significant differences in expression in the root and leaf of each genotype. An array of key isoform target genes showed protein kinase or phosphorylation functions, and their protein interaction network contained most of the circadian oscillator proteins. Spliced isoforms from the GIGANTEA (GI) protien were differentially regulated in each genotype and might be expected to regulate translational activities, including the sequence and function of target proteins. In addition, these spliced isoforms generate diurnal expression profiles in cotton leaves, which may alter the transcriptional regulatory network of seedling growth. Silencing of the novel spliced GI isoform Gh_A02G0645_N17 significantly affected biomass traits, contributed to variable growth, and increased transcription of the early flowering pathway gene ELF in cotton. Our high-throughput hybrid sequencing results will be useful to dissect functional differences among spliced isoforms in the polyploid cotton genome.

Polyploidy and environmental stress response: a comparative study of fern gametophytes

Climate change is rapidly altering natural habitats and generating complex patterns of environmental stress. Ferns are major components of many forest understories and, given their independent gametophyte generation, may experience unique pressures in emerging temperature and drought regimes. Polyploidy is widespread in ferns and may provide a selective advantage in these rapidly changing environments. This work aimed to understand whether the gametophytes of allopolyploid ferns respond differently to climate-related physiological stress than their diploid parents. The experimental approach involved a multifactorial design with 27 treatment combinations including exposure to multiple levels of drought and temperature over three treatment durations, with recovery measured at multiple timepoints. We measured Chl fluorescence from over 2000 gametophytes to evaluate stress avoidance and tolerance in diploid and polyploid species. Polyploids generally showed a greater ability to avoid and/or tolerate a range of stress conditions compared with their diploid counterparts, suggesting that polyploidy may confer enhanced flexibility and resilience under climate stress. Overall, these results suggest that polyploidy may provide some resilience to climate change in mixed ploidy populations. However, all species remain susceptible to the impacts of extreme drought and heat stress.

Harnessing cold adaptation for postglacial colonisation: Galactinol synthase expression and raffinose accumulation in a polyploid and its progenitors

Allotetraploid white clover (Trifolium repens) formed during the last glaciation through hybridisation of two European diploid progenitors from restricted niches: one coastal, the other alpine. Here, we examine which hybridisation-derived molecular events may have underpinned white clover's postglacial niche expansion. We compared the transcriptomic frost responses of white clovers (an inbred line and an alpine-adapted ecotype), extant descendants of its progenitor species and a resynthesised white clover neopolyploid to identify genes that were exclusively frost-induced in the alpine progenitor and its derived subgenomes. From these analyses we identified galactinol synthase, the rate-limiting enzyme in biosynthesis of the cryoprotectant raffinose, and found that the extant descendants of the alpine progenitor as well as the neopolyploid white clover rapidly accumulated significantly more galactinol and raffinose than the coastal progenitor under cold stress. The frost-induced galactinol synthase expression and rapid raffinose accumulation derived from the alpine progenitor likely provided an advantage during early postglacial colonisation for white clover compared to its coastal progenitor.

Major chromosome rearrangements in intergeneric wheat × rye hybrids in compatible and incompatible crosses detected by GBS read coverage analysis

The presence of incompatibility alleles in primary amphidiploids constitutes a reproductive barrier in newly synthesized wheat-rye hybrids. To overcome this barrier, the genome stabilization process includes large-scale chromosome rearrangements. In incompatible crosses resulting in fertile amphidiploids, the elimination of one of the incompatible alleles Eml-A1 or Eml-R1b can occur already in the somatic tissue of the wheat × rye hybrid embryo. We observed that the interaction of incompatible loci Eml-A1 of wheat and Eml-R1b of rye after overcoming embryo lethality leads to hybrid sterility in primary triticale. During subsequent seed reproductions (R1, R2 or R3) most of the chromosomes of A, B, D and R subgenomes undergo rearrangement or eliminations to increase the fertility of the amphidiploid by natural selection. Genotyping-by-sequencing (GBS) coverage analysis showed that improved fertility is associated with the elimination of entire and partial chromosomes carrying factors that either cause the disruption of plant development in hybrid plants or lead to the restoration of the euploid number of chromosomes (2n = 56) in the absence of one of the incompatible alleles. Highly fertile offspring obtained in compatible and incompatible crosses can be successfully adapted for the production of triticale pre-breeding stocks.

Polyploidization: A Biological Force That Enhances Stress Resistance

Organisms with three or more complete sets of chromosomes are designated as polyploids. Polyploidy serves as a crucial pathway in biological evolution and enriches species diversity, which is demonstrated to have significant advantages in coping with both biotic stressors (such as diseases and pests) and abiotic stressors (like extreme temperatures, drought, and salinity), particularly in the context of ongoing global climate deterioration, increased agrochemical use, and industrialization. Polyploid cultivars have been developed to achieve higher yields and improved product quality. Numerous studies have shown that polyploids exhibit substantial enhancements in cell size and structure, physiological and biochemical traits, gene expression, and epigenetic modifications compared to their diploid counterparts. However, some research also suggested that increased stress tolerance might not always be associated with polyploidy. Therefore, a more comprehensive and detailed investigation is essential to complete the underlying stress tolerance mechanisms of polyploids. Thus, this review summarizes the mechanism of polyploid formation, the polyploid biochemical tolerance mechanism of abiotic and biotic stressors, and molecular regulatory networks that confer polyploidy stress tolerance, which can shed light on the theoretical foundation for future research.

Effects of Allopolyploidization and Homoeologous Chromosomal Segment Exchange on Homoeolog Expression in a Synthetic Allotetraploid Wheat under Variable Environmental Conditions

Allopolyploidy through the combination of divergent genomes into a common nucleus at doubled dosage is known as a potent genetic and evolutionary force. As a macromutation, a striking feature of allopolyploidy in comparison with other mutational processes is that 'genome shock' can be evoked, thereby generating rapid and saltational biological consequences. A major manifestation of genome shock is genome-wide gene expression rewiring, which previously remained to be fully elucidated. Here, using a large set of RNAseq-based transcriptomic data of a synthetic allotetraploid wheat (genome AADD) and its parental species, we performed in-depth analyses of changes in the genome-wide gene expression under diverse environmental conditions at the subgenome (homoeolog) level and investigated the additional effects of homoeologous chromosomal segment exchanges (abbreviated HEs). We show that allopolyploidy caused large-scale changes in gene expression that were variable across the conditions and exacerbated by both stresses and HEs. Moreover, although both subgenomes (A and D) showed clear commonality in the changes, they responded differentially under variable conditions. The subgenome- and condition-dependent differentially expressed genes were enriched for different gene ontology terms implicating different biological functions. Our results provide new insights into the direct impacts of allopolyploidy on condition-dependent changes in subgenome expression and the additional effects of HEs in nascent allopolyploidy.

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.

Subgenome-dominant expression and alternative splicing in response to Sclerotinia infection in polyploid Brassica napus and progenitors

Polyploidy has played an extensive role in the evolution of flowering plants. Allopolyploids, with subgenomes containing duplicated gene pairs called homeologs, can show rapid transcriptome changes including novel alternative splicing (AS) patterns. The extent to which abiotic stress modulates AS of homeologs is a nascent topic in polyploidy research. We subjected both resynthesized and natural lines of polyploid Brassica napus, along with the progenitors Brassica rapa and Brassica oleracea, to infection with the fungal pathogen Sclerotinia sclerotiorum. RNA-sequencing analyses revealed widespread divergence between polyploid subgenomes in both gene expression and AS patterns. Resynthesized B. napus displayed significantly more A and C subgenome biased homeologs under pathogen infection than during uninfected growth. Differential AS (DAS) in response to infection was highest in natural B. napus (12 709 DAS events) and lower in resynthesized B. napus (8863 DAS events). Natural B. napus had more upregulated events and fewer downregulated events. There was a global expression bias towards the B. oleracea-derived (C) subgenome in both resynthesized and natural B. napus, enhanced by widespread non-parental downregulation of the B. rapa-derived (A) homeolog. In the resynthesized B. napus, this resulted in a disproportionate C subgenome contribution to the pathogen defense response, characterized by biases in both transcript expression levels and the proportion of induced genes. Our results elucidate the complex ways in which Sclerotinia infection affects expression and AS of homeologous genes in resynthesized and natural B. napus.

Genome-Wide Investigation and Co-Expression Network Analysis of SBT Family Gene in Gossypium

Subtilases (SBTs), which belong to the serine peptidases, control plant development by regulating cell wall properties and the activity of extracellular signaling molecules, and affect all stages of the life cycle, such as seed development and germination, and responses to biotic and abiotic environments. In this study, 146 Gossypium hirsutum, 138 Gossypium barbadense, 89 Gossypium arboreum and 84 Gossypium raimondii SBTs were identified and divided into six subfamilies. Cotton SBTs are unevenly distributed on chromosomes. Synteny analysis showed that the members of SBT1 and SBT4 were expanded in cotton compared to Arabidopsis thaliana. Co-expression network analysis showed that six Gossypium arboreum SBT gene family members were in a network, among which five SBT1 genes and their Gossypium hirsutum and Arabidopsis thaliana direct homologues were down-regulated by salt treatment, indicating that the co-expression network might share conserved functions. Through co-expression network and annotation analysis, these SBTs may be involved in the biological processes of auxin transport, ABA signal transduction, cell wall repair and root tissue development. In summary, this study provides valuable information for the study of SBT genes in cotton and excavates SBT genes in response to salt stress, which provides ideas for cotton breeding for salinity resistance.

Genome-wide characterization of R2R3-MYB gene family in Santalum album and their expression analysis under cold stress

Sandalwood (Santalum album) is a high-value multifunctional tree species that is rich in aromatic substances and is used in medicine and global cosmetics. Due to the scarcity of land resources in tropical and subtropical regions, land in temperate regions is a potential resource for the development of S. album plantations in order to meet the needs of S. album production and medicine. The R2R3-MYB transcription factor family is one of the largest in plants and plays an important role in the response to various abiotic stresses. However, the R2R3-MYB gene family of S. album has not been studied. In this study, 144 R2R3-MYB genes were successfully identified in the assembly genome sequence, and their characteristics and expression patterns were investigated under various durations of low temperature stress. According to the findings, 31 of the 114 R2R3-MYB genes showed significant differences in expression after cold treatment. Combining transcriptome and weighted gene co-expression network analysis (WGCNA) revealed three key candidate genes (SaMYB098, SaMYB015, and SaMYB068) to be significantly involved in the regulation of cold resistance in S. album. The structural characteristics, evolution, and expression pattern of the R2R3-MYB gene in S. album were systematically examined at the whole genome level for the first time in this study. It will provide important information for future research into the function of the R2R3-MYB genes and the mechanism of cold stress response in S. album.

Molecular responses to salinity stress in Salix matsudana (Koidz) females and males

Sexual dimorphism has commonly been found in many species. The phenotypes of Salix matsudana females and males are different under salinity stress. An F₁ population was selected to compare the differences between males and females. As a result, males showed stronger roots and heavier dry weights than females. The unique molecular mechanisms of males and females under salinity stress were further analyzed based on the root transcriptome of males and females. Both males and females up-regulated systemic acquired resistance genes, such as ADH and oxygenase-related genes, to resist salt. Moreover, many other abiotic stress response genes were up-regulated in males to adjust to salinity stress, while females showed more down-regulation of nitrogen metabolism-related genes to decrease the harm from salinity stress. The research on salinity tolerance in Salix matsudana males and females would help to further understand sexual dimorphism under selection pressure and provide benefits to the ecological environment.

Genome-wide analysis of polyamine biosynthesis genes in wheat reveals gene expression specificity and involvement of STRE and MYB-elements in regulating polyamines under drought

BACKGROUND

Polyamines (PAs) are considered promising biostimulants that have diverse key roles during growth and stress responses in plants. Nevertheless, the molecular basis of these roles by PAs has not been completely realized even now, and unfortunately, the transcriptional analyses of the biosynthesis pathway in various wheat tissues have not been investigated under normal or stress conditions. In this research, the findings of genome-wide analyses of genes implicated in the PAs biosynthesis in wheat (ADC, Arginine decarboxylase; ODC, ornithine decarboxylase; AIH, agmatine iminohydrolase; NPL1, Nitrlase like protein 1; SAMDC, S-adenosylmethionine decarboxylase; SPDS, spermidine synthase; SPMS, spermine synthase and ACL5, thermospermine synthase) are shown.

RESULTS

In total, thirty PAs biosynthesis genes were identified. Analysis of gene structure, subcellular compartmentation and promoters were discussed. Furthermore, experimental gene expression analyses in roots, shoot axis, leaves, and spike tissues were investigated in adult wheat plants under control and drought conditions. Results revealed structural similarity within each gene family and revealed the identity of two new motifs that were conserved in SPDS, SPMS and ACL5. Analysis of the promoter elements revealed the incidence of conserved elements (STRE, CAAT-box, TATA-box, and MYB TF) in all promoters and highly conserved CREs in >80% of promoters (G-Box, ABRE, TGACG-motif, CGTCA-motif, as1, and MYC). The results of the quantification of PAs revealed higher levels of putrescine (Put) in the leaves and higher spermidine (Spd) in the other tissues. However, no spermine (Spm) was detected in the roots. Drought stress elevated Put level in the roots and the Spm in the leaves, shoots and roots, while decreased Put in spikes and elevated the total PAs levels in all tissues. Interestingly, PA biosynthesis genes showed tissue-specificity and some homoeologs of the same gene family showed differential gene expression during wheat development. Additionally, gene expression analysis showed that ODC is the Put biosynthesis path under drought stress in roots.

CONCLUSION

The information gained by this research offers important insights into the transcriptional regulation of PA biosynthesis in wheat that would result in more successful and consistent plant production.

Parental legacy versus regulatory innovation in salt stress responsiveness of allopolyploid cotton (Gossypium) species

Polyploidy provides an opportunity for evolutionary innovation and species diversification, especially under stressful conditions. In allopolyploids, the conditional dynamics of homoeologous gene expression can be either inherited from ancestral states pre-existing in the parental diploids or novel upon polyploidization, the latter potentially permitting a wider range of phenotypic responses to stresses. To gain insight into regulatory mechanisms underlying the diversity of salt resistance in Gossypium species, we compared global transcriptomic responses to modest salinity stress in two allotetraploid (AD-genome) cotton species, Gossypium hirsutum and G. mustelinum, relative to their model diploid progenitors (A-genome and D-genome). Multivariate and pairwise analyses of salt-responsive changes revealed a profound alteration of gene expression for about one third of the transcriptome. Transcriptional responses and associated functional implications of salt acclimation varied across species, as did species-specific coexpression modules among species and ploidy levels. Salt responsiveness in both allopolyploids was strongly biased toward the D-genome progenitor. A much lower level of transgressive downregulation was observed in the more salt-tolerant G. mustelinum than in the less tolerant G. hirsutum. By disentangling inherited effects from evolved responses, we show that expression biases that are not conditional upon salt stress approximately equally reflect parental legacy and regulatory novelty upon allopolyploidization, whereas stress-responsive biases are predominantly novel, or evolved, in allopolyploids. Overall, our work suggests that allopolyploid cottons acquired a wide range of stress response flexibility relative to their diploid ancestors, most likely mediated by complex suites of duplicated genes and regulatory factors.

How do suboptimal temperatures affect polyploid sterlet Acipenser ruthenus during early development?

Sturgeons are ancient fish exhibiting unique genome plasticity and a high tendency to produce spontaneously autopolyploid genome states. The temperature profiles of the rivers in which sturgeon live and reproduce have been severely altered by human intervention, and the effect of global warming is expected to cause further temperature shifts, which may be detrimental for early developmental stages with narrow windows of thermal tolerance. The comparison of the performance of diploid and autopolyploid sturgeon kept at unfavourable temperatures contributes to scientific knowledge of the effects of polyploid genome states on organisms and can shed light on the ability of polyploids to cope with human-induced alterations to natural conditions. Using the sterlet Acipenser ruthenus as a model species, we carried out conventional artificial fertilization, as well as the induction of the second polar body retention (SPBR), of the first mitotic division suppression (FMDS) and of the second polar body retention followed by the first mitotic division suppression (SPBR+FMDS). Two experiments were conducted to evaluate the effect of polyploidy on two basic performance parameters, survival and growth. In Experiment 1, fish belonging to untreated, SPBR-, FMDS- and SPBR+FMDS-induced groups were kept at 10, 16 and 20°C from the neurula stage until the end of endogenous feeding. In Experiment 2, larvae from the untreated and SPBR-induced groups were reared at 10, 16 and 20°C after their endogenous feeding transition for 3 weeks. Based on our findings, we report that the embryos, prelarvae and larvae of triploid A. ruthenus do not differ from diploids in their ability to survive, grow and develop under suboptimal temperature conditions, while the survival of tetraploids was significantly reduced even at the optimal temperature and even more so at temperatures far from the optimum. This was also the case in the 2n/4n mosaics observed in FMDS-induced group. Thus, we assume that in tetraploid and 2n/4n individuals, the limits of thermal tolerance are closer to the optimum than in diploids. We also conclude that the hexaploid genome state is probably lethal in A. ruthenus since none of the hexaploids or 3n/6n mosaics arising from the SPBR+FMDS induction survived the prelarval period.

Expression Level Dominance and Homeolog Expression Bias Upon Cold Stress in the F1 Hybrid Between the Invasive Sphagneticola trilobata and the Native S. calendulacea in South China, and Implications for Its Invasiveness

The role of hybridization is significant in biological invasion, and thermotolerance is a trait critical to range expansions. The South American Sphagneticola trilobata is now widespread in South China, threatening the native S. calendulacea by competition and hybridization. Furthermore, upon formation, their F1 hybrid can quickly replace both parents. In this study, the three taxa were used as a model to investigate the consequences of hybridization on cold tolerance, particularly the effect of subgenome dominance in the hybrid. Upon chilling treatments, physiological responses and transcriptome profiles were compared across different temperature points to understand their differential responses to cold. While both parents showed divergent responses, the hybrid’s responses showed an overall resemblance to S. calendulacea, but the contribution of homeolog expression bias to cold stress was not readily evident in the F1 hybrid possibly due to inherent bias that comes with the sampling location. Our findings provided insights into the role of gene expression in differential cold tolerance, and further contribute to predicting the invasive potential of other hybrids between S. trilobata and its congeners around the world.

Impact of polyploidy on plant tolerance to abiotic and biotic stresses

Polyploidy, defined as the coexistence of three or more complete sets of chromosomes in an organism's cells, is considered as a pivotal moving force in the evolutionary history of vascular plants and has played a major role in the domestication of several crops. In the last decades, improved cultivars of economically important species have been developed artificially by inducing autopolyploidy with chemical agents. Studies on diverse species have shown that the anatomical and physiological changes generated by either natural or artificial polyploidization can increase tolerance to abiotic and biotic stresses as well as disease resistance, which may positively impact on plant growth and net production. The aim of this work is to review the current literature regarding the link between plant ploidy level and tolerance to abiotic and biotic stressors, with an emphasis on the physiological and molecular mechanisms responsible for these effects, as well as their impact on the growth and development of both natural and artificially generated polyploids, during exposure to adverse environmental conditions. We focused on the analysis of those types of stressors in which more progress has been made in the knowledge of the putative morpho-physiological and/or molecular mechanisms involved, revealing both the factors in common, as well as those that need to be addressed in future research.

A complex interplay of Gβ and Gγ proteins regulates plant growth and defence traits in the allotetraploid Brassica juncea

Gene expression analysis coupled with in-planta studies showed that specific Gβγ combination regulates plant growth and defence traits in the allotetraploid Brassica juncea. Plant heterotrimeric G-proteins regulate a wide range of responses despite their limited repertoire of core components. The roles and functional interactions between different G-protein subunits are quite perplexing, which get further complicated with polyploidy. Here, we show that the allotetraploid Brassica juncea comprises multiple homologs of G-protein genes, encoding six BjuGβ and ten highly divergent BjuGγ subunit proteins, later being classified into type-A1, type-A2 and type-C Gγ proteins. The encoded BjuGβ and BjuGγ proteins shared close evolutionary relationship and have retained distinct spatio-temporal expression patterns during plant developmental stages and in response to the necrotrophic pathogen, Sclerotinia sclerotiorum. RNAi based suppression of BjuGβ and BjuGγ genes suggested functional overlap and selectivity of BjuGβs with three distinct BjuGγ type subunits, to regulate plant height (BjuGβγA2 and BjuGβγC), seed weight (BjuGβGγA1 and BjuGβGγC), silique size (BjuGβGγC) and pathogen response (BjuGβGγA1 and BjuGβGγC). Further, the triplicated BjuGβ genes, formed due to Brassica specific whole-genome-triplication event, showed differential involvement during pathogen response, wherein overexpression of BjuGβ2 displayed higher resistance to Sclerotinia infection. Taken together, our study demonstrates that multiple BjuGβ and BjuGγ proteins have retained distinct spatio-temporal expression and functional selectivity to regulate specific plant growth and defence traits in the oilseed B. juncea.

Comparative transcriptome profiling of a resistant vs susceptible bread wheat (Triticum aestivum L.) cultivar in response to water deficit and cold stress

Bread wheat (Triticum aestivum L.) is one of the most important agricultural plants wearing abiotic stresses, such as water deficit and cold, that cause its productivity reduction. Since resistance to abiotic factors is a multigenic trait, therefore modern genome-wide approaches can help to involve various genetic material in breeding. One technique is full transcriptome analysis that reveals groups of stress response genes serving marker-assisted selection markers. Comparing transcriptome profiles of the same genetic material under several stresses is essential and makes the whole picture. Here, we addressed this by studying the transcriptomic response to water deficit and cold stress for two evolutionarily distant bread wheat varieties: stress-resistant cv. Saratovskaya 29 (S29) and stress-sensitive cv. Yanetzkis Probat (YP). For the first time, transcriptomes for these cultivars grown under abiotic stress conditions were obtained using Illumina based MACE technology. We identified groups of genes involved in response to cold and water deficiency stresses, including responses to each stress factor and both factors simultaneously that may be candidates for resistance genes. We discovered a core group of genes that have a similar pattern of stress-induced expression changes. The particular expression pattern was revealed not only for the studied varieties but also for the published transcriptomic data on cv. Jing 411 and cv. Fielder. Comparative transcriptome profiling of cv. S29 and cv. YP in response to water deficit and cold stress confirmed the hypothesis that stress-induced expression change is unequal within a homeologous gene group. As a rule, at least one changed significantly while the others had a relatively lower expression. Also, we found several SNPs distributed throughout the genomes of cv. S29 and cv. YP and distinguished the studied varieties from each other and the reference cv. Chinese Spring. Our results provide new data for genomics-assisted breeding of stress-tolerant wheat cultivars.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Global insights into duplicated gene expression and alternative splicing in polyploid Brassica napus under heat, cold, and drought stress

Polyploidy has been a prevalent process during plant evolution and it has made a major impact on the structure and evolution of plant genomes. Many important crop plants are polyploid. There is considerable interest in expression patterns of duplicated genes in polyploids. Alternative splicing (AS) is a fundamental aspect of gene expression that produces multiple final transcript types from a single type of mRNAs. The effects of abiotic stress conditions on AS in polyploids has received little attention. We conducted a global transcriptome analysis of Brassica napus, an allotetraploid derived from B. rapa (A_T ) and B. oleracea (C_T ), by RNA-Seq of plants subjected to cold, heat, and drought stress treatments. Analyses of 27,360 pairs of duplicated genes revealed overall A_T subgenome biases in gene expression and C_T subgenome biases in the extent of alternative splicing under all three stress treatments. More genes increased in expression than decreased in response to the stresses. Negative correlations were found between expression levels and AS frequency for each type of AS. Cold stress produced the greatest changes in gene expression and AS. Cold-induced AS changes were more likely to be shared with those generated by drought than by heat stress. We used homeologs of FLC and CCA1 as case studies to show the dynamics of how duplicates in a polyploid respond to cold stress. Our results suggest that divergence in gene expression and AS patterns between duplicated genes may increase the flexibility of polyploids when responding to abiotic stressors.

Transcriptomic and anatomic profiling reveal the germination process of different wheat varieties in response to waterlogging stress

BACKGROUND

Waterlogging is one of the most serious abiotic stresses affecting wheat-growing regions in China. Considerable differences in waterlogging tolerance have been found among different wheat varieties, and the mechanisms governing the waterlogging tolerance of wheat seeds during germination have not been elucidated.

RESULTS

The results showed no significant difference between the germination rate of 'Bainong 207' (BN207) (after 72 h of waterlogging treatment) and that of the control seeds. However, the degree of emulsification and the degradation rate of endosperm cells under waterlogging stress were higher than those obtained with the control treatment, and the number of amyloplasts in the endosperm was significantly reduced by waterlogging. Transcriptomic data were obtained from seed samples (a total of 18 samples) of three wheat varieties, 'Zhoumai 22' (ZM22), BN207 and 'Bainong 607' (BN607), subjected to the waterlogging and control treatments. A comprehensive analysis identified a total of 2775 differentially expressed genes (DEGs). In addition, an analysis of the correlations among the expression difference levels of DEGs and the seed germination rates of the three wheat varieties under waterlogging stress revealed that the relative expression levels of 563 and 398 genes were positively and negatively correlated with the germination rate of the wheat seeds, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the difference in the waterlogging tolerance among the three wheat varieties was related to the abundance of key genes involved in the glycolysis pathway, the starch and sucrose metabolism pathway, and the lactose metabolism pathway. The alcohol dehydrogenase (ADH) gene in the endosperm of BN607 was induced immediately after short-term waterlogging, and the energy provided by the glycolysis pathway enabled the BN607 seeds to germinate as early as possible; in addition, the expression of the AP2/ERF transcription factor was upregulated to further enhance the waterlogging tolerance of this cultivar.

CONCLUSIONS

Taken together, the results of this study help elucidate the mechanisms through which different wheat varieties respond to waterlogging stress during germination.

Genome-wide identification and evolution of Dof transcription factor family in cultivated and ancestral cotton species

The DNA-binding with one finger (Dof) proteins are transcription factors involved in many biological processes in plants. To predict the evolutionary pattern, a genome-wide in-silico analysis of Dof TFs family in diploid (Gossypium arboreum and Gossypium raimondii) and allotetraploid (Gossypium hirsutum and Gossypium barbadense) cotton species were carried out. In G. arboreum, we have identified 58 non-redundant genes encoding Dof proteins renamed as GaDof (G. arboreum Dof), 55 Dof genes were identified in G. raimondii (GrDof), 89 were predicted ffrom G. hirsutum (GhDof) and the highest, 110 Dof genes were identified in G. barbadense (GbDof). The phylogenetic analysis, physical location, gene structure, conserved domain analyses were also investigated for G. arboreum, G. raimondii, and G. hirsutum. The gene expression pattern in G. hirsutum, at different growth stages, revealing the probable involvement of some GhDof genes in growth and development. These genes may improve seed germination and growth in cotton.

Knock-out of TERMINAL FLOWER 1 genes altered flowering time and plant architecture in Brassica napus

Background

TERMINAL FLOWER 1 (TFL1) is a member of phosphatidylethanolamine-binding protein (PEBP) family, which plays an important role in the determination of floral meristem identity and regulates flowering time in higher plants.

Results

Five BnaTFL1 gene copies were identified in the genome of Brassica napus. The phylogenetic analysis indicated that all five BnaTFL1 gene copies were clustered with their corresponding homologous copies in the ancestral species, B. rapa and B. oleracea. The expression of the BnaTFL1s were confined to flower buds, flowers, seeds, siliques and stem tissues and displayed distinct expression profiles. Knockout mutants of BnaC03.TFL1 generated by CRISPR/Cas9 exhibited early flowering phenotype, while the knockout mutants of the other gene copies had similar flowering time as the wild type. Furthermore, knock-out mutants of individual BnaTFL1 gene copy displayed altered plant architecture. The plant height, branch initiation height, branch number, silique number, number of seeds per silique and number of siliques on the main inflorescence were significantly reduced in the BnaTFL1 mutants.

Conclusions

Our results indicated that BnaC03.TFL1 negatively regulates flowering time in B. napus. BnaC03.TFL1 together with the other BnaTFL1 paralogues are essential for controlling the plant architecture.

Salt-tolerance diversity in diploid and polyploid cotton (Gossypium) species

The development of salt-tolerant genotypes is pivotal for the effective utilization of salinized land and to increase global crop productivity. Several cotton species comprise the most important source of textile fibers globally, and these are increasingly grown on marginal or increasingly saline agroecosystems. The allopolyploid cotton species also provide a model system for polyploid research, of relevance here because polyploidy was suggested to be associated with increased adaptation to stress. To evaluate the genetic variation of salt tolerance among cotton species, 17 diverse accessions of allopolyploid (AD-genome) and diploid (A- and D-genome) Gossypium were evaluated for a total of 29 morphological and physiological traits associated with salt tolerance. For most morphological and physiological traits, cotton accessions showed highly variable responses to 2 weeks of exposure to moderate (50 mm NaCl) and high (100 mm NaCl) hydroponic salinity treatments. Our results showed that the most salt-tolerant species were the allopolyploid Gossypium mustelinum from north-east Brazil, the D-genome diploid Gossypium klotzschianum from the Galapagos Islands, followed by the A-genome diploids of Africa and Asia. Generally, A-genome accessions outperformed D-genome cottons under salinity conditions. Allopolyploid accessions from either diploid genomic group did not show significant differences in salt tolerance, but they were more similar to one of the two progenitor lineages. Our findings demonstrate that allopolyploidy in itself need not be associated with increased salinity stress tolerance and provide information for using the secondary Gossypium gene pool to breed for improved salt tolerance.

BnaMPK6 is a determinant of quantitative disease resistance against Sclerotinia sclerotiorum in oilseed rape

Sclerotinia sclerotiorum causes a devastating disease in oilseed rape (Brassica napus), resulting in major economic losses. Resistance response of B. napus against S. sclerotiorum exhibits a typical quantitative disease resistance (QDR) characteristic, but the molecular determinants of this QDR are largely unknown. In this study, we isolated a B. napus mitogen-activated protein kinase gene, BnaMPK6, and found that BnaMPK6 expression is highly responsive to infection by S. sclerotiorum and treatment with salicylic acid (SA) or jasmonic acid (JA). Moreover, overexpression (OE) of BnaMPK6 significantly enhances resistance to S. sclerotiorum, whereas RNAi in BnaMPK6 significantly reduces this resistance. These results showed that BnaMPK6 plays an important role in defense to S. sclerotiorum. Furthermore, expression of defense genes associated with SA-, JA- and ethylene (ET)-mediated signaling was investigated in BnaMPK6-RNAi, WT and BnaMPK6-OE plants after S. sclerotiorum infection, and consequently, it was indicated that the activation of ET signaling by BnaMPK6 may play a role in the defense. Further, four BnaMPK6-encoding homologous loci were mapped in the B. napus genome. Using the allele analysis and expression analysis on the four loci, we demonstrated that the locus BnaA03.MPK6 makes an important contribution to QDR against S. sclerotiorum. Our data indicated that BnaMPK6 is a previously unknown determinant of QDR against S. sclerotiorum in B. napus.

Assembly of the 373k gene space of the polyploid sugarcane genome reveals reservoirs of functional diversity in the world's leading biomass crop

BACKGROUND

Sugarcane cultivars are polyploid interspecific hybrids of giant genomes, typically with 10-13 sets of chromosomes from 2 Saccharum species. The ploidy, hybridity, and size of the genome, estimated to have >10 Gb, pose a challenge for sequencing.

RESULTS

Here we present a gene space assembly of SP80-3280, including 373,869 putative genes and their potential regulatory regions. The alignment of single-copy genes in diploid grasses to the putative genes indicates that we could resolve 2-6 (up to 15) putative homo(eo)logs that are 99.1% identical within their coding sequences. Dissimilarities increase in their regulatory regions, and gene promoter analysis shows differences in regulatory elements within gene families that are expressed in a species-specific manner. We exemplify these differences for sucrose synthase (SuSy) and phenylalanine ammonia-lyase (PAL), 2 gene families central to carbon partitioning. SP80-3280 has particular regulatory elements involved in sucrose synthesis not found in the ancestor Saccharum spontaneum. PAL regulatory elements are found in co-expressed genes related to fiber synthesis within gene networks defined during plant growth and maturation. Comparison with sorghum reveals predominantly bi-allelic variations in sugarcane, consistent with the formation of 2 "subgenomes" after their divergence ∼3.8-4.6 million years ago and reveals single-nucleotide variants that may underlie their differences.

CONCLUSIONS

This assembly represents a large step towards a whole-genome assembly of a commercial sugarcane cultivar. It includes a rich diversity of genes and homo(eo)logous resolution for a representative fraction of the gene space, relevant to improve biomass and food production.

Homeologous regulation of Frigida-like genes provides insights on reproductive development and somatic embryogenesis in the allotetraploid Coffea arabica

Coffea arabica is an allotetraploid of high economic importance. C. arabica transcriptome is a combination of the transcripts of two parental genomes (C. eugenioides and C. canephora) that gave rise to the homeologous genes of the species. Previous studies have reported the transcriptional dynamics of C. arabica. In these reports, the ancestry of homeologous genes was identified and the overall regulation of homeologous differential expression (HDE) was explored. One of these genes is part of the FRIGIDA-like family (FRL), which includes the Arabidopsis thaliana flowering-time regulation protein, FRIGIDA (FRI). As nonfunctional FRI proteins give rise to rapid-cycling summer annual ecotypes instead of vernalization-responsive winter-annuals, allelic variation in FRI can modulate flowering time in A. thaliana. Using bioinformatics, genomic analysis, and the evaluation of gene expression of homeologs, we characterized the FRL gene family in C. arabica. Our findings indicate that C. arabica expresses 10 FRL homeologs, and that, throughout flower and fruit development, these genes are differentially transcribed. Strikingly, in addition to confirming the expression of FRL genes during zygotic embryogenesis, we detected FRL expression during direct somatic embryogenesis, a novel finding regarding the FRL gene family. The HDE profile of FRL genes suggests an intertwined homeologous gene regulation. Furthermore, we observed that FLC gene of C. arabica has an expression profile similar to that of CaFRL genes.

Next-generation transcriptome assembly and analysis: Impact of ploidy

Whole genome duplications (WGD) occur widely in plants, but the effects of these events impact all branches of life. WGD events have major evolutionary impacts, often leading to major structural changes within the chromosomes and massive changes in gene expression that facilitate rapid speciation and gene diversification. Even for species that currently have diploid genomes, the impact of ancestral duplication events is still present in the genomes, especially in the context of highly similar gene families that are retained from WGD. However, the impact of these ploidies on various bioinformatics workflows has not been studied well. In this review, we overview biological significance of polyploidy in different organisms. We describe the impact of having polyploid transcriptomes on bioinformatics analyses, especially focusing on transcriptome assembly and transcript quantification. We discuss the benefits of using simulated benchmarking data when we examine the performance of various methods. We also present an example strategy to generate simulated allopolyploid transcriptomes and RNAseq datasets and how these benchmark datasets can be used to assess the performance of transcript assembly and quantification methods. Our benchmarking study shows that all transcriptome assembly methods are affected by having polyploid genomes. Quantification accuracy is also impacted by polyploidy depending on the method. These simulated datasets can be adapted for testing, such as, read mapping, variant calling, and differential expression using biologically realistic conditions.

Evolutionary Conservation and Divergence of Genes Encoding 3-Hydroxy-3-methylglutaryl Coenzyme A Synthase in the Allotetraploid Cotton Species Gossypium hirsutum

Polyploidization is important for the speciation and subsequent evolution of many plant species. Analyses of the duplicated genes produced via polyploidization events may clarify the origin and evolution of gene families. During terpene biosynthesis, 3-hydroxy-3-methylglutaryl coenzyme A synthase (HMGS) functions as a key enzyme in the mevalonate pathway. In this study, we first identified a total of 53 HMGS genes in 23 land plant species, while no HMGS genes were detected in three green algae species. The phylogenetic analysis suggested that plant HMGS genes may have originated from a common ancestral gene before clustering in different branches during the divergence of plant lineages. Then, we detected six HMGS genes in the allotetraploid cotton species (Gossypium hirsutum), which was twice that of the two diploid cotton species (Gossypium raimondii and Gossypium arboreum). The comparison of gene structures and phylogenetic analysis of HMGS genes revealed conserved evolution during polyploidization in Gossypium. Moreover, the expression patterns indicated that six GhHMGS genes were expressed in all tested tissues, with most genes considerably expressed in the roots, and they were responsive to various phytohormone treatments and abiotic stresses. The sequence and expression divergence of duplicated genes in G. hirsutum implied the sub-functionalization of GhHMGS1A and GhHMGS1D as well as GhHMGS3A and GhHMGS3D, whereas it implied the pseudogenization of GhHMGS2A and GhHMGS2D. Collectively, our study unraveled the evolutionary history of HMGS genes in green plants and from diploid to allotetraploid in cotton and illustrated the different evolutionary fates of duplicated HMGS genes resulting from polyploidization.

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.

Homeologous Epistasis in Wheat: The Search for an Immortal Hybrid

Hybridization between related species results in the formation of an allopolyploid with multiple subgenomes. These subgenomes will each contain complete, yet evolutionarily divergent, sets of genes. Like a diploid hybrid, allopolyploids will have two versions, or homeoalleles, for every gene. Partial functional redundancy between homeologous genes should result in a deviation from additivity. These epistatic interactions between homeoalleles are analogous to dominance effects, but are fixed across subgenomes through self pollination. An allopolyploid can be viewed as an immortalized hybrid, with the opportunity to select and fix favorable homeoallelic interactions within inbred varieties. We present a subfunctionalization epistasis model to estimate the degree of functional redundancy between homeoallelic loci and a statistical framework to determine their importance within a population. We provide an example using the homeologous dwarfing genes of allohexaploid wheat, Rht-1, and search for genome-wide patterns indicative of homeoallelic subfunctionalization in a breeding population. Using the IWGSC RefSeq v1.0 sequence, 23,796 homeoallelic gene sets were identified and anchored to the nearest DNA marker to form 10,172 homeologous marker sets. Interaction predictors constructed from products of marker scores were used to fit the homeologous main and interaction effects, as well as estimate whole genome genetic values. Some traits displayed a pattern indicative of homeoallelic subfunctionalization, while other traits showed a less clear pattern or were not affected. Using genomic prediction accuracy to evaluate importance of marker interactions, we show that homeologous interactions explain a portion of the nonadditive genetic signal, but are less important than other epistatic interactions.

Functional trait divergence and trait plasticity confer polyploid advantage in heterogeneous environments

Polyploidy, or whole-genome duplication often with hybridization, is common in eukaryotes and is thought to drive ecological and evolutionary success, especially in plants. The mechanisms of polyploid success in ecologically relevant contexts, however, remain largely unknown. We conducted an extensive test of functional trait divergence and plasticity in conferring polyploid fitness advantage in heterogeneous environments, by growing clonal replicates of a worldwide genotype collection of six allopolyploid and five diploid wild strawberry (Fragaria) taxa in three climatically different common gardens. Among leaf functional traits, we detected divergence in trait means but not plasticities between polyploids and diploids, suggesting that increased genomic redundancy in polyploids does not necessarily translate into greater trait plasticity in response to environmental change. Across the heterogeneous garden environments, however, polyploids exhibited fitness advantage, which was conferred by both trait means and adaptive trait plasticities, supporting a 'jack-and-master' hypothesis for polyploids. Our findings elucidate essential ecological mechanisms underlying polyploid adaptation to heterogeneous environments, and provide an important insight into the prevalence and persistence of polyploid plants.

A Robust Methodology for Assessing Differential Homeolog Contributions to the Transcriptomes of Allopolyploids

Polyploidy has played a pivotal and recurring role in angiosperm evolution. Allotetraploids arise from hybridization between species and possess duplicated gene copies (homeologs) that serve redundant roles immediately after polyploidization. Although polyploidization is a major contributor to plant evolution, it remains poorly understood. We describe an analytical approach for assessing homeolog-specific expression that begins with de novo assembly of parental transcriptomes and effectively (i) reduces redundancy in de novo assemblies, (ii) identifies putative orthologs, (iii) isolates common regions between orthologs, and (iv) assesses homeolog-specific expression using a robust Bayesian Poisson-Gamma model to account for sequence bias when mapping polyploid reads back to parental references. Using this novel methodology, we examine differential homeolog contributions to the transcriptome in the recently formed allopolyploids Tragopogon mirus and T. miscellus (Compositae). Notably, we assess a larger Tragopogon gene set than previous studies of this system. Using carefully identified orthologous regions and filtering biased orthologs, we find in both allopolyploids largely balanced expression with no strong parental bias. These new methods can be used to examine homeolog expression in any tetrapolyploid system without requiring a reference genome.

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.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

Background

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

Results

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

Conclusion

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Electronic supplementary material

The online version of this article (10.1186/s12863-017-0596-1) contains supplementary material, which is available to authorized users.

Characterization of the late embryogenesis abundant (LEA) proteins family and their role in drought stress tolerance in upland cotton

BACKGROUND

Late embryogenesis abundant (LEA) proteins are large groups of hydrophilic proteins with major role in drought and other abiotic stresses tolerance in plants. In-depth study and characterization of LEA protein families have been carried out in other plants, but not in upland cotton. The main aim of this research work was to characterize the late embryogenesis abundant (LEA) protein families and to carry out gene expression analysis to determine their potential role in drought stress tolerance in upland cotton. Increased cotton production in the face of declining precipitation and availability of fresh water for agriculture use is the focus for breeders, cotton being the backbone of textile industries and a cash crop for many countries globally.

RESULTS

In this work, a total of 242, 136 and 142 LEA genes were identified in G. hirsutum, G. arboreum and G. raimondii respectively. The identified genes were classified into eight groups based on their conserved domain and phylogenetic tree analysis. LEA 2 were the most abundant, this could be attributed to their hydrophobic character. Upland cotton LEA genes have fewer introns and are distributed in all chromosomes. Majority of the duplicated LEA genes were segmental. Syntenic analysis showed that greater percentages of LEA genes are conserved. Segmental gene duplication played a key role in the expansion of LEA genes. Sixty three miRNAs were found to target 89 genes, such as miR164, ghr-miR394 among others. Gene ontology analysis revealed that LEA genes are involved in desiccation and defense responses. Almost all the LEA genes in their promoters contained ABRE, MBS, W-Box and TAC-elements, functionally known to be involved in drought stress and other stress responses. Majority of the LEA genes were involved in secretory pathways. Expression profile analysis indicated that most of the LEA genes were highly expressed in drought tolerant cultivars Gossypium tomentosum as opposed to drought susceptible, G. hirsutum. The tolerant genotypes have a greater ability to modulate genes under drought stress than the more susceptible upland cotton cultivars.

CONCLUSION

The finding provides comprehensive information on LEA genes in upland cotton, G. hirsutum and possible function in plants under drought stress.

Gene Co-expression Network Analysis Suggests the Existence of Transcriptional Modules Containing a High Proportion of Transcriptionally Differentiated Homoeologs in Hexaploid Wheat

Genome duplications aid in the formation of novel molecular networks through regulatory differentiation of the duplicated genes and facilitate adaptation to environmental change. Hexaploid wheat, Triticum aestivum, contains three homoeologous chromosome sets, the A-, B-, and D-subgenomes, which evolved through interspecific hybridization and subsequent whole-genome duplication. The divergent expression patterns of the homoeologs in hexaploid wheat suggest that they have undergone transcriptional and/or functional differentiation during wheat evolution. However, the distribution of transcriptionally differentiated homoeologs in gene regulatory networks and their related biological functions in hexaploid wheat are still largely unexplored. Therefore, we retrieved 727 publicly available wheat RNA-sequencing (RNA-seq) datasets from various tissues, developmental stages, and conditions, and identified 10,415 expressed homoeologous triplets. Examining the co-expression modules in the wheat transcriptome, we found that 66% of the expressed homoeologous triplets possess all three homoeologs grouped in the same co-expression modules. Among these, 15 triplets contain co-expressed homoeologs with differential expression levels between homoeoalleles across ≥ 95% of the 727 RNA-seq datasets, suggesting a consistent trend of homoeolog expression bias. In addition, we identified 2,831 differentiated homoeologs that showed gene expression patterns that deviated from those of the other two homoeologs. We found that seven co-expression modules contained a high proportion of such differentiated homoeologs, which accounted for ≥ 20% of the genes in each module. We also found that five of the co-expression modules are abundantly composed of genes involved in biological processes such as chloroplast biogenesis, RNA metabolism, putative defense response, putative posttranscriptional modification, and lipid metabolism, thereby suggesting that, the differentiated homoeologs might highly contribute to these biological functions in the gene network of hexaploid wheat.

BnDGAT1s Function Similarly in Oil Deposition and Are Expressed with Uniform Patterns in Tissues of Brassica napus

As an allotetraploid oilcrop, Brassica napus contains four duplicated Acyl-CoA:diacylglycerol acyltransferase 1 (DGAT1) genes, which catalyze one of the rate-limiting steps in triacylglycerol (TAG) biosynthesis in plants. While all four BnDGAT1s have been expressed functionally in yeast, their expression patterns in different germplasms and tissues and also consequent contribution to seed oil accumulation in planta remain to be elucidated. In this study, the coding regions of the four BnDGAT1s were expressed in an Arabidopsis dgat1 mutant. All four BnDGAT1s showed similar effects on oil content and fatty acid composition, a result which is different from that observed in previous studies of their expression in yeast. Expression patterns of BnDGAT1s were analyzed in developing seeds of 34 B. napus inbred lines and in different tissues of 14 lines. Different expression patterns were observed for the four BnDGAT1s, which suggests that they express independently or randomly in different germplasm sources. Higher expression of BnDGAT1s was correlated with higher seed oil content lines. Tissue-specific analyses showed that the BnDGAT1s were expressed in a uniform pattern in different tissues. Our results suggest that it is important to maintain expression of the four BnDGAT1s for maximum return on oil content.

Empirical evidence of fixed and homeostatic patterns of polyploid advantage in a keystone grass exposed to drought and heat stress

A long-standing hypothesis in evolutionary biology is that polyploid plants have a fitness advantage over diploids in climatically variable or extreme habitats. Here we provide the first empirical evidence that polyploid advantage in these environments is caused by two distinct processes: homeostatic maintenance of reproductive output under elevated abiotic stress, and fixed differences in seed development. In an outdoor climate manipulation experiment using coastal to inland Australian populations of the perennial grass Themeda triandra Forssk., we found that total output of viable seed in drought- and heat-stressed tetraploid plants was over four times higher than in diploids, despite being equal under more favourable growing conditions. Tetraploids also consistently produced heavier seeds with longer hygroscopic awns, traits which increase propagule fitness in extreme environments. These differences add to fitness benefits associated with broader-scale local adaptation of inland T. triandra populations to drought stress. Our study provides evidence that nucleotypic effects of genome size and increased reproductive flexibility can jointly underlie polyploid advantage in plants in stressful environments, and argue that ploidy can be an important criterion for selecting plant populations for use in genetic rescue, restoration and revegetation projects, including in habitats affected by climate change.

Redox and Ionic Homeostasis Regulations against Oxidative, Salinity and Drought Stress in Wheat (A Systems Biology Approach)

Systems biology and omics has provided a comprehensive understanding about the dynamics of the genome, metabolome, transcriptome, and proteome under stress. In wheat, abiotic stresses trigger specific networks of pathways involved in redox and ionic homeostasis as well as osmotic balance. These networks are considerably more complicated than those in model plants, and therefore, counter models are proposed by unifying the approaches of omics and stress systems biology. Furthermore, crosstalk among these pathways is monitored by the regulation and streaming of transcripts and genes. In this review, we discuss systems biology and omics as a promising tool to study responses to oxidative, salinity, and drought stress in wheat.

Segmental allotetraploidy generates extensive homoeologous expression rewiring and phenotypic diversity at the population level in rice

Allopolyploidization, that is, concomitant merging and doubling of two or more divergent genomes in a common nucleus/cytoplasm, is known to instantly alter genomewide transcriptome dynamics, a phenomenon referred to as "transcriptomic shock." However, the immediate effects of transcriptomic alteration in generating phenotypic diversity at the population level remain underinvestigated. Here, we employed the MassARRAY-based Sequenom platform to assess and compare orthologous, allelic and homoeologous gene expression status in two tissues (leaf and root) of a set of randomly chosen individuals from populations of parental rice subspecies (indica and japonica), in vitro "hybrids" (parental mixes), reciprocal F1 hybrids and reciprocal tetraploids at the 5th-selfed generation (S5). We show that hybridization and whole genome duplication (WGD) have opposing effects on allelic and homoeologous expression in the F1 hybrids and tetraploids, respectively. Whereas hybridization exerts strong attenuating effects on allelic expression differences in diploid hybrids, WGD augments the intrinsic parental differences and generates extensive and variable homoeolog content which triggers diversification in expression patterning among the tetraploid plants. Coupled with the vast phenotypic diversity observed among the tetraploid individuals, our results provide experimental evidence in support of the notion that allopolyploidy catalyses rapid phenotypic diversification in higher plants. Our data further suggest that largely stochastic homoeolog content reshuffling rather than alteration in total expression level may be an important feature of evolution in young segmental allopolyploids, which underlies rapid expression diversity at the population level.

Down-regulation of BnDA1, whose gene locus is associated with the seeds weight, improves the seeds weight and organ size in Brassica napus

Brassica napus L. is an important oil crop worldwide and is the main raw material for biofuel. Seed weight and seed size are the main contributors to seed yield. DA1 (DA means big in Chinese) is an ubiquitin receptor and negatively regulates seed size. Down-regulation of AtDA1 in Arabidopsis leads to larger seeds and organs by increasing cell proliferation in integuments. In this study, BnDA1 was down-regulated in B. napus by over expressed of AtDA1^R358K , which is a functional deficiency of DA1 with an arginine-to-lysine mutation at the 358th amino acid. The results showed that the biomass and size of the seeds, cotyledons, leaves, flowers and siliques of transgenic plants all increased significantly. In particular, the 1000 seed weight increased 21.23% and the seed yield per plant increased 13.22% in field condition. The transgenic plants had no negative traits related to yield. The candidate gene association analysis demonstrated that the BnDA1 locus was contributed to the seeds weight. Therefore, our study showed that regulation of DA1 in B. napus can increase the seed yield and biomass, and DA1 is a promising target for crop improvement.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Differential transcript profiling alters regulatory gene expression during the development of Gossypium arboreum, G.stocksii and somatic hybrids

Polyploidy or genome doubling (i.e., the presence of two or more diploid parental genome sets within an organism) are very important in higher plants. Of particular interest are the mechanisms in the new microenvironment of the common nucleus, where doubled regulatory networks interact to generate a viable genetic system capable of regulating growth, development and responses to the environment. To determine the effects of whole genome merging and doubling on the global gene expression architecture of a new polyploid, derived from protoplast fusion of the A₁A₁ genome of Gossypium arboreum and the E₁E₁ genome of Gossypium stocksii, we monitored gene expression through cDNA-AFLP in the somatic hybrids (G. arboreum + G. stocksii). The genomic expression patterns of the somatic hybrids revealed that changes in expression levels mainly involved regulatory genes (31.8% of the gene expression profiles), and the AA and EE genomes contributed equally to genome-wide expression in the newly formed AAEE genome from additivity and dominance effects. These results provide a novel perspective on polyploid gene regulation and hint at the underlying genetic basis of allopolyploid adaption in the new microenvironmental nucleus.

The defence-associated transcriptome of hexaploid wheat displays homoeolog expression and induction bias

Bread wheat (Triticum aestivum L.) is an allopolyploid species containing three ancestral genomes. Therefore, three homoeologous copies exist for the majority of genes in the wheat genome. Whether different homoeologs are differentially expressed (homoeolog expression bias) in response to biotic and abiotic stresses is poorly understood. In this study, we applied a RNA-seq approach to analyse homoeolog-specific global gene expression patterns in wheat during infection by the fungal pathogen Fusarium pseudograminearum, which causes crown rot disease in cereals. To ensure specific detection of homoeologs, we first optimized read alignment methods and validated the results experimentally on genes with known patterns of subgenome-specific expression. Our global analysis identified widespread patterns of differential expression among homoeologs, indicating homoeolog expression bias underpins a large proportion of the wheat transcriptome. In particular, genes differentially expressed in response to Fusarium infection were found to be disproportionately contributed from B and D subgenomes. In addition, we found differences in the degree of responsiveness to pathogen infection among homoeologous genes with B and D homoeologs exhibiting stronger responses to pathogen infection than A genome copies. We call this latter phenomenon as 'homoeolog induction bias'. Understanding how homoeolog expression and induction biases operate may assist the improvement of biotic stress tolerance in wheat and other polyploid crop species.

Phytophthora megakarya and P. palmivora, closely related causal agents of cacao black pod rot, underwent increases in genome sizes and gene numbers by different mechanisms

Phytophthora megakarya (Pmeg) and Phytophthora palmivora (Ppal) are closely related species causing cacao black pod rot. Although Ppal is a cosmopolitan pathogen, cacao is the only known host of economic importance for Pmeg. Pmeg is more virulent on cacao than Ppal. We sequenced and compared the Pmeg and Ppal genomes and identified virulence-related putative gene models (PGeneM) that may be responsible for their differences in host specificities and virulence. Pmeg and Ppal have estimated genome sizes of 126.88 and 151.23 Mb and PGeneM numbers of 42,036 and 44,327, respectively. The evolutionary histories of Pmeg and Ppal appear quite different. Postspeciation, Ppal underwent whole-genome duplication whereas Pmeg has undergone selective increases in PGeneM numbers, likely through accelerated transposable element-driven duplications. Many PGeneMs in both species failed to match transcripts and may represent pseudogenes or cryptic genetic reservoirs. Pmeg appears to have amplified specific gene families, some of which are virulence-related. Analysis of mycelium, zoospore, and in planta transcriptome expression profiles using neural network self-organizing map analysis generated 24 multivariate and nonlinear self-organizing map classes. Many members of the RxLR, necrosis-inducing phytophthora protein, and pectinase genes families were specifically induced in planta . Pmeg displays a diverse virulence-related gene complement similar in size to and potentially of greater diversity than Ppal but it remains likely that the specific functions of the genes determine each species’ unique characteristics as pathogens.

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.