The SR Splicing Factors: Providing Perspectives on Their Evolution, Expression, Alternative Splicing, and Function in Populus trichocarpa

Differences in alternative splicing and their potential underlying factors between animals and plants

BACKGROUND

Alternative splicing (AS), a posttranscriptional process, contributes to the complexity of transcripts from a limited number of genes in a genome, and AS is considered a great source of genetic and phenotypic diversity in eukaryotes. In animals, AS is tightly regulated during the processes of cell growth and differentiation, and its dysregulation is involved in many diseases, including cancers. Likewise, in plants, AS occurs in all stages of plant growth and development, and it seems to play important roles in the rapid reprogramming of genes in response to environmental stressors. To date, the prevalence and functional roles of AS have been extensively reviewed in animals and plants. However, AS differences between animals and plants, especially their underlying molecular mechanisms and impact factors, are anecdotal and rarely reviewed.

AIM OF REVIEW

This review aims to broaden our understanding of AS roles in a variety of biological processes and provide insights into the underlying mechanisms and impact factors likely leading to AS differences between animals and plants.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We briefly summarize the roles of AS regulation in physiological and biochemical activities in animals and plants. Then, we underline the differences in the process of AS between plants and animals and especially analyze the potential impact factors, such as gene exon/intron architecture, 5'/3' untranslated regions (UTRs), spliceosome components, chromatin dynamics and transcription speeds, splicing factors [serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs)], noncoding RNAs, and environmental stimuli, which might lead to the differences. Moreover, we compare the nonsense-mediated mRNA decay (NMD)-mediated turnover of the transcripts with a premature termination codon (PTC) in animals and plants. Finally, we summarize the current AS knowledge published in animals versus plants and discuss the potential development of disease therapies and superior crops in the future.

Pan-transcriptomic analysis reveals alternative splicing control of cold tolerance in rice

Plants have evolved complex mechanisms to adapt to harsh environmental conditions. Rice (Oryza sativa) is a staple food crop that is sensitive to low temperatures. However, its cold stress responses remain poorly understood, thus limiting possibilities for crop engineering to achieve greater cold tolerance. In this study, we constructed a rice pan-transcriptome and characterized its transcriptional regulatory landscape in response to cold stress. We performed Iso-Seq and RNA-Seq of 11 rice cultivars subjected to a time-course cold treatment. Our analyses revealed that alternative splicing-regulated gene expression plays a significant role in the cold stress response. Moreover, we identified CATALASE C (OsCATC) and Os03g0701200 as candidate genes for engineering enhanced cold tolerance. Importantly, we uncovered central roles for the 2 serine-arginine-rich proteins OsRS33 and OsRS2Z38 in cold tolerance. Our analysis of cold tolerance and resequencing data from a diverse collection of 165 rice cultivars suggested that OsRS2Z38 may be a key selection gene in japonica domestication for cold adaptation, associated with the adaptive evolution of rice. This study systematically investigated the distribution, dynamic changes, and regulatory mechanisms of alternative splicing in rice under cold stress. Overall, our work generates a rich resource with broad implications for understanding the genetic basis of cold response mechanisms in plants.

Evidence Supporting a Role of Alternative Splicing Participates in Melon (Cucumis melo L.) Fruit Ripening

One key post-transcriptional modification mechanism that dynamically controls a number of physiological processes in plants is alternative splicing (AS). However, the functional impacts of AS on fruit ripening remain unclear. In this research, we used RNA-seq data from climacteric (VED, Harukei 3) and non-climacteric (PI, PS) melon cultivars to explore alternative splicing (AS) in immature and mature fruit. The results revealed dramatic changes in differential AS genes (DAG) between the young and mature fruit stages, particularly in genes involved in fruit development/ripening, carotenoid and capsaicinoid biosynthesis, and starch and sucrose metabolism. Serine/arginine-rich (SR) family proteins are known as important splicing factors in AS events. From the melon genome, a total of 17 SR members were discovered in this study. These genes could be classified into eight distinct subfamilies based on gene structure and conserved motifs. Promoter analysis detected various cis-acting regulatory elements involved in hormone pathways and fruit development. Interestingly, these SR genes exhibited specific expression patterns in reproductive organs such as flowers and ovaries. Additionally, concurrent with the increase in AS levels in ripening fruit, the transcripts of these SR genes were activated during fruit maturation in both climacteric and non-climacteric melon varieties. We also found that most SR genes were under selection during domestication. These results represent a novel finding of increased AS levels and SR gene expression during fruit ripening, indicating that alternative splicing may play a role in fruit maturation.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

