Identification and characterization of cadmium stress-related LncRNAs from Betula platyphylla

Unraveling cadmium tolerance mechanisms in Betula platyphylla through a hierarchical gene regulatory network in hormone signaling

Cadmium (Cd), a toxic heavy metal, is a significant pollutant that impacts plant productivity. While some studies have been conducted, the underlying mechanisms by which plants respond to Cd stress remain largely unclear. Here, we performed RNA-seq analysis of Betula platyphylla (birch) under CdCl2 treatment. The findings revealed a substantial enrichment of differentially expressed genes (DEGs) in pathways associated with plant hormones. A gene regulatory network (GRN) was constructed, and the regulatory relationships between genes were determined using a partial correlation coefficient algorithm. The GRN comprises 2151 regulatory interactions, including 7 transcription factors (TFs) from the first layer, 25 TFs from the second layer, and 168 structural genes from the third layer, all of which are linked to ten enriched biological processes. ChIP-PCR and qRT-PCR assays validated approximately 85.2 % of the predicted interactions between the first and second layers, along with 88.3 % of the interactions between the second and third layers, supporting the validity of the GRN. Eighteen genes were selected from the third layer of multiple biological pathways to analyze their functions, and the results indicated that these genes can enhance Cd tolerance in birch plants. Additionally, two TFs in the first layer, BpHD-zip7 and BpRAV1, were successfully introduced into birch plants, confirming their role in improving Cd tolerance. Our findings elucidate the regulatory mechanisms and key determinants that function in the adaptation of B. platyphylla to Cd stress.

Identification and Network Construction of mRNAs, miRNAs, lncRNAs, and circRNAs in Sweetpotato (Ipomoea batatas L.) Adventitious Roots Under Salt Stress via Whole-Transcriptome RNA Sequencing

Sweetpotato is the seventh largest crop worldwide, and soil salinization is a major environmental stress limiting its yield. Recent studies have shown that noncoding RNAs (ncRNAs) play important regulatory roles in plant responses to abiotic stress. However, ncRNAs in sweetpotato remain largely unexplored. This study analyzed the characteristics of salt-responsive ncRNAs in sweetpotato adventitious roots under salt stress via whole-transcriptome RNA sequencing. The results revealed that 3175 messenger RNAs (mRNAs), 458 microRNAs (miRNAs), 544 long-chain ncRNAs (lncRNAs), and 23 circular RNAs (circRNAs) were differentially expressed. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that most differentially expressed mRNAs (DEmRNAs) and miRNAs (DEmiRNAs) were enriched primarily in phenylpropanoid biosynthesis, starch and sucrose metabolism, the Mitogen-Activated Protein Kinase (MAPK) signaling pathway, plant hormone signal transduction, the mRNA surveillance pathway, and ATP-binding cassette (ABC) transporters. Gene Ontology (GO) enrichment analysis revealed that the majority of DEmRNAs, their target DEmiRNAs, and differentially expressed lncRNAs (DElncRNAs) were associated with the cell wall, oxidation-reduction, the plasma membrane, protein phosphorylation, metabolic processes, transcription factor activity, and the regulation of transcription. Additionally, based on the competitive endogenous RNA (ceRNA) hypothesis, we predicted interactions among different RNAs and constructed a salt-responsive ceRNA network comprising 22 DEmiRNAs, 42 DEmRNAs, 27 DElncRNAs, and 10 differentially expressed circRNAs (DEcircRNAs). Some miRNAs, such as miR408, miR169, miR160, miR5139, miR5368, and miR6179, were central to the network, suggesting their crucial roles in the sweetpotato salt response. Our findings provide a foundation for further research into the potential functions of ncRNAs and offer new targets for salt stress resistance improvement through the manipulation of ncRNAs.

Long non-coding RNAs modulate glutathione metabolism gene expression and tolerance to Pb stress in root tissue of endophyte-infected rice seedling

Endophyte can improve plant resistance to Pb stress, and long non-coding RNAs (lncRNAs) have been proved to play a vital role in response to environmental stress. However, there are few studies on the role of lncRNAs induced by endophyte in host plants under Pb stress. Therefore, we conducted high-throughput sequencing on root tissue of endophyte-infected and -uninfected rice seedlings under Pb stress, and analyzed the target genes of differentially expressed lncRNAs (DElncRNAs). The results showed that endophyte infection significantly increased plant height, above-ground fresh weight and dry weight, but significantly decreased root length and under-ground dry weight after 5 d of treatment. A total of 413 DElncRNAs (167 down-regulated and 246 up-regulated) were screened. Their target differentially expressed mRNAs (DEmRNAs) were significantly enriched in glutathione metabolism, plasma membrane and mineral elements transfer etc. DEmRNAs were most significantly enriched in glutathione metabolism, thereinto detected total glutathione, reduced and oxidized glutathione contents, glutathione S-transferase and glutathione reductase activities were significantly increased after 5 d of treatment. DElncRNAs and DEmRNAs were combined with miRNA database to construct ceRNA network. DEmRNAs in ceRNA network were mainly participated in carbohydrate metabolic process, peroxidase activity and phenylpropanoid biosynthesis. This study will help to understand the molecular mechanism elicited by endophytic infection within rice seedlings under Pb stress from the perspective of lncRNA.

