RNA-seq analysis reveals alternative splicing under salt stress in cotton, Gossypium davidsonii

Alternative splicing analysis of stress tolerance to Al and flg22 in Vitis quinquangularis

MAIN CONCLUSION

Alternative splicing of transcriptomes after Al and flg22 treatment for 12 h in response to plant defense of Chinese wild Vitis quinquangularis: genes related to stress resistance and splicing factors were identified in response to Al and flg22 treatment. Alternative splicing (AS) is one of the major post-transcriptional regulation processes that potentially regulates the response to biotic and abiotic stresses in plants. So far, the insight into potential roles of AS in grapevine response to aluminium (Al) and flagellin 22 (flg22) stresses remains poorly understood. We performed transcriptome sequencing of grape leaves before and after Al treatment and flg22 treatment, respectively, to identify AS genes. In this study, a total of 11,805 AS events were identified in Al treatment, of which the skipped exon (SE; 88.72%) type was the most frequent. 9156 AS events were identified under flg22 treatment, of which the SE (88.52%) type was the most frequent. Compared with Al-treated and flg22-treated 0 h, there were 42 and 147 differential alternative splicing (DAS) genes differentially expressed (DASEGs) in Al-treated and flg22-treated 12 h, respectively. Functional analysis showed that DASEGs after Al treatment were mainly enriched in glutathione metabolism pathway; DASEGs after flg22 treatment were enriched in MAPK signaling and plant hormone signal transduction. We further verified seven resistance-related DASEGs with up-regulated expression in Al-treated 12 h, including beta-glucosidase, calcineurin B-like protein, synaptotagmin-3, cysteine synthase and glutathione reductase. Several genes function as leucine-rich repeats receptor-like serine/threonine protein kinase, BRI1 associated receptor kinase 1 and receptor-like protein kinase were also verified by RT-qPCR. We also verified four serine/arginine (SR)-rich proteins SCL30A, SCL28, RS2Z32 and SR45A, which were up-regulated in both Al and flg22 stresses. In conclusion, this study provides an in-depth analysis of the correlation between alternative splicing and grapevine stress tolerance, which helps to identify potential candidate genes for useful traits, provides a theoretical basis for grapevine breeding in plant stress tolerance, and offers new perspectives for understanding grapevine environmental adaptation strategies.

The Regulatory Roles of RNA-Binding Proteins in Plant Salt Stress Response

Salt stress is one of the most prominent abiotic stresses. Behind the intricate adaptive responses of plants to salt stress, the regulation of gene expression assumes a pivotal role. Complementing transcriptional mechanisms, post-transcriptional regulation performed by RNA-binding proteins provides an additional layer of control through sophisticated molecular machinery. RBPs interact with both RNA molecules and protein partners to coordinate RNA metabolism and, thus, fine-tune the expression of salt-responsive genes, enabling plants to rapidly adapt to ionic challenges. This review systematically evaluates the functional roles of RBPs localized in distinct subcellular compartments, including nuclear, cytoplasmic, chloroplastic, and mitochondrial systems, in mediating post-transcriptional regulatory networks under salinity challenges. Specific classes of RBPs are discussed in detail, including glycine-rich RNA-binding proteins (GR-RBPs), serine/arginine-rich splicing factors (SR proteins), zinc finger domain-containing proteins, DEAD-box RNA helicases (DBRHs), KH domain-containing proteins, Pumilio domain-containing proteins (PUMs), pentatricopeptide repeat proteins (PPRs), and RBPs involved in cytoplasmic RNA granule formation. By integrating their subcellular localization and current mechanistic insights, this review concludes by summarizing the current knowledge and highlighting potential future research directions, aiming to inspire further investigations into the complex network of RBPs in modulating plant responses to salt stress and facilitating the development of strategies to enhance plant salt tolerance.

The alternative splicing of HvLHCA4.2 enhances drought tolerance in barley by regulating ROS scavenging and stomatal closure

Alternative splicing (AS) provides flexible strategies for plants to cope with abiotic stress, but how AS regulates drought response in barley (Hordeum vulgare) remains poorly understood. This study investigates the genome-wide AS patterns in drought-sensitive (Baudin) and tolerant barley genotypes (Tadmor, EC_S1) during drought and recovery stages. Baudin showed maximal AS events during drought, while tolerant genotypes activated AS predominantly during recovery. Drought-tolerant EC_S1 exhibited differential AS events enriched in metabolic pathways, amino acid biosynthesis, and chromatin remodeling. WGCNA analysis identified eight modules linked to stomatal dynamics, biomass, and ROS scavenging, with the 'MEmidnightblue' module (related to stoma and ROS regulation) being drought-responsive. A drought-inducible splice variant, LHCA4.2b, was exclusively identified in tolerant genotypes. Silencing HvLHCA4.2a (lhca4.2a) and dual-silencing HvLHCA4.2a + b (lhca4.2a + b) drastically impaired drought tolerance, manifesting as wilting, accelerated water loss, reduced biomass, and elevated electrolyte leakage, coupled with chlorophyll degradation, ROS overaccumulation, and malondialdehyde overproduction. Enhanced stress sensitivity in lhca4.2a + b versus lhca4.2a underscores HvLHCA4.2b's critical role. Additionally, lhca4.2a + b displayed ABA-insensitive stomata with unaltered stomatal conductance under drought or exogenous ABA, implicating LHCA4.2b in ABA signaling. Mechanistically, AS-generated LHCA4.2b isoform enhances PSI-LHCI super-complex stability through strengthened interaction interfaces with LHCA1, thereby improving energy transfer efficiency and reducing ROS generation. This isoform simultaneously coordinates ABA signaling by elevating ABF1/ABF3 transcription factors while suppressing ABI3/ABI4/ABI5 repressors, collectively modulating stomatal closure and ROS homeostasis. Our findings elucidate AS-mediated drought adaptation mechanisms in barley and highlight HvLHCA4.2b as a potential target for breeding drought-tolerant cultivars.

Exploring Alternative Splicing in Response to Salinity: A Tissue-Level Comparative Analysis Using Arabidopsis thaliana Public Transcriptomic Data

Increased soil salinity is a major threat to global agriculture and food security, caused mainly by anthropogenic activities and changing climatic cycles. Plants responses to salinity involve multiple regulatory layers, from transcriptome reprogramming to proteomic and metabolomic changes. Alternative splicing (AS) plays a role in coordinating the response to salinity, yet its extent, tissue, and condition specificity, remain poorly understood aspects. In this study, we used 52 publicly available RNA-seq datasets of salinity treatment to identify differential alternative splicing (DAS) events and genes participating in the response to this stimulus. Our findings reveal that either independently or coordinately, AS can regulate up to 20% of the transcriptome detected in Arabidopsis, with treatment intensity being the most determining factor. Moreover, we show that AS regulation was highly tissue-specific, with roots displaying strong AS-mediated stress responses. Furthermore, cross-stress comparisons showed that roots have a core set of AS-regulated genes associated with stress response and development, with functionally distinct sets of genes when comparing salt with other stresses, while also conserving a relevant condition-specific response. We demonstrate the need to integrate AS analysis to better understand plant adaptation mechanisms and highlight the key role of AS in salinity responses, revealing shared AS regulation between salt, heat, and drought responses.

SMRT sequencing of a full-length transcriptome reveals cold induced alternative splicing in Vitis amurensis root

Alternative splicing enhances diversity at the transcriptional and protein levels that widely involved in plant response to biotic and abiotic stresses. V. amurensis is an extremely cold-tolerant wild grape variety, however, studies on alternative splicing (AS) in amur grape at low temperatures are currently poorly understood. In this study, we analyzed full-length transcriptome and RNA seq data at 0, 2, and 24 h after cold stress in V. amurensis roots. Following quality control and correction, 221,170 high-quality full-length non-concatemer (FLNC) reads were identified. A total of 16,181 loci and 30,733 isoforms were identified. These included 22,868 novel isoforms from annotated genes and 2815 isoforms from 2389 novel genes. Among the distinguished novel isoforms, 673 Long non-coding RNAs (LncRNAs) and 18,164 novel isoforms open reading frame (ORF) region were found. A total of 2958 genes produced 8797 AS events, of which 189 genes were involved in the low-temperature response. Twelve transcription factors show AS during cold treatment and VaMYB108 was selected for initial exploration. Two transcripts, Chr05.63.1 (VaMYB108short) and Chr05.63.2 (VaMYB108normal) of VaMYB108, display up-regulated expression after cold treatment in amur grape roots and are both localized in the nucleus. Only VaMYB108normal exhibits transcriptional activation activity. Overexpression of either VaMYB108short or VaMYB108normal in grape roots leads to increased expression of the other transcript and both increased chilling resistance of amur grape roots. The results improve and supplement the genome annotations and provide insights for further investigation into AS mechanisms during cold stress in V. amurensis.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

RNA-seq analysis reveals transcriptome reprogramming and alternative splicing during early response to salt stress in tomato root