Plant RNA-binding proteins: Phase separation dynamics and functional mechanisms underlying plant development and stress responses

RNA-binding proteins (RBPs) accompany RNA from synthesis to decay, mediating every aspect of RNA metabolism and impacting diverse cellular and developmental processes in eukaryotes. Many RBPs undergo phase separation along with their bound RNA to form and function in dynamic membraneless biomolecular condensates for spatiotemporal coordination or regulation of RNA metabolism. Increasing evidence suggests that phase-separating RBPs with RNA-binding domains and intrinsically disordered regions play important roles in plant development and stress adaptation. Here, we summarize the current knowledge about how dynamic partitioning of RBPs into condensates controls plant development and enables sensing of experimental changes to confer growth plasticity under stress conditions, with a focus on the dynamics and functional mechanisms of RBP-rich nuclear condensates and cytoplasmic granules in mediating RNA metabolism. We also discuss roles of multiple factors, such as environmental signals, protein modifications, and N6-methyladenosine RNA methylation, in modulating the phase separation behaviors of RBPs, and highlight the prospects and challenges for future research on phase-separating RBPs in crops.

Full-length transcriptome sequencing provides insights into alternative splicing under cold stress in peanut

Introduction

Peanut (Arachis hypogaea L.), also called groundnut is an important oil and cash crop grown widely in the world. The annual global production of groundnuts has increased to approximately 50 million tons, which provides a rich source of vegetable oils and proteins for humans. Low temperature (non-freezing) is one of the major factors restricting peanut growth, yield, and geographic distribution. Since the complexity of cold-resistance trait, the molecular mechanism of cold tolerance and related gene networks were largely unknown in peanut.

Methods

In this study, comparative transcriptomic analysis of two peanut cultivars (SLH vs. ZH12) with differential cold tolerance under low temperature (10°C) was performed using Oxford Nanopore Technology (ONT) platform.

Results and discussion

As a result, we identified 8,949 novel gene loci and 95,291 new/novel isoforms compared with the reference database. More differentially expressed genes (DEGs) were discovered in cold-sensitive cultivar (ZH12) than cold-tolerant cultivar (SLH), while more alternative splicing events were found in SLH compared to ZH12. Gene Ontology (GO) analyses of the common DEGs showed that the "response to stress", "chloroplast part", and "transcription factor activity" were the most enriched GO terms, indicating that photosynthesis process and transcription factors play crucial roles in cold stress response in peanut. We also detected a total of 708 differential alternative splicing genes (DASGs) under cold stress compared to normal condition. Intron retention (IR) and exon skipping (ES) were the most prevalent alternative splicing (AS) events. In total, 4,993 transcription factors and 292 splicing factors were detected, many of them had differential expression levels and/or underwent AS events in response to cold stress. Overexpression of two candidate genes (encoding trehalose-6-phosphatephosphatases, AhTPPs) in yeast improves cold tolerance. This study not only provides valuable resources for the study of cold resistance in peanut but also lay a foundation for genetic modification of cold regulators to enhance stress tolerance in crops.

Unveiling unique alternative splicing responses to low temperature in Zoysia japonica through ZjRTD1.0, a high-quality reference transcript dataset

Inadequate reference databases in RNA-seq analysis can hinder data utilization and interpretation. In this study, we have successfully constructed a high-quality reference transcript dataset, ZjRTD1.0, for Zoysia japonica, a widely-used turfgrass with exceptional tolerance to various abiotic stress, including low temperatures and salinity. This dataset comprises 113,089 transcripts from 57,143 genes. BUSCO analysis demonstrates exceptional completeness (92.4%) in ZjRTD1.0, with reduced proportions of fragmented (3.3%) and missing (4.3%) orthologs compared to prior datasets. ZjRTD1.0 enables more precise analyses, including transcript quantification and alternative splicing assessments using public datasets, which identified a substantial number of differentially expressed transcripts (DETs) and differential alternative splicing (DAS) events, leading to several novel findings on Z. japonica's responses to abiotic stresses. First, spliceosome gene expression influenced alternative splicing significantly under abiotic stress, with a greater impact observed during low-temperature stress. Then, a significant positive correlation was found between the number of differentially expressed genes (DEGs) encoding protein kinases and the frequency of DAS events, suggesting the role of protein phosphorylation in regulating alternative splicing. Additionally, our results suggest possible involvement of serine/arginine-rich (SR) proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs) in generating inclusion/exclusion isoforms under low-temperature stress. Furthermore, our investigation revealed a significantly enhanced overlap between DEGs and differentially alternatively spliced genes (DASGs) in response to low-temperature stress, suggesting a unique co-regulatory mechanism governing transcription and splicing in the context of low-temperature response. In conclusion, we have proven that ZjRTD1.0 will serve as a reliable and useful resource for future transcriptomic analyses in Z. japonica.

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.

