Identification of circularRNAs and their targets in Gossypium under Verticillium wilt stress based on RNA-seq

The circular RNA circANK suppresses rice resistance to bacterial blight by inhibiting microRNA398b-mediated defense

Circular RNAs (circRNAs) are prevalent in eukaryotic cells and have been linked to disease progressions. Their unique circular structure and stability make them potential biomarkers and therapeutic targets. Compared with animal models, plant circRNA research is still in its infancy. The lack of effective tools to specifically knock down circRNAs without affecting host gene expression has slowed the progress of plant circRNA research. Here, we have developed a CRISPR-Cas13d tool that can specifically knock down circRNAs in plant systems, successfully achieving the targeted knockdown of circRNAs in rice (Oryza sativa). We further focused on Os-circANK (a circRNA derived from Ankyrin repeat-containing protein), a circRNA differentially expressed in rice upon pathogen infection. Physiological and biochemical experiments revealed that Os-circANK functions as a sponge for miR398b, suppressing the cleavage of Cu/Zn-superoxidase dismutase (CSD)1/CSD2/copper chaperone for superoxide dismutase/superoxidase dismutaseX through competing endogenous RNA, leading to reduced reactive oxygen species levels following Xanthomonas oryzae pv. oryzae (Xoo) infection and a negative regulation of rice resistance to bacterial blight. Our findings indicate Os-circANK inhibits rice resistance to bacterial blight via the microRNA398b(miR398b)/CSD/SOD pathway.

Viroids and Retrozymes: Plant Circular RNAs Capable of Autonomous Replication

Among the long non-coding RNAs that are currently recognized as important regulatory molecules influencing a plethora of processes in eukaryotic cells, circular RNAs (circRNAs) represent a distinct class of RNAs that are predominantly produced by back-splicing of pre-mRNA. The most studied regulatory mechanisms involving circRNAs are acting as miRNA sponges, forming R-loops with genomic DNA, and encoding functional proteins. In addition to circRNAs generated by back-splicing, two types of circRNAs capable of autonomous RNA-RNA replication and systemic transport have been described in plants: viroids, which are infectious RNAs that cause a number of plant diseases, and retrozymes, which are transcripts of retrotransposon genomic loci that are capable of circularization due to ribozymes. Based on a number of common features, viroids and retrozymes are considered to be evolutionarily related. Here, we provide an overview of the biogenesis mechanisms and regulatory functions of non-replicating circRNAs produced by back-splicing and further discuss in detail the currently available data on viroids and retrozymes, focusing on their structural features, replication mechanisms, interaction with cellular components, and transport in plants. In addition, biotechnological approaches involving replication-capable plant circRNAs are discussed, as well as their potential applications in research and agriculture.

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.

Comparative analysis, diversification, and functional validation of plant nucleotide-binding site domain genes

Nucleotide-binding site (NBS) domain genes are one of the superfamily of resistance genes involved in plant responses to pathogens. The current study identified 12,820 NBS-domain-containing genes across 34 species covering from mosses to monocots and dicots. These identified genes are classified into 168 classes with several novel domain architecture patterns encompassing significant diversity among plant species. Several classical (NBS, NBS-LRR, TIR-NBS, TIR-NBS-LRR, etc.) and species-specific structural patterns (TIR-NBS-TIR-Cupin_1-Cupin_1, TIR-NBS-Prenyltransf, Sugar_tr-NBS etc.) were discovered. We observed 603 orthogroups (OGs) with some core (most common orthogroups; OG0, OG1, OG2, etc.) and unique (highly specific to species; OG80, OG82, etc.) OGs with tandem duplications. The expression profiling presented the putative upregulation of OG2, OG6, and OG15 in different tissues under various biotic and abiotic stresses in susceptible and tolerant plants to cotton leaf curl disease (CLCuD). The genetic variation between susceptible (Coker 312) and tolerant (Mac7) Gossypium hirsutum accessions identified several unique variants in NBS genes of Mac7 (6583 variants) and Coker312 (5173 variants). The protein-ligand and proteins-protein interaction showed a strong interaction of some putative NBS proteins with ADP/ATP and different core proteins of the cotton leaf curl disease virus. The silencing of GaNBS (OG2) in resistant cotton through virus-induced gene silencing (VIGS) demonstrated its putative role in virus tittering. The presented study will be further helpful in understanding the plant adaptation mechanism.