Salt stress is one of the dominant abiotic stress conditions that cause severe damage to plant growth and, in turn, limiting crop productivity. It is therefore crucial to understand the molecular mechanism underlying plant root responses to high salinity as such knowledge will aid in efforts to develop salt-tolerant crops. Alternative splicing (AS) of precursor RNA is one of the important RNA processing steps that regulate gene expression and proteome diversity, and, consequently, many physiological and biochemical processes in plants, including responses to abiotic stresses like salt stress. In the current study, we utilized high-throughput RNA-sequencing to analyze the changes in the transcriptome and characterize AS landscape during the early response of tomato root to salt stress. Under salt stress conditions, 10,588 genes were found to be differentially expressed, including those involved in hormone signaling transduction, amino acid metabolism, and cell cycle regulation. More than 700 transcription factors (TFs), including members of the MYB, bHLH, and WRKY families, potentially regulated tomato root response to salt stress. AS events were found to be greatly enhanced under salt stress, where exon skipping was the most prevalent event. There were 3709 genes identified as differentially alternatively spliced (DAS), the most prominent of which were serine/threonine protein kinase, pentatricopeptide repeat (PPR)-containing protein, E3 ubiquitin-protein ligase. More than 100 DEGs were implicated in splicing and spliceosome assembly, which may regulate salt-responsive AS events in tomato roots. This study uncovers the stimulation of AS during tomato root response to salt stress and provides a valuable resource of salt-responsive genes for future studies to improve tomato salt tolerance.

Identification of Key Genes and Pathways for Anaerobic Germination Tolerance in Rice Using Weighted Gene Co-Expression Network Analysis (WGCNA) in Association with Quantitative Trait Locus (QTL) Mapping

BACKGROUND

Rice is one of the most important food crops in the world, and with the development of direct seeding methods for rice, exposure to anaerobic stress has become a major factor limiting its growth.

RESULTS

In this experiment, we tested the tolerance to anaerobic germination of rice varieties NIP and HD84, and they were used as parents to construct a DH (doubled-haploid) population. The transcriptomes of NIP (highly tolerant) and HD86 (intolerant), and their progeny HR (highly tolerant) and NHR (intolerant) were sequenced from normal and anaerobic environments. The differentially-expressed genes (DEGs) were subjected to GO (Gene ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes), and WGCNA analyses. QTL mapping of the DH population identified tolerance to anaerobic germination-related chromosomal segments. The transcriptome results from 24 samples were combined with the anaerobic stress QTL results for 159 DH population lines to construct a metabolic network to identify key pathways and a gene interaction network to study the key genes. Essential genes were initially subjected to rigorous functional validation, followed by a comprehensive analysis aimed at elucidating their potential utility in domestication and breeding practices, particularly focusing on the exploitation of dominant haplotypes.

CONCLUSION

The results show that pyruvate decarboxylase (PDC) and alcohol dehydrogenase (ADH) are the starting signals of energy metabolism for coleoptile length growth, the auxin transporter EXPA is the determining signal for coleoptile length growth. The pivotal genes Os05g0498700 and Os01g0866100 exert a negative regulatory influence on coleoptile length, ultimately enhancing tolerance to anaerobic germination in rice. Analyses of breeding potential underscore the additional value of Os05g0498700-hyp2 and Os01g0866100-hyp2, highlighting their potential utility in further improving rice through breeding programs. The results of our study will provide a theoretical basis for breeding anaerobic-tolerant rice varieties.

Exploring salt tolerance mechanisms using machine learning for transcriptomic insights: case study in Spartina alterniflora

Salt stress poses a significant threat to global cereal crop production, emphasizing the need for a comprehensive understanding of salt tolerance mechanisms. Accurate functional annotations of differentially expressed genes are crucial for gaining insights into the salt tolerance mechanism. The challenge of predicting gene functions in under-studied species, especially when excluding infrequent GO terms, persists. Therefore, we proposed the use of NetGO 3.0, a machine learning-based annotation method that does not rely on homology information between species, to predict the functions of differentially expressed genes under salt stress. Spartina alterniflora, a halophyte with salt glands, exhibits remarkable salt tolerance, making it an excellent candidate for in-depth transcriptomic analysis. However, current research on the S. alterniflora transcriptome under salt stress is limited. In this study we used S. alterniflora as an example to investigate its transcriptional responses to various salt concentrations, with a focus on understanding its salt tolerance mechanisms. Transcriptomic analysis revealed substantial changes impacting key pathways, such as gene transcription, ion transport, and ROS metabolism. Notably, we identified a member of the SWEET gene family in S. alterniflora, SA_12G129900.m1, showing convergent selection with the rice ortholog SWEET15. Additionally, our genome-wide analyses explored alternative splicing responses to salt stress, providing insights into the parallel functions of alternative splicing and transcriptional regulation in enhancing salt tolerance in S. alterniflora. Surprisingly, there was minimal overlap between differentially expressed and differentially spliced genes following salt exposure. This innovative approach, combining transcriptomic analysis with machine learning-based annotation, avoids the reliance on homology information and facilitates the discovery of unknown gene functions, and is applicable across all sequenced species.

Unraveling the evolutionary origin of the P5CS gene: a story of gene fusion and horizontal transfer

The accumulation of proline in response to the most diverse types of stress is a widespread defense mechanism. In prokaryotes, fungi, and certain unicellular eukaryotes (green algae), the first two reactions of proline biosynthesis occur through two distinct enzymes, γ-glutamyl kinase (GK E.C. 2.7.2.11) and γ-glutamyl phosphate reductase (GPR E.C. 1.2.1.41), encoded by two different genes, ProB and ProA, respectively. Plants, animals, and a few unicellular eukaryotes carry out these reactions through a single bifunctional enzyme, the Δ1-pyrroline-5-carboxylate synthase (P5CS), which has the GK and GPR domains fused. To better understand the origin and diversification of the P5CS gene, we use a robust phylogenetic approach with a broad sampling of the P5CS, ProB and ProA genes, including species from all three domains of life. Our results suggest that the collected P5CS genes have arisen from a single fusion event between the ProA and ProB gene paralogs. A peculiar fusion event occurred in an ancestral eukaryotic lineage and was spread to other lineages through horizontal gene transfer. As for the diversification of this gene family, the phylogeny of the P5CS gene in plants shows that there have been multiple independent processes of duplication and loss of this gene, with the duplications being related to old polyploidy events.

DEAD-Box RNA Helicase Family in Physic Nut (Jatropha curcas L.): Structural Characterization and Response to Salinity

Helicases, motor proteins present in both prokaryotes and eukaryotes, play a direct role in various steps of RNA metabolism. Specifically, SF2 RNA helicases, a subset of the DEAD-box family, are essential players in plant developmental processes and responses to biotic and abiotic stresses. Despite this, information on this family in the physic nut (Jatropha curcas L.) remains limited, spanning from structural patterns to stress responses. We identified 79 genes encoding DEAD-box RNA helicases (JcDHX) in the J. curcas genome. These genes were further categorized into three subfamilies: DEAD (42 genes), DEAH (30 genes), and DExH/D (seven genes). Characterization of the encoded proteins revealed a remarkable diversity, with observed patterns in domains, motifs, and exon-intron structures suggesting that the DEAH and DExH/D subfamilies in J. curcas likely contribute to the overall versatility of the family. Three-dimensional modeling of the candidates showed characteristic hallmarks, highlighting the expected functional performance of these enzymes. The promoter regions of the JcDHX genes revealed potential cis-elements such as Dof-type, BBR-BPC, and AP2-ERF, indicating their potential involvement in the response to abiotic stresses. Analysis of RNA-Seq data from the roots of physic nut accessions exposed to 150 mM of NaCl for 3 h showed most of the JcDHX candidates repressed. The protein-protein interaction network indicated that JcDHX proteins occupy central positions, connecting events associated with RNA metabolism. Quantitative PCR analysis validated the expression of nine DEAD-box RNA helicase transcripts, showing significant associations with key components of the stress response, including RNA turnover, ribosome biogenesis, DNA repair, clathrin-mediated vesicular transport, phosphatidyl 3,5-inositol synthesis, and mitochondrial translation. Furthermore, the induced expression of one transcript (JcDHX44) was confirmed, suggesting that it is a potential candidate for future functional analyses to better understand its role in salinity stress tolerance. This study represents the first global report on the DEAD-box family of RNA helicases in physic nuts and displays structural characteristics compatible with their functions, likely serving as a critical component of the plant's response pathways.

Responsive Alternative Splicing Events of Opisthopappus Species against Salt Stress

Salt stress profoundly affects plant growth, prompting intricate molecular responses, such as alternative splicing (AS), for environmental adaptation. However, the response of AS events to salt stress in Opisthopappus (Opisthopappus taihangensis and Opisthopappus longilobus) remains unclear, which is a Taihang Mountain cliff-dwelling species. Using RNA-seq data, differentially expressed genes (DEGs) were identified under time and concentration gradients of salt stress. Two types of AS, skipped exon (SE) and mutually exclusive exons (MXE), were found. Differentially alternative splicing (DAS) genes in both species were significantly enriched in "protein phosphorylation", "starch and sucrose metabolism", and "plant hormone signal transduction" pathways. Meanwhile, distinct GO terms and KEGG pathways of DAS occurred between two species. Only a small subset of DAS genes overlapped with DEGs under salt stress. Although both species likely adopted protein phosphorylation to enhance salt stress tolerance, they exhibited distinct responses. The results indicated that the salt stress mechanisms of both Opisthopappus species exhibited similarities and differences in response to salt stress, which suggested that adaptive divergence might have occurred between them. This study initially provides a comprehensive description of salt responsive AS events in Opisthopappus and conveys some insights into the molecular mechanisms behind species tolerance on the Taihang Mountains.

