Identification of cucumber circular RNAs responsive to salt stress

PlantCircRNA: a comprehensive database for plant circular RNAs

Circular RNAs (circRNAs) represent recently discovered novel regulatory non-coding RNAs. While they are present in many eukaryotes, there has been limited research on plant circRNAs. We developed PlantCircRNA (https://plant.deepbiology.cn/PlantCircRNA/) to fill this gap. The two most important features of PlantCircRNA are (i) it incorporates circRNAs from 94 plant species based on 39 245 RNA-sequencing samples and (ii) it imports the original AtCircDB and CropCircDB databases. We manually curated all circRNAs from published articles, and imported them into the database. Furthermore, we added detailed information of tissue as well as abiotic stresses to the database. To help users understand these circRNAs, the database includes a detection score to measure their consistency and a naming system following the guidelines recently proposed for eukaryotes. Finally, we developed a comprehensive platform for users to visualize, analyze, and download data regarding specific circRNAs. This resource will serve as a home for plant circRNAs and provide the community with unprecedented insights into these mysterious molecule.

Advances in CircRNAs in the Past Decade: Review of CircRNAs Biogenesis, Regulatory Mechanisms, and Functions in Plants

Circular RNA (circRNA) is a type of non-coding RNA with multiple biological functions. Whole circRNA genomes in plants have been identified, and circRNAs have been demonstrated to be widely present and highly expressed in various plant tissues and organs. CircRNAs are highly stable and conserved in plants, and exhibit tissue specificity and developmental stage specificity. CircRNAs often interact with other biomolecules, such as miRNAs and proteins, thereby regulating gene expression, interfering with gene function, and affecting plant growth and development or response to environmental stress. CircRNAs are less studied in plants than in animals, and their regulatory mechanisms of biogenesis and molecular functions are not fully understood. A variety of circRNAs in plants are involved in regulating growth and development and responding to environmental stress. This review focuses on the biogenesis and regulatory mechanisms of circRNAs, as well as their biological functions during growth, development, and stress responses in plants, including a discussion of plant circRNA research prospects. Understanding the generation and regulatory mechanisms of circRNAs is a challenging but important topic in the field of circRNAs in plants, as it can provide insights into plant life activities and their response mechanisms to biotic or abiotic stresses as well as new strategies for plant molecular breeding and pest control.

CircRNA: An emerging star in plant research: A review

CircRNAs are a class of covalently closed non-coding RNA formed by linking the 5' terminus and the 3' terminus after reverse splicing. CircRNAs are widely found in eukaryotes, and they are highly conserved, with spatio-temporal expression specificity and stability. CircRNAs can act as miRNA sponges to regulate the expression of downstream target genes, regulating the transcription of parental genes and some can even be translated into peptides or proteins. Research on circRNAs in plants is still in its infancy compared to that in animals. With the deepening of research, the results of a variety of plant circRNAs suggest that they play an important role in growth and development, and tolerance towards abiotic stresses such as salt, drought, low temperature, high temperature and other adverse environments. In this review paper, we elaborated the molecular characteristics, mechanism of action, function and bioinformatics databases of plant circRNAs, combined with the progress of circRNA research in animals, discussed the potential mechanism of action of plant circRNAs, and proposed the unsolved problems and prospects for future application of plant circRNAs.

Whole-transcriptome profiling and identification of cold tolerance-related ceRNA networks in japonica rice varieties

Introduction

Low-temperature stress negatively impacts rice yield, posing a significant risk to food security. While previous studies have explored the physiological and linear gene expression alterations in rice under low-temperature conditions, the changes in competing endogenous RNA (ceRNA) networks remain largely unexamined.

Methods

We conducted RNA sequencing on two japonica rice varieties with differing cold-tolerance capabilities to establish ceRNA networks. This enabled us to investigate the transcriptional regulatory network and molecular mechanisms that rice employs in response to low-temperature stress.

Results

We identified 364 differentially expressed circular RNAs (circRNAs), 224 differentially expressed microRNAs (miRNAs), and 12,183 differentially expressed messenger RNAs (mRNAs). WRKY family was the most prominent transcription factor family involved in cold tolerance. Based on the expression patterns and targeted relationships of these differentially expressed RNAs, we discerned five potential ceRNA networks related to low-temperature stress in rice: osa-miR166j-5p from the miR166 family was associated with cold tolerance; osa-miR528-3p and osa-miR156j-3p were linked to stress response; and osa-miR156j-3p was involved in the antioxidant system. In addition, Os03g0152000 in the antioxidant system, as well as Os12g0491800 and Os05g0381400, correlated with the corresponding stress response and circRNAs in the network. A gene sequence difference analysis and phenotypic validation of Os11g0685700 (OsWRKY61) within the WRKY family suggested its potential role in regulating cold tolerance in rice.

Discussion and conclusion

We identified Os11g0685700 (OsWRKY61) as a promising candidate gene for enhancing cold tolerance in japonica rice. The candidate miRNAs, mRNAs, and circRNAs uncovered in this study are valuable targets for researchers and breeders. Our findings will facilitate the development of cold-tolerant rice varieties from multiple angles and provide critical directions for future research into the functions of cold-tolerance-related miRNAs, mRNAs, and circRNAs in rice.

Vv-circSIZ1 mediated by pre-mRNA processing machinery contributes to salt tolerance

CircRNAs exist widely in plants, but the regulatory mechanisms for the biogenesis and function of plant circRNAs remain largely unknown. Using extensive mutagenesis of expression plasmids and genetic transformation methods, we analyzed the biogenesis and anti-salt functions of a new grape circRNA Vv-circSIZ1. We identified Vv-circSIZ1 that is mainly expressed in the cytoplasm of xylem. CircSIZ1 is species-specific, and genomic circSIZ1-forming region of seven tested species could be backspliced in Nicotiana benthamiana, but not in Arabidopsis. The retention length of Vv-circSIZ1 flanking introns was significantly positively correlated with its generation efficiency. The precise splicing of Vv-circSIZ1 does not depend on its mature exon sequence or internal intron sequences, but on the AG/GT splicing signal sites and branch site of the flanking introns. The spliceosome activity was inversely proportional to the expression level of Vv-circSIZ1. Furthermore, RNA-binding proteins can regulate the expression of Vv-circSIZ1. The overexpression of Vv-circSIZ1 improved salt tolerance of grape and N. benthamiana. Additionally, Vv-circSIZ1 could relieve the repressive effect of VvmiR3631 on its target VvVHAc1. Vv-circSIZ1 also promoted transcription of its parental gene. Overall, these results broaden our understanding of circRNAs in plants.

Role of non-coding RNAs against salinity stress in Oryza species: Strategies and challenges in analyzing miRNAs, tRFs and circRNAs

Salinity is an imbalanced concentration of mineral salts in the soil or water that causes yield loss in salt-sensitive crops. Rice plant is vulnerable to soil salinity stress at seedling and reproductive stages. Different non-coding RNAs (ncRNAs) post-transcriptionally regulate different sets of genes during different developmental stages under varying salinity tolerance levels. While microRNAs (miRNAs) are well known small endogenous ncRNAs, tRNA-derived RNA fragments (tRFs) are an emerging class of small ncRNAs derived from tRNA genes with a demonstrated regulatory role, like miRNAs, in humans but unexplored in plants. Circular RNA (circRNA), another ncRNA produced by back-splicing events, acts as target mimics by preventing miRNAs from binding with their target mRNAs, thereby reducing the miRNA's action upon its target. Same may hold true between circRNAs and tRFs. Hence, the work done on these ncRNAs was reviewed and no reports were found for circRNAs and tRFs under salinity stress in rice, either at seedling or reproductive stages. Even the reports on miRNAs are restricted to seedling stage only, in spite of severe effects on rice crop production due to salt stress during reproductive stage. Moreover, this review sheds light on strategies to predict and analyze these ncRNAs in an effective manner.

Genome-wide identification and characterization of tissue-specific non-coding RNAs in black pepper (Piper nigrum L.)

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are the two classes of non-coding RNAs (ncRNAs) present predominantly in plant cells and have various gene regulatory functions at pre- and post-transcriptional levels. Previously deemed as "junk", these ncRNAs have now been reported to be an important player in gene expression regulation, especially in stress conditions in many plant species. Black pepper, scientifically known as Piper nigrum L., despite being one of the most economically important spice crops, lacks studies related to these ncRNAs. From a panel of 53 RNA-Seq datasets of black pepper from six tissues, namely, flower, fruit, leaf, panicle, root, and stem of six black pepper cultivars, covering eight BioProjects across four countries, we identified and characterized a total of 6406 lncRNAs. Further downstream analysis inferred that these lncRNAs regulated 781 black pepper genes/gene products via miRNA-lncRNA-mRNA network interactions, thus working as competitive endogenous RNAs (ceRNAs). The interactions may be various mechanisms like miRNA-mediated gene silencing or lncRNAs acting as endogenous target mimics (eTMs) of the miRNAs. A total of 35 lncRNAs were also identified to be potential precursors of 94 miRNAs after being acted upon by endonucleases like Drosha and Dicer. Tissue-wise transcriptome analysis revealed 4621 circRNAs. Further, miRNA-circRNA-mRNA network analysis showed 432 circRNAs combining with 619 miRNAs and competing for the binding sites on 744 mRNAs in different black pepper tissues. These findings can help researchers to get a better insight to the yield regulation and responses to stress in black pepper in endeavor of higher production and improved breeding programs in black pepper varieties.

Cellobiose elicits immunity in lettuce conferring resistance to Botrytis cinerea

Cellobiose is the primary product of cellulose hydrolysis and is expected to function as a type of pathogen/damage-associated molecular pattern in evoking plant innate immunity. In this study, cellobiose was demonstrated to be a positive regulator in the immune response of lettuce, but halted autoimmunity when lettuce was exposed to concentrations of cellobiose >60 mg l-1. When lettuce plants were infected by Botrytis cinerea, cellobiose endowed plants with enhanced pre-invasion resistance by activating high β-1,3-glucanase and antioxidative enzyme activities at the initial stage of pathogen infection. Cellobiose-activated core regulatory factors such as EDS1, PTI6, and WRKY70, as well as salicylic acid signaling, played an indispensable role in modulating plant growth-defense trade-offs. Transcriptomics data further suggested that the cellobiose-activated plant-pathogen pathways are involved in microbe/pathogen-associated molecular pattern-triggered immune responses. Genes encoding receptor-like kinases, transcription factors, and redox homeostasis, phytohormone signal transduction, and pathogenesis-related proteins were also up- or down-regulated by cellobiose. Taken together, the findings of this study demonstrated that cellobiose serves as an elicitor to directly activate disease-resistance-related cellular functions. In addition, multiple genes have been identified as potential modulators of the cellobiose-induced immune response, which could aid understanding of underlying molecular events.

Identification, biogenesis, function, and mechanism of action of circular RNAs in plants

Circular RNAs (circRNAs) are a class of single-stranded, closed RNA molecules with unique functions that are ubiquitously expressed in all eukaryotes. The biogenesis of circRNAs is regulated by specific cis-acting elements and trans-acting factors in humans and animals. circRNAs mainly exert their biological functions by acting as microRNA sponges, forming R-loops, interacting with RNA-binding proteins, or being translated into polypeptides or proteins in human and animal cells. Genome-wide identification of circRNAs has been performed in multiple plant species, and the results suggest that circRNAs are abundant and ubiquitously expressed in plants. There is emerging compelling evidence to suggest that circRNAs play essential roles during plant growth and development as well as in the responses to biotic and abiotic stress. However, compared with recent advances in human and animal systems, the roles of most circRNAs in plants are unclear at present. Here we review the identification, biogenesis, function, and mechanism of action of plant circRNAs, which will provide a fundamental understanding of the characteristics and complexity of circRNAs in plants.

CircRNA: a rising star in plant biology

Circular RNAs (circRNAs) are covalently closed single-stranded RNA molecules, which are widespread in eukaryotic cells. As regulatory molecules, circRNAs have various functions, such as regulating gene expression, binding miRNAs or proteins, and being translated into proteins, which are important for cell proliferation and cell differentiation, individual growth and development, as well as many other biological processes. However, compared with that in animal models, studies of circRNAs in plants lags behind and, particularly, the regulatory mechanisms of biogenesis and molecular functions of plant circRNAs remain elusive. Recent studies have shown that circRNAs are wide spread in plants with tissue- or development-specific expression patterns and are responsive to a variety of environmental stresses. In this review, we summarize these advances, focusing on the regulatory mechanisms of biogenesis, molecular and biological functions of circRNAs, and the methods for investigating circRNAs. We also discuss the challenges and the prospects of plant circRNA studies.

Regulatory non-coding RNA: The core defense mechanism against plant pathogens

Plant pathogens damage crops and threaten global food security. Plants have evolved complex defense networks against pathogens, using crosstalk among various signaling pathways. Key regulators conferring plant immunity through signaling pathways include protein-coding genes and non-coding RNAs (ncRNAs). The discovery of ncRNAs in plant transcriptomes was first considered "transcriptional noise". Recent reviews have highlighted the importance of non-coding RNAs. However, understanding interactions among different types of noncoding RNAs requires additional research. This review attempts to consider how long-ncRNAs, small-ncRNAs and circular RNAs interact in response to pathogenic diseases within different plant species. Developments within genomics and bioinformatics could lead to the further discovery of plant ncRNAs, knowledge of their biological roles, as well as an understanding of their importance in exploiting the recent molecular-based technologies for crop protection.

Identification and Characterization of Circular RNAs Involved in the Flower Development and Senescence of Rhododendron delavayi Franch

Floral development and senescence are a crucial determinant for economic and ornamental value. CircRNAs play an essential role in regulating plant growth and development; however, there is no systematic identification of circRNAs during the lifespan of flowers. This study aims to explore the expression profile and functional role of circRNAs in the full flowering stages of Rhododendron delavayi Franch. We carried out transcriptome sequencing of the six stages of Rhododendron delavayi Franch flowers to identify the circular RNA expression profile. In addition, using bioinformatics methods, we explored the functions of circRNAs, including analysis of the circRNA-miRNA-mRNA network, short time-series expression miner (STEM), and so on. We identified 146 circRNAs, of which 79 were differentially expressed from the budding to fading stages. Furthermore, using STEM analysis, one of the 42 circRNA expression model profiles was significantly upregulated during the senescence stage, including 16 circRNAs. Additionally, 7 circRNA-miRNA-mRNA networks were constructed with 10 differentially expressed circRNAs, in which some target mRNA may regulate the development and senescence of the Rhododendron flowers. Finally, by analyzing the correlation between circRNAs and mRNA, combined with existing reports, we proposed that circRNAs play a regulatory role during flower development and senescence by mediating the jasmonate signaling pathway. Overall, these results provide new clues to the potential mechanism of circRNAs acting as novel post-transcriptional regulators in the development and senescence process of flowers.

Mining the Roles of Cucumber DUF966 Genes in Fruit Development and Stress Response

DUF966 genes are widely found in monocotyledons, dicotyledons, mosses, and other species. Current evidence strongly suggests that they are involved in growth regulation and stress tolerance in crops. However, their functions in cucumbers remain unexplored. Here, cucumber CsDUF966 was systemically identified and characterized using bioinformatics. Eight CsDUF966 genes were identified in the cucumber genome. These were phylogenetically separated into three groups. All CsDUF966 proteins were hydrophilic and localized to the nucleus. Moreover, three acidic and five basic proteins were identified. Evolutionary analysis of DUF966 between cucumber and 33 other Cucurbitaceae species/cultivars revealed that most CsDUF966 genes were conserved, whereas CsDUF966_4.c and CsDUF966_7.c were positively selected among the five cucumber cultivars. Expression profiling analysis showed that CsDUF966 had variable expression patterns, and that miRNA164, miRNA166, and Csa-novel-35 were involved in the post-transcriptional regulation of CsDUF966_4.c and CsDUF966_7.c. The expression of CsDUF966_4.c and CsDUF966_7.c, which were under strong neofunctionalization selection, was strictly regulated in fruit and tissues, including seeds, pericarps, peels, and spines, suggesting that these genes are fruit growth regulators and were strongly selected during the cucumber breeding program. In conclusion, the results reveal the roles of CsDUF966s in regulating cucumber fruit development and lay the foundation for further functional studies.

Genome-wide circular RNA profiling and competing endogenous RNA regulatory network analysis provide new insights into the molecular mechanisms underlying early somatic embryogenesis in Dimocarpus longan Lour

Circular RNAs (circRNAs) are widely involved in plant growth and development. However, the function of circRNAs in plant somatic embryogenesis (SE) remains elusive. Here, by using high-throughput sequencing, a total of 5029 circRNAs were identified in the three stages of longan (Dimocarpus longan Lour.) early SE. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed (DE) circRNA host genes were enriched in the 'non-homologous end-joining' (NHEJ) and 'butanoate metabolism' pathways. In addition, the reactive oxygen species (ROS) content during longan early SE was determined. The results indicated that ROS-induced DNA double-strand breaks may not depend on the NHEJ repair pathway. Correlation analyses of the levels of related metabolites (glutamate, γ-aminobutyrate and pyruvate) and the expression levels of circRNAs and their host genes involved in butanoate metabolism were performed. The results suggested that circRNAs may act as regulators of the expression of cognate mRNAs, thereby affecting the accumulation of related compounds. A competing endogenous RNA (ceRNA) network of DE circRNAs, DE mRNAs, DE long noncoding RNAs (lncRNAs) and DE microRNAs (miRNAs) was constructed. The results showed that the putative targets of the noncoding RNA (ncRNAs) were significantly enriched in the KEGG pathways 'mitogen-activated protein kinase signaling' and 'nitrogen metabolism'. Furthermore, the expression patterns of the candidate circRNAs, lncRNAs, miRNAs and mRNAs confirmed the negative correlation between miRNAs and ceRNAs. In addition, two circRNA overexpression vectors were constructed to further verify the ceRNA network correlations in longan early SE. Our study revealed the potential role of circRNAs in longan early SE, providing new insights into the intricate regulatory mechanism underlying plant SE.

Transcriptome Profiles of Circular RNAs in Common Wheat during Fusarium Head Blight Disease

Circular RNAs (circRNAs) are covalently closed RNA molecules, and have been identified in many crops. However, there are few datasets for circRNA junctions from common wheat during Fusarium head blight disease. In the present study, we used RNA-seq to determine the changes in circRNAs among the control (CK) and 1, 3, and 5 days post-Fusarium graminearum inoculation (dpi) samples. More than one billion reads were produced from 12 libraries, and 99.99% of the reads were successfully mapped to a wheat reference genome. In total, 2091 high-confidence circRNAs—which had two or more junction reads and were supported by at least two circRNA identification algorithms—were detected. The completed expression profiling revealed a distinct expression pattern of circRNAs among the CK, 1dpi, 3dpi and 5dpi samples. This study provides a valuable resource for identifying F. graminearum infection-responsive circRNAs in wheat and for further functional characterization of circRNAs that participated in the Fusarium head blight disease response of wheat.

A Combination of a Genome-Wide Association Study and a Transcriptome Analysis Reveals circRNAs as New Regulators Involved in the Response to Salt Stress in Maize

Salinization seriously threatens the normal growth of maize, especially at the seedling stage. Recent studies have demonstrated that circular RNAs (circRNAs) play vital roles in the regulation of plant stress resistance. Here, we performed a genome-wide association study (GWAS) on the survival rate of 300 maize accessions under a salt stress treatment. A total of 5 trait-associated SNPs and 86 candidate genes were obtained by the GWAS. We performed RNA sequencing for 28 transcriptome libraries derived from 2 maize lines with contrasting salt tolerance under normal and salt treatment conditions. A total of 1217 highly expressed circRNAs were identified, of which 371 were responsive to a salt treatment. Using PCR and Sanger sequencing, we verified the reliability of these differentially expressed circRNAs. An integration of the GWAS and RNA-Seq analyses uncovered two differentially expressed hub genes (Zm00001eb013650 and Zm00001eb198930), which were regulated by four circRNAs. Based on these results, we constructed a regulation model of circRNA/miRNA/mRNA that mediated salt stress tolerance in maize. By conducting hub gene-based association analyses, we detected a favorable haplotype in Zm00001eb198930, which was responsible for high salt tolerance. These results help to clarify the regulatory relationship between circRNAs and their target genes as well as to develop salt-tolerant lines for maize breeding.

Advances of circRNA-miRNA-mRNA regulatory network in cerebral ischemia/reperfusion injury

Ischemic stroke (IS) is the most common type of stroke, and its pathophysiological process is more complex. In recent years, the key regulatory roles of non-coding RNA (miRNA, circRNA) and mRNA in the development of IS have attracted more attention. In the process of cerebral ischemia/reperfusion injury, circRNA can regulate nerves, blood vessels and immune system through miRNA/mRNA axis, so as to affect the neurovascular unit of IS. The combination of these noncoding RNAs and mRNAs can be used as non-invasive biomarkers and therapeutic tools for IS diagnosis, prognosis and brain injury. Therefore, it is very important to study the potential molecular mechanism, activation pathway and treatment methods of circRNA/miRNA/mRNA network in IS. This review will focus on the latest progress of circRNA/miRNA/mRNA regulatory network, we have also included some circRNAs, which does not mediate through a miRNA, so we also include circRNA -mRNA network. And explore the application prospect of these RNAs as potential therapeutic targets in the prevention and treatment of IS.

Identification of Reference Genes for Reverse Transcription-Quantitative PCR Analysis of Ginger Under Abiotic Stress and for Postharvest Biology Studies

Gene expression analysis largely improves our understanding of the molecular basis underpinning various plant biological processes. Stable reference genes play a foundational role during the normalization of gene expression levels. However, until now, there have been few reference genes suitable for ginger reverse transcription-quantitative PCR (RT-qPCR) research. In this study, 29 candidate reference genes with stable expression patterns across multiple ginger tissues and 13 commonly used reference genes were selected to design RT-qPCR primers. After amplification specificity validation, 32 candidates were selected and further evaluated by RT-qPCR using samples from various organs subjected to NaCl, drought, heat, waterlogging, and chilling stress. Four strategies, including delta-CT, BestKeeper, geNorm, and NormFinder, were used to rank the stability of reference genes, and the ranks produced by these four strategies were comprehensively evaluated by RefFinder to determine the final rank. Overall, the top three stability reference genes indicated by RefFinder were RBP > ATPase > 40S_S3. Their expression pattern correlation analysis showed that the coefficients among each pair of RBP, ATPase, and 40S_S3 were larger than 0.96, revealing consistent and stable expression patterns under various treatments. Then, the expression of three pathogenesis-related (PR) genes and seven MYB genes in rhizomes during postharvest storage and subjected to pathogen infection was normalized by RBP, ATPase, 40S_S3, RBP and ATPase, ATPase and 40S-S3, and RBP and 40S-S3. The results showed that PR and MYB genes were induced by postharvest deterioration and pathogen infection. The correlation coefficients of RBP/ATPase, RBP/40S_S3, ATPase/40S_S3, RBP and ATPase/ATPase and 40S-S3, RBP and ATPase/RBP and 40S-S3, and ATPase and 40S-S3/RBP and 40S-S3 were 0.99, 0.96, 0.99, 0.99, 1.00, and 1.00, respectively, which confirmed the stability of these three reference genes in postharvest biology studies of ginger. In summary, this study identified appropriate reference genes for RT-qPCR in ginger and facilitated gene expression studies under biotic and abiotic stress conditions.

Drought tolerance improvement in Solanum lycopersicum: an insight into "OMICS" approaches and genome editing

Solanum lycopersicum (tomato) is an internationally acclaimed vegetable crop that is grown worldwide. However, drought stress is one of the most critical challenges for tomato production, and it is a crucial task for agricultural biotechnology to produce drought-resistant cultivars. Although breeders have done a lot of work on the tomato to boost quality and quantity of production and enhance resistance to biotic and abiotic stresses, conventional tomato breeding approaches have been limited to improving drought tolerance because of the intricacy of drought traits. Many efforts have been made to better understand the mechanisms involved in adaptation and tolerance to drought stress in tomatoes throughout the years. "Omics" techniques, such as genomics, transcriptomics, proteomics, and metabolomics in combination with modern sequencing technologies, have tremendously aided the discovery of drought-responsive genes. In addition, the availability of biotechnological tools, such as plant transformation and the recently developed genome editing system for tomatoes, has opened up wider opportunities for validating the function of drought-responsive genes and the generation of drought-tolerant varieties. This review highlighted the recent progresses for tomatoes improvement against drought stress through "omics" and "multi-omics" technologies including genetic engineering. We have also discussed the roles of non-coding RNAs and genome editing techniques for drought stress tolerance improvement in tomatoes.

Identification and evolutionary analysis of the metal-tolerance protein family in eight Cucurbitaceae species

Metal-tolerance proteins (MTPs) are divalent cation transporters and play fundamental roles in plant metal tolerance and ion homeostasis. Despite that, a systematic investigation of MTPs in Cucurbitacea is still lacking. In this study, 142 MTPs were identified from 11 released genomes of 8 Cucurbitaceae species. They were phylogenetically separated into three clusters (Zn-cation diffusion facilitator proteins [CDFs], Fe/Zn-CDFs, and Mn-CDFs) and further subdivided into seven groups (G1, G5, G6, G7, G8, G9, and G12). Characterization analysis revealed that most MTPs were plasma membrane-located hydrophobic proteins. Motif and exon/intron analysis showed that members in the same group contained similar conserved motifs and gene structures. Moreover, 98 pairs of segmental-like duplication events were found. The nonsynonymous/synonymous substitution ratios between each pair were less than 1, implying that Cucurbitaceae MTPs were under purification selection. Expression profiling suggested that several MTP genes, such as CsCLMTP1, CmeMTP3, LsMTP3, and Cl97103MTP3, were constitutively expressed in corresponding Cucurbitaceae species, and their expression levels were not significantly altered by NaCl, drought, or pathogen infection. The expression patterns of cucumber MTP genes under Zn²⁺ , Cu²⁺ , Mn²⁺ , and Cd²⁺ stress were studied by quantitative real-time polymerase chain reaction and the results showed that these MTPs were induced by at least one metal ion, suggesting their involvement in metal tolerance or transportation. The identification and comprehensive investigation of MTP family members will provide a basis for the analysis of ion transport functions and ion tolerance mechanisms of Cucurbitaceae MTPs.

Identification and Characterization of Rice Circular RNAs Responding to Xanthomonas oryzae pv. oryzae Invasion

Emerging roles of circular RNAs (circRNAs) in various biological processes have advanced our knowledge of transcriptional and posttranscriptional gene regulation. To date, no research has been conducted to explore their roles in the rice-Xanthomonas oryzae pv. oryzae interaction. Therefore, we identified 3,517 circRNAs from rice leaves infected with the highly virulent X. oryzae pv. oryzae strain PXO99^A by using rRNA depleted RNA sequencing technique coupled with the CIRI2 and CIRCexplorer2 pipeline. Characterization analyses showed that these circRNAs were distributed across the whole genome of rice, and most circRNAs arose from exons (85.13%), ranged from 200 to 1,000 bp, and were with a noncanonical GT/AG (including CT/AC equivalent) splicing signal. Functional annotation and enrichment analysis of the host genes that produced the differentially expressed circRNAs (DEcircRNAs) suggested that these identified circRNAs might play an important role in reprogramming rice responses to PXO99^A invasion, mainly by mediating photorespiration and chloroplast, peroxisome, and diterpenoid biosynthesis. Moreover, 31 DEcircRNAs were predicted to act as microRNA decoys in rice. The expression profile of four DEcircRNAs were validated by quantitative real-time PCR with divergent primers, and the back-splicing sites of seven DEcircRNAs were verified by PCR analysis and Sanger sequencing. Collectively, these results inferred a potential functional role of circRNAs in the regulation of rice immunity and provide novel clues about the molecular mechanisms of rice-PXO99^A interaction.

Identification and functional deciphering suggested the regulatory roles of long intergenic ncRNAs (lincRNAs) in increasing grafting pepper resistance to Phytophthora capsici

BACKGROUND

As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown.

RESULTS

In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks.

CONCLUSIONS

This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Comparative genomic analysis of expansin superfamily gene members in zucchini and cucumber and their expression profiles under different abiotic stresses

Zucchini and cucumber belong to the Cucurbitaceae family, a group of economical and nutritious food plants that is consumed worldwide. Expansin superfamily proteins are generally localized in the cell wall of plants and are known to possess an effect on cell wall modification by causing the expansion of this region. Although the whole genome sequences of cucumber and zucchini plants have been resolved, the determination and characterization of expansin superfamily members in these plants using whole genomic data have not been implemented yet. In the current study, a genome-wide analysis of zucchini (Cucurbita pepo) and cucumber (Cucumis sativus) genomes was performed to determine the expansin superfamily genes. In total, 49 and 41 expansin genes were identified in zucchini and cucumber genomes, respectively. All expansin superfamily members were subjected to further bioinformatics analysis including gene and protein structure, ontology of the proteins, phylogenetic relations and conserved motifs, orthologous relations with other plants, targeting miRNAs of those genes and in silico gene expression profiles. In addition, various abiotic stress responses of zucchini and cucumber expansin genes were examined to determine their roles in stress tolerance. CsEXPB-04 and CsEXPA-11 from cucumber and CpEXPA-20 and CpEXPLA-14 from zucchini can be candidate genes for abiotic stress response and tolerance in addition to their roles in the normal developmental processes, which are supported by the gene expression analysis. This work can provide new perspectives for the roles of expansin superfamily genes and offers comprehensive knowledge for future studies investigating the modes of action of expansin proteins.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01108-w.

Combined Stress Conditions in Melon Induce Non-additive Effects in the Core miRNA Regulatory Network

Climate change has been associated with a higher incidence of combined adverse environmental conditions that can promote a significant decrease in crop productivity. However, knowledge on how a combination of stresses might affect plant development is still scarce. MicroRNAs (miRNAs) have been proposed as potential targets for improving crop productivity. Here, we have combined deep-sequencing, computational characterization of responsive miRNAs and validation of their regulatory role in a comprehensive analysis of response of melon to several combinations of four stresses (cold, salinity, short day, and infection with a fungus). Twenty-two miRNA families responding to double and/or triple stresses were identified. The regulatory role of the differentially expressed miRNAs was validated by quantitative measurements of the expression of the corresponding target genes. A high proportion (ca. 60%) of these families (mainly highly conserved miRNAs targeting transcription factors) showed a non-additive response to multiple stresses in comparison with that observed under each one of the stresses individually. Among those miRNAs showing non-additive response to stress combinations, most interactions were negative, suggesting the existence of functional convergence in the miRNA-mediated response to combined stresses. Taken together, our results provide compelling pieces of evidence that the response to combined stresses cannot be easily predicted from the study individual stresses.

Genome-Wide Characterization of Salt-Responsive miRNAs, circRNAs and Associated ceRNA Networks in Tomatoes

Soil salinization is a major environmental stress that causes crop yield reductions worldwide. Therefore, the cultivation of salt-tolerant crops is an effective way to sustain crop yield. Tomatoes are one of the vegetable crops that are moderately sensitive to salt stress. Global market demand for tomatoes is huge and growing. In recent years, the mechanisms of salt tolerance in tomatoes have been extensively investigated; however, the molecular mechanism through which non-coding RNAs (ncRNAs) respond to salt stress is not well understood. In this study, we utilized small RNA sequencing and whole transcriptome sequencing technology to identify salt-responsive microRNAs (miRNAs), messenger RNAs (mRNAs), and circular RNAs (circRNAs) in roots of M82 cultivated tomato and Solanum pennellii (S. pennellii) wild tomato under salt stress. Based on the theory of competitive endogenous RNA (ceRNA), we also established several salt-responsive ceRNA networks. The results showed that circRNAs could act as miRNA sponges in the regulation of target mRNAs of miRNAs, thus participating in the response to salt stress. This study provides insights into the mechanisms of salt tolerance in tomatoes and serves as an effective reference for improving the salt tolerance of salt-sensitive cultivars.

Might exogenous circular RNAs act as protein-coding transcripts in plants?

Circular RNAs (circRNAs) are regulatory molecules involved in the modulation of gene expression. Although originally assumed as non-coding RNAs, recent studies have evidenced that animal circRNAs can act as translatable transcripts. The study of plant-circRNAs is incipient, and no autonomous coding plant-circRNA has been described yet. Viroids are the smallest plant-pathogenic circRNAs known to date. Since their discovery 50 years ago, viroids have been considered valuable systems for the study of the structure-function relationships in RNA, essentially because they have not been shown to have coding capacity. We used two pathogenic circRNAs (Hop stunt viroid and Eggplant latent viroid) as experimental tools to explore the coding potential of plant-circRNAs. Our work supports that the analysed viroids contain putative ORFs able to encode peptides carrying subcellular localization signals coincident with the corresponding replication-specific organelle. Bioassays in well-established hosts revealed that mutations in these ORFs diminish their biological efficiency. Interestingly, circular forms of HSVd and ELVd were found to co-sediment with polysomes, revealing their physical interaction with the translational machinery of the plant cell. Based on this evidence we hypothesize about the possibility that plant circRNAs in general, and viroids in particular, can act, under certain cellular conditions, as non-canonical translatable transcripts.

Genome-wide identification of ZF-HD gene family in Triticum aestivum: Molecular evolution mechanism and function analysis

ZF-HD family genes play important roles in plant growth and development. Studies about the whole genome analysis of ZF-HD gene family have been reported in some plant species. In this study, the whole genome identification and expression profile of the ZF-HD gene family were analyzed for the first time in wheat. A total of 37 TaZF-HD genes were identified and divided into TaMIF and TaZHD subfamilies according to the conserved domain. The phylogeny tree of the TaZF-HD proteins was further divided into six groups based on the phylogenetic relationship. The 37 TaZF-HDs were distributed on 18 of 21 chromosomes, and almost all the genes had no introns. Gene duplication and Ka/Ks analysis showed that the gene family may have experienced powerful purification selection pressure during wheat evolution. The qRT-PCR analysis showed that TaZF-HD genes had significant expression patterns in different biotic stress and abiotic stress. Through subcellular localization experiments, we found that TaZHD6-3B was located in the nucleus, while TaMIF4-5D was located in the cell membrane and nucleus. Our research contributes to a comprehensive understanding of the TaZF-HD family, provides a new perspective for further research on the biological functions of TaZF-HD genes in wheat.

Pan-transcriptome identifying master genes and regulation network in response to drought and salt stresses in Alfalfa (Medicago sativa L.)

Alfalfa is an important legume forage grown worldwide and its productivity is affected by environmental stresses such as drought and high salinity. In this work, three alfalfa germplasms with contrasting tolerances to drought and high salinity were used for unraveling the transcriptomic responses to drought and salt stresses. Twenty-one different RNA samples from different germplasm, stress conditions or tissue sources (leaf, stem and root) were extracted and sequenced using the PacBio (Iso-Seq) and the Illumina platforms to obtain full-length transcriptomic profiles. A total of 1,124,275 and 91,378 unique isoforms and genes were obtained, respectively. Comparative analysis of transcriptomes identified differentially expressed genes and isoforms as well as transcriptional and post-transcriptional modifications such as alternative splicing events, fusion genes and nonsense-mediated mRNA decay events and non-coding RNA such as circRNA and lncRNA. This is the first time to identify the diversity of circRNA and lncRNA in response to drought and high salinity in alfalfa. The analysis of weighted gene co-expression network allowed to identify master genes and isoforms that may play important roles on drought and salt stress tolerance in alfalfa. This work provides insight for understanding the mechanisms by which drought and salt stresses affect alfalfa growth at the whole genome level.

NGS Methodologies and Computational Algorithms for the Prediction and Analysis of Plant Circular RNAs

Circular RNAs (circRNAs) are a class of single-stranded RNAs derived from exonic, intronic, and intergenic regions from precursor messenger RNAs (pre-mRNA), where a noncanonical back-splicing event occurs, in which the 5' and 3' ends are attached by covalent bond. CircRNAs participate in the regulation of gene expression at the transcriptional and posttranscriptional level primarily as miRNA and RNA-binding protein (RBP) sponges, but also involved in the regulation of alternative RNA splicing and transcription. CircRNAs are widespread and abundant in plants where they have been involved in stress responses and development. Through the analysis of all publications in this field in the last five years, we can summarize that the identification of these molecules is carried out through next generation sequencing studies, where samples have been previously treated to eliminate DNA, rRNA, and linear RNAs as a means to enrich circRNAs. Once libraries are prepared, they are sequenced and subsequently studied from a bioinformatics point of view. Among the different tools for identifying circRNAs, we can highlight CIRI as the most used (in 60% of the published studies), as well as CIRCExplorer (20%) and find_circ (20%). Although it is recommended to use more than one program in combination, and preferably developed specifically to treat with plant samples, this is not always the case. It should also be noted that after identifying these circular RNAs, most of the authors validate their findings in the laboratory in order to obtain bona fide results.

Genome-wide profiling of circular RNAs, alternative splicing, and R-loops in stem-differentiating xylem of Populus trichocarpa

Circular RNAs (circRNAs) are a recently discovered type of non-coding RNA derived from pre-mRNAs. R-loops consist of a DNA:RNA hybrid and the associated single-stranded DNA. In Arabidopsis thaliana, circRNA:DNA R-loops regulate alternative splicing (AS) of SEPALLATA3 (SEP3). However, the occurrence and functions of circRNAs and R-loops in Populus trichocarpa are largely unexplored. Here, we performed circRNA-enriched sequencing in the stem-differentiating xylem (SDX) of P. trichocarpa and identified 2,742 distinct circRNAs, including circ-CESA4, circ-IRX7, and circ-GUX1, which are generated from genes involved in cellulose, and hemicellulose biosynthesis, respectively. To investigate the roles of circRNAs in modulating alternative splicing (AS), we detected 7,836 AS events using PacBio Iso-Seq and identified 634 circRNAs that overlapped with 699 AS events. Furthermore, using DNA:RNA hybrid immunoprecipitation followed by sequencing (DRIP-seq), we identified 8,932 R-loop peaks that overlapped with 181 circRNAs and 672 AS events. Notably, several SDX-related circRNAs overlapped with R-loop peaks, pointing to their possible roles in modulating AS in SDX. Indeed, overexpressing circ-IRX7 increased the levels of R-loop structures and decreased the frequency of intron retention in linear IRX7 transcripts. This study provides a valuable R-loop atlas resource and uncovers the interplay between circRNAs and AS in SDX of P. trichocarpa.

Identification and characterization of melon circular RNAs involved in powdery mildew responses through comparative transcriptome analysis

Circular RNAs (circRNAs) are a class of newly discovered non-coding RNAs that are typically derived from a genome's exonic, intronic, and intergenic regions. Recent studies of circRNAs in animals and plants have shown that circRNAs are vital in response to various abiotic and biotic stresses. Powdery mildew disease (PM) is a serious fungal disease threatening the melon industry. We performed whole transcriptome sequencing using the leaves of a PM-resistant (M1) and a PM-susceptible (B29) melon to identify circRNAs and determine their molecular functions. A total of 303 circRNAs were identified and >50% circRNAs were derived from exonic regions. Expression levels were significantly altered in 17 and 23 circRNAs after PM infections in B29 and M1, respectively. Melon circRNAs may participate in the response to biotic stimuli, oxidation reduction, metabolic processes, and the regulation of gene expression based on the functional annotation of circRNA parental genes. Furthermore, 27 circRNAs were predicted to be potential targets or 'sponges' for 18 microRNAs (miRNAs). Our results are the first to identify and characterize circRNA functions in melon and may contribute to a better understanding of the role and regulatory mechanisms of circRNAs in resisting PM.

Genome-wide analysis of growth-regulating factors (GRFs) in Triticum aestivum

The Growth-Regulating Factor (GRF) family encodes a type of plant-specific transcription factor (TF). GRF members play vital roles in plant development and stress response. Although GRF family genes have been investigated in a variety of plants, they remain largely unstudied in bread wheat (Triticum aestivum L.). The present study was conducted to comprehensively identify and characterize the T. aestivum GRF (TaGRF) gene family members. We identified 30 TaGRF genes, which were divided into four groups based on phylogenetic relationship. TaGRF members within the same subgroup shared similar motif composition and gene structure. Synteny analysis suggested that duplication was the dominant reason for family member expansion. Expression pattern profiling showed that most TaGRF genes were highly expressed in growing tissues, including shoot tip meristems, stigmas and ovaries, suggesting their key roles in wheat growth and development. Further qRT-PCR analysis revealed that all 14 tested TaGRFs were significantly differentially expressed in responding to drought or salt stresses, implying their additional involvement in stress tolerance of wheat. Our research lays a foundation for functional determination of TaGRFs, and will help to promote further scrutiny of their regulatory network in wheat development and stress response.

Systematic identification and characterization of circular RNAs involved in flag leaf senescence of rice

Circular RNAs (circRNAs) identification, expression profiles, and construction of circRNA-parental gene relationships and circRNA-miRNA-mRNA ceRNA networks indicate that circRNAs are involved in flag leaf senescence of rice. Circular RNAs (circRNAs) are a class of 3'-5' head-to-tail covalently closed non-coding RNAs which have been proved to play important roles in various biological processes. However, no systematic identification of circRNAs associated with leaf senescence in rice has been studied. In this study, a genome-wide high-throughput sequencing analysis was performed using rice flag leaves developing from normal to senescence. Here, a total of 6612 circRNAs were identified, among which, 113 circRNAs were differentially expressed (DE) during the leaf senescence process. Moreover, 4601 (69.59%) circRNAs were derived from the exons or introns of their parental genes, while 2110 (71%) of the parental genes produced only one circRNA. The sequence alignment analysis showed that hundreds of rice circRNAs were conserved among different plant species. Gene Ontology (GO) enrichment analysis revealed that parental genes of DE circRNAs were enriched in many biological processes closely related to leaf senescence. Through weighted gene co-expression network analysis (WGCNA), six continuously down-expressed circRNAs, 18 continuously up-expressed circRNAs and 15 turn-point high-expressed circRNAs were considered to be highly associated with leaf senescence. Additionally, a total of 17 senescence-associated circRNAs were predicted to have parental genes, in which, regulations of three circRNAs to their parental genes were validated by qRT-PCR. The competing endogenous RNA (ceRNA) networks were also constructed. And a total of 11 senescence-associated circRNAs were predicted to act as miRNA sponges to regulate mRNAs, in which, regulation of two circRNAs to eight mRNAs was validated by qRT-PCR. It is discussed that senescence-associated circRNAs were involved in flag leaf senescence probably through mediating their parental genes and ceRNA networks, to participate in several well-studied senescence-associated processes, mainly including the processes of transcription, translation, and posttranslational modification (especially protein glycosylation), oxidation-reduction process, involvement of senescence-associated genes, hormone signaling pathway, proteolysis, and DNA damage repair. This study not only showed the systematic identification of circRNAs involved in leaf senescence of rice, but also laid a foundation for functional research on candidate circRNAs.

Identification of Circular RNAs by Multiple Displacement Amplification and Their Involvement in Plant Development

With the innovative knowledge and bioinformatics tools in the identification and characterization of noncoding RNAs, circular RNA (circRNA) is added as a new member to the noncoding RNAs family. CircRNA enrichment by rRNA depletion/RNase R or poly-A removal/RNase R treatment followed by NGS analysis is the most frequently adopted method for circular RNA identification and characterization. In this chapter, we describe the multiple displacement amplification (MDA) as a convenient method to augment the identification of even the abysmally expressed circular RNAs at low sequencing depth. Total RNA, extracted at three different developmental stages of rice, is subjected to RiboMinus and RNase R treatment to deplete the linear RNAs. The enriched circular RNAs are reverse transcribed with random hexamers. The resulting cDNA is subjected to phi29 DNA polymerase amplification using exo-resistant random pentamers to yield high molecular weight dsDNA product, followed by Illumina sequencing at ten million paired end reads per sample. The sequence analysis yielded a promising number of circRNAs with the appreciable inclusion of differentially regulated and minimally expressed circRNAs at a comparatively reduced cost.

Generation of Transgenic Rice Expressing CircRNA and Its Functional Characterization

Circular RNA (CircRNA) is yet another vital addition to the noncoding RNA family. They are mainly derived by fusion of downstream 3' splice donor with upstream 5' splice acceptor by a noncanonical form of alternative splicing mechanism called backsplicing. An array of functional aspects of these circRNAs has been reported in animal systems. However, functional investigation of circRNA in plants is very limited. In this chapter, we described a methodological outline to study the circRNA biogenesis and to characterize its function(s). Sequence of a newly identified Oryza sativa Indica circRNA flanked by complementary repeat sequences of a rice intron was assembled to yield a circRNA expression cassette. This cassette can be cloned into any plant expression vector which has a suitable promoter (CaMV 35S or ubiquitin promoter) and terminator, and can be used for any circRNA-mediated functional studies. Subsequent agroinfection of rice calli with this cassette yielded circRNA expressing transgenic plants. These transgenic plants were used to establish a correlation between the expressing circRNA, parental gene, and interacting miRNAs. Moreover, effect of circRNA overexpression on plant phenotype under various stress conditions can be studied using these transgenic plants. Also, RNA pull-down assay can be performed to identify the circRNA interacting proteins and the expression of these RBPs can also be studied from these transgenic plants.

Computational methods for annotation of plant regulatory non-coding RNAs using RNA-seq

Plant transcriptome encompasses numerous endogenous, regulatory non-coding RNAs (ncRNAs) that play a major biological role in regulating key physiological mechanisms. While studies have shown that ncRNAs are extremely diverse and ubiquitous, the functions of the vast majority of ncRNAs are still unknown. With ever-increasing ncRNAs under study, it is essential to identify, categorize and annotate these ncRNAs on a genome-wide scale. The use of high-throughput RNA sequencing (RNA-seq) technologies provides a broader picture of the non-coding component of transcriptome, enabling the comprehensive identification and annotation of all major ncRNAs across samples. However, the detection of known and emerging class of ncRNAs from RNA-seq data demands complex computational methods owing to their unique as well as similar characteristics. Here, we discuss major plant endogenous, regulatory ncRNAs in an RNA sample followed by computational strategies applied to discover each class of ncRNAs using RNA-seq. We also provide a collection of relevant software packages and databases to present a comprehensive bioinformatics toolbox for plant ncRNA researchers. We assume that the discussions in this review will provide a rationale for the discovery of all major categories of plant ncRNAs.

Interactions and links among the noncoding RNAs in plants under stresses

The complex interplay among sRNAs, lncRNAs and circRNAs has been implicated in plants under biotic and abiotic stresses. Here, we review current advances in our understanding of ncRNA interactions and links, which have considerable potential for improving the agronomic traits and the environmental adaptability of plants. Plants can respond to biotic or abiotic stresses. To cope with various conditions, numerous intricate molecular regulatory mechanisms have evolved in plants. Noncoding RNAs (ncRNAs) can be divided into small noncoding RNAs (sRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs). Emerging evidence has demonstrated that interplay among the ncRNAs acts as a novel layer in the regulatory mechanisms, which has attracted substantial interest. Links between sRNAs can affect plant immune responses and development in synergistic or antagonistic manners. Additionally, multiple interactions between lncRNAs and sRNAs are involved in crop breeding, disease resistance and high tolerance to environmental stresses. Here, we summarize current knowledge of the interactions and links among the ncRNAs in plant responses to stresses and the methods for identifying ncRNA interactions. Furthermore, challenges and prospects for further progress in elucidating ncRNA interactions and links are highlighted.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.

Silicon confers cucumber resistance to salinity stress through regulation of proline and cytokinins

Salt stress is a continuous threat to global crop production. Here, we studied the alleviation role of exogenous silicon (Si) in NaCl-stressed cucumber, with special emphasis on plant growth, proline (Pro) and hormone metabolisms. The results showed that Si supplementation ameliorated the adverse effects of NaCl on plants growth, biomass, and oxidative stress. Salt stress greatly increased the content of Pro throughout the experiment, while Si regulated Pro content in two distinct ways. Si promoted the salt-induced Pro levels after 3 and 6 days of treatment, but decreased it after 9 and 12 days of treatment. Moreover, P5CS and ProDH activities and P5CS gene play important roles in Si and salt-regulated Pro levels in different stress phase. Under stress condition, Si addition tend to revert the content of ABA, IAA, cytokinin and SA to the control levels in most cases. Further correlation analysis revealed a negative correlation between the root cytokinin and Pro content after 3 days of treatment, suggesting the interaction between cytokinin and Pro metabolism. Exogenous application of Pro and ProDH competitive inhibitor D-Lactate confirmed the possible interplay between Pro and cytokinin metabolism. Further study identified several CKX (Csa4G647490 and Csa1G589070) and IPT (Csa7G392940 and Csa3G150100) genes that may be responsible for the regulation of cytokinin accumulation by Si and/or Pro after short-term of treatment. The results suggested that Pro is a key factor in Si-induced salt tolerance, and Si-increased Pro content may participate in the regulation of cytokinin metabolism under short-term of salt stress.

Genome-wide identification and characterization of SRO gene family in wheat: Molecular evolution and expression profiles during different stresses

SRO (SIMILAR TO RCD ONE), a type of plant-specific small protein family, play important roles in plant growth and development, as well as in response to biotic/abiotic stresses. Although characterization of SROs have been performed in model plants, little is known about their function in wheat, especially under stress conditions. In this study, 30 SRO genes were identified from the wheat genome (TaSRO). They were phylogenetically separated into two groups with distinct structures. The cis-regulatory elements in the promoter region of TaSROs were analyzed and numerous elements functionally associated with stress responding and hormones were interpreted, implying the reason for induction expression patterns of TaSROs during abiotic and biotic stresses in wheat. Whole-genome replication events in the SRO gene family of wheat and seven other species (Arabidopsis thaliana, rice, maize, barley, soybean, upland cotton, and cucumber) were analyzed, resulting in 1, 12, 9, 23, 6, 5, and 3 of gene pairs, respectively. The tissue-specific expression pattern profiling revealed that most TaSROs are highly expressed in one or more tissues and may play an important role in wheat growth and development. In addition, qRT-PCR results further confirmed that these TaSRO genes are involved in wheat stress response. In summary, our study laid a theoretical basis for molecular function deciphering of TaSROs, especially in plant hormones and biotic/abiotic stress responses.

Non-coding RNAs as emerging targets for crop improvement

Food security is affected by climate change, population growth, as well as abiotic and biotic stresses. Conventional and molecular marker assisted breeding and genetic engineering techniques have been employed extensively for improving resistance to biotic stress in crop plants. Advances in next-generation sequencing technologies have permitted the exploration and identification of parts of the genome that extend beyond the regions with protein coding potential. These non-coding regions of the genome are transcribed to generate many types of non-coding RNAs (ncRNAs). These ncRNAs are involved in the regulation of growth, development, and response to stresses at transcriptional and translational levels. ncRNAs, including long ncRNAs (lncRNAs), small RNAs and circular RNAs have been recognized as important regulators of gene expression in plants and have been suggested to play important roles in plant immunity and adaptation to abiotic and biotic stresses. In this article, we have reviewed the current state of knowledge with respect to lncRNAs and their mechanism(s) of action as well as their regulatory functions, specifically within the context of biotic stresses. Additionally, we have provided insights into how our increased knowledge about lncRNAs may be used to improve crop tolerance to these devastating biotic stresses.

Genome-wide characterization of the C2H2 zinc-finger genes in Cucumis sativus and functional analyses of four CsZFPs in response to stresses

BACKGROUNDS

C2H2-type zinc finger protein (ZFPs) form a relatively large family of transcriptional regulators in plants, and play many roles in plant growth, development, and stress response. However, the comprehensive analysis of C2H2 ZFPs in cucumber (CsZFPs) and their regulation function in cucumber are still lacking.

RESULTS

In the current study, the whole genome identification and characterization of CsZFPs, including the gene structure, genome localization, phylogenetic relationship, and gene expression were performed. Functional analysis of 4 selected genes by transient transformation were also conducted. A total of 129 full-length CsZFPs were identified, which could be classified into four groups according to the phylogenetic analysis. The 129 CsZFPs unequally distributed on 7 chromosomes. Promoter cis-element analysis showed that the CsZFPs might involve in the regulation of phytohormone and/or abiotic stress response, and 93 CsZFPs were predicted to be targeted by one to 20 miRNAs. Moreover, the subcellular localization analysis indicated that 10 tested CsZFPs located in the nucleus and the transcriptome profiling analysis of CsZFPs demonstrated that these genes are involved in root and floral development, pollination and fruit spine. Furthermore, the transient overexpression of Csa1G085390 and Csa7G071440 into Nicotiana benthamiana plants revealed that they could decrease and induce leave necrosis in response to pathogen attack, respectively, and they could enhance salt and drought stresses through the initial induction of H2O2. In addition, Csa4G642460 and Csa6G303740 could induce cell death after 5 days transformation.

CONCLUSIONS

The identification and function analysis of CsZFPs demonstrated that some key individual CsZFPs might play essential roles in response to biotic and abiotic stresses. These results could lay the foundation for understanding the role of CsZFPs in cucumber development for future genetic engineering studies.

Characterization and Function of Circular RNAs in Plants

CircRNAs are covalently closed-loop single-stranded RNA molecules ubiquitously expressing in eukaryotes. As an important member of the endogenous ncRNA family, circRNAs are associated with diverse biological processes and can regulate transcription, modulate alternative splicing, and interact with miRNAs or proteins. Compared to abundant advances in animals, studies of circRNAs in plants are rapidly emerging. The databases and analysis tools for plant circRNAs are constantly being developed. Large numbers of circRNAs have been identified and characterized in plants and proved to play regulatory roles in plant growth, development, and stress responses. Here, we review the biogenesis, characteristics, bioinformatics resources, and biological functions of plant circRNAs, and summarize the distinct circularization features and differentially expression patterns comparison with animal-related results.

Genome-wide identification of long non-coding RNAs and circular RNAs reveal their ceRNA networks in response to cucumber green mottle mosaic virus infection in watermelon

Long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) play vital roles in plant defense responses against viral infections. However, there is no systematic understanding of lncRNAs and circRNAs and their competing endogenous RNA (ceRNA) networks in watermelon under cucumber green mottle mosaic virus (CGMMV) stress. Here, we present the characterization and expression profiles of lncRNAs and circRNAs in watermelon leaves 48-h post-inoculation (48 hpi) with CGMMV, with mock inoculation as a control. Deep sequencing analysis revealed 2373 lncRNAs and 606 circRNAs in the two libraries. Among them, 67 lncRNAs (40 upregulated and 27 downregulated) and 548 circRNAs (277 upregulated and 271 downregulated) were differentially expressed (DE) in the 48 hpi library compared with the control library. Furthermore, 263 cis-acting matched lncRNA-mRNA pairs were detected for 49 of the DE-lncRNAs. KEGG pathway analysis of the cis target genes of the DE-lncRNAs revealed significant associations with phenylalanine metabolism, the citrate cycle (TCA cycle), and endocytosis. Additionally, 30 DE-lncRNAs were identified as putative target mimics of 33 microRNAs (miRNAs), and 153 DE-circRNAs were identified as putative target mimics of 88 miRNAs. Furthermore, ceRNA networks of lncRNA/circRNA-miRNA-mRNA in response to CGMMV infection are described, with 12 DE-lncRNAs and 65 DE-circRNAs combining with 22 miRNAs and competing for the miRNA binding sites on 29 mRNAs. The qRT-PCR validation of selected lncRNAs and circRNAs showed a general correlation with the high-throughput sequencing results. This study provides a valuable resource of lncRNAs and circRNAs involved in the response to CGMMV infection in watermelon.