Identification and Characterisation of the CircRNAs Involved in the Regulation of Leaf Colour in Quercus mongolica

Circular RNAs (circRNAs) are important regulatory molecules involved in various biological processes. However, the potential function of circRNAs in the turning red process of Quercus mongolica leaves is unclear. This study used RNA-seq data to identify 6228 circRNAs in leaf samples from four different developmental stages and showed that 88 circRNAs were differentially expressed. A correlation analysis was performed between anthocyanins and the circRNAs. A total of 16 circRNAs that may be involved in regulating the colour of Mongolian oak leaves were identified. CircRNAs may affect the colour of Q. mongolica leaves by regulating auxin, cytokinin, gibberellin, ethylene, and abscisic acid. This study revealed the potential role of circRNAs in the colour change of Q. mongolica leaves.

An accretive detection method for in silico identification and validation of circular RNAs in wheat (Triticum aestivum L.) using RT-qPCR

BACKGROUND

Circular RNAs (circRNAs) are novel class of non-coding RNAs, which are involved in various functions at the transcriptional and post-transcriptional level in response to a fungal pathogen (Puccinia triticina), including microRNA (miRNA) sponge, RNA binding proteins sponge, regulation of parental gene and biomarkers. Detailed analysis of wheat circRNAs is essential to accelerate the regulated expression of fungal miRNAs. Therefore, we suggest a protocol to aid circRNA identification through RNA-Seq data using various algorithms based on perl script followed by validation through divergent primer designing, standard PCR, and RT-qPCR assays.

METHODS AND RESULT

The divergent primer has been widely used to detect, validate, and quantify back-spliced junction (BSJ) of circRNAs. The procedure covers index file formation, circRNA identification and BSJ detections. However, the laboratory validation of circRNA includes wheat genomic DNA isolation, RNA isolation and its cDNA conversion upto validation. In this study, we identified 28 circRNAs from RNA-Seq of S0 and R0, wherein six circRNAs are commonly present and 75% of the identified circRNAs were belongs to inter-genic, 14% were exonic and intronic category were 11%. Divergent primer designing method successfully validated the two circRNAs via RT-qPCR assay, where circRNA_2 showed less relative expression pattern than circRNA_1 in contrast with housekeeping genes.

CONCLUSION

Thus, our results of identified and validated circRNAs showed that, this protocol is quite helpful, relatively easy, reliable, and accurate for large datasets as other algorithms need various dependencies and have complex scripts with high chances of error occurrence. Additionally, analysis time will vary depending on the expertise level and the number of RNA-Seq data. This proposed protocol can also be used for a wide range of monocotyledons belonging to the Poaceae plant family.

Identification and characterization of circular RNAs involved in the fertility stability of cotton CMS-D2 restorer line under heat stress

BACKGROUND

As a vital type of noncoding RNAs, circular RNAs (circRNAs) play important roles in plant growth and development and stress response. However, little is known about the biological roles of circRNAs in regulating the stability of male fertility restoration for cytoplasmic male sterility (CMS) conditioned by Gossypium harknessii cytoplasm (CMS-D2) cotton under high-temperature (HT) stress.

RESULTS