Orthogroup and phylotranscriptomic analyses identify transcription factors involved in the plant cold response: A case study of Arabidopsis BBX29

C-repeat binding factors (CBFs) are well-known transcription factors (TFs) that regulate plant cold acclimation. RNA sequencing (RNA-seq) data from diverse plant species provide opportunities to identify other TFs involved in the cold response. However, this task is challenging because gene gain and loss has led to an intertwined community of co-orthologs and in-paralogs between and within species. Using orthogroup (closely related homologs) analysis, we identified 10,549 orthogroups in five representative eudicots. A phylotranscriptomic analysis of cold-treated seedlings from eudicots identified 35 high-confidence conserved cold-responsive transcription factor orthogroups (CoCoFos). These 35 CoCoFos included the well-known cold-responsive regulators CBFs, HSFC1, ZAT6/10, and CZF1 among others. We used Arabidopsis BBX29 for experimental validation. Expression and genetic analyses showed that cold-induction of BBX29 is CBF- and abscisic acid-independent, and BBX29 is a negative regulator of cold tolerance. Integrative RNA-seq and Cleavage Under Targets and Tagmentation followed by sequencing analyses revealed that BBX29 represses a set of cold-induced TFs (ZAT12, PRR9, RVE1, MYB96, etc.). Altogether, our analysis yielded a library of eudicot CoCoFos and demonstrated that BBX29 is a negative regulator of cold tolerance in Arabidopsis.

Interplays between cis- and trans-Acting Factors for Alternative Splicing in Response to Environmental Changes during Biological Invasions of Ascidians

Alternative splicing (AS), a pivotal biological process contributing to phenotypic plasticity, creates a bridge linking genotypes with phenotypes. Despite its importance, the AS mechanisms underlying environmental response and adaptation have not been well studied, and more importantly, the cis- and trans-acting factors influencing AS variation remain unclear. Using the model invasive congeneric ascidians, Ciona robusta, and Ciona savignyi, we compared their AS responses to environmental changes and explored the potential determinants. Our findings unveiled swift and dynamic AS changes in response to environmental challenges, and differentially alternative spliced genes (DASGs) were functionally enriched in transmembrane transport processes. Interestingly, both the prevalence and level of AS in C. robusta were lower than those observed in C. savignyi. Furthermore, these two indices were higher under temperature stresses compared to salinity stresses in C. savignyi. All the observed patterns underscore the species-specific and environmental context-dependent AS responses to environmental challenges. The dissimilarities in genomic structure and exon/intron size distributions between these two species likely contributed to the observed AS variation. Moreover, we identified a total of 11 and 9 serine/arginine-rich splicing factors (SRSFs) with conserved domains and gene structures in the genomes of C. robusta and C. savignyi, respectively. Intriguingly, our analysis revealed that all detected SRSFs did not exhibit prevalent AS regulations. Instead, we observed AS control over a set of genes related to splicing factors and spliceosome components. Altogether, our results elucidate species-specific and environmental challenge-dependent AS response patterns in closely related invasive ascidians. The identified splicing factors and spliceosome components under AS control offer promising candidates for further investigations into AS-mediated rapid responses to environmental challenges complementary to SRSFs.

A phylotranscriptomic dataset of angiosperm species under cold stress

