Integrative lncRNA landscape reveals lncRNA-coding gene networks in the secondary cell wall biosynthesis pathway of moso bamboo (Phyllostachys edulis)

Expression and function analysis of phenylalanine ammonia-lyase genes involved in Bamboo lignin biosynthesis

Bamboo, renowned as the fastest-growing plant globally, matures within an astonishingly short period of 40-50 days from shoots, reaching heights of 10-20 meters. Moreover, it can be harvested for various uses within 3-5 years. Bamboo exhibits exceptional mechanical properties, characterized by high hardness and flexibility, largely attributed to its lignin content. Phenylalanine ammonia-lyase (PAL) catalyzes the crucial initial step in lignin biosynthesis, but its precise role in bamboo lignification processes remains elusive. Thus, elucidating the functions of PAL genes in bamboo lignification processes is imperative for understanding its rapid growth and mechanical strength. Here, we systematically identified and classified PAL genes in Moso bamboo, ensuring nomenclature consistency across prior studies. Subsequently, we evaluated PAL gene expression profiles using publicly available transcriptome data. The downregulation of PePALs expression in Moso bamboo through in planta gene editing resulted in a decrease in PAL activity and a subsequent reduction in lignin content. In contrast, overexpression of PePAL led to enhanced PAL activity and an increase in lignin content. These findings highlight the critical role of PAL in the lignin biosynthesis process of Moso bamboo, which will help to unravel the mechanism underpinning bamboo's rapid growth and mechanical strength, with a specific emphasis on elucidating the functions of PAL genes.

Noncoding RNAs in regulation of plant secondary metabolism

Plant secondary metabolites (PSMs) are a large class of structurally diverse molecules, mainly consisting of terpenoids, phenolic compounds, and nitrogen-containing compounds, which play active roles in plant development and stress responses. The biosynthetic processes of PSMs are governed by a sophisticated regulatory network at multiple levels. Noncoding RNAs (ncRNAs) such as microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs) may serve as post-transcriptional regulators for plant secondary metabolism through acting on genes encoding either transcription factors or participating enzymes in relevant metabolic pathways. High-throughput sequencing technologies have facilitated the large-scale identifications of ncRNAs potentially involved in plant secondary metabolism in model plant species as well as certain species with enriched production of specific types of PSMs. Moreover, a series of miRNA-target modules have been functionally characterized to be responsible for regulating PSM biosynthesis and accumulation in plants under abiotic or biotic stresses. In this review, we will provide an overview of current findings on the ncRNA-mediated regulation of plant secondary metabolism with special attention to its participation in plant stress responses, and discuss possible issues to be addressed in future fundamental research and breeding practice.

Jasmonic acid's impact on Sedum alfredii growth and cadmium tolerance: A physiological and transcriptomic study

Soil cadmium (Cd) pollution is escalating, necessitating effective remediation strategies. This study investigated the effects of exogenous jasmonic acid (JA) on Sedum alfredii Hance under Cd stress, aiming to enhance its phytoextraction efficiency. Initially, experiments were conducted to assess the impact of various concentrations of JA added to environments with Cd concentrations of 100, 300, and 500 μmol/L. The results determined that a concentration of 1 μmol/L JA was optimal. This concentration effectively mitigated the level of ROS products by enhancing the activity of antioxidant enzymes. Additionally, JA fostered Cd absorption and accumulation, while markedly improving plant biomass and photosynthetic performance. In further experiments, treatment with 1 μmol/L JA under 300 μmol/L Cd stress was performed and transcriptomic analysis unveiled a series of differentially expressed genes (DEGs) instrumental in the JA-mediated Cd stress response. These DEGs encompass not only pathways of JA biosynthesis and signaling but also genes encoding functions that influence antioxidant systems and photosynthesis, alongside genes pertinent to cell wall synthesis, and metal chelation and transport. This study highlights that JA treatment significantly enhances S. alfredii's Cd tolerance and accumulation, offering a promising strategy for plant remediation and deepening our understanding of plant responses to heavy metal stress.

Deciphering shared attributes of plant long non-coding RNAs through a comparative computational approach