In this study, RNA-sequencing and bioinformatics analysis were performed on pollen grains of isonuclear alloplasmic near-isogenic restorer lines NH [N(Rf1rf1)] and SH [S(Rf1rf1)] with obvious differences in fertility stability under HT stress at two environments. A total of 967 circRNAs were identified, with 250 differentially expressed under HT stress. We confirmed the back-splicing sites of eight selected circRNAs using divergent primers and Sanger sequencing. Tissue-specific expression patterns of five differentially expressed circRNAs (DECs) were also verified by RT-PCR and qRT-PCR. Functional enrichment and metabolic pathway analysis revealed that the parental genes of DECs were significantly enriched in fertility-related biological processes such as pollen tube guidance and cell wall organization, as well as the Pentose and glucuronate interconversions, Steroid biosynthesis, and N-Glycan biosynthesis pathways. Moreover, we also constructed a putative circRNA-mediated competing endogenous RNA (ceRNA) network consisting of 21 DECs, eight predicted circRNA-binding miRNAs, and their corresponding 22 mRNA targets, especially the two ceRNA modules circRNA346-miR159a-MYB33 and circRNA484-miR319e-MYB33, which might play important biological roles in regulating pollen fertility stability of cotton CMS-D2 restorer line under HT stress.

CONCLUSIONS

Through systematic analysis of the abundance, characteristics and expression patterns of circRNAs, as well as the potential functions of their parent genes, our findings suggested that circRNAs and their mediated ceRNA networks acted vital biological roles in cotton pollen development, and might be also essential regulators for fertility stability of CMS-D2 restorer line under heat stress. This study will open a new door for further unlocking complex regulatory mechanisms underpinning the fertility restoration stability for CMS-D2 in cotton.

Comprehensive non-coding RNA analysis reveals specific lncRNA/circRNA-miRNA-mRNA regulatory networks in the cotton response to drought stress

Root and leaf are essential organs of plants in sensing and responding to drought stress. However, comparative knowledge of non-coding RNAs (ncRNAs) of root and leaf tissues in the regulation of drought response in cotton is limited. Here, we used deep sequencing data of leaf and root tissues of drought-resistant and drought-sensitive cotton varieties for identifying miRNAs, lncRNAs and circRNAs. A total of 1531 differentially expressed (DE) ncRNAs was identified, including 77 DE miRNAs, 1393 DE lncRNAs and 61 DE circRNAs. The tissue-specific and variety-specific competing endogenous RNA (ceRNA) networks of DE lncRNA-miRNA-mRNA response to drought were constructed. Furthermore, the novel drought-responsive lncRNA 1 (DRL1), specifically and differentially expressed in root, was verified to positively affect phenotypes of cotton seedlings under drought stress, competitively binding to miR477b with GhNAC1 and GhSCL3. In addition, we also constructed another ceRNA network consisting of 18 DE circRNAs, 26 DE miRNAs and 368 DE mRNAs. Fourteen circRNA were characterized, and a novel molecular regulatory system of circ125- miR7484b/miR7450b was proposed under drought stress. Our findings revealed the specificity of ncRNA expression in tissue- and variety-specific patterns involved in the response to drought stress, and uncovered novel regulatory pathways and potentially effective molecules in genetic improvement for crop drought resistance.

Insights to Gossypium defense response against Verticillium dahliae: the Cotton Cancer

The soil-borne pathogen Verticillium dahliae, also referred as "The Cotton Cancer," is responsible for causing Verticillium wilt in cotton crops, a destructive disease with a global impact. To infect cotton plants, the pathogen employs multiple virulence mechanisms such as releasing enzymes that degrade cell walls, activating genes that contribute to virulence, and using protein effectors. Conversely, cotton plants have developed numerous defense mechanisms to combat the impact of V. dahliae. These include strengthening the cell wall by producing lignin and depositing callose, discharging reactive oxygen species, and amassing hormones related to defense. Despite the efforts to develop resistant cultivars, there is still no permanent solution to Verticillium wilt due to a limited understanding of the underlying molecular mechanisms that drive both resistance and pathogenesis is currently prevalent. To address this challenge, cutting-edge technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), host-induced gene silencing (HIGS), and gene delivery via nano-carriers could be employed as effective alternatives to control the disease. This article intends to present an overview of V. dahliae virulence mechanisms and discuss the different cotton defense mechanisms against Verticillium wilt, including morphophysiological and biochemical responses and signaling pathways including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and strigolactones (SLs). Additionally, the article highlights the significance of microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) in gene expression regulation, as well as the different methods employed to identify and functionally validate genes to achieve resistance against this disease. Gaining a more profound understanding of these mechanisms could potentially result in the creation of more efficient strategies for combating Verticillium wilt in cotton crops.