Angiosperms are one of the most diverse and abundant plant groups that are widely distributed on Earth, from tropical to temperate and polar zones. The wide distribution of angiosperms may be attributed to the evolution of sophisticated mechanisms of environmental adaptability, including cold tolerance. Since the development of high-throughput sequencing, transcriptome has been widely utilized to gain insights into the molecular mechanisms of plants in response to cold stress. However, previous studies generally focused on single or two species, and comparative transcriptome analyses for multispecies responding to cold stress were limited. In this study, we selected 11 representative angiosperm species, performed phylotranscriptome experiments at four time points before and after cold stress, and presented a profile of cold-induced transcriptome changes in angiosperms. Our multispecies cold-responsive RNA-seq datasets provide valuable references for exploring conserved and evolutionary mechanisms of angiosperms in adaptation to cold stress.

Plant serine/arginine-rich proteins: versatile players in RNA processing

MAIN CONCLUSION

Serine/arginine-rich (SR) proteins participate in RNA processing by interacting with precursor mRNAs or other splicing factors to maintain plant growth and stress responses. Alternative splicing is an important mechanism involved in mRNA processing and regulation of gene expression at the posttranscriptional level, which is the main reason for the diversity of genes and proteins. The process of alternative splicing requires the participation of many specific splicing factors. The SR protein family is a splicing factor in eukaryotes. The vast majority of SR proteins' existence is an essential survival factor. Through its RS domain and other unique domains, SR proteins can interact with specific sequences of precursor mRNA or other splicing factors and cooperate to complete the correct selection of splicing sites or promote the formation of spliceosomes. They play essential roles in the composition and alternative splicing of precursor mRNAs, providing pivotal functions to maintain growth and stress responses in animals and plants. Although SR proteins have been identified in plants for three decades, their evolutionary trajectory, molecular function, and regulatory network remain largely unknown compared to their animal counterparts. This article reviews the current understanding of this gene family in eukaryotes and proposes potential key research priorities for future functional studies.

Integrative multi-omics analysis of three early diverged rosid species reveals an ancient hierarchical cold-responsive regulatory network

Elucidating regulators, including transcription factors (TFs) and RNA-binding proteins (RBPs), underlying gene transcriptional and post-transcriptional co-regulatory network is key to understand plant cold responses. Previous studies were mainly conducted on single species, and whether the regulators are conserved across different species remains elusive. Here, we selected three species that diverged at the early evolution of rosids (~99-113 million years ago), performed cold-responsive phylotranscriptome experiments, and integrated chromatin immunoprecipitation- and DNA affinity purification-sequencing (ChIP/DAP-seq) analysis to explore cold-responsive regulators and their regulatory networks. First, we detected over 10,000 cold-induced differentially expressed genes (DEGs) and alternative splicing genes (DASGs) in each species. Among the DEGs, a set of TFs and RBPs were conserved in rosid cold response. Compared to TFs, RBPs displayed a delayed cold-responsive pattern, implying a hierarchical regulation of DEGs and DASGs. By integrating DEGs and DASGs, we identified 259 overlapping DE-DASG orthogroups (closely-related homologs) that were cold-regulated at both transcriptional and post-transcriptional levels in all three studied species. Notably, pathway analysis on each of the DEGs, DASGs, and DE-DASGs in the three species showed a common enrichment connected to the circadian rhythm. Evidently, 226 cold-responsive genes were directly targeted by at least two circadian rhythm components (CCA1, LHY, RV4, RVE7, and RVE8). Finally, we revealed an ancient hierarchy of cold-responsive regulatory networks at transcriptional and post-transcriptional levels launched by circadian components in rosids. Altogether, this study sheds light on conserved regulators underlying cold-responsive regulatory networks across rosid species, despite a long evolutionary history after their divergence.

NTR1 is involved in heat stress tolerance through mediating expression regulation and alternative splicing of heat stress genes in Arabidopsis

As a common adverse environmental factor, heat stress (HS) not only drastically changes the plant transcriptome at the transcription level but also increases alternative splicing (AS), especially intron retention (IR) events. However, the exact mechanisms are not yet well understood. Here, we reported that NTC-related protein 1 (NTR1), which acts as an accessory component for spliceosome disassembly, is necessary for this process. The mutants of NTR1, both the T-DNA insertion and the point mutation identified through ethyl methanesulfonate (EMS) mutagenesis screening, are vulnerable to HS, indicating that NTR1 is essential for plant HS tolerance. At the molecular level, genes of response to heat and response to temperature stimulus are highly enriched among those of heat-induced but less-expressed ntr1 mutants. Moreover, a large portion of HS response (HSR) genes such as heat shock transcription factors (HSFs) and heat shock proteins (HSPs) are less induced by heat treatment, and more AS events, especially IR events, were found in heat-treated ntr1 mutants. Furthermore, HS suppressed the expression of NTR1 and NTR1-associated complex components. Thus, it is very likely that upon HS, the plant reduces the expression of the NTR1-associated complex to fulfill the fast demands for transcription of HSR genes such as HSFs and HSPs, which in turn results in the accumulation of improperly spliced especially IR products and eventually causes harm to plants.