Alternative Splicing under Cold Stress in Paper Mulberry

The paper mulberry is a commonly found tree species with a long history of cultivation. It also serves as a crucial case study for understanding how woody plants adapt to low temperatures. Under cold treatment, we observed a substantial number of alternative splicing (AS) genes, showcasing the intricate landscape of AS events. We have detected all seven types of AS events, with the alternative 3' splice site (A3) having the most. We observed that many genes that underwent differential AS were significantly enriched in starch and sucrose metabolism and circadian rhythm pathways. Moreover, a considerable proportion of differentially spliced genes (DSGs) also showed differential expression, with 20.38% and 25.65% under 12 h and 24 h cold treatments, respectively. This suggests a coordinated regulation between gene AS and expression, playing a pivotal role in the paper mulberry's adaptation to cold stress. We further investigated the regulatory mechanisms of AS, identifying 41 serine/arginine-rich (SR) splicing factors, among which 11 showed differential expression under cold treatment, while 29 underwent alternative splicing. Additionally, genes undergoing AS displayed significantly higher DNA methylation levels under cold stress, while normal splicing (non-AS) genes exhibited relatively lower methylation levels. These findings suggest that methylation may play an important role in governing gene AS. Finally, our research will provide useful information on the role of AS in the cold acclimation tolerance of the paper mulberry.

Alternative Splicing Variation: Accessing and Exploiting in Crop Improvement Programs

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology, bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization, and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside a gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

Alternative Splicing during Fiber Development in G. hirsutum

Cotton is a valuable cash crop in many countries. Cotton fiber is a trichome that develops from a single epidermal cell and serves as an excellent model for understanding cell differentiation and other life processes. Alternative splicing (AS) of genes is a common post-transcriptional regulatory process in plants that is essential for plant growth and development. The process of AS during cotton fiber formation, on the other hand, is mainly unknown. A substantial number of multi-exon genes were discovered to be alternatively spliced during cotton fiber formation in this study, accounting for 23.31% of the total number of genes in Gossypium hirsutum. Retention intron (RI) is not necessarily the most common AS type, indicating that AS genes and processes during fiber development are very temporal and tissue-specific. When compared to fiber samples, AS is more prevalent at the fiber initiation stages and in the ovule, indicating that development stages and tissues use different AS strategies. Genes involved in fiber development have gone through stage-specific AS, demonstrating that AS regulates cotton fiber development. Furthermore, AS can be regulated by trans-regulation elements such as splicing factor and cis-regulation elements such as gene length, exon numbers, and GC content, particularly at exon-intron junction sites. Our findings also suggest that increased DNA methylation may aid in the efficiency of AS, and that gene body methylation is key in AS control. Finally, our research will provide useful information about the roles of AS during the cotton fiber development process.

Mapping intron retention events contributing to complex traits using splice quantitative trait locus

BACKGROUND

Alternative splicing (AS) of mRNA plays an important roles in transcriptome diversity, involving regulation of plant growth and stress response. Understanding the variation of AS events underlying GWAS loci in a crop population can provide insight into the molecular mechanisms of complex agronomic traits. To date, genome-wide association studies relating AS events to agronomic traits have rarely been conducted at the population level in crops.

RESULTS

Here, a pipeline was constructed to identify candidate AS events related to complex traits. Firstly, ovule transcriptome data were used to characterize intron retention (IR), the predominant type of AS in plants, on a genome-wide scale. This was done in a natural population consisting of 279 upland cotton lines. Secondly, splice quantitative trait locus (sQTL) analysis was carried out, which yielded a total of 2295 sQTLs involving 1607 genes. Of these, 14.25% (n = 427) were cis-regulatory loci. Integration with expression quantitative trait loci (eQTL) revealed that 53 (21.4%) cis-sGenes were regulated by both cis-sQTLs and cis-eQTLs. Finally, co-localization analysis integrated with GWAS loci in this population showed 32 cis-QTLs to be co-located with genetic regulatory loci related to fiber yield and quality traits, indicating that sQTLs are likely to participate in regulating cotton fiber yield and quality. An in-depth evaluation confirmed that differences in the IR rates of sQTL-regulated candidate genes such as GhLRRK1 and GhGC1 are associated with lint percentage (LP), which has potential in breeding applications.

CONCLUSION

This study provides a clue that AS of mRNA has an impact on crop yield, along with functional sQTLs are new genetic resources for cotton precision breeding.

Deciphering the role of alternative splicing as modulators of defense response in the MYMIV- Vigna mungo pathosystem

Alternative splicing (AS) is a crucial regulatory mechanism that impacts transcriptome and proteome complexity under stressful situations. Although its role in abiotic stresses is somewhat understood, our understanding of the mechanistic regulation of pre-mRNA splicing in plant-pathogen interaction is meagre. To comprehend this unexplored immune reprogramming mechanism, transcriptome profiles of Mungbean Yellow Mosaic India Virus (MYMIV)-resistant and susceptible Vigna mungo genotypes were analysed for AS genes that may underlie the resistance mechanism. Results revealed a repertoire of AS-isoforms accumulated during pathogenic infestation, with intron retention being the most common AS mechanism. Identification of 688 differential alternatively spliced (DAS) genes in the resistant host elucidates its robust antiviral response, whereas 322 DAS genes were identified in the susceptible host. Enrichment analyses confirmed DAS transcripts pertaining to stress, signalling, and immune system pathways have undergone maximal perturbations. Additionally, a strong regulation of the splicing factors has been observed both at transcriptional and post-transcriptional levels. qPCR validation of candidate DAS transcripts with induced expression upon MYMIV-infection demonstrated a competent immune response in the resistant background. The AS-impacted genes resulted either in partial/complete loss of functional domains or altered sensitivity to miRNA-mediated gene silencing. A complex regulatory module, miR7517-ATAF2, has been identified in an aberrantly spliced ATAF2 isoform that exposes an intronic miR7517 binding site, thereby suppressing the negative regulator to enhance defense reaction. The present study establishes AS as a non-canonical immune reprogramming mechanism that operates in parallel, thereby offering an alternative strategy for developing yellow mosaic-resistant V. mungo cultivars. This article is protected by copyright. All rights reserved.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses

BACKGROUND

Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions.

RESULTS

We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes.

CONCLUSION

Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

Transcriptome, proteome and functional characterization reveals salt stress tolerance mechanisms in upland cotton (Gossypium hirsutum L.)

Uncovering the underlying mechanism of salt tolerance is important to breed cotton varieties with improved salt tolerance. In this study, transcriptome and proteome sequencing were performed on upland cotton (Gossypium hirsutum L.) variety under salt stress, and integrated analysis was carried out to exploit salt-tolerance genes in cotton. Enrichment analysis using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was performed on differentially expressed genes (DEGs) obtained from transcriptome and proteome sequencing. GO enrichment was carried out mainly in the cell membrane, organelle, cellular process, metabolic process, and stress response. The expression of 23,981 genes was changed in physiological and biochemical processes such as cell metabolism. The metabolic pathways obtained by KEGG enrichment included glycerolipid metabolism, sesquiterpene and triterpenoid biosynthesis, flavonoid production, and plant hormone signal transduction. Combined transcriptome and proteome analysis to screen and annotate DEGs yielded 24 candidate genes with significant differential expression. The quantitative real-time polymerase chain reaction (qRT-PCR) validation of the candidate genes showed that two genes (Gh_D11G0978 and Gh_D10G0907) responded significantly to the induction of NaCl, and these two genes were further selected as target genes for gene cloning and functional validation through virus-induced gene silencing (VIGS). The silenced plants exhibited early wilting with a greater degree of salt damage under salt treatment. Moreover, they showed higher levels of reactive oxygen species (ROS) than the control. Therefore, we can infer that these two genes have a pivotal role in the response to salt stress in upland cotton. The findings in this research will facilitate the breeding of salt tolerance cotton varieties that can be grown on saline alkaline lands.

Comparative transcriptome responses of leaf and root tissues to salt stress in wheat strains with different salinity tolerances

Background: Salinity stress is a major adverse environmental factor that can limit crop yield and restrict normal land use. The selection of salt-tolerant strains and elucidation of the underlying mechanisms by plant breeding scientists are urgently needed to increase agricultural production in arid and semi-arid regions. Results: In this study, we selected the salt-tolerant wheat (Triticum aestivum) strain ST9644 as a model to study differences in expression patterns between salt-tolerant and salt-sensitive strains. High-throughput RNA sequencing resulted in more than 359.10 Gb of clean data from 54 samples, with an average of 6.65 Gb per sample. Compared to the IWGSC reference annotation, we identified 50,096 new genes, 32,923 of which have functional annotations. Comparisons of abundances between salt-tolerant and salt-sensitive strains revealed 3,755, 5,504, and 4,344 genes that were differentially expressed at 0, 6, and 24 h, respectively, in root tissue under salt stress. KEGG pathway analysis of these genes showed that they were enriched for phenylpropanoid biosynthesis (ko00940), cysteine and methionine metabolism (ko00270), and glutathione metabolism (ko00480). We also applied weighted gene co-expression network analysis (WGCNA) analysis to determine the time course of root tissue response to salt stress and found that the acute response lasts >6 h and ends before 12 h. We also identified key alternative splicing factors showing different splicing patterns in salt-sensitive and salt-tolerant strains; however, only few of them were differentially expressed in the two groups. Conclusion: Our results offer a better understanding of wheat salt tolerance and improve wheat breeding.