Identification and characterization of CircRNA-associated CeRNA networks in moso bamboo under nitrogen stress

BACKGROUND

Nitrogen is a macronutrient element for plant growth and development. Circular RNAs (circRNAs) serve as pivotal regulators for the coordination between nutrient supply and plant demand. Moso bamboo (Phyllostachys edulis) is an excellent plant with fast growth, and the mechanism of the circRNA-target module in response to nitrogen remains unclear.

RESULTS

Deep small RNA sequencing results of moso bamboo seedlings under different concentrations of KNO₃ (N0 = 0 mM, N6 = 6 mM, N18 = 18 mM) were used to identify circRNAs. A total of 549 circRNAs were obtained, of which 309 were generated from corresponding parental coding genes including 66 new ones. A total of 536 circRNA-parent genes were unevenly distributed in 24 scaffolds and were associated with root growth and development. Furthermore, 52 differentially expressed circRNAs (DECs) were obtained, including 24, 33 and 15 DECs from three comparisons of N0 vs. N6, N0 vs. N18 and N6 vs. N18, respectively. Based on integrative analyses of the identified DECs, differentially expressed mRNAs (DEGs), and miRNAs (DEMs), a competitive endogenous RNA (ceRNA) network was constructed, including five DECs, eight DEMs and 32 DEGs. A regulatory module of PeSca_6:12,316,320|12,372,905-novel_miR156-PH02Gene35622 was further verified by qPCR and dual-luciferase reporter assays.

CONCLUSION

The results indicated that circRNAs could participate in multiple biological processes as miRNA sponges, including organ nitrogen compound biosynthesis and metabolic process regulation in moso bamboo. Our results provide valuable information for further study of circRNAs in moso bamboo under fluctuating nitrogen conditions.

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 in Brassica rapa in Response to Plasmodiophora brassicae

Plasmodiophora brassicae is a soil-borne pathogen that attacks the roots of cruciferous plants and causes clubroot disease. CircRNAs are noncoding RNAs, widely existing in plant and animal species. Although knowledge of circRNAs has been updated continuously and rapidly, information about circRNAs in the regulation of clubroot disease resistance is extremely limited in Brassica rapa. Here, Chinese cabbage (BJN 222) containing clubroot resistance genes (CRa) against P. brassicae Pb4 was susceptible to PbE. To investigate the mechanism of cicRNAs responsible for clubroot disease resistance in B. rapa, circRNA-seq was performed with roots of 'BJN 222' at 0, 8, and 23 days post-inoculated (dpi) with Pb4 and PbE. A total of 231 differentially expressed circRNAs were identified between the groups. Based on the differentially expressed circRNAs, the circRNA-miRNA-mRNA network was constructed using the target genes directly or indirectly related to plant resistance. Upregulated novel_circ_000495 suppressed the expression of miR5656-y, leading to the upregulation of Bra026508, which might cause plant resistance. Our results provide new insights into clubroot resistance mechanisms and lay a foundation for further studies exploring complex gene regulation networks in B. rapa.

lncRNA7 and lncRNA2 modulate cell wall defense genes to regulate cotton resistance to Verticillium wilt