Genome-wide identification of CNGC genes in Chinese jujube (Ziziphus jujuba Mill.) and ZjCNGC2 mediated signalling cascades in response to cold stress

BACKGROUNDS

Cyclic nucleotide gated channels (CNGCs) play multifaceted roles in plant physiological processes, especially with respect to signalling processes, plant development, and responses to environmental stresses. However, little information is known about the CNGC family in the large cosmopolitan family Rhamnaceae, which has strong tolerance to biotic and abiotic stresses.

RESULTS

In the current study, a total of 15 ZjCNGCs which located on 7 chromosomes were firstly identified in Chinese jujube (Ziziphus jujuba Mill.), the most important species of Rhamnaceae in terms of economic and ecological values. Phylogenetic analysis showed that these ZjCNGCs could be classified into four groups, ZjCNGC12 belonged to group IVA, and ZjCNGC13, 14, 15 belonged to group IVB. In addition, the paralogous and orthologous homology duplication of ZjCNGC15 occurred during the evolutionary process. The characteristics of ZjCNGCs regarding to exon-intron numbers and post-translational modifications showed diversified structures and functions. Motif composition and protein sequence analysis revealed that the phosphate-binding cassette and hinge regions were conserved among ZjCNGCs. Prediction of the cis-acting regulatory elements and expression profiles by real-time quantitative PCR analysis showed that some of the ZjCNGCs responded to environmental changes, especially ZjCNGC2, which was significantly downregulated in response to cold stress, and ZjCNGC4 was highly induced in response to cold, salt and alkaline stresses. ZjCNGC13 and 14 were highly induced in the phytoplasma-resistant cultivar and downregulated in the susceptible cultivar. Furthermore, ZjCNGC2 could be regulated by cAMP treatment, microtubule changes and interact with ZjMAPKK4, which suggested that cAMP and microtubule might play important roles in ZjCNGC2 mediated ZjMAPKK4 signalling transduction involved in cold stress.