Identification and expression analysis of lncRNAs in rice roots (Oryza sativa L.) under elevated CO2 concentration and/or cadmium stress

The gradual rise of CO2 is one of the global climate changes, Cd stress is also a major abiotic stress factor that affects rice (Oryza sativa L.). The rice seedlings were treated under two CO2 concentrations and two CdCl2 concentrations for 7 days (treatments names: 400 ± 20 μmol mol-1 CO2 and 0 μmol L-1 CdCl2 concentrations, AC; 400 ± 20 μmol mol-1 CO2 and 150 μmol L-1 CdCl2 concentrations, Cd; 800 ± 20 μmol mol-1 CO2 and 0 μmol L-1 CdCl2 concentrations, EC; 800 ± 20 μmol mol-1 CO2 and 150 μmol L-1 CdCl2 concentrations, EC + Cd). The lncRNAs informations were analyzed and excavated using high-throughput sequencing, target genes annotation, and qRT-PCR analysis techniques so as to reveal the regulatory mechanism of lncRNAs in rice roots under high CO2 concentrations and/or Cd stress. The results show that: (1) 326 (AC vs Cd), 331 (AC vs EC), 343 (AC vs EC + Cd), 112 (Cd vs EC + Cd) DE-lncRNAs were identified. (2) MAPK signaling pathway-plant (relevant genes Os04g0534166, Os05g0399800 regulated by MSTRG.18576.11, MSTRG.20864.1) and diterpenoid biosynthesis (relevant genes Os12g0491800, Os02g0570400 regulated by MSTRG.8965.1, MSTRG.11509.1) were annotated in AC vs Cd; Under EC relative to AC, DE-lncRNAs were annotated significantly to the flavonoid biosynthesis (relevant genes Os10g0196100, Os10g0320100, Os11g0116300, Os03g0819600 regulated by MSTRG.4612.1, MSTRG.4668.1, MSTRG.6051.1, MSTRG.16669.1); Under composite treatments, relative to AC, DE-lncRNAs were mainly annotated in the plant hormone signal transduction pathway (relevant genes Os03g0180800, Os03g0180900, Os03g0181100 regulated by MSTRG.13776.1). Under combined treatment, elevated CO2 alleviates Cd stress damage by regulating phenylpropanoid biosynthesis through DE-lncRNAs (relevant genes Os09g0419200 regulated by MSTRG. 29,573.1).

Nano-Fe3O4: Enhancing the tolerance of Elymus nutans to Cd stress through regulating programmed cell death

Cadmium (Cd) poses a significant threat to plant growth and the environment. Nano-Fe3O4 is effective in alleviating Cd stress in plants. Elymus nutans Griseb. is an important fodder crop on the Qinghai-Tibetan Plateau (QTP). However, the potential mechanism by which nano-Fe3O4 alleviates Cd stress in E. nutans is not well understood. E. nutans were subjected to single Cd, single nano-Fe3O4, and co-treatment with nano-Fe3O4 and Cd, and the effects on morphology, Cd uptake, antioxidant enzyme activity, reactive oxygen species (ROS) levels and programmed cell death (PCD) were studied to clarify the regulatory mechanism of nano-Fe3O4. The results showed that Cd stress significantly decreased the germination percentage and biomass of E. nutans. The photosynthetic pigment content decreased significantly under Cd stress. Cd stress also caused oxidative stress and lipid peroxidation, accumulation of excessive ROS, resulting in PCD, but the effect of nano-Fe3O4 was different. Seed germination, seedling growth, and physiological processes were analyzed to elucidate the regulatory role of nano-Fe3O4 nanoparticles in promoting photosynthesis, reducing Cd accumulation, scavenging ROS, and regulating PCD, to promote seed germination and seedling growth in E. nutans. This report provides a scientific basis for improving the tolerance of Elymus to Cd stress by using nano-Fe3O4.

Comparative Long Non-Coding Transcriptome Analysis of Three Contrasting Barley Varieties in Response to Aluminum Stress