In plants, long noncoding RNAs (lncRNAs) regulate disease resistance against fungi and other pathogens. However, the specific mechanism behind this regulation remains unclear. In this study, we identified disease resistance-related lncRNAs as well as their regulating genes and assessed their functions by infection of cotton (Gossypium) chromosome segment substitution lines with Verticillium dahliae. Our results demonstrated that lncRNA7 and its regulating gene Pectin methylesterase inhibitor 13 (GbPMEI13) positively regulated disease resistance via the silencing approach, while ectopic overexpression of GbPMEI13 in Arabidopsis (Arabidopsis thaliana) promoted growth and enhanced resistance to V. dahliae. In contrast, lncRNA2 and its regulating gene Polygalacturonase 12 (GbPG12) negatively regulated resistance to V. dahliae. We further found that fungal disease-related agents, including the pectin-derived oligogalacturonide (OG), could downregulate the expression of lncRNA2 and GbPG12, leading to pectin accumulation. Conversely, OG upregulated the expression of lncRNA7, which encodes a plant peptide phytosulfokine (PSK-α), which was confirmed by lncRNA7 overexpression and Ultra Performance Liquid Chromatography Tandem Mass Spectrometry (UPLC-MS) experiments. We showed that PSK-α promoted 3-Indoleacetic acid (IAA) accumulation and activated GbPMEI13 expression through Auxin Response Factor 5. Since it is an inhibitor of pectin methylesterase (PME), GbPMEI13 promotes pectin methylation and therefore increases the resistance to V. dahliae. Consistently, we also demonstrated that GbPMEI13 inhibits the mycelial growth and spore germination of V. dahliae in vitro. In this study, we demonstrated that lncRNA7, lncRNA2, and their regulating genes modulate cell wall defense against V. dahliae via auxin-mediated signaling, providing a strategy for cotton breeding.

Transcriptome Analysis of a Cotton Cultivar Provides Insights into the Differentially Expressed Genes Underlying Heightened Resistance to the Devastating Verticillium Wilt

Cotton is an important economic crop worldwide. Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is a serious disease in cotton, resulting in massive yield losses and decline of fiber quality. Breeding resistant cotton cultivars is an efficient but elaborate method to improve the resistance of cotton against V. dahliae infection. However, the functional mechanism of several excellent VW resistant cotton cultivars is poorly understood at present. In our current study, we carried out RNA-seq to discover the differentially expressed genes (DEGs) in the roots of susceptible cotton Gossypium hirsutum cultivar Junmian 1 (J1) and resistant cotton G.hirsutum cultivar Liaomian 38 (L38) upon Vd991 inoculation at two time points compared with the mock inoculated control plants. The potential function of DEGs uniquely expressed in J1 and L38 was also analyzed by GO enrichment and KEGG pathway associations. Most DEGs were assigned to resistance-related functions. In addition, resistance gene analogues (RGAs) were identified and analyzed for their role in the heightened resistance of the L38 cultivar against the devastating Vd991. In summary, we analyzed the regulatory network of genes in the resistant cotton cultivar L38 during V. dahliae infection, providing a novel and comprehensive insight into VW resistance in cotton.

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.

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.

Identification, characterization and expression profiling of circular RNAs in the early cotton fiber developmental stages

Circular RNA is one of the endogenous non-coding RNAs with a covalently closed loop structure and largely involved in regulating gene expression. However, the abundance of circular RNAs and their regulatory functions during the early stages of fiber development are still not known. In this work, we conducted high-throughput sequencing of the Ligonlintless-1 and its wild-type at 0 DPA, 8 DPA and stem. A total of 2811 circular RNAs were identified and unevenly distributed across cotton chromosomes. We found 34, 142 and 27 circular RNAs were differentially expressed between Ligonlintless-1 and wild-type at 0 DPA, 8 DPA and stem, respectively. Both circular RNA-microRNA-mRNA network and MeJA treatment results in Ligonlintless-1 and wild-type might provide a strong indication of four circular RNAs and ghr_miR169b being important biological molecular associating with fiber development. The results provide new insight into the putative molecular function of circular RNAs in the regulation of fiber development.

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.

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.

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.

Epigenomic landscape and epigenetic regulation in maize

Epigenetic regulation has been implicated in the control of multiple agronomic traits in maize. Here, we review current advances in our understanding of epigenetic regulation, which has great potential for improving agronomic traits and the environmental adaptability of crops. Epigenetic regulation plays vital role in the control of complex agronomic traits. Epigenetic variation could contribute to phenotypic diversity and can be used to improve the quality and productivity of crops. Maize (Zea mays L.), one of the most widely cultivated crops for human food, animal feed, and ethanol biofuel, is a model plant for genetic studies. Recent advances in high-throughput sequencing technology have made possible the study of epigenetic regulation in maize on a genome-wide scale. In this review, we discuss recent epigenetic studies in maize many achieved by Chinese research groups. These studies have explored the roles of DNA methylation, posttranslational modifications of histones, chromatin remodeling, and noncoding RNAs in the regulation of gene expression in plant development and environment response. We also provide our future prospects for manipulating epigenetic regulation to improve crops.

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.