CONCLUSIONS

The identification and classification analysis of ZjCNGCs were firstly reported, and some key individual ZjCNGCs might play essential roles in the response to biotic and abiotic stresses, especially ZjCNGC2 mediated ZjMAPKK4 signalling transduction involved in cold stress. This systematic analysis could provide important information for further functional characterization of ZjCNGCs with the aim of breeding stress-resistant cultivars.

Genome-Wide Identification, Structure Characterization, and Expression Profiling of Dof Transcription Factor Gene Family in Wheat (Triticum aestivum L.)

DNA binding with one finger (Dof) proteins are plant-specific transcription factors with crucial roles in plant growth and stress response. Even so, little is known about them in wheat. In this study, 108 wheat Dof (TaDof) genes across 21 chromosomes were detected. Although variable in sequence length, molecular weight, and isoelectric point, all TaDof proteins contained conserved zinc-finger structures and were phylogenetically divided into 7 sub-groups. Exon/intron and motif analyses suggested that TaDof structures and conserved motifs were similar within sub-groups but diverse among sub-groups. Many segmental duplications were identified and Ka/Ks and inter-species synthetic analyses indicated that polyploidization was main reason for increased number of TaDofs. Prediction and experimental confirmation revealed that TaDofs functioned as transcription factors in the nucleus. Expression pattern profiling showed that TaDofs specifically affected growth and development, and biotic and abiotic stress responses. Wheat miRNAs and cis-regulator were predicted as essential players in molding TaDofs expression patterns. qRT-PCR analysis revealed that TaDofs were induced by salt and drought stresses. Customized annotation revealed that TaDofs were widely involved in phytohormone response, defense, growth and development, and metabolism. Our study provided a comprehensive understanding to wheat TaDofs.