Aluminum toxicity is a major abiotic stress on acidic soils, leading to restricted root growth and reduced plant yield. Long non-coding RNAs are crucial signaling molecules regulating the expression of downstream genes, particularly under abiotic stress conditions. However, the extent to which lncRNAs participate in the response to aluminum (Al) stress in barley remains largely unknown. Here, we conducted RNA sequencing of root samples under aluminum stress and compared the lncRNA transcriptomes of two Tibetan wild barley genotypes, XZ16 (Al-tolerant) and XZ61 (Al-sensitive), as well as the aluminum-tolerant cultivar Dayton. In total, 268 lncRNAs were identified as aluminum-responsive genes on the basis of their differential expression profiles under aluminum treatment. Through target gene prediction analysis, we identified 938 candidate lncRNA-messenger RNA (mRNA) pairs that function in a cis-acting manner. Subsequently, enrichment analysis showed that the genes targeted by aluminum-responsive lncRNAs were involved in diterpenoid biosynthesis, peroxisome function, and starch/sucrose metabolism. Further analysis of genotype differences in the transcriptome led to the identification of 15 aluminum-responsive lncRNAs specifically altered by aluminum stress in XZ16. The RNA sequencing data were further validated by RT-qPCR. The functional roles of lncRNA-mRNA interactions demonstrated that these lncRNAs are involved in the signal transduction of secondary messengers, and a disease resistance protein, such as RPP13-like protein 4, is probably involved in aluminum tolerance in XZ16. The current findings significantly contribute to our understanding of the regulatory roles of lncRNAs in aluminum tolerance and extend our knowledge of their importance in plant responses to aluminum stress.

Underground communication: Long non-coding RNA signaling in the plant rhizosphere

Long non-coding RNAs (lncRNAs) have emerged as integral gene-expression regulators underlying plant growth, development, and adaptation. To adapt to the heterogeneous and dynamic rhizosphere, plants use interconnected regulatory mechanisms to optimally fine-tune gene-expression-governing interactions with soil biota, as well as nutrient acquisition and heavy metal tolerance. Recently, high-throughput sequencing has enabled the identification of plant lncRNAs responsive to rhizosphere biotic and abiotic cues. Here, we examine lncRNA biogenesis, classification, and mode of action, highlighting the functions of lncRNAs in mediating plant adaptation to diverse rhizosphere factors. We then discuss studies that reveal the significance and target genes of lncRNAs during developmental plasticity and stress responses at the rhizobium interface. A comprehensive understanding of specific lncRNAs, their regulatory targets, and the intricacies of their functional interaction networks will provide crucial insights into how these transcriptomic switches fine-tune responses to shifting rhizosphere signals. Looking ahead, we foresee that single-cell dissection of cell-type-specific lncRNA regulatory dynamics will enhance our understanding of the precise developmental modulation mechanisms that enable plant rhizosphere adaptation. Overcoming future challenges through multi-omics and genetic approaches will more fully reveal the integral roles of lncRNAs in governing plant adaptation to the belowground environment.

Dynamic transcriptome profiling revealed a key gene ZmJMJ20 and pathways associated with cadmium stress in maize

Cadmium (Cd) pollution in soil poses a global concern due to its serious impacts on human health and ecological security. In plants, tremendous efforts have been made to identify some key genes and pathways in Cd stress responses. However, studies on the roles of epigenetic factors in response to Cd stress were still limited. In the study, we first gain insight into the gene expression dynamics for maize seedlings under 0 h, 12 h, and 72 h Cd stress. As a result, six distinct groups of genes were identified by hierarchical clustering and principal component analysis. The key pathways associated with 12 h Cd stress were protein modifications including protein ubiquitination, signal transduction by protein phosphorylation, and histone modification. Whereas, under 72 h stress, main pathways were involved in biological processes including phenylalanine metabolism, response to oxygen-containing compounds and metal ions. Then to be noted, one of the most highly expressed genes at 12 h under Cd treatment is annotated as histone demethylases (ZmJMJ20). The evolutionary tree analysis and domain analysis showed that ZmJMJ20 belonged to the JmjC-only subfamily of the Jumonji-C (JmjC) family, and ZmJMJ20 was conserved in rice and Arabidopsis. After 72 h of Cd treatment, the zmjmj20 mutant created by EMS treatment manifested less severe chlorosis/leaf yellowing symptoms compared with wild-type plants, and there was no significant difference in Fv/Fm and φPSII value before and after Cd treatment. Moreover, the expression levels of several photosynthesis-related down-regulated genes in EMS mutant plants were dramatically increased compared with those in wild-type plants at 12 h under Cd treatment. Our results suggested that ZmJMJ20 plays an important role in the Cd tolerance response pathway and will facilitate the development of cultivars with improved Cd stress tolerance.

A long non-coding RNA associated with H3K7me3 methylation negatively regulates OsZIP16 transcription under cadmium stress