circRNAs Are Involved in the Rice-Magnaporthe oryzae Interaction

Circular RNAs (circRNAs) play roles in various biological processes, but their functions in the rice (Oryza sativa) response to Magnaporthe oryzae remain elusive. Here, we demonstrate that circRNAs are involved in the rice-M. oryzae interaction using comparative circRNA-sequencing and transgenic approaches. We identified 2932 high-confidence circRNAs from young leaves of the blast-resistant accession International Rice Blast Line Pyricularia-Kanto51-m-Tsuyuake (IR25) and the blast-susceptible accession Lijiangxin Tuan Heigu (LTH) under M oryzae-infected or uninfected conditions; 636 were detected specifically upon M oryzae infection. The circRNAs in IR25 were significantly more diverse than those in LTH, especially under M. oryzae infection. Particularly, the number of circRNAs generated per parent gene was much higher in IR25 than in LTH and increased in IR25 but decreased in LTH upon M. oryzae infection. The higher diversity of circRNAs in IR25 was further associated with more frequent 3' and 5' alternative back-splicing and usage of complex splice sites. Moreover, a subset of circRNAs was differentially responsive to M oryzae in IR25 and LTH. We further confirmed that circR5g05160 promotes rice immunity against M oryzae Therefore, our data indicate that circRNA diversity is associated with different responses to M oryzae infection in rice and provide a starting point to investigate a new layer of regulation in the rice-M oryzae interaction.

Plant Noncoding RNAs: Hidden Players in Development and Stress Responses

A large and significant portion of eukaryotic transcriptomes consists of noncoding RNAs (ncRNAs) that have minimal or no protein-coding capacity but are functional. Diverse ncRNAs, including both small RNAs and long ncRNAs (lncRNAs), play essential regulatory roles in almost all biological processes by modulating gene expression at the transcriptional and posttranscriptional levels. In this review, we summarize the current knowledge of plant small RNAs and lncRNAs, with a focus on their biogenesis, modes of action, local and systemic movement, and functions at the nexus of plant development and environmental responses. The complex connections among small RNAs, lncRNAs, and small peptides in plants are also discussed, along with the challenges of identifying and investigating new classes of ncRNAs.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

Identification and characterization of CircRNAs involved in the regulation of wheat root length

BACKGROUND

Recent studies indicate that circular RNAs (circRNAs) may play important roles in the regulation of plant growth and development. Plant roots are the main organs of nutrient and water uptake. However, whether circRNAs involved in the regulation of plant root growth remains to be elucidated.

METHODS

LH9, XN979 and YN29 are three Chinese wheat varieties with contrasting root lengths. Here, the root circRNA expression profiles of LH9, XN979 and YN29 were examined by using high-throughput sequencing technology.

RESULTS

Thirty-three and twenty-two differentially expressed circRNAs (DECs) were identified in the YN29-LH9 comparison and YN29-XN979 comparison, respectively. Among them, ten DECs coexisted in both comparisons. As the roots of both LH9 and XN979 were significantly larger and deeper than YN29, the ten DECs coexisting in the two comparisons were highly likely to be involved in the regulation of wheat root length. Moreover, three of the ten DECs have potential miRNA binding sites. Real-time PCR analysis showed that the expression levels of the potential binding miRNAs exhibited significant differences between the long root plants and the short root plants.

CONCLUSIONS

The expression levels of some circRNAs exhibited significant differences in wheat varieties with contrasting root phenotypes. Ten DECs involved in the regulation of wheat root length were successfully identified in which three of them have potential miRNAs binding sites. The expression levels of putative circRNA-binding miRNAs were correlated with their corresponding circRNAs. Our results provide new clues for studying the potential roles of circRNAs in the regulation of wheat root length.