Transcriptome-wide identification of novel circular RNAs in soybean in response to low-phosphorus stress

Low-phosphorus (LP) stress is a major factor limiting the growth and yield of soybean. Circular RNAs (circRNAs) are novel noncoding RNAs that play a crucial role in plant responses to abiotic stress. However, how LP stress mediates the biogenesis of circRNAs in soybean remains unclear. Here, to explore the response mechanisms of circRNAs to LP stress, the roots of two representative soybean genotypes with different P-use efficiency, Bogao (a LP-sensitive genotype) and Nannong 94156 (a LP-tolerant genotype), were used for the construction of RNA sequencing (RNA-seq) libraries and circRNA identification. In total, 371 novel circRNA candidates, including 120 significantly differentially expressed (DE) circRNAs, were identified across different P levels and genotypes. More DE circRNAs were significantly regulated by LP stress in Bogao than in NN94156, suggesting that the tolerant genotype was less affected by LP stress than the sensitive genotype was; in other words, NN94156 may have a better ability to maintain P homeostasis under LP stress. Moreover, a positive correlation was observed between the expression patterns of P stress-induced circRNAs and their circRNA-host genes. Gene Ontology (GO) enrichment analysis of these circRNA-host genes and microRNA (miRNA)-targeted genes indicated that these DE circRNAs were involved mainly in defense responses, ADP binding, nucleoside binding, organic substance catabolic processes, oxidoreductase activity, and signal transduction. Together, our results revealed that LP stress can significantly alter the genome-wide profiles of circRNAs and indicated that the regulation of circRNAs was both genotype and environment specific in response to LP stress. LP-induced circRNAs might provide a rich resource for LP-responsive circRNA candidates for future studies.