Comparative Analysis of Alternative Splicing in Two Contrasting Apple Cultivars Defense against Alternaria alternata Apple Pathotype Infection

Alternaria blotch disease, caused by the Alternaria alternata apple pathotype (A. alternata AP), is one of the most serious fungal diseases in apples. Alternative splicing (AS), one of the pivotal post-transcriptional regulatory mechanisms, plays essential roles in various disease resistance responses. Here, we performed RNA-Seq for two apple cultivars (resistant cultivar 'Jonathan' (J) and susceptible cultivar 'Starking Delicious' (SD)) infected by A. alternata AP to further investigate their AS divergence. In total, 1454, 1780, 1367 and 1698 specifically regulated differential alternative splicing (DAS) events were detected in J36, J72, SD36 and SD72 groups, respectively. Retained intron (RI) was the dominant AS pattern. Conformably, 642, 764, 585 and 742 uniquely regulated differentially spliced genes (DSGs) were found during A. alternata AP infection. Comparative analysis of AS genes in differential splicing and expression levels suggested that only a small proportion of DSGs overlapped with differentially expressed genes (DEGs). Gene ontology (GO) enrichment analysis demonstrated that the DSGs were significantly enriched at multiple levels of gene expression regulation. Briefly, the specific AS was triggered in apple defense against A. alternata AP. Therefore, this study facilitates our understanding on the roles of AS regulation in response to A. alternata AP infection in apples.

Full-length fruit transcriptomes of southern highbush (Vaccinium sp.) and rabbiteye (V. virgatum Ait.) blueberry

Background

Blueberries (Vaccinium sp.) are native to North America and breeding efforts to improve blueberry fruit quality are focused on improving traits such as increased firmness, enhanced flavor and greater shelf-life. Such efforts require additional genomic resources, especially in southern highbush and rabbiteye blueberries.

Results

We generated the first full-length fruit transcriptome for the southern highbush and rabbiteye blueberry using the cultivars, Suziblue and Powderblue, respectively. The transcriptome was generated using the Pacific Biosciences single-molecule long-read isoform sequencing platform with cDNA pooled from seven stages during fruit development and postharvest storage. Raw reads were processed through the Isoseq pipeline and full-length transcripts were mapped to the ‘Draper’ genome with unmapped reads collapsed using Cogent. Finally, we identified 16,299 and 15,882 non-redundant transcripts in ‘Suziblue’ and ‘Powderblue’ respectively by combining the reads mapped to Northern Highbush blueberry ‘Draper’ genome and Cogent analysis. In both cultivars, > 80% of sequences were longer than 1,000 nt, with the median transcript length around 1,700 nt. Functionally annotated transcripts using Blast2GO were > 92% in both ‘Suziblue’ and ‘Powderblue’ with overall equal distribution of gene ontology (GO) terms in the two cultivars. Analyses of alternative splicing events indicated that around 40% non-redundant sequences exhibited more than one isoform. Additionally, long non-coding RNAs were predicted to represent 5.6% and 7% of the transcriptomes in ‘Suziblue’ and ‘Powderblue’, respectively. Fruit ripening is regulated by several hormone-related genes and transcription factors. Among transcripts associated with phytohormone metabolism/signaling, the highest number of transcripts were related to abscisic acid (ABA) and auxin metabolism followed by those for brassinosteroid, jasmonic acid and ethylene metabolism. Among transcription factor-associated transcripts, those belonging to ripening-related APETALA2/ethylene-responsive element-binding factor (AP2/ERF), NAC (NAM, ATAF1/2 and CUC2), leucine zipper (HB-zip), basic helix-loop-helix (bHLH), MYB (v-MYB, discovered in avian myeloblastosis virus genome) and MADS-Box gene families, were abundant.

Further we measured three fruit ripening quality traits and indicators [ABA, and anthocyanin concentration, and texture] during fruit development and ripening. ABA concentration increased during the initial stages of fruit ripening and then declined at the Ripe stage, whereas anthocyanin content increased during the final stages of fruit ripening in both cultivars. Fruit firmness declined during ripening in ‘Powderblue’. Genes associated with the above parameters were identified using the full-length transcriptome. Transcript abundance patterns of these genes were consistent with changes in the fruit ripening and quality-related characteristics.

Conclusions

A full-length, well-annotated fruit transcriptome was generated for two blueberry species commonly cultivated in the southeastern United States. The robustness of the transcriptome was verified by the identification and expression analyses of multiple fruit ripening and quality–regulating genes. The full-length transcriptome is a valuable addition to the blueberry genomic resources and will aid in further improving the annotation. It will also provide a useful resource for the investigation of molecular aspects of ripening and postharvest processes.

Cytokinin and abiotic stress tolerance -What has been accomplished and the way forward?

More than a half-century has passed since it was discovered that phytohormone cytokinin (CK) is essential to drive cytokinesis and proliferation in plant tissue culture. Thereafter, cytokinin has emerged as the primary regulator of the plant cell cycle and numerous developmental processes. Lately, a growing body of evidence suggests that cytokinin has a role in mitigating both abiotic and biotic stress. Cytokinin is essential to defend plants against excessive light exposure and a unique kind of abiotic stress generated by an altered photoperiod. Secondly, cytokinin also exhibits multi-stress resilience under changing environments. Furthermore, cytokinin homeostasis is also affected by several forms of stress. Therefore, the diverse roles of cytokinin in reaction to stress, as well as its interactions with other hormones, are discussed in detail. When it comes to agriculture, understanding the functioning processes of cytokinins under changing environmental conditions can assist in utilizing the phytohormone, to increase productivity. Through this review, we briefly describe the biological role of cytokinin in enhancing the performance of plants growth under abiotic challenges as well as the probable mechanisms underpinning cytokinin-induced stress tolerance. In addition, the article lays forth a strategy for using biotechnological tools to modify genes in the cytokinin pathway to engineer abiotic stress tolerance in plants. The information presented here will assist in better understanding the function of cytokinin in plants and their effective investigation in the cropping system.

Alternative Splicing and Its Roles in Plant Metabolism

Plant metabolism, including primary metabolism such as tricarboxylic acid cycle, glycolysis, shikimate and amino acid pathways as well as specialized metabolism such as biosynthesis of phenolics, alkaloids and saponins, contributes to plant survival, growth, development and interactions with the environment. To this end, these metabolic processes are tightly and finely regulated transcriptionally, post-transcriptionally, translationally and post-translationally in response to different growth and developmental stages as well as the constantly changing environment. In this review, we summarize and describe the current knowledge of the regulation of plant metabolism by alternative splicing, a post-transcriptional regulatory mechanism that generates multiple protein isoforms from a single gene by using alternative splice sites during splicing. Numerous genes in plant metabolism have been shown to be alternatively spliced under different developmental stages and stress conditions. In particular, alternative splicing serves as a regulatory mechanism to fine-tune plant metabolism by altering biochemical activities, interaction and subcellular localization of proteins encoded by splice isoforms of various genes.

Hypoxia-Induced Aquaporins and Regulation of Redox Homeostasis by a Trans-Plasma Membrane Electron Transport System in Maize Roots

In plants, flooding-induced oxygen deficiency causes severe stress, leading to growth reduction and yield loss. It is therefore important to understand the molecular mechanisms for adaptation to hypoxia. Aquaporins at the plasma membrane play a crucial role in water uptake. However, their role during hypoxia and membrane redox changes is still not fully understood. The influence of 24 h hypoxia induction on hydroponically grown maize (Zea mays L.) was investigated using an oil-based setup. Analyses of physiological parameters revealed typical flooding symptoms such as increased ethylene and H₂O₂ levels, an increased alcohol dehydrogenase activity, and an increased redox activity at the plasma membrane along with decreased oxygen of the medium. Transcriptomic analysis and shotgun proteomics of plasma membranes and soluble fractions were performed to determine alterations in maize roots. RNA-sequencing data confirmed the upregulation of genes involved in anaerobic metabolism, biosynthesis of the phytohormone ethylene, and its receptors. Transcripts of several antioxidative systems and other oxidoreductases were regulated. Mass spectrometry analysis of the plasma membrane proteome revealed alterations in redox systems and an increased abundance of aquaporins. Here, we discuss the importance of plasma membrane aquaporins and redox systems in hypoxia stress response, including the regulation of plant growth and redox homeostasis.

Analysis of Genomic Alternative Splicing Patterns in Rat under Heat Stress Based on RNA-Seq Data

Heat stress is one of the most severe challenges faced in livestock production in summer. Alternative splicing as an important post-transcriptional regulation is rarely studied in heat-stressed animals. Here, we performed and analyzed RNA-sequencing assays on the liver of Sprague-Dawley rats in control (22 °C, n = 5) and heat stress (4 °C for 120 min, H120; n = 5) groups, resulting in the identification of 636 differentially expressed genes. Identification analysis of the alternative splicing events revealed that heat stress-induced alternative splicing events increased by 20.18%. Compared with other types of alternative splicing events, the alternative start increased the most (43.40%) after heat stress. Twenty-eight genes were differentially alternatively spliced (DAS) between the control and H120 groups, among which Acly, Hnrnpd and mir3064 were also differentially expressed. For DAS genes, Srebf1, Shc1, Srsf5 and Ensa were associated with insulin, while Cast, Srebf1, Tmem33, Tor1aip2, Slc39a7 and Sqstm1 were enriched in the composition of the endoplasmic reticulum. In summary, our study conducts a comprehensive profile of alternative splicing in heat-stressed rats, indicating that alternative splicing is one of the molecular mechanisms of heat stress response in mammals and providing reference data for research on heat tolerance in mammalian livestock.