Cadmium (Cd) is a toxic environmental pollutant, posing a high risk to crop production and human health. However, the genetic mechanisms for regulation of Cd accumulation in crops are poorly understood. In this study, we functionally identified a novel long non-coding RNA (lncRNA, TCONS_00502780) that repressed a locus encoding an uncharacterized metal transporter ZIP16 (ZRT/IRT-like Protein) in rice. LncRNA-OsZIP16 (L16) is resident in the antisense strand of OsZIP16. Both L16 and OsZIP16 were transcriptionally expressed during the life cycle, but under Cd stress the L16 transcription was repressed, whereas the OsZIP16 expression was upregulated. OsZIP16 is localized to the endoplasmic reticulum. Knocking out OsZIP16 by CRISPR-Cas9 (C16) resulted in Cd sensitivity, manifested by reduced plant growth and intense cellular damage with a slightly higher Cd translocation from roots to shoots, suggesting that OsZIP16 expression is required for rice growth and development under Cd stress. Conversely, OsZIP16 constitutive overexpression (OE16) lines displayed a growth phenotype compatible to the wide-type with lower Cd translocation ratio from roots to shoots. L16 knock-down lines by RNA interference (L16-R) showed a similar phenotype to the OE16 lines, while the L16 overexpression (L16-OE) lines were phenotypically similar to the C16 lines. The OsZIP16 transcription was upregulated in the L16-R lines but downregulated in the L16-OE lines. Using an antibody against the trimethylation of histone H3 lysine 27 (H3K27me3) followed by chromatin immunoprecipitation (ChIP), we found the reduced H3K27me3 methylation marks surrounding the OsZIP16 gene under Cd stress. Further examination of H3K27me3 in the L16-R lines revealed that the methylation levels were also significantly lower than those in WT. Taken together, these data suggest that the L16 could negatively regulate the OsZIP16 transcriptional expression through an epigenetic mechanism for rice adaption to Cd stress.

Melatonin enhances cadmium tolerance in rice via long non-coding RNA-mediated modulation of cell wall and photosynthesis

In plants, melatonin (MLT) is a versatile signaling molecule involved in promoting plant development and mitigating the damage caused by heavy metal exposure. Long non-coding RNAs (lncRNAs) are essential components in the plant's response to various abiotic stress, functioning within the gene regulatory network. Here, a hydroponic experiment was performed to explore the involvement of lncRNAs in MLT-mediated amelioration of cadmium (Cd) toxicity in rice plants. The results demonstrated that applying 250 mg L-1 MLT in a solution containing 10 μM Cd leads to an effective reduction of 30.0% in shoot Cd concentration. Remarkably, the treatment resulted in a 21.2% improvement in potassium and calcium uptake, a 164.5% enhancement in net photosynthetic rate, and a 33.2% decrease in malondialdehyde accumulation, resulting increases in plant height, root length, and biomass accumulation. Moreover, a transcriptome analysis revealed 2510 differentially expressed transcripts, including the Cd transporters (-3.82-fold downregulated) and the Cd tolerance-associated genes (1.24-fold upregulated). Notably, regulatory network prediction uncovered 6 differentially expressed lncRNAs that act as competitive endogenous RNA or in RNA complex interactions. These key lncRNAs regulate the expression of target genes that are involved in pectin and cellulose metabolism, scavenging of reactive oxygen species, salicylic acid-mediated defense response, and biosynthesis of brassinosteroids, which ultimately modify the cell wall for Cd adsorption, safeguard photosynthesis, and control hormone signaling to reduce Cd toxicity. Our results unveiled a crucial lncRNA-mediated mechanism underlying MLT's role in Cd detoxification in rice plants, providing potential applications in agricultural practices and environmental remediation.

Long noncoding RNA from Betula platyphylla, BplncSIR1, confers salt tolerance by regulating BpNAC2 to mediate reactive oxygen species scavenging and stomatal movement

Long noncoding RNAs (lncRNAs) play an important role in abiotic stress tolerance. However, their function in conferring abiotic stress tolerance is still unclear. Herein, we characterized the function of a salt-responsive nuclear lncRNA (BplncSIR1) from Betula platyphylla (birch). Birch plants overexpressing and knocking out for BplncSIR1 were generated. BplncSIR1 was found to improve salt tolerance by inducing antioxidant activity and stomatal closure, and also accelerate plant growth. Chromatin isolation by RNA purification (ChIRP) combined with RNA sequencing indicated that BplncSIR1 binds to the promoter of BpNAC2 (encoding NAC domain-containing protein 2) to activate its expression. Plants overexpressing and knocking out for BpNAC2 were generated. Consistent with that of BplncSIR1, overexpression of BpNAC2 also accelerated plant growth and conferred salt tolerance. In addition, BpNAC2 binds to different cis-acting elements, such as G-box and 'CCAAT' sequences, to regulate the genes involved in salt tolerance, resulting in reduced ROS accumulation and decreased water loss rate by stomatal closure. Taken together, BplncSIR1 serves as the regulator of BpNAC2 to induce its expression in response to salt stress, and activated BpNAC2 accelerates plant growth and improves salt tolerance. Therefore, BplncSIR1 might be a candidate gene for molecular breeding to cultivate plants with both a high growth rate and improved salt tolerance.

Integrated physio-biochemical and RNA sequencing analysis revealed mechanisms of long non-coding RNA-mediated response to cadmium toxicity in wheat