Genome-wide identification, structure characterization, and expression pattern profiling of aquaporin gene family in cucumber

BACKGROUND

Aquaporin (AQP) proteins comprise a group of membrane intrinsic proteins (MIPs) that are responsible for transporting water and other small molecules, which is crucial for plant survival under stress conditions including salt stress. Despite the vital role of AQPs, little is known about them in cucumber (Cucumis sativus L.).

RESULTS

In this study, we identified 39 aquaporin-encoding genes in cucumber that were separated by phylogenetic analysis into five sub-families (PIP, TIP, NIP, SIP, and XIP). Their substrate specificity was then assessed based on key amino acid residues such as the aromatic/Arginine (ar/R) selectivity filter, Froger's positions, and specificity-determining positions. The putative cis-regulatory motifs available in the promoter region of each AQP gene were analyzed and results revealed that their promoter regions contain many abiotic related cis-regulatory elements. Furthermore, analysis of previously released RNA-seq data revealed tissue- and treatment-specific expression patterns of cucumber AQP genes (CsAQPs). Three aquaporins (CsTIP1;1, CsPIP2;4, and CsPIP1;2) were the most transcript abundance genes, with CsTIP1;1 showing the highest expression levels among all aquaporins. Subcellular localization analysis in Nicotiana benthamiana epidermal cells revealed the diverse and broad array of sub-cellular localizations of CsAQPs. We then performed RNA-seq to identify the expression pattern of CsAQPs under salt stress and found a general decreased expression level of root CsAQPs. Moreover, qRT-PCR revealed rapid changes in the expression levels of CsAQPs in response to diverse abiotic stresses including salt, polyethylene glycol (PEG)-6000, heat, and chilling stresses. Additionally, transient expression of AQPs in N. benthamiana increased leaf water loss rate, suggesting their potential roles in the regulation of plant water status under stress conditions.