Early responses to salt stress in quinoa genotypes with opposite behavior

Soil salinity is among the major abiotic stresses that plants must cope with, mainly in arid and semiarid regions. The tolerance to high salinity is an important agronomic trait to sustain food production. Quinoa is a halophytic annual pseudo-cereal species with high nutritional value that can secrete salt out of young leaves in external non-glandular cells called epidermal bladder cells (EBC). Previous work showed high salt tolerance, but low EBC density was associated with an improved response in the early phases of salinity stress, mediated by tissue-tolerance traits mainly in roots. We compared the transcript profiling of two quinoa genotypes with contrasting salt tolerance patterning to identify the candidate genes involved in the differentially early response among genotypes. The transcriptome profiling, supported by in vitro physiological analyses, provided insights into the early-stage molecular mechanisms, both at the shoot and root level, based on the sensitive/tolerance traits. Results showed the presence of numerous differentially expressed genes among genotypes, tissues, and treatments, with genes involved in hormonal and stress response upregulated mainly in the sensitive genotype, suggesting that tolerance may be correlated to restricted changes in gene expression, at least after a short salt stress. These data, showing constitutive differences between the two genotypes, represent a solid basis for further studies to characterize the salt tolerance traits. Additionally, new information provided by this work might be useful for the development of plant breeding or genome engineering programs in quinoa.

Gene regulation at transcriptional and post-transcriptional levels to combat salt stress in plants

Soil salinity is a major challenge that will be faced more and more by human population in the near future. Higher salt concentrations in the soil limit the growth and production of crops, which poses serious threats to global food production. Various plant breeding approaches have been followed in the past which are reported to reduce the effect of salt stress by inducing the level of protective metabolites like osmolytes and antioxidants. Conventional breeding approaches are time-consuming and not cost-effective. In recent times, genetic engineering has been largely followed to confer salt tolerance through introgressions of single transgenes or stacking multiple transgenes. However, most of such works are limited only at the laboratory level and field trials are still awaited to prove the long-term efficacy of such transgenics. In this review, we attempt to present a broad overview of the current strategies undertaken to develop halophytic and salt-tolerant crops. The salt-induced damages in the plants are highlighted, followed by representing the novel traits, associated with salt stress, which can be used for engineering salt tolerance in glycophytic crops. Additionally, the role of transcriptional and epigenetic regulation in plants for amelioration of salt-induced damages has been reviewed. The role of post-transcriptional mechanisms such as microRNA regulation, genome editing and alternative splicing, during salt stress, and their implications in the development of salt-tolerant crops are also discussed. Finally, we present a short overview about the role of ion transporters and rhizobacteria in the engineering of salt tolerance in crop species.

The Landscape of Alternative Splicing Regulating Potassium Use Efficiency in Nicotiana tabacum

Alternative splicing (AS) occurs extensively in eukaryotes as an essential mechanism for regulating transcriptome complexity and diversity, but the AS landscape regulating potassium (K) use efficiency in plants is unclear. In this study, we performed high-throughput transcriptome sequencing of roots and shoots from allopolyploid Nicotiana tabacum under K⁺ deficiency. Preliminary physiological analysis showed that root system architecture was dramatically changed due to potassium deficiency and that IAA content was significantly reduced in root and shoot. AS analysis showed that a total of 28,179 genes exhibited 54,457 AS events, and 1,510 and 1,732 differentially alternatively spliced (DAS) events were identified in shoots and roots under low K⁺ stress. Nevertheless, only 120 DAS events occurred in both shoots and roots, implying that most DAS events were tissue-specific. Both in shoot and the root, the proportion of DAS genes in differentially expressed (DE) genes equaled that in non-DE genes, which indicated that AS might play a unique regulatory role in response to low potassium. Gene ontology analysis further indicated that transcription regulation and AS modulation worked independently in response to low K⁺ stress in tobacco, as their target biological processes were different. Totally 45 DAS transcription factors (TFs) were found, which were involved in 18 TF families. Five Auxin response factor (ARF) TFs were significantly DAS in root, suggesting that response to auxin was probably subject to AS regulation in the tobacco root. Our study shows that AS variation occurs extensively and has a particular regulatory mechanism under K⁺ deficiency in tobacco. The study also links changes in root system architecture with the changes in AS of ARF TFs, which implied the functional significance of these AS events for root growth and architecture.

PacBio and Illumina RNA Sequencing Identify Alternative Splicing Events in Response to Cold Stress in Two Poplar Species

In eukaryotes, alternative splicing (AS) is a crucial regulatory mechanism that modulates mRNA diversity and stability. The contribution of AS to stress is known in many species related to stress, but the posttranscriptional mechanism in poplar under cold stress is still unclear. Recent studies have utilized the advantages of single molecular real-time (SMRT) sequencing technology from Pacific Bioscience (PacBio) to identify full-length transcripts. We, therefore, used a combination of single-molecule long-read sequencing and Illumina RNA sequencing (RNA-Seq) for a global analysis of AS in two poplar species (Populus trichocarpa and P. ussuriensis) under cold stress. We further identified 1,261 AS events in P. trichocarpa and 2,101 in P. ussuriensis among which intron retention, with a frequency of more than 30%, was the most prominent type under cold stress. RNA-Seq data analysis and annotation revealed the importance of calcium, abscisic acid, and reactive oxygen species signaling in cold stress response. Besides, the low temperature rapidly induced multiple splicing factors, transcription factors, and differentially expressed genes through AS. In P. ussuriensis, there was a rapid occurrence of AS events, which provided a new insight into the complexity and regulation of AS during cold stress response in different poplar species for the first time.

Full-length transcriptome analysis and identification of transcript structures in Eimeria necatrix from different developmental stages by single-molecule real-time sequencing

Background

Eimeria necatrix is one of the most pathogenic parasites, causing high mortality in chickens. Although its genome sequence has been published, the sequences and complete structures of its mRNA transcripts remain unclear, limiting exploration of novel biomarkers, drug targets and genetic functions in E. necatrix.

Methods

Second-generation merozoites (MZ-2) of E. necatrix were collected using Percoll density gradients, and high-quality RNA was extracted from them. Single-molecule real-time (SMRT) sequencing and Illumina sequencing were combined to generate the transcripts of MZ-2. Combined with the SMRT sequencing data of sporozoites (SZ) collected in our previous study, the transcriptome and transcript structures of E. necatrix were studied.

Results

SMRT sequencing yielded 21,923 consensus isoforms in MZ-2. A total of 17,151 novel isoforms of known genes and 3918 isoforms of novel genes were successfully identified. We also identified 2752 (SZ) and 3255 (MZ-2) alternative splicing (AS) events, 1705 (SZ) and 1874 (MZ-2) genes with alternative polyadenylation (APA) sites, 4019 (SZ) and 2588 (MZ-2) fusion transcripts, 159 (SZ) and 84 (MZ-2) putative transcription factors (TFs) and 3581 (SZ) and 2039 (MZ-2) long non-coding RNAs (lncRNAs). To validate fusion transcripts, reverse transcription-PCR was performed on 16 candidates, with an accuracy reaching up to 87.5%. Sanger sequencing of the PCR products further confirmed the authenticity of chimeric transcripts. Comparative analysis of transcript structures revealed a total of 3710 consensus isoforms, 815 AS events, 1139 genes with APA sites, 20 putative TFs and 352 lncRNAs in both SZ and MZ-2.

Conclusions

We obtained many long-read isoforms in E. necatrix SZ and MZ-2, from which a series of lncRNAs, AS events, APA events and fusion transcripts were identified. Information on TFs will improve understanding of transcriptional regulation, and fusion event data will greatly improve draft versions of gene models in E. necatrix. This information offers insights into the mechanisms governing the development of E. necatrix and will aid in the development of novel strategies for coccidiosis control.

Graphical Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-05015-7.

Comparative Proteomic Analysis of Tolerant and Sensitive Varieties Reveals That Phenylpropanoid Biosynthesis Contributes to Salt Tolerance in Mulberry

Mulberry, an important woody tree, has strong tolerance to environmental stresses, including salinity, drought, and heavy metal stress. However, the current research on mulberry resistance focuses mainly on the selection of resistant resources and the determination of physiological indicators. In order to clarify the molecular mechanism of salt tolerance in mulberry, the physiological changes and proteomic profiles were comprehensively analyzed in salt-tolerant (Jisang3) and salt-sensitive (Guisangyou12) mulberry varieties. After salt treatment, the malondialdehyde (MDA) content and proline content were significantly increased compared to control, and the MDA and proline content in G12 was significantly lower than in Jisang3 under salt stress. The calcium content was significantly reduced in the salt-sensitive mulberry varieties Guisangyou12 (G12), while sodium content was significantly increased in both mulberry varieties. Although the Jisang3 is salt-tolerant, salt stress caused more reductions of photosynthetic rate in Jisang3 than Guisangyou12. Using tandem mass tags (TMT)-based proteomics, the changes of mulberry proteome levels were analyzed in salt-tolerant and salt-sensitive mulberry varieties under salt stress. Combined with GO and KEGG databases, the differentially expressed proteins were significantly enriched in the GO terms of amino acid transport and metabolism and posttranslational modification, protein turnover up-classified in Guisangyou12 while down-classified in Jisang3. Through the comparison of proteomic level, we identified the phenylpropanoid biosynthesis may play an important role in salt tolerance of mulberry. We clarified the molecular mechanism of mulberry salt tolerance, which is of great significance for the selection of excellent candidate genes for saline-alkali soil management and mulberry stress resistance genetic engineering.