The yield and quality of wheat (Triticum aestivum L.) is seriously affected by soil cadmium (Cd), a hazardous material to plant and human health. Long non-coding RNAs (lncRNAs) of plants are shown actively involved in response to various biotic and abiotic stresses by mediating the gene regulatory networks. However, the functions of lncRNAs in wheat against Cd stress are still obscure. Using deep strand-specific RNA sequencing, 10,044 confident novel lncRNAs in wheat roots response to Cd stress were identified. It was found that 69 lncRNA-target pairs referred to cis-acting regulation and impacted the expressions of their neighboring genes involving in Cd transport and detoxification, photosynthesis, and antioxidant defense. These findings were positively corelated with the physio-biochemical results, i.e. Cd stress affected Cd accumulation, photosynthesis system and ROS in wheat. Overexpression of lncRNA37228 (targeted to a photosystem II protein D1 coding gene), resulted in enhancing Arabidopsis thaliana resistance against Cd stress. By genome-wide identification and characterization, the possible functions of photosystem II protein gene family in wheat under Cd condition were illustrated. Our findings provide novel knowledge into the molecular mechanisms of lncRNAs-regulated wheat tolerance to Cd toxicity and lay foundations for the further studies concerning lncRNAs in food safety production.

Genome-wide identification and characterization of long non-coding RNA in barley roots in response to Piriformospora indica colonization

Currently, there is very limited information about long noncoding RNAs (lncRNAs) found in barley. It remains unclear whether barley lncRNAs are responsive to Piriformospora indica (P. indica) colonization.We found that barley roots exhibited fast development and that large roots branched after P. indica colonization. Genome-wide high-throughput RNA-seq and bioinformatic analysis showed that 4356 and 5154 differentially expressed LncRNAs (DELs) were found in response to P. indica at 3 and 7 days after colonization (dai), respectively, and 2456 DELs were found at 7 dai compared to 3 dai. Based on the coexpression correlation of lncRNAmRNA, we found that 98.6% of lncRNAs were positively correlated with 3430 mRNAs at 3 dai and 7 dai. Further GO analysis showed that 30 lncRNAs might be involved in the regulation of gene transcription; 23 lncRNAs might participate in cell cycle regulation. Moreover, the metabolite analysis indicated that chlorophyll a, sucrose, protein, gibberellin, and auxin were in accordance with the results of the transcriptome, and the respective lncRNAs were positively correlated with these target RNAs. Gene silencing suggested that lncRNA TCONS_00262342 is probably a key regulator of GA₃ synthesis pathway, which participates in P. indica and barley interactions. We concluded that acting as a molecular material basis and resource, lncRNAs respond to P. indica colonization by regulating metabolite content in barley and coordinate the complex regulatory process of higher life by constructing highly positive correlations with their target mRNAs.

ASLncR: a novel computational tool for prediction of abiotic stress-responsive long non-coding RNAs in plants

Abiotic stresses are detrimental to plant growth and development and have a major negative impact on crop yields. A growing body of evidence indicates that a large number of long non-coding RNAs (lncRNAs) are key to many abiotic stress responses. Thus, identifying abiotic stress-responsive lncRNAs is essential in crop breeding programs in order to develop crop cultivars resistant to abiotic stresses. In this study, we have developed the first machine learning-based computational model for predicting abiotic stress-responsive lncRNAs. The lncRNA sequences which were responsive and non-responsive to abiotic stresses served as the two classes of the dataset for binary classification using the machine learning algorithms. The training dataset was created using 263 stress-responsive and 263 non-stress-responsive sequences, whereas the independent test set consists of 101 sequences from both classes. As the machine learning model can adopt only the numeric data, the Kmer features ranging from sizes 1 to 6 were utilized to represent lncRNAs in numeric form. To select important features, four different feature selection strategies were utilized. Among the seven learning algorithms, the support vector machine (SVM) achieved the highest cross-validation accuracy with the selected feature sets. The observed 5-fold cross-validation accuracy, AU-ROC, and AU-PRC were found to be 68.84, 72.78, and 75.86%, respectively. Furthermore, the robustness of the developed model (SVM with the selected feature) was evaluated using an independent test dataset, where the overall accuracy, AU-ROC, and AU-PRC were found to be 76.23, 87.71, and 88.49%, respectively. The developed computational approach was also implemented in an online prediction tool ASLncR accessible at https://iasri-sg.icar.gov.in/aslncr/ . The proposed computational model and the developed prediction tool are believed to supplement the existing effort for the identification of abiotic stress-responsive lncRNAs in plants.

Long Non-Coding RNAs of Plants in Response to Abiotic Stresses and Their Regulating Roles in Promoting Environmental Adaption