Alternative splicing: transcriptional regulatory network in agroforestry

Alternative splicing (AS) in plants plays a key role in regulating the expression of numerous transcripts from a single gene in a regulatory pathway. Variable concentrations of growth regulatory hormones and external stimuli trigger alternative splicing to switch among different growth stages and adapt to environmental stresses. In the AS phenomenon, a spliceosome causes differential transcriptional modifications in messenger RNA (mRNAs), resulting in partial or complete retention of one or more introns as compared to fully spliced mRNA. Differentially expressed proteins translated from intron-retaining messenger RNA (mRNA^ir) perform vital functions in the feedback mechanism. At the post-transcriptional level, AS causes the remodeling of transcription factors (TFs) by the addition or deletion of binding domains to activate and/or repress transcription. In this study, we have summarized the specific role of AS in the regulation of gene expression through repression and activation of the transcriptional regulatory network under external stimuli and switch among developmental stages.

Innovations and stepwise evolution of CBFs/DREB1s and their regulatory networks in angiosperms

The C-repeat binding factors/dehydration-responsive element binding protein 1s (CBFs/DREB1s) have been identified as major regulators of cold acclimation in many angiosperm plants. However, their origin and evolutionary process associated to cold responsiveness are still lacking. By integrating multi-omics data of genomes, transcriptomes, and CBFs/DREB1s genome-wide binding profiles, we unveil the origin and evolution of CBFs/DREB1s and their regulatory network. Gene collinearity and phylogeny analyses show that CBF/DREB1 is an innovation evolved from tandem duplication-derived DREB III gene. A subsequent event of ε-whole genome duplication led to two CBF/DREB1 archetypes (Clades I and II) in ancient angiosperms. In contrast to cold-insensitivity of Clade I and their parent DREB III genes, Clade II evolved a further innovation in cold-sensitive response and was stepwise expanded in eudicots and monocots by independent duplications. In geological time, the duplication events were mainly enriched around the Cretaceous-Paleogene (K-Pg) boundary and/or in the Late Cenozoic Ice Age, when the global average temperature significantly decreased. Consequently, the duplicated CBF/DREB1 genes contributed to the rewiring of CBFs/DREB1s-regulatory network for cold tolerance. Altogether, our results highlight an origin and convergent evolution of CBFs/DREB1s and their regulatory network probably for angiosperms adaptation to global cooling.

Regulatory Network of Serine/Arginine-Rich (SR) Proteins: The Molecular Mechanism and Physiological Function in Plants

Serine/arginine-rich (SR) proteins are a type of splicing factor. They play significant roles in constitutive and alternative pre-mRNA splicing, and are involved in post-splicing activities, such as mRNA nuclear export, nonsense-mediated mRNA decay, mRNA translation, and miRNA biogenesis. In plants, SR proteins function under a complex regulatory network by protein–protein and RNA–protein interactions between SR proteins, other splicing factors, other proteins, or even RNAs. The regulatory networks of SR proteins are complex—they are regulated by the SR proteins themselves, they are phosphorylated and dephosphorylated through interactions with kinase, and they participate in signal transduction pathways, whereby signaling cascades can link the splicing machinery to the exterior environment. In a complex network, SR proteins are involved in plant growth and development, signal transduction, responses to abiotic and biotic stresses, and metabolism. Here, I review the current status of research on plant SR proteins, construct a model of SR proteins function, and ask many questions about SR proteins in plants.

OsSCL30 overexpression reduces the tolerance of rice seedlings to low temperature, drought and salt