Decoding co-/post-transcriptional complexities of plant transcriptomes and epitranscriptome using next-generation sequencing technologies

Next-generation sequencing (NGS) technologies - Illumina RNA-seq, Pacific Biosciences isoform sequencing (PacBio Iso-seq), and Oxford Nanopore direct RNA sequencing (DRS) - have revealed the complexity of plant transcriptomes and their regulation at the co-/post-transcriptional level. Global analysis of mature mRNAs, transcripts from nuclear run-on assays, and nascent chromatin-bound mRNAs using short as well as full-length and single-molecule DRS reads have uncovered potential roles of different forms of RNA polymerase II during the transcription process, and the extent of co-transcriptional pre-mRNA splicing and polyadenylation. These tools have also allowed mapping of transcriptome-wide start sites in cap-containing RNAs, poly(A) site choice, poly(A) tail length, and RNA base modifications. The emerging theme from recent studies is that reprogramming of gene expression in response to developmental cues and stresses at the co-/post-transcriptional level likely plays a crucial role in eliciting appropriate responses for optimal growth and plant survival under adverse conditions. Although the mechanisms by which developmental cues and different stresses regulate co-/post-transcriptional splicing are largely unknown, a few recent studies indicate that the external cues target spliceosomal and splicing regulatory proteins to modulate alternative splicing. In this review, we provide an overview of recent discoveries on the dynamics and complexities of plant transcriptomes, mechanistic insights into splicing regulation, and discuss critical gaps in co-/post-transcriptional research that need to be addressed using diverse genomic and biochemical approaches.

Adding Some "Splice" to Stress Eating: Autophagy, ESCRT and Alternative Splicing Orchestrate the Cellular Stress Response

Autophagy is a widely studied self-renewal pathway that is essential for degrading damaged cellular organelles or recycling biomolecules to maintain cellular homeostasis, particularly under cellular stress. This pathway initiates with formation of an autophagosome, which is a double-membrane structure that envelopes cytosolic components and fuses with a lysosome to facilitate degradation of the contents. The endosomal sorting complexes required for transport (ESCRT) proteins play an integral role in controlling autophagosome fusion events and disruption to this machinery leads to autophagosome accumulation. Given the central role of autophagy in maintaining cellular health, it is unsurprising that dysfunction of this process is associated with many human maladies including cancer and neurodegenerative diseases. The cell can also rapidly respond to cellular stress through alternative pre-mRNA splicing that enables adaptive changes to the cell's proteome in response to stress. Thus, alternative pre-mRNA splicing of genes that are involved in autophagy adds another layer of complexity to the cell's stress response. Consequently, the dysregulation of alternative splicing of genes associated with autophagy and ESCRT may also precipitate disease states by either reducing the ability of the cell to respond to stress or triggering a maladaptive response that is pathogenic. In this review, we summarize the diverse roles of the ESCRT machinery and alternative splicing in regulating autophagy and how their dysfunction can have implications for human disease.

Cotton D genome assemblies built with long-read data unveil mechanisms of centromere evolution and stress tolerance divergence

BACKGROUND

Many of genome features which could help unravel the often complex post-speciation evolution of closely related species are obscured because of their location in chromosomal regions difficult to accurately characterize using standard genome analysis methods, including centromeres and repeat regions.

RESULTS

Here, we analyze the genome evolution and diversification of two recently diverged sister cotton species based on nanopore long-read sequence assemblies and Hi-C 3D genome data. Although D genomes are conserved in gene content, they have diversified in gene order, gene structure, gene family diversification, 3D chromatin structure, long-range regulation, and stress-related traits. Inversions predominate among D genome rearrangements. Our results support roles for 5mC and 6mA in gene activation, and 3D chromatin analysis showed that diversification in proximal-vs-distal regulatory-region interactions shape the regulation of defense-related-gene expression. Using a newly developed method, we accurately positioned cotton centromeres and found that these regions have undergone obviously more rapid evolution relative to chromosome arms. We also discovered a cotton-specific LTR class that clarifies evolutionary trajectories among diverse cotton species and identified genetic networks underlying the Verticillium tolerance of Gossypium thurberi (e.g., SA signaling) and salt-stress tolerance of Gossypium davidsonii (e.g., ethylene biosynthesis). Finally, overexpression of G. thurberi genes in upland cotton demonstrated how wild cottons can be exploited for crop improvement.

CONCLUSIONS

Our study substantially deepens understanding about how centromeres have developed and evolutionarily impacted the divergence among closely related cotton species and reveals genes and 3D genome structures which can guide basic investigations and applied efforts to improve crops.

A manipulative interplay between positive and negative regulators of phytohormones: A way forward for improving drought tolerance in plants

Among different abiotic stresses, drought stress is the leading cause of impaired plant growth and low productivity worldwide. It is therefore essential to understand the process of drought tolerance in plants and thus to enhance drought resistance. Accumulating evidence indicates that phytohormones are essential signaling molecules that regulate diverse processes of plant growth and development under drought stress. Plants can often respond to drought stress through a cascade of phytohormones signaling as a means of plant growth regulation. Understanding biosynthesis pathways and regulatory crosstalk involved in these vital compounds could pave the way for improving plant drought tolerance while maintaining overall plant health. In recent years, the identification of phytohormones related key regulatory genes and their manipulation through state-of-the-art genome engineering tools have helped to improve drought tolerance plants. To date, several genes linked to phytohormones signaling networks, biosynthesis, and metabolism have been described as a promising contender for engineering drought tolerance. Recent advances in functional genomics have shown that enhanced expression of positive regulators involved in hormone biosynthesis could better equip plants against drought stress. Similarly, knocking down negative regulators of phytohormone biosynthesis can also be very effective to negate the negative effects of drought on plants. This review explained how manipulating positive and negative regulators of phytohormone signaling could be improvised to develop future crop varieties exhibiting higher drought tolerance. In addition, we also discuss the role of a promising genome editing tool, CRISPR/Cas9, on phytohormone mediated plant growth regulation for tackling drought stress.

Transcriptome analysis of the impact of exogenous methyl jasmonate on the opening of sorghum florets

Background

Methyl Jasmonate (MeJA) could promote the opening of sorghum florets, but the molecular mechanism remains unclear.

Objective

We aimed to investigate the molecular mechanism of exogenous MeJA in promoting the opening of sorghum florets.

Methods

Hybrid sorghum Aikang-8 was selected as the test material in this study. Sorghum plants of uniform growth with approximately 20%-25% florets open were selected and treated with 0, 0.5 and 2.0 mmol/L of MeJA. Totally there were 27 samples with lodicules removed were obtained at different time points and used for the transcriptome analysis using the BGISEQ_500RS platform.

Results

The results showed the sorghum florets opened earlier than the control after the treatment with exogenous MeJA, and the promotive effect increased along with the increase of exogenous MeJA concentration. The number of differentially expressed genes (DEGs) in plasma cells increased with the increase of MeJA concentration, whether up- or down-regulated, after the exogenous MeJA treatment. Besides, the number of metabolic pathways was also positively correlated with the concentration of MeJA. GO and KEGG analysis suggested the DEGs were mainly enriched in starch and sucrose metabolism-related pathways (i.e., LOC8063704, LOC8083539 and LOC8056206), plant hormone signal transduction pathways (i.e., LOC8084842, LOC8072010, and LOC8057408), energy metabolic pathway (i.e., LOC8076139) and the α-linolenic acid metabolic pathway (i.e., LOC8055636, LOC8057399, LOC8063048 and LOC110430730). Functional analysis of target genes showed that two genes named LOC-1 (LOC8063704) and LOC-2 (LOC8076139) could induce the earlier flowering of Arabidopsis thaliana.

Conclusion

The results of this study suggest that exogenous MeJA treatments could induce the up- or down- regulation of genes related to starch and sucrose metabolism, -linolenic acid metabolism and plant hormone signal transduction pathways in the plasma cells of sorghum florets, thereby promoting the opening of sorghum florets.

Physiological changes and transcriptome profiling in Saccharum spontaneum L. leaf under water stress and re-watering conditions

As the polyploidy progenitor of modern sugarcane, Saccharum spontaneum is considered to be a valuable resistance source to various biotic and abiotic stresses. However, little has been reported on the mechanism of drought tolerance in S. spontaneum. Herein, the physiological changes of S. spontaneum GXS87-16 at three water-deficit levels (mild, moderate, and severe) and after re-watering during the elongation stage were investigated. RNA sequencing was utilized for global transcriptome profiling of GXS87-16 under severe drought and re-watered conditions. There were significant alterations in the physiological parameters of GXS87-16 in response to drought stress and then recovered differently after re-watering. A total of 1569 differentially expressed genes (DEGs) associated with water stress and re-watering were identified. Notably, the majority of the DEGs were induced by stress. GO functional annotations and KEGG pathway analysis assigned the DEGs to 47 GO categories and 93 pathway categories. The pathway categories were involved in various processes, such as RNA transport, mRNA surveillance, plant hormone signal transduction, and plant-pathogen interaction. The reliability of the RNA-seq results was confirmed by qRT-PCR. This study shed light on the regulatory processes of drought tolerance in S. spontaneum and identifies useful genes for genetic improvement of drought tolerance in sugarcane.