Abiotic stresses triggered by climate change and human activity cause substantial agricultural and environmental problems which hamper plant growth. Plants have evolved sophisticated mechanisms in response to abiotic stresses, such as stress perception, epigenetic modification, and regulation of transcription and translation. Over the past decade, a large body of literature has revealed the various regulatory roles of long non-coding RNAs (lncRNAs) in the plant response to abiotic stresses and their irreplaceable functions in environmental adaptation. LncRNAs are recognized as a class of ncRNAs that are longer than 200 nucleotides, influencing a variety of biological processes. In this review, we mainly focused on the recent progress of plant lncRNAs, outlining their features, evolution, and functions of plant lncRNAs in response to drought, low or high temperature, salt, and heavy metal stress. The approaches to characterize the function of lncRNAs and the mechanisms of how they regulate plant responses to abiotic stresses were further reviewed. Moreover, we discuss the accumulating discoveries regarding the biological functions of lncRNAs on plant stress memory as well. The present review provides updated information and directions for us to characterize the potential functions of lncRNAs in abiotic stresses in the future.

Transcriptional Regulatory Network of Plant Cadmium Stress Response

Cadmium (Cd) is a non-essential heavy metal with high toxicity to plants. Plants have acquired specialized mechanisms to sense, transport, and detoxify Cd. Recent studies have identified many transporters involved in Cd uptake, transport, and detoxification. However, the complex transcriptional regulatory networks involved in Cd response remain to be elucidated. Here, we provide an overview of current knowledge regarding transcriptional regulatory networks and post-translational regulation of the transcription factors involved in Cd response. An increasing number of reports indicate that epigenetic regulation and long non-coding and small RNAs are important in Cd-induced transcriptional responses. Several kinases play important roles in Cd signaling that activate transcriptional cascades. We also discuss the perspectives to reduce grain Cd content and improve crop tolerance to Cd stress, which provides a theoretical reference for food safety and the future research of plant varieties with low Cd accumulation.

Genome-wide characterization and functional identification of MYB genes in Malus sieversii infected by Valsa mali

Among the most important transcription factors in plants, the v-myb avian myeloblastosis viral oncogene homolog (MYB) regulates the expression network of response genes under stresses such as fungal infection. In China, the canker disease Valsa mali threatens the survival of Malus sieversii, an ancestor of cultivated apples. Using the M. sieversii genome, we identified 457 MsMYB and 128 R2R3-MsMYB genes that were randomly distributed across 17 chromosomes. Based on protein sequence and structure, the R2R3-MsMYB genes were phylogenetically divided into 29 categories, and 26 conserved motifs were identified. We further predicted cis-elements in the 2000-kb promoter region of R2R3-MsMYBs based on the genome. Transcriptome analysis of M. sieversii under V. mali infection showed that 27 R2R3-MsMYBs were significantly differentially expressed, indicating their key role in the response to V. mali infection. Using transient transformation, MsMYB14, MsMYB24, MsMYB39, MsMYB78, and MsMYB108, which were strongly induced by V. mali infection, were functionally identified. Among the five MsMYBs, MsMYB14 and MsMYB78 were both important in enhancing resistance to diseases, whereas MsMYB24 inhibited resistance. Based on the results of this study, we gained a better understanding of the MsMYB transcription factor family and laid the foundation for a future research program on disease prevention strategies in M. sieversii.

A novel long noncoding RNA CIL1 enhances cold stress tolerance in Arabidopsis

With the intensification of global warming, extreme weather events have occurred more frequently, among which cold stress has become one of the major environmental factors that restrict global crop yield and production. Multiple long noncoding RNAs (lncRNAs) have been predicted or recognized in the plant response to cold stress, however, the molecular biological functions of most of these RNAs are still poorly understood. Here, we identified a novel lncRNA, COLD INDUCED lncRNA 1 (CIL1), as a positive regulator of the plant response to cold stress in Arabidopsis. CIL1 was significantly induced when the plant was exposed to cold stress. Moreover, knockdown mutants showed more sensitivity to cold stress than the wild type did, accompanied by an increased content of endogenous ROS (reactive oxygen species) and reduced osmoregulatory substances. Genome-wide transcriptome analysis indicated that 256 genes were downregulated and 34 genes were upregulated in cil1 mutants under cold stress, which were mainly involved in hormone signal transduction, ROS homeostasis and glucose metabolism. Our study implies that CIL1 has a positive effect on the plant response to cold stress by regulating the expression of multiple stress-related genes during the seedling stage.

LncRNA PMAT-PtoMYB46 module represses PtoMATE and PtoARF2 promoting Pb2+ uptake and plant growth in poplar