Rice is one of the main food crops for the world population. Various abiotic stresses, such as low temperature, drought, and high salinity, affect rice during the entire growth period, determining its yield and quality, and even leading to plant death. In this study, by constructing overexpression vectors D-163 + 1300:OsSCL30 and D-163 + 1300-AcGFP:OsSCL30-GFP, the mechanism of action of OsSCL30 in various abiotic stresses was explored. Bioinformatics analysis showed that OsSCL30 was located on the chromosome 12 of rice Nipponbare, belonging to the plant-specific SCL subfamily of the SR protein family. The 1500 bp section upstream of the open reading frame start site contains stress-related cis-acting elements such as ABRE, MYC, and MYB. Under normal conditions, the expression of OsSCL30 was higher in leaves and leaf sheaths. The results of reverse transcription polymerase chain reaction showed that the expression of OsSCL30 decreased after low temperature, drought and salt treatment. In root cells OsSCL30 was localized in the nuclei. The results of the rice seedling tolerance and recovery tests showed that overexpression of OsSCL30 diminished the resistance to low temperature, drought and salt stresses in transgenic rice and resulted in larger accumulation of reactive oxygen species. This study is of great significance for exploring the response mechanisms of SR proteins under abiotic stresses.

Role of lncRNAs in cis- and trans-regulatory responses to salt in Populus trichocarpa

Long non-coding RNAs (lncRNAs) are emerging as versatile regulators in diverse biological processes. However, little is known about their cis- and trans-regulatory contributions in gene expression under salt stress. Using 27 RNA-seq data sets from Populus trichocarpa leaves, stems and roots, we identified 2988 high-confidence lncRNAs, including 1183 salt-induced differentially expressed lncRNAs. Among them, 301 lncRNAs have potential for positively affecting their neighboring genes, predominantly in a cis-regulatory manner rather than by co-transcription. Additionally, a co-expression network identified six striking salt-associated modules with a total of 5639 genes, including 426 lncRNAs, and in these lncRNA sequences, the DNA/RNA binding motifs are enriched. This suggests that lncRNAs might contribute to distant gene expression of the salt-associated modules in a trans-regulatory manner. Moreover, we found 30 lncRNAs that have potential to simultaneously cis- and trans-regulate salt-responsive homologous genes, and Ptlinc-NAC72, significantly induced under long-term salt stress, was selected for validating its regulation of the expression and functional roles of the homologs PtNAC72.A and PtNAC72.B (PtNAC72.A/B). The transient transformation of Ptlinc-NAC72 and a dual-luciferase assay of Ptlinc-NAC72 and PtNAC72.A/B promoters confirmed that Ptlinc-NAC72 can directly upregulate PtNAC72.A/B expression, and a presence/absence assay was further conducted to show that the regulation is probably mediated by Ptlinc-NAC72 recognizing the tandem elements (GAAAAA) in the PtNAC72.A/B 5' untranslated region (5'-UTR). Finally, the overexpression of Ptlinc-NAC72 produces a hypersensitive phenotype under salt stress. Altogether, our results shed light on the cis- and trans-regulation of gene expression by lncRNAs in Populus and provides an example of long-term salt-induced Ptlinc-NAC72 that could be used to mitigate growth costs by conferring plant resilience to salt stress.

A Comprehensive Identification and Function Analysis of Serine/Arginine-Rich (SR) Proteins in Cotton (Gossypium spp.)

As one of the most important factors in alternative splicing (AS) events, serine/arginine-rich (SR) proteins not only participate in the growth and development of plants but also play pivotal roles in abiotic stresses. However, the research about SR proteins in cotton is still lacking. In this study, we performed an extensive comparative analysis of SR proteins and determined their phylogeny in the plant lineage. A total of 169 SR family members were identified from four Gossypium species, and these genes could be divided into eight distinct subfamilies. The domain, motif distribution and gene structure of cotton SR proteins are conserved within each subfamily. The expansion of SR genes is mainly contributed by WGD and allopolyploidization events in cotton. The selection pressure analysis showed that all the paralogous gene pairs were under purifying selection pressure. Many cis-elements responding to abiotic stress and phytohormones were identified in the upstream sequences of the GhSR genes. Expression profiling suggested that some GhSR genes may involve in the pathways of plant resistance to abiotic stresses. The WGCNA analysis showed that GhSCL-8 co-expressed with many abiotic responding related genes in a salt-responding network. The Y2H assays showed that GhSCL-8 could interact with GhSRs in other subfamilies. The subcellular location analysis showed that GhSCL-8 is expressed in the nucleus. The further VIGS assays showed that the silencing of GhSCL-8 could decrease salt tolerance in cotton. These results expand our knowledge of the evolution of the SR gene family in plants, and they will also contribute to the elucidation of the biological functions of SR genes in the future.