Combined Transcriptomic and Metabolomic Analysis Reveals the Role of Phenylpropanoid Biosynthesis Pathway in the Salt Tolerance Process of Sophora alopecuroides

Salt stress is the main abiotic stress that limits crop yield and agricultural development. Therefore, it is imperative to study the effects of salt stress on plants and the mechanisms through which plants respond to salt stress. In this study, we used transcriptomics and metabolomics to explore the effects of salt stress on Sophora alopecuroides. We found that salt stress incurred significant gene expression and metabolite changes at 0, 4, 24, 48, and 72 h. The integrated transcriptomic and metabolomic analysis revealed that the differentially expressed genes (DEGs) and differential metabolites (DMs) obtained in the phenylpropanoid biosynthesis pathway were significantly correlated under salt stress. Of these, 28 DEGs and seven DMs were involved in lignin synthesis and 23 DEGs and seven DMs were involved in flavonoid synthesis. Under salt stress, the expression of genes and metabolites related to lignin and flavonoid synthesis changed significantly. Lignin and flavonoids may participate in the removal of reactive oxygen species (ROS) in the root tissue of S. alopecuroides and reduced the damage caused under salt stress. Our research provides new ideas and genetic resources to study the mechanism of plant responses to salt stress and further improve the salt tolerance of plants.

Temporal salt stress-induced transcriptome alterations and regulatory mechanisms revealed by PacBio long-reads RNA sequencing in Gossypium hirsutum

Background

Cotton (Gossypium hirsutum) is considered a fairly salt tolerant crop however, salinity can still cause significant economic losses by affecting the yield and deteriorating the fiber quality. We studied a salt-tolerant upland cotton cultivar under temporal salt stress to unfold the salt tolerance molecular mechanisms. Biochemical response to salt stress (400 mM) was measured at 0 h, 3 h, 12 h, 24 h and 48 h post stress intervals and single-molecule long-read sequencing technology from Pacific Biosciences (PacBio) combined with the unique molecular identifiers approach was used to identify differentially expressed genes (DEG).

Results

Antioxidant enzymes including, catalase (CAT), peroxidase (POD), superoxide dismutase (SOD) were found significantly induced under temporal salt stress, suggesting that reactive oxygen species scavenging antioxidant machinery is an essential component of salt tolerance mechanism in cotton. We identified a wealth of novel transcripts based on the PacBio long reads sequencing approach. Prolonged salt stress duration induces high number of DEGs. Significant numbers of DEGs were found under key terms related to stress pathways such as “response to oxidative stress”, “response to salt stress”, “response to water deprivation”, “cation transport”, “metal ion transport”, “superoxide dismutase”, and “reductase”. Key DEGs related to hormone (abscisic acid, ethylene and jasmonic acid) biosynthesis, ion homeostasis (CBL-interacting serine/threonine-protein kinase genes, calcium-binding proteins, potassium transporter genes, potassium channel genes, sodium/hydrogen exchanger or antiporter genes), antioxidant activity (POD, SOD, CAT, glutathione reductase), transcription factors (myeloblastosis, WRKY, Apetala 2) and cell wall modification were found highly active in response to salt stress in cotton. Expression fold change of these DEGs showed both positive and negative responses, highlighting the complex nature of salt stress tolerance mechanisms in cotton.

Conclusion

Collectively, this study provides a good insight into the regulatory mechanism under salt stress in cotton and lays the foundation for further improvement of salt stress tolerance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07260-z.

Photosynthetic and transcriptomic responses of two C4 grass species with different NaCl tolerance

This report reveals the effects of salt on the photosynthetic electron transport and transcriptome of the glycophyte Setaria viridis (S. viridis) and its salt-tolerant close relative halophyte Spartina alterniflora (S. alterniflora). S. viridis was unable to survive exposed to sodium chloride (NaCl) levels higher than 100 mM, in contrast, S. alterniflora could tolerate NaCl up to 550 mM, with negligible effect on gas exchange related parameters and conductance of electrons transport chain (g_ETC). Under salt, the prompt fluorescence (OJIP-curves) exhibits an increase in the O- and J-steps in S. viridis and much less for S. alterniflora. Flowing NaCl stress, a dramatic decline in the photosystem II (PSII) primary photochemistry was observed for S. viridis, as reflected by the drastic drop in F_v/F_m, F_v/F_o and Φ_PSII; however, no substantial change was recorded for these parameters in S. alterniflora. Interestingly, we found an increase in the primary PSII photochemistry (Φ_PSII) for S. alterniflora with increasing either NaCl concentration or NaCl treatment duration. The NPQ magnitude was strongly enhanced for S. viridis even at a low NaCl (50 mM); however, it remains unchangeable or slightly increased for S. alterniflora at NaCl levels above 400 mM. After NaCl treatment, we found an increase in both the proportion of oxidized P₇₀₀ and the amount of active P₇₀₀ in S. viridis and almost no change for S. alterniflora. Under salt, the net photosynthetic rate (A) and stomatal conductance (gs) measurements demonstrate that A decreases earlier in S. viridis, even after one week exposure to only 50 mM NaCl; in contrast, in S. alterniflora, the effect of NaCl on A and g_s was minor even after exposure for two weeks to high NaCl levels. For S. viridis exposed to 50 mM NaCl for 12 d, carbon dioxide (CO₂) at a concentration of 2000 μL L^-1 could not fully restore A to the control (Ctrl) level. Conversely, in S. alterniflora, high CO₂ can fully restore A for all NaCl treatments except at 550 mM. RNA-seq data shows a major impact of NaCl on metabolic pathways in S. viridis and we found a number of transcription factors potentially related to NaCl responses. For S. alterniflora, no major changes in the transcriptomic levels were recorded under NaCl stress. To confirm our data analysis of RNA-seq, we performed quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis for randomly selected four genes for each species (8 genes in total) and we found that our results (up- and/or down-regulated genes) are fully consistent and match well our RNA-seq data. Overall, this study showed drastically different photosynthetic and transcriptomic responses of a salt-tolerant C₄ grass species and one salt-sensitive C₄ grass species to NaCl stress, which suggests that S. alterniflora could be used as a promising model species to study salt tolerance in C₄ or monocot species.

A histone deacetylase, GhHDT4D, is positively involved in cotton response to drought stress

Acetylation and deacetylation of histones are important for regulating a series of biological processes in plants. Histone deacetylases (HDACs) control the histone deacetylation that plays an important role in plant response to abiotic stress. In our study, we show the evidence that GhHDT4D (a member of the HD2 subfamily of HDACs) is involved in cotton (Gossypium hirsutum) response to drought stress. Overexpression of GhHDT4D in Arabidopsis increased plant tolerance to drought, whereas silencing GhHDT4D in cotton resulted in plant sensitivity to drought. Simultaneously, the H3K9 acetylation level was altered in the GhHDT4D silenced cotton, compared with the controls. Further study revealed that GhHDT4D suppressed the transcription of GhWRKY33, which plays a negative role in cotton defense to drought, by reducing its H3K9 acetylation level. The expressions of the stress-related genes, such as GhDREB2A, GhDREB2C, GhSOS2, GhRD20-1, GhRD20-2 and GhRD29A, were significantly decreased in the GhHDT4D silenced cotton, but increased in the GhWRKY33 silenced cotton. Given these data together, our findings suggested that GhHDT4D may enhance drought tolerance by suppressing the expression of GhWRKY33, thereby activating the downstream drought response genes in cotton.

Tissue-specific changes in the RNA structurome mediate salinity response in Arabidopsis

Little is known concerning the effects of abiotic factors on in vivo RNA structures. We applied Structure-seq to assess the in vivo mRNA structuromes of Arabidopsis thaliana under salinity stress, which negatively impacts agriculture. Structure-seq utilizes dimethyl sulfate reactivity to identify As and Cs that lack base-pairing or protection. Salt stress refolded transcripts differentially in root versus shoot, evincing tissue specificity of the structurome. Both tissues exhibited an inverse correlation between salt stress-induced changes in transcript reactivity and changes in abundance, with stress-related mRNAs showing particular structural dynamism. This inverse correlation is more pronounced in mRNAs wherein the mean reactivity of the 5'UTR, CDS, and 3'UTR concertedly change under salinity stress, suggesting increased susceptibility to abundance control mechanisms in transcripts exhibiting this phenomenon, which we name "concordancy." Concordant salinity-induced increases in reactivity were notably observed in photosynthesis genes, thereby implicating mRNA structural loss in the well-known depression of photosynthesis by salt stress. Overall, changes in secondary structure appear to impact mRNA abundance, molding the functional specificity of the transcriptome under stress.

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.

Genome-wide association reveals genetic variation of lint yield components under salty field conditions in cotton (Gossypium hirsutum L.)

BACKGROUND

Salinity is one of the most significant environmental factors limiting the productivity of cotton. However, the key genetic components responsible for the reduction in cotton yield in saline-alkali soils are still unclear.

RESULTS