Over the past decade, long non-coding RNA (lncRNA), which lacks protein-coding potential, has emerged as an essential regulator of the genome. The present study examined 13,599 lncRNAs in Arabidopsis thaliana, 11,565 in Oryza sativa, and 32,397 in Zea mays for their characteristic features and explored the associated genomic and epigenomic features. We found lncRNAs were distributed throughout the chromosomes and the Helitron family of transposable elements (TEs) enriched, while the terminal inverted repeat depleted in lncRNA transcribing regions. Our analyses determined that lncRNA transcribing regions show rare or weak signals for most epigenetic marks except for H3K9me2 and cytosine methylation in all three plant species. LncRNAs showed preferential localization in the nucleus and cytoplasm; however, the distribution ratio in the cytoplasm and nucleus varies among the studied plant species. We identified several conserved endogenous target mimic sites in the lncRNAs among the studied plants. We found 233, 301, and 273 unique miRNAs, potentially targeting the lncRNAs of A. thaliana, O. sativa, and Z. mays, respectively. Our study has revealed that miRNAs, which interact with lncRNAs, target genes that are involved in a diverse array of biological and molecular processes. The miRNA-targeted lncRNAs displayed a strong affinity for several transcription factors, including ERF and BBR-BPC, mutually present in all three plants, advocating their conserved functions. Overall, the present study showed that plant lncRNAs exhibit conserved genomic and epigenomic characteristics and potentially govern the growth and development of plants.

Inhibition of DNA and RNA methylation disturbs root development of moso bamboo

DNA methylation (5mC) and N6-methyladenosine (m6A) are two important epigenetics regulators, which have a profound impact on plant growth development. Phyllostachys edulis (P. edulis) is one of the fastest spreading plants due to its well-developed root system. However, the association between 5mC and m6A has seldom been reported in P. edulis. In particular, the connection between m6A and several post-transcriptional regulators remains uncharacterized in P. edulis. Here, our morphological and electron microscope observations showed the phenotype of increased lateral root under RNA methylation inhibitor (DZnepA) and DNA methylation inhibitor (5-azaC) treatment. RNA epitranscriptome based on Nanopore direct RNA sequencing revealed that DZnepA treatment exhibits significantly decreased m6A level in the 3'-untranslated region (3'-UTR), which was accompanied by increased gene expression, full-length ratio, higher proximal poly(A) site usage and shorter poly(A) tail length. DNA methylation levels of CG and CHG were reduced in both coding sequencing and transposable element upon 5-azaC treatment. Cell wall synthesis was impaired under methylation inhibition. In particular, differentially expressed genes showed a high percentage of overlap between DZnepA and 5-azaC treatment, which suggested a potential correlation between two methylations. This study provides preliminary information for a better understanding of the link between m6A and 5mC in root development of moso bamboo.

Identification and characterization of long noncoding RNAs in two contrasting olive (Olea europaea L.) genotypes subjected to aluminum toxicity

Aluminum (Al) toxcity is considered to be the primary factor limiting crop productivity in acidic soil. Many studies indicate that long non-coding RNAs (lncRNAs) fulfil a crucial role in plant growth and responses to different abiotic stress. However, identification and characterization of lncRNAs responsive to Al stress at a genome-wide level in olive tree is still lacking. Here, we performed comparative analysis on lncRNA transcriptome between Zhonglan (an Al-tolerant genotype) and Frantoio selezione (Al-sensitive) responding to Al exposure. A total of 19,498 novel lncRNAs were identified from both genotypes, and 6900 lncRNA-target pairs were identified as cis-acting and 2311 supposed to be trans-acting. Among them, 2076 lncRNAs were appraised as Al tolerance-associated lncRNAs due to their distinctly genotype-specific expression profiles under Al exposure. Target prediction and functional analyses revealed several key lncRNAs are related to genes encoding pectinesterases, xyloglucan endotransglucosylase/hydrolase, WRKY and MYB transcription factors, which mainly participate in the modification of cell wall for Al tolerance. Furthermore, gene co-expression network analysis showed 8 lncRNA-mRNA-miRNA modules participate in transcriptional regulation of downstream Al resistant genes. Our findings increased our understanding about the function of lncRNAs in responding to Al stress in olive and identified potential promising lncRNAs for further investigation.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

A long non-coding RNA PelncRNA1 is involved in Phyllostachys edulis response to UV-B stress

Phyllostachys edulis (moso bamboo) is China's most widespread bamboo species, with significant economic and ecological values. Long non-coding RNA (lncRNA) is a type of regulatory RNA that is longer than 200 nucleotides and incapable of encoding proteins, and is frequently involved in regulating biotic and abiotic stress and plant development. However, the biological functions of lncRNA in moso bamboo are unknown. In this study, a lncRNA (named PelncRNA1) differentially expressed following UV-B treatment was discovered in the whole transcriptome sequencing database of moso bamboo. The target genes were filtered and defined by correlation analysis of PelncRNA1 and gene expression pattern. The expression levels of PelncRNA1 and its target genes were verified using qRT-PCR. The results demonstrated that the expression levels of PelncRNA1 and its target genes increased during UV-B treatment. In Arabidopsis transgenic seedlings and moso bamboo protoplasts, PelncRNA1 was discovered to influence the expression of its target genes when overexpressed. In addition, transgenic Arabidopsis showed higher tolerance to UV-B stress. These results suggest that PelncRNA1 and its target genes are involved in the response of moso bamboo to UV-B stress. The novel findings would contribute to our understanding of how lncRNAs regulate the response to abiotic stresses in moso bamboo.