Lead (Pb²⁺) is one of the most toxic heavy-metal contaminants. Fast-growing woody plants with substantial biomass are ideal for bioremediation. However, the transcriptional regulation of Pb²⁺ uptake in woody plants remains unclear. Here, we identified 226 Pb²⁺-induced, differentially expressed long non-coding RNAs (DELs) in Populus tomentosa. Functional annotation revealed that these DELs mainly regulate carbon metabolism, biosynthesis of secondary metabolites, energy metabolism, and signal transduction through their potential target genes. Association and epistasis analysis showed that the lncRNA PMAT (Pb²⁺-induced multidrug and toxic compound extrusion (MATE) antisense lncRNA) interacts epistatically with PtoMYB46 to regulate leaf dry weight, photosynthesis rate, and transketolase activity. Genetic transformation and molecular assays showed that PtoMYB46 reduces the expression of PtoMATE directly or indirectly through PMAT, thereby reducing the secretion of citric acid (CA) and ultimately promoting Pb²⁺ uptake. Meanwhile, PtoMYB46 targets auxin response factor 2 (ARF2) and reduces its expression, thus positively regulating plant growth. We concluded that the PMAT-PtoMYB46-PtoMATE-PtoARF2 regulatory module control Pb²⁺ tolerance, uptake, and plant growth. This study demonstrates the involvement of lncRNAs in response to Pb²⁺ in poplar, yielding new insight into the potential for developing genetically improved woody plant varieties for phytoremediating lead-contaminated soils.

The Intersection of Non-Coding RNAs Contributes to Forest Trees' Response to Abiotic Stress

Non-coding RNAs (ncRNAs) play essential roles in plants by modulating the expression of genes at the transcriptional or post-transcriptional level. In recent years, ncRNAs have been recognized as crucial regulators for growth and development in forest trees, and ncRNAs that respond to various abiotic stresses are now under intense study. In this review, we summarized recent advances in the understanding of abiotic stress-responsive microRNAs (miRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs) in forest trees. Furthermore, we analyzed the intersection of miRNAs, and epigenetic modified ncRNAs of forest trees in response to abiotic stress. In particular, the abiotic stress-related lncRNA/circRNA-miRNA-mRNA regulatory network of forest trees was explored.

Transcriptomic Analysis Reveals LncRNAs Associated with Flowering of Angelica sinensis during Vernalization

Angelica sinensis is a “low-temperature and long-day” perennial plant that produces bioactive compounds such as phthalides, organic acids, and polysaccharides for various types of clinical agents, including those with cardio-cerebrovascular, hepatoprotective, and immunomodulatory effects. To date, the regulatory mechanism of flowering under the photoperiod has been revealed, while the regulatory network of flowering genes during vernalization, especially in the role of lncRNAs, has yet to be identified. Here, lncRNAs associated with flowering were identified based on the full-length transcriptomic analysis of A. sinensis at vernalization and freezing temperatures, and the coexpressed mRNAs of lncRNAs were validated by qRT-PCR. We obtained a total of 2327 lncRNAs after assessing the protein-coding potential of coexpressed mRNAs, with 607 lncRNAs aligned against the TAIR database of model plant Arabidopsis, 345 lncRNAs identified, and 272 lncRNAs characterized on the SwissProt database. Based on the biological functions of coexpressed mRNAs, the 272 lncRNAs were divided into six categories: (1) chromatin, DNA/RNA and protein modification; (2) flowering; (3) stress response; (4) metabolism; (5) bio-signaling; and (6) energy and transport. The differential expression levels of representatively coexpressed mRNAs were almost consistent with the flowering of A. sinensis. It can be concluded that the flowering of A. sinensis is positively or negatively regulated by lncRNAs, which provides new insights into the regulation mechanism of the flowering of A. sinensis.

The genome and preliminary single-nuclei transcriptome of Lemna minuta reveals mechanisms of invasiveness

The ability to trace every cell in some model organisms has led to the fundamental understanding of development and cellular function. However, in plants the complexity of cell number, organ size, and developmental time makes this a challenge even in the diminutive model plant Arabidopsis (Arabidopsis thaliana). Duckweed, basal nongrass aquatic monocots, provide an opportunity to follow every cell of an entire plant due to their small size, reduced body plan, and fast clonal growth habit. Here we present a chromosome-resolved genome for the highly invasive Lesser Duckweed (Lemna minuta) and generate a preliminary cell atlas leveraging low cell coverage single-nuclei sequencing. We resolved the 360 megabase genome into 21 chromosomes, revealing a core nonredundant gene set with only the ancient tau whole-genome duplication shared with all monocots, and paralog expansion as a result of tandem duplications related to phytoremediation. Leveraging SMARTseq2 single-nuclei sequencing, which provided higher gene coverage yet lower cell count, we profiled 269 nuclei covering 36.9% (8,457) of the L. minuta transcriptome. Since molecular validation was not possible in this nonmodel plant, we leveraged gene orthology with model organism single-cell expression datasets, gene ontology, and cell trajectory analysis to define putative cell types. We found that the tissue that we computationally defined as mesophyll expressed high levels of elemental transport genes consistent with this tissue playing a role in L. minuta wastewater detoxification. The L. minuta genome and preliminary cell map provide a paradigm to decipher developmental genes and pathways for an entire plant.

Integrated analysis of long non-coding RNAs and mRNAs reveals the regulatory network of maize seedling root responding to salt stress

Background