Here, we evaluated three main components of lint yield, single boll weight (SBW), lint percentage (LP) and boll number per plant (BNPP), across 316 G. hirsutum accessions under four salt conditions over two years. Phenotypic analysis indicated that LP was unchanged under different salt conditions, however BNPP decreased significantly and SBW increased slightly under high salt conditions. Based on 57,413 high-quality single nucleotide polymorphisms (SNPs) and genome-wide association study (GWAS) analysis, a total of 42, 91 and 25 stable quantitative trait loci (QTLs) were identified for SBW, LP and BNPP, respectively. Phenotypic and QTL analysis suggested that there was little correlation among the three traits. For LP, 8 stable QTLs were detected simultaneously in four different salt conditions, while fewer repeated QTLs for SBW or BNPP were identified. Gene Ontology (GO) analysis indicated that their regulatory mechanisms were also quite different. Via transcriptome profile data, we detected that 10 genes from the 8 stable LP QTLs were predominantly expressed during fiber development. Further, haplotype analyses found that a MYB gene (GhMYB103), with the two SNP variations in cis-regulatory and coding regions, was significantly correlated with lint percentage, implying a crucial role in lint yield. We also identified that 40 candidate genes from BNPP QTLs were salt-inducible. Genes related to carbohydrate metabolism and cell structure maintenance were rich in plants grown in high salt conditions, while genes related to ion transport were active in plants grown in low salt conditions, implying different regulatory mechanisms for BNPP at high and low salt conditions.

CONCLUSIONS

This study provides a foundation for elucidating cotton salt tolerance mechanisms and contributes gene resources for developing upland cotton varieties with high yields and salt stress tolerance.

A perspective on crassulacean acid metabolism photosynthesis evolution of orchids on different continents: Dendrobium as a case study

Members of the Orchidaceae, one of the largest families of flowering plants, evolved the crassulacean acid metabolism (CAM) photosynthesis strategy. It is thought that CAM triggers adaptive radiation into new niche spaces, yet very little is known about its origin and diversification on different continents. Here, we assess the prevalence of CAM in Dendrobium, which is one of the largest genera of flowering plants and found in a wide range of environments, from the high altitudes of the Himalayas to relatively arid habitats in Australia. Based on phylogenetic time trees, we estimated that CAM, as determined by δ 13C values less negative than -20.0‰, evolved independently at least eight times in Dendrobium. The oldest lineage appeared in the Asian clade during the middle Miocene, indicating the origin of CAM was associated with a pronounced climatic cooling that followed a period of aridity. Divergence of the four CAM lineages in the Asian clade appeared to be earlier than divergence of those in the Australasian clade. However, CAM species in the Asian clade are much less diverse (25.6%) than those in the Australasian clade (57.9%). These findings shed new light on CAM evolutionary history and the aridity levels of the paleoclimate on different continents.

Transcriptomic and Co-Expression Network Profiling of Shoot Apical Meristem Reveal Contrasting Response to Nitrogen Rate between Indica and Japonica Rice Subspecies

Reducing nitrogen (N) input is a key measure to achieve a sustainable rice production in China, especially in Jiangsu Province. Tiller is the basis for achieving panicle number that plays as a major factor in the yield determination. In actual production, excessive N is often applied in order to produce enough tillers in the early stages. Understanding how N regulates tillering in rice plants is critical to generate an integrative management to reduce N use and reaching tiller number target. Aiming at this objective, we utilized RNA sequencing and weighted gene co-expression network analysis (WGCNA) to compare the transcriptomes surrounding the shoot apical meristem of indica (Yangdao6, YD6) and japonica (Nipponbare, NPB) rice subspecies. Our results showed that N rate influenced tiller number in a different pattern between the two varieties, with NPB being more sensitive to N enrichment, and YD6 being more tolerant to high N rate. Tiller number was positively related to N content in leaf, culm and root tissue, but negatively related to the soluble carbohydrate content, regardless of variety. Transcriptomic comparisons revealed that for YD6 when N rate enrichment from low (LN) to medium (MN), it caused 115 DEGs (LN vs. MN), from MN to high level (HN) triggered 162 DEGs (MN vs. HN), but direct comparison of low with high N rate showed a 511 DEGs (LN vs. HN). These numbers of DEG in NPB were 87 (LN vs. MN), 40 (MN vs. HN), and 148 (LN vs. HN). These differences indicate that continual N enrichment led to a bumpy change at the transcription level. For the reported sixty-five genes which affect tillering, thirty-six showed decent expression in SAM at tiller starting phase, among them only nineteen being significantly influenced by N level, and two genes showed significant interaction between N rate and variety. Gene ontology analysis revealed that the majority of the common DEGs are involved in general stress responses, stimulus responses, and hormonal signaling process. WGCNA network identified twenty-two co-expressing gene modules and ten candidate hubgenes for each module. Several genes associated with tillering and N rate fall on the related modules. These indicate that there are more genes participating in tillering regulation in response to N enrichment.

Genome-wide association and differential expression analysis of salt tolerance in Gossypium hirsutum L at the germination stage

BACKGROUND

Salinity is a major abiotic stress seriously hindering crop yield. Development and utilization of tolerant varieties is the most economical way to address soil salinity. Upland cotton is a major fiber crop and pioneer plant on saline soil and thus its genetic architecture underlying salt tolerance should be extensively explored.

RESULTS

In this study, genome-wide association analysis and RNA sequencing were employed to detect salt-tolerant qualitative-trait loci (QTLs) and candidate genes in 196 upland cotton genotypes at the germination stage. Using comprehensive evaluation values of salt tolerance in four environments, we identified 33 significant single-nucleotide polymorphisms (SNPs), including 17 and 7 SNPs under at least two and four environments, respectively. The 17 stable SNPs were located within or near 98 candidate genes in 13 QTLs, including 35 genes that were functionally annotated to be involved in salt stress responses. RNA-seq analysis indicated that among the 98 candidate genes, 13 were stably differentially expressed. Furthermore, 12 of the 13 candidate genes were verified by qRT-PCR. RNA-seq analysis detected 6640, 3878, and 6462 differentially expressed genes at three sampling time points, of which 869 were shared.

CONCLUSIONS

These results, including the elite cotton accessions with accurate salt tolerance evaluation, the significant SNP markers, the candidate genes, and the salt-tolerant pathways, could improve our understanding of the molecular regulatory mechanisms under salt stress tolerance and genetic manipulation for cotton improvement.

Reconstruction of the full-length transcriptome atlas using PacBio Iso-Seq provides insight into the alternative splicing in Gossypium australe

BACKGROUND

Gossypium australe F. Mueller (2n = 2x = 26, G₂ genome) possesses valuable characteristics. For example, the delayed gland morphogenesis trait causes cottonseed protein and oil to be edible while retaining resistance to biotic stress. However, the lack of gene sequences and their alternative splicing (AS) in G. australe remain unclear, hindering to explore species-specific biological morphogenesis.

RESULTS

Here, we report the first sequencing of the full-length transcriptome of the Australian wild cotton species, G. australe, using Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq) from the pooled cDNA of ten tissues to identify transcript loci and splice isoforms. We reconstructed the G. australe full-length transcriptome and identified 25,246 genes, 86 pre-miRNAs and 1468 lncRNAs. Most genes (12,832, 50.83%) exhibited two or more isoforms, suggesting a high degree of transcriptome complexity in G. australe. A total of 31,448 AS events in five major types were found among the 9944 gene loci. Among these five major types, intron retention was the most frequent, accounting for 68.85% of AS events. 29,718 polyadenylation sites were detected from 14,536 genes, 7900 of which have alternative polyadenylation sites (APA). In addition, based on our AS events annotations, RNA-Seq short reads from germinating seeds showed that differential expression of these events occurred during seed germination. Ten AS events that were randomly selected were further confirmed by RT-PCR amplification in leaf and germinating seeds.

CONCLUSIONS

The reconstructed gene sequences and their AS in G. australe would provide information for exploring beneficial characteristics in G. australe.

Transcriptome profiling of Gossypium arboreum during fiber initiation and the genome-wide identification of trihelix transcription factors

Cotton fiber initiation is the first step in fiber development, and it determines the yield. Here, genome-wide transcriptome profiling of Gossypium arboreum was performed to determine the molecular basis of cotton fiber initiation. A comparison of the transcriptomes of fiber-bearing ovules at -0.5, 0, 0.5, 1, 1.5, 2, 2.5 and 3 d post-anthesis detected 12,049 differentially expressed genes that mainly participated in ribosome, carbon metabolism and amino acid biosynthesis pathways. Genes encoding alcohol dehydrogenase 1 and hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase, involving in fatty acid degradation and flavonoid biosynthesis, were enriched. Furthermore, 1049 differentially expressed transcription factors were identified. Among these, 17 were trihelix family transcription factors, which play important roles in plant development and responses to biotic and abiotic stresses. In total, 52 full-length trihelix genes, named as GaGTs, were identified in G. arboreum and located in 12 of the 13 cotton chromosomes. Transcriptomic data and a quantitative real-time PCR analysis indicated that several GaGTs were significantly induced during fiber initiation in G. arboreum. Thus, the genome-wide comprehensive analysis of gene expression in G. arboreum fiber initiation will serve as a useful resource for unraveling the functions of specific genes. The phylogenetic relationships and expression analyses of the G. arboreum trihelix genes established a solid foundation for future comprehensive functional analyses of the GaGTs.