Identification and differential analysis of noncoding RNAs in response to drought in Phyllostachys aureosulcata f. spectabilis

The role of noncoding RNAs (ncRNAs) in plant resistance to abiotic stresses is increasingly being discovered. Drought stress is one of the most common stresses that affecting plant growth, and high intensity drought has a significant impact on the normal growth of plants. In this study, a high-throughput sequencing was performed on plant tissue samples of Phyllostachys aureosulcata f. spectabilis C. D. Chu et C. S. Chao by drought treatment for 0, 2, 4 and 6 days. The sequencing results were analysed bioinformatically. We detected 336,946 RNAs among all 12 samples, including 192,098 message RNAs (mRNAs), 142,761 long noncoding RNAs (lncRNAs), 1,670 circular RNAs (circRNAs), and 417 microRNAs (miRNAs). We detected 2,419 differentially expressed (DE) ncRNAs, including 213 DE circRNAs, 2,088 DE lncRNAs and 118 DE miRNAs. Then, we used Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) to functionally predict DE ncRNAs. The results showed that most DE ncRNAs are involved in the response to drought stress, mainly in biochemical reactions involved in some metabolites, as well as in organelle activities. In addition, we validated two random circRNAs and demonstrated their circularity. We also found a stable internal reference gene available for Phyllostachys aureosulcata f. spectabilis and validated the accuracy of this experiment by quantitative real-time polymerase chain reaction (qRT-PCR).

Systematic Investigation of Immune-Related lncRNA Landscape Reveals a Potential Long Non-Coding RNA Signature for Predicting Prognosis in Renal Cell Carcinoma

Background: Renal cell carcinoma (RCC) is the predominant type of malignant tumor in kidney cancer. Finding effective biomarkers, particularly those based on the tumor immune microenvironments (TIME), is critical for the prognosis and diagnosis of RCC. Increasing evidence has revealed that long non-coding RNAs (lncRNAs) play a crucial role in cancer immunity. However, the comprehensive landscape of immune infiltration-associated lncRNAs and their potential roles in the prognosis and diagnosis of RCC remain largely unexplored.

Methods: Based on transcriptomic data of 261 RCC samples, novel lncRNAs were identified using a custom pipeline. RCC patients were classified into different immune groups using unsupervised clustering algorithms. Immune-related lncRNAs were obtained according to the immune status of RCC. Competing endogenous RNAs (ceRNA) regulation network was constructed to reveal their functions. Expression patterns and several tools such as miRanda, RNAhybrid, miRWalk were used to define lncRNAs-miRNAs-mRNAs interactions. Univariate Cox, LASSO, and multivariate Cox regression analyses were performed on the training set to construct a tumorigenesis-immune-infiltration-related (TIR)-lncRNA signature for predicting the prognosis of RCC. Independent datasets involving 531 RCC samples were used to validate the TIR-lncRNA signature.

Results: Tens of thousands of novel lncRNAs were identified in RCC samples. Comparing tumors with controls, 1,400 tumorigenesis-related (TR)-lncRNAs, 1269 TR-mRNAs, and 192 TR-miRNAs were obtained. Based on the infiltration of immune cells, RCC patients were classified into three immune clusters. By comparing immune-high with immune-low groups, 241 TIR-lncRNAs were identified, many of which were detected in urinary samples. Based on lncRNA-miRNA-mRNA interactions, we constructed a ceRNA network, which included 25 TR-miRNAs, 28 TIR-lncRNAs, and 66 TIR-mRNAs. Three TIR lncRNAs were identified as a prognostic signature for RCC. RCC patients in the high-risk group exhibited worse OS than those in the low-risk group in the training and testing sets (p < 0.01). The AUC was 0.9 in the training set. Univariate and multivariate Cox analyses confirmed that the TIR-lncRNA signature was an independent prognostic factor in the training and testing sets.

Conclusion: Based on the constructed immune-related lncRNA landscape, 241 TIR-lncRNAs were functionally characterized, three of which were identified as a novel TIR-lncRNA signature for predicting the prognosis of RCC.