Long non-coding RNAs (lncRNAs) play important roles in response to abiotic stresses in plants, by acting as cis- or trans-acting regulators of protein-coding genes. As a widely cultivated crop worldwide, maize is sensitive to salt stress particularly at the seedling stage. However, it is unclear how the expressions of protein-coding genes are affected by non-coding RNAs in maize responding to salt tolerance.

Results

The whole transcriptome sequencing was employed to investigate the differential lncRNAs and target transcripts responding to salt stress between two maize inbred lines with contrasting salt tolerance. We developed a flexible, user-friendly, and modular RNA analysis workflow, which facilitated the identification of lncRNAs and novel mRNAs from whole transcriptome data. Using the workflow, 12,817 lncRNAs and 8,320 novel mRNAs in maize seedling roots were identified and characterized. A total of 742 lncRNAs and 7,835 mRNAs were identified as salt stress-responsive transcripts. Moreover, we obtained 41 cis- and 81 trans-target mRNA for 88 of the lncRNAs. Among these target transcripts, 11 belonged to 7 transcription factor (TF) families including bHLH, C2H2, Hap3/NF-YB, HAS, MYB, WD40, and WRKY. The above 8,577 salt stress-responsive transcripts were further classified into 28 modules by weighted gene co-expression network analysis. In the salt-tolerant module, we constructed an interaction network containing 79 nodes and 3081 edges, which included 5 lncRNAs, 18 TFs and 56 functional transcripts (FTs). As a trans-acting regulator, the lncRNA MSTRG.8888.1 affected the expressions of some salt tolerance-relative FTs, including protein-serine/threonine phosphatase 2C and galactinol synthase 1, by regulating the expression of the bHLH TF.

Conclusions

The contrasting genetic backgrounds of the two inbred lines generated considerable variations in the expression abundance of lncRNAs and protein-coding transcripts. In the co-expression networks responding to salt stress, some TFs were targeted by the lncRNAs, which further regulated the salt tolerance-related functional transcripts. We constructed a regulatory pathway of maize seedlings to salt stress, which was mediated by the hub lncRNA MSTRG.8888.1 and participated by the bHLH TF and its downstream target transcripts. Future work will be focused on the functional revelation of the regulatory pathway.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08286-7.

Biologia Futura: progress and future perspectives of long non-coding RNAs in forest trees

Forest trees are affected by climate change, anthropogenic pressure, as well as abiotic and biotic stresses. Conventional tree breeding has so far been limited to enhance overall productivity, and our understanding of the genetic basis of quantitative traits is still inadequate. Quantum leaps in next-generation sequencing technologies and bioinformatics have permitted the exploration and identification of various non-coding regions of the genome other than protein coding genes. These genomic regions produce various types of non-coding RNAs and regulate myriads of biological functions at epigenetic, transcriptional and translational levels. Recently, long non-coding RNAs (lncRNAs) which act as molecular switch have been identified to be pivotal molecules in forest trees. This review focuses on progress made in regulatory mechanisms in various developmental phases like wood formation, adventitious rooting and flowering and stress responses. It was predicted that complex regulatory interactions among lncRNA, miRNA and gene exist. LncRNAs can function as a sponge for miRNAs, reducing the suppressive effect of miRNAs on target mRNAs and perhaps adding a new layer of regulatory interactions among non-coding RNA classes in trees. Furthermore, network analysis revealed the interactions of lncRNA and genes during the expression of several important genes. The insights generated about lncRNAs in forest trees would enable improvement of economically important traits including the devastating abiotic and biotic stresses. In addition, solid understanding on the wide range of regulatory functions of lncRNAs on traits influencing biomass productivity and adaptation would aid the applications of biotechnology in genetic improvement of forest trees.

LncRNAs as Therapeutic Targets and Potential Biomarkers for Lipid-Related Diseases

Lipid metabolism is an essential biological process involved in nutrient adjustment, hormone regulation, and lipid homeostasis. An irregular lifestyle and long-term nutrient overload can cause lipid-related diseases, including atherosclerosis, myocardial infarction (MI), obesity, and fatty liver diseases. Thus, novel tools for efficient diagnosis and treatment of dysfunctional lipid metabolism are urgently required. Furthermore, it is known that lncRNAs based regulation like sponging microRNAs (miRNAs) or serving as a reservoir for microRNAs play an essential role in the progression of lipid-related diseases. Accordingly, a better understanding of the regulatory roles of lncRNAs in lipid-related diseases would provide the basis for identifying potential biomarkers and therapeutic targets for lipid-related diseases. This review highlighted the latest advances on the potential biomarkers of lncRNAs in lipid-related diseases and summarised current knowledge on dysregulated lncRNAs and their potential molecular mechanisms. We have also provided novel insights into the underlying mechanisms of lncRNAs which might serve as potential biomarkers and therapeutic targets for lipid-related diseases. The information presented here may be useful for designing future studies and advancing investigations of lncRNAs as biomarkers for diagnosis, prognosis, and therapy of lipid-related diseases.