CONCLUSIONS

Our results indicated that CsAQPs play important roles in response to salt stress. The genome-wide identification and primary function characterization of cucumber aquaporins provides insight to elucidate the complexity of the AQP gene family and their biological functions in cucumber.

Role of Silicon in Mediating Salt Tolerance in Plants: A Review

Salt stress is a major threat for plant growth worldwide. The regulatory mechanisms of silicon in alleviating salt stress have been widely studied using physiological, molecular genetics, and genomic approaches. Recently, progresses have been made in elucidating the alleviative effects of silicon in salt-induced osmotic stress, Na toxicity, and oxidative stress. In this review, we highlight recent development on the impact of silicon application on salt stress responses. Emphasis will be given to the following aspects. (1) Silicon transporters have been experimentally identified in different plant species and their structure feature could be an important molecular basis for silicon permeability. (2) Silicon could mediate salt-induced ion imbalance by (i) regulating Na⁺ uptake, transport, and distribution and (ii) regulating polyamine levels. (3) Si-mediated upregulation of aquaporin gene expression and osmotic adjustment play important roles in alleviating salinity-induced osmotic stress. (4) Silicon application direct/indirectly mitigates oxidative stress via regulating the antioxidant defense and polyamine metabolism. (5) Omics studies reveal that silicon could regulate plants' response to salt stress by modulating the expression of various genes including transcription factors and hormone-related genes. Finally, research areas that require further investigation to provide a deeper understanding of the role of silicon in plants are highlighted.