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

Regulation of seed germination: ROS, epigenetic, and hormonal aspects

BACKGROUND

The whole life of a plant is regulated by complex environmental or hormonal signaling networks that control genomic stability, environmental signal transduction, and gene expression affecting plant development and viability. Seed germination, responsible for the transformation from seed to seedling, is a key initiation step in plant growth and is controlled by unique physiological and biochemical processes. It is continuously modulated by various factors including epigenetic modifications, hormone transport, ROS signaling, and interaction among them. ROS showed versatile crucial functions in seed germination including various physiological oxidations to nucleic acid, protein, lipid, or chromatin in the cytoplasm, cell wall, and nucleus.

AIM

of review: This review intends to provide novel insights into underlying mechanisms of seed germination especially associated with the ROS, and considers how these versatile regulatory mechanisms can be developed as useful tools for crop improvement.

KEY SCIENTIFIC CONCEPTS OF REVIEW

We have summarized the generation and elimination of ROS during seed germination, with a specific focus on uncovering and understanding the mechanisms of seed germination at the level of phytohormones, ROS, and epigenetic switches, as well as the close connections between them. The findings exhibit that ROS plays multiple roles in regulating the ethylene, ABA, and GA homeostasis as well as the Ca2+ signaling, NO signaling, and MAPK cascade in seed germination via either the signal trigger or the oxidative modifier agent. Further, ROS shows the potential in the nuclear genome remodeling and some epigenetic modifiers function, although the detailed mechanisms are unclear in seed germination. We propose that ROS functions as a hub in the complex network regulating seed germination.

The roles of non-coding RNAs in male reproductive development and abiotic stress responses during this unique process in flowering plants

Successful male reproductive development is the guarantee for sexual reproduction of flowering plants. Male reproductive development is a complicated and multi-stage process that integrates physiological processes and adaptation and tolerance to a myriad of environmental stresses. This well-coordinated process is governed by genetic and epigenetic machineries. Non-coding RNAs (ncRNAs) play pleiotropic roles in the plant growth and development. The identification, characterization and functional analysis of ncRNAs and their target genes have opened a new avenue for comprehensively revealing the regulatory network of male reproductive development and its response to environmental stresses in plants. This review briefly addresses the types, origin, biogenesis and mechanisms of ncRNAs in plants, highlights important updates on the roles of ncRNAs in regulating male reproductive development and emphasizes the contribution of ncRNAs, especially miRNAs and lncRNAs, in responses to abiotic stresses during this unique process in flowering plants.

Transposable elements: multifunctional players in the plant genome

Transposable elements (TEs) are indispensable components of eukaryotic genomes that play diverse roles in gene regulation, recombination, and environmental adaptation. Their ability to mobilize within the genome leads to gene expression and DNA structure changes. TEs serve as valuable markers for genetic and evolutionary studies and facilitate genetic mapping and phylogenetic analysis. They also provide insight into how organisms adapt to a changing environment by promoting gene rearrangements that lead to new gene combinations. These repetitive sequences significantly impact genome structure, function and evolution. This review takes a comprehensive look at TEs and their applications in biotechnology, particularly in the context of plant biology, where they are now considered "genomic gold" due to their extensive functionalities. The article addresses various aspects of TEs in plant development, including their structure, epigenetic regulation, evolutionary patterns, and their use in gene editing and plant molecular markers. The goal is to systematically understand TEs and shed light on their diverse roles in plant biology.

Integrated Transcriptomics and Metabolomics Analysis Promotes the Understanding of Adventitious Root Formation in Eucommia ulmoides Oliver

As a primary approach to nutrient propagation for many woody plants, cutting roots is essential for the breeding and production of Eucommia ulmoides Oliver. In this study, hormone level, transcriptomics, and metabolomics analyses were performed on two E. ulmoides varieties with different adventitious root (AR) formation abilities. The higher JA level on the 0th day and the lower JA level on the 18th day promoted superior AR development. Several hub genes executed crucial roles in the crosstalk regulation of JA and other hormones, including F-box protein (EU012075), SAUR-like protein (EU0125382), LOB protein (EU0124232), AP2/ERF transcription factor (EU0128499), and CYP450 protein (EU0127354). Differentially expressed genes (DEGs) and metabolites of AR formation were enriched in phenylpropanoid biosynthesis, flavonoid biosynthesis, and isoflavonoid biosynthesis pathways. The up-regulated expression of PAL, CCR, CAD, DFR, and HIDH genes on the 18th day could contribute to AR formation. The 130 cis-acting lncRNAs had potential regulatory functions on hub genes in the module that significantly correlated with JA and DEGs in three metabolism pathways. These revealed key molecules, and vital pathways provided more comprehensive insight for the AR formation mechanism of E. ulmoides and other plants.

Conserved and non-conserved RNA-target modules in plants: lessons for a better understanding of Marchantia development

A wide variety of functional regulatory non-coding RNAs (ncRNAs) have been identified as essential regulators of plant growth and development. Depending on their category, ncRNAs are not only involved in modulating target gene expression at the transcriptional and post-transcriptional levels but also are involved in processes like RNA splicing and RNA-directed DNA methylation. To fulfill their molecular roles properly, ncRNAs must be precisely processed by multiprotein complexes. In the case of small RNAs, DICER-LIKE (DCL) proteins play critical roles in the production of mature molecules. Land plant genomes contain at least four distinct classes of DCL family proteins (DCL1-DCL4), of which DCL1, DCL3 and DCL4 are also present in the genomes of bryophytes, indicating the early divergence of these genes. The liverwort Marchantia polymorpha has become an attractive model species for investigating the evolutionary history of regulatory ncRNAs and proteins that are responsible for ncRNA biogenesis. Recent studies on Marchantia have started to uncover the similarities and differences in ncRNA production and function between the basal lineage of bryophytes and other land plants. In this review, we summarize findings on the essential role of regulatory ncRNAs in Marchantia development. We provide a comprehensive overview of conserved ncRNA-target modules among M. polymorpha, the moss Physcomitrium patens and the dicot Arabidopsis thaliana, as well as Marchantia-specific modules. Based on functional studies and data from the literature, we propose new connections between regulatory pathways involved in Marchantia's vegetative and reproductive development and emphasize the need for further functional studies to understand the molecular mechanisms that control ncRNA-directed developmental processes.

Evidence of a novel silencing effect on transgenes in the Arabidopsis thaliana sperm cell

We encountered unexpected transgene silencing in Arabidopsis thaliana sperm cells; transgenes encoding proteins with no specific intracellular localization (cytoplasmic proteins) were silenced transcriptionally or posttranscriptionally. The mRNA of cytoplasmic protein transgenes tagged with a fluorescent protein gene was significantly reduced, resulting in undetectable fluorescent protein signals in the sperm cell. Silencing of the cytoplasmic protein transgenes in the sperm cell did not affect the expression of either its endogenous homologous genes or cotransformed transgenes encoding a protein with targeted intracellular localization. This transgene silencing in the sperm cell persisted in mutants of the major gene silencing machinery including DNA methylation. The incomprehensible, yet real, transgene silencing phenotypes occurring in the sperm cell could mislead the interpretation of experimental results in plant reproduction, and this Commentary calls attention to that risk and highlights details of this novel cytoplasmic protein transgene silencing.

Potential role of the regulatory miR1119-MYC2 module in wheat (Triticum aestivum L.) drought tolerance

MicroRNA (miRNA)-target gene modules are essential components of plants' abiotic stress signalling pathways Little is known about the drought-responsive miRNA-target modules in wheat, but systems biology approaches have enabled the prediction of these regulatory modules and systematic study of their roles in responses to abiotic stresses. Using such an approach, we sought miRNA-target module(s) that may be differentially expressed under drought and non-stressed conditions by mining Expressed Sequence Tag (EST) libraries of wheat roots and identified a strong candidate (miR1119-MYC2). We then assessed molecular and physiochemical differences between two wheat genotypes with contrasting drought tolerance in a controlled drought experiment and assessed possible relationships between their tolerance and evaluated traits. We found that the miR1119-MYC2 module significantly responds to drought stress in wheat roots. It is differentially expressed between the contrasting wheat genotypes and under drought versus non-stressed conditions. We also found significant associations between the module's expression profiles and ABA hormone content, water relations, photosynthetic activities, H₂O₂ levels, plasma membrane damage, and antioxidant enzyme activities in wheat. Collectively, our results suggest that a regulatory module consisting of miR1119 and MYC2 may play an important role in wheat's drought tolerance.

De novo full-length transcriptome analysis of two ecotypes of Phragmites australis (swamp reed and dune reed) provides new insights into the transcriptomic complexity of dune reed and its long-term adaptation to desert environments

BACKGROUND

The extremely harsh environment of the desert is changing dramatically every moment, and the rapid adaptive stress response in the short term requires enormous energy expenditure to mobilize widespread regulatory networks, which is all the more detrimental to the survival of the desert plants themselves. The dune reed, which has adapted to desert environments with complex and variable ecological factors, is an ideal type of plant for studying the molecular mechanisms by which Gramineae plants respond to combinatorial stress of the desert in their natural state. But so far, the data on the genetic resources of reeds is still scarce, therefore most of their research has focused on ecological and physiological studies.

RESULTS

In this study, we obtained the first De novo non-redundant Full-Length Non-Chimeric (FLNC) transcriptome databases for swamp reeds (SR), dune reeds (DR) and the All of Phragmites australis (merged of iso-seq data from SR and DR), using PacBio Iso-Seq technology and combining tools such as Iso-Seq3 and Cogent. We then identified and described long non-coding RNAs (LncRNA), transcription factor (TF) and alternative splicing (AS) events in reeds based on a transcriptome database. Meanwhile, we have identified and developed for the first time a large number of candidates expressed sequence tag-SSR (EST-SSRs) markers in reeds based on UniTransModels. In addition, through differential gene expression analysis of wild-type and homogenous cultures, we found a large number of transcription factors that may be associated with desert stress tolerance in the dune reed, and revealed that members of the Lhc family have an important role in the long-term adaptation of dune reeds to desert environments.

CONCLUSIONS

Our results provide a positive and usable genetic resource for Phragmites australis with a widespread adaptability and resistance, and provide a genetic database for subsequent reeds genome annotation and functional genomic studies.

Regulatory Networks of lncRNAs, miRNAs, and mRNAs in Response to Heat Stress in Wheat (Triticum Aestivum L.): An Integrated Analysis

Climate change has become a major source of concern, particularly in agriculture, because it has a significant impact on the production of economically important crops such as wheat, rice, and maize. In the present study, an attempt has been made to identify differentially expressed heat stress-responsive long non-coding RNAs (lncRNAs) in the wheat genome using publicly available wheat transcriptome data (24 SRAs) representing two conditions, namely, control and heat-stressed. A total of 10,965 lncRNAs have been identified and, among them, 153, 143, and 211 differentially expressed transcripts have been found under 0 DAT, 1 DAT, and 4 DAT heat-stress conditions, respectively. Target prediction analysis revealed that 4098 lncRNAs were targeted by 119 different miRNA responses to a plethora of environmental stresses, including heat stress. A total of 171 hub genes had 204 SSRs (simple sequence repeats), and a set of target sequences had SNP potential as well. Furthermore, gene ontology analysis revealed that the majority of the discovered lncRNAs are engaged in a variety of cellular and biological processes related to heat stress responses. Furthermore, the modeled three-dimensional (3D) structures of hub genes encoding proteins, which had an appropriate range of similarity with solved structures, provided information on their structural roles. The current study reveals many elements of gene expression regulation in wheat under heat stress, paving the way for the development of improved climate-resilient wheat cultivars.

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

Flooding and herbivory: the effect of concurrent stress factors on plant volatile emissions and gene expression in two heirloom tomato varieties

BACKGROUND

In nature and in cultivated fields, plants encounter multiple stress factors. Nonetheless, our understanding of how plants actively respond to combinatorial stress remains limited. Among the least studied stress combination is that of flooding and herbivory, despite the growing importance of these stressors in the context of climate change. We investigated plant chemistry and gene expression changes in two heirloom tomato varieties: Cherokee Purple (CP) and Striped German (SG) in response to flooding, herbivory by Spodoptera exigua, and their combination.

RESULTS

Volatile organic compounds (VOCs) identified in tomato plants subjected to flooding and/or herbivory included several mono- and sesquiterpenes. Flooding was the main factor altering VOCs emission rates, and impacting plant biomass accumulation, while different varieties had quantitative differences in their VOC emissions. At the gene expression levels, there were 335 differentially expressed genes between the two tomato plant varieties, these included genes encoding for phenylalanine ammonia-lyase (PAL), cinnamoyl-CoA-reductase-like, and phytoene synthase (Psy1). Flooding and variety effects together influenced abscisic acid (ABA) signaling genes with the SG variety showing higher levels of ABA production and ABA-dependent signaling upon flooding. Flooding downregulated genes associated with cytokinin catabolism and general defense response and upregulated genes associated with ethylene biosynthesis, anthocyanin biosynthesis, and gibberellin biosynthesis. Combining flooding and herbivory induced the upregulation of genes including chalcone synthase (CHS), PAL, and genes encoding BAHD acyltransferase and UDP-glucose iridoid glucosyltransferase-like genes in one of the tomato varieties (CP) and a disproportionate number of heat-shock proteins in SG. Only the SG variety had measurable changes in gene expression due to herbivory alone, upregulating zeatin, and O-glucosyltransferase and thioredoxin among others.

CONCLUSION

Our results suggest that both heirloom tomato plant varieties differ in their production of secondary metabolites including phenylpropanoids and terpenoids and their regulation and activation of ABA signaling upon stress associated with flooding. Herbivory and flooding together had interacting effects that were evident at the level of plant chemistry (VOCs production), gene expression and biomass markers. Results from our study highlight the complex nature of plant responses to combinatorial stresses and point at specific genes and pathways that are affected by flooding and herbivory combined.

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.

Functional Characterization and Molecular Marker Development of the Proenkephalin as Biomarker of Food Addiction in Food Habit Domestication of Mandarin Fish (Siniperca Chuatsi)

Proenkephalin (PENK), as the precursor of endogenous opioid enkephalin (ENK), is widely present in the nervous system and plays an important role in animal food addiction and rewarding behavior. In our study, we intend to study the functional characterization and molecular marker development of the penk gene related to food habit domestication of mandarin fish. We found that the penk gene of mandarin fish had three types of endogenous opioid peptide sequences. Compared with other tissues, penk mRNA was highly expressed in the whole brain. Intracerebroventricular (ICV) injection of lysine or methionine significantly increased the expression of penk mRNA. The expression of penk mRNA in the brain of mandarin fish that could be easily domesticated from eating live prey fish to artificial diets was significantly higher than those that could not. After feeding with high-carbohydrate artificial diets, the expression of penk mRNA showed no significant difference between mandarin fish with hypophagia and those that still ate normally. A total of four single nucleotide polymorphisms (SNP) loci related to easy domestication toward eating artificial diets were screened from the mandarin fish population. Additionally, the TT genotype at one of the loci was significantly correlated with the food habit domestication of mandarin fish.

Moving Beyond DNA Sequence to Improve Plant Stress Responses

Plants offer a habitat for a range of interactions to occur among different stress factors. Epigenetics has become the most promising functional genomics tool, with huge potential for improving plant adaptation to biotic and abiotic stresses. Advances in plant molecular biology have dramatically changed our understanding of the molecular mechanisms that control these interactions, and plant epigenetics has attracted great interest in this context. Accumulating literature substantiates the crucial role of epigenetics in the diversity of plant responses that can be harnessed to accelerate the progress of crop improvement. However, harnessing epigenetics to its full potential will require a thorough understanding of the epigenetic modifications and assessing the functional relevance of these variants. The modern technologies of profiling and engineering plants at genome-wide scale provide new horizons to elucidate how epigenetic modifications occur in plants in response to stress conditions. This review summarizes recent progress on understanding the epigenetic regulation of plant stress responses, methods to detect genome-wide epigenetic modifications, and disentangling their contributions to plant phenotypes from other sources of variations. Key epigenetic mechanisms underlying stress memory are highlighted. Linking plant response with the patterns of epigenetic variations would help devise breeding strategies for improving crop performance under stressed scenarios.

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

Transposable Elements and Human Diseases: Mechanisms and Implication in the Response to Environmental Pollutants

Transposable elements (TEs) are recognized as major players in genome plasticity and evolution. The high abundance of TEs in the human genome, especially the Alu and Long Interspersed Nuclear Element-1 (LINE-1) repeats, makes them responsible for the molecular origin of several diseases. This involves several molecular mechanisms that are presented in this review: insertional mutation, DNA recombination and chromosomal rearrangements, modification of gene expression, as well as alteration of epigenetic regulations. This literature review also presents some of the more recent and/or more classical examples of human diseases in which TEs are involved. Whether through insertion of LINE-1 or Alu elements that cause chromosomal rearrangements, or through epigenetic modifications, TEs are widely implicated in the origin of human cancers. Many other human diseases can have a molecular origin in TE-mediated chromosomal recombination or alteration of gene structure and/or expression. These diseases are very diverse and include hemoglobinopathies, metabolic and neurological diseases, and common diseases. Moreover, TEs can also have an impact on aging. Finally, the exposure of individuals to stresses and environmental contaminants seems to have a non-negligible impact on the epigenetic derepression and mobility of TEs, which can lead to the development of diseases. Thus, improving our knowledge of TEs may lead to new potential diagnostic markers of diseases.

An Epigenetic Alphabet of Crop Adaptation to Climate Change

Crop adaptation to climate change is in a part attributed to epigenetic mechanisms which are related to response to abiotic and biotic stresses. Although recent studies increased our knowledge on the nature of these mechanisms, epigenetics remains under-investigated and still poorly understood in many, especially non-model, plants, Epigenetic modifications are traditionally divided into two main groups, DNA methylation and histone modifications that lead to chromatin remodeling and the regulation of genome functioning. In this review, we outline the most recent and interesting findings on crop epigenetic responses to the environmental cues that are most relevant to climate change. In addition, we discuss a speculative point of view, in which we try to decipher the “epigenetic alphabet” that underlies crop adaptation mechanisms to climate change. The understanding of these mechanisms will pave the way to new strategies to design and implement the next generation of cultivars with a broad range of tolerance/resistance to stresses as well as balanced agronomic traits, with a limited loss of (epi)genetic variability.

Computational Analysis of Transposable Elements and CircRNAs in Plants

This chapter provides two main contributions: (1) a description of computational tools and databases used to identify and analyze transposable elements (TEs) and circRNAs in plants; and (2) data analysis on public TE and circRNA data. Our goal is to highlight the primary information available in the literature on circular noncoding RNAs and transposable elements in plants. The exploratory analysis performed on publicly available circRNA and TEs data help discuss four sequence features. Finally, we investigate the association on circRNAs:TE in plants in the model organism Arabidopsis thaliana.

RAD gene family analysis in cotton provides some key genes for flowering and stress tolerance in upland cotton G. hirsutum

Background

RADIALIS (RAD), belongs to the MYB gene family and regulates a variety of functions including floral dorsoventral asymmetry in Antirrhinum majus and development of fruit proteins in Solanum lycopersicum. RAD genes contain an SNF2_N superfamily domain. Here, we comprehensively identified 68 RAD genes from six different species including Arabidopsis and five species of cotton.

Results

Phylogenetic analysis classified RAD genes into five groups. Gene structure, protein motifs and conserved amino acid residues indicated that GhRAD genes were highly conserved during the evolutionary process. Chromosomal location information showed that GhRAD genes were distributed unevenly on different chromosomes. Collinearity and selection pressure analysis indicated RAD gene family expansion in G. hirsutum and G. barbadense with purifying selection pressure. Further, various growth and stress related promotor cis-acting elements were observed. Tissue specific expression level indicated that most GhRAD genes were highly expressed in roots and flowers (GhRAD2, GhRAD3, GhRAD4 and GhRAD11). Next, GhRAD genes were regulated by phytohormonal stresses (JA, BL and IAA). Moreover, Ghi-miRN1496, Ghi-miR1440, Ghi-miR2111b, Ghi-miR2950a, Ghi-miR390a, Ghi-miR390b and Ghi-miR7495 were the miRNAs targeting most of GhRAD genes.

Conclusions

Our study revealed that RAD genes are evolutionary conserved and might be involved in different developmental processes and hormonal stress response. Data presented in our study could be used as the basis for future studies of RAD genes in cotton.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08248-z.

Dissection of full-length transcriptome and metabolome of Dichocarpum (Ranunculaceae): implications in evolution of specialized metabolism of Ranunculales medicinal plants

Several main families of Ranunculales are rich in alkaloids and other medicinal compounds; many species of these families are used in traditional and folk medicine. Dichocarpum is a representative medicinal genus of Ranunculaceae, but the genetic basis of its metabolic phenotype has not been investigated, which hinders its sustainable conservation and utilization. We use the third-generation high-throughput sequencing and metabolomic techniques to decipher the full-length transcriptomes and metabolomes of five Dichocarpum species endemic in China, and 71,598 non-redundant full-length transcripts were obtained, many of which are involved in defense, stress response and immunity, especially those participating in the biosynthesis of specialized metabolites such as benzylisoquinoline alkaloids (BIAs). Twenty-seven orthologs extracted from trancriptome datasets were concatenated to reconstruct the phylogenetic tree, which was verified by the clustering analysis based on the metabolomic profile and agreed with the Pearson correlation between gene expression patterns of Dichocarpum species. The phylogenomic analysis of phytometabolite biosynthesis genes, e.g., (S)-norcoclaurine synthase, methyltransferases, cytochrome p450 monooxygenases, berberine bridge enzyme and (S)-tetrahydroprotoberberine oxidase, revealed the evolutionary trajectories leading to the chemodiversity, especially that of protoberberine type, aporphine type and bis-BIA abundant in Dichocarpum and related genera. The biosynthesis pathways of these BIAs are proposed based on full-length transcriptomes and metabolomes of Dichocarpum. Within Ranunculales, the gene duplications are common, and a unique whole genome duplication is possible in Dichocarpum. The extensive correlations between metabolite content and gene expression support the co-evolution of various genes essential for the production of different specialized metabolites. Our study provides insights into the transcriptomic and metabolomic landscapes of Dichocarpum, which will assist further studies on genomics and application of Ranunculales plants.

Intergenic spaces: a new frontier to improving plant health

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

The Evolutionary Volte-Face of Transposable Elements: From Harmful Jumping Genes to Major Drivers of Genetic Innovation

Transposable elements (TEs) are self-replicating DNA elements that constitute major fractions of eukaryote genomes. Their ability to transpose can modify the genome structure with potentially deleterious effects. To repress TE activity, host cells have developed numerous strategies, including epigenetic pathways, such as DNA methylation or histone modifications. Although TE neo-insertions are mostly deleterious or neutral, they can become advantageous for the host under specific circumstances. The phenomenon leading to the appropriation of TE-derived sequences by the host is known as TE exaptation or co-option. TE exaptation can be of different natures, through the production of coding or non-coding DNA sequences with ultimately an adaptive benefit for the host. In this review, we first give new insights into the silencing pathways controlling TE activity. We then discuss a model to explain how, under specific environmental conditions, TEs are unleashed, leading to a TE burst and neo-insertions, with potential benefits for the host. Finally, we review our current knowledge of coding and non-coding TE exaptation by providing several examples in various organisms and describing a method to identify TE co-option events.

Function identification of miR394 in tomato resistance to Phytophthora infestans

MiR394 plays a negative role in tomato resistance to late blight. The lncRNA40787 severing as an eTM for miR394 to regulate LCR and exerting functions in tomato resistance. Tomato (Solanum lycopersicum), which was used as model species for studying the mechanism of plant disease defense, is susceptible to multiple pathogens. Non-coding RNA (ncRNA) has a pivotal role in plants response to biological stresses. It has previously been observed that the expression level of miR394 changed significantly after the infection of various pathogens. However, there has been no detailed investigation of the accumulated or suppressed mechanism of miR394. Our previous study predicted three lncRNAs (lncRNA40787, lncRNA27177, and lncRNA42566) that contain miR394 endogenous target mimics (eTM), which may exist as the competitive endogenous RNAs (ceRNAs) of miR394. In our study, the transcription levels of these three lncRNAs were strongly up-regulated in tomato upon infection with P. infestans. In contrast with the three lncRNAs, the accumulation of miR394 was significantly suppressed. Based on the expression pattern, and value of minimum free energy (mfes) that represents the binding ability between lncRNA and miRNA, lncRNA40787 was chosen for further investigation. Results showed that overexpression of lncRNA40787 reduced the expression of miR394 along with decreased lesion area and enhanced disease resistance. Overexpression of miR394, however, decreased the expression of its target gene Leaf Curling Responsiveness (LCR), and suppressed the synthesis components genes of jasmonic acid (JA), depressing the resistance of tomato to P. infestans infection. Taken together, our findings indicated that miR394 can be decoyed by lncRNA40787, and negatively regulated the expression of LCR to enhance tomato susceptibility under P. infestans infection. Our study provided detailed information on the lncRNA40787-miR394-LCR regulatory network and serves as a reference for future research.

Transition from Seeds to Seedlings: Hormonal and Epigenetic Aspects

Transition from seed to seedling is one of the critical developmental steps, dramatically affecting plant growth and viability. Before plants enter the vegetative phase of their ontogenesis, massive rearrangements of signaling pathways and switching of gene expression programs are required. This results in suppression of the genes controlling seed maturation and activation of those involved in regulation of vegetative growth. At the level of hormonal regulation, these events are controlled by the balance of abscisic acid and gibberellins, although ethylene, auxins, brassinosteroids, cytokinins, and jasmonates are also involved. The key players include the members of the LAFL network-the transcription factors LEAFY COTYLEDON1 and 2 (LEC 1 and 2), ABSCISIC ACID INSENSITIVE3 (ABI3), and FUSCA3 (FUS3), as well as DELAY OF GERMINATION1 (DOG1). They are the negative regulators of seed germination and need to be suppressed before seedling development can be initiated. This repressive signal is mediated by chromatin remodeling complexes-POLYCOMB REPRESSIVE COMPLEX 1 and 2 (PRC1 and PRC2), as well as PICKLE (PKL) and PICKLE-RELATED2 (PKR2) proteins. Finally, epigenetic methylation of cytosine residues in DNA, histone post-translational modifications, and post-transcriptional downregulation of seed maturation genes with miRNA are discussed. Here, we summarize recent updates in the study of hormonal and epigenetic switches involved in regulation of the transition from seed germination to the post-germination stage.

Genomic reconfiguration in parasitic plants involves considerable gene losses alongside global genome size inflation and gene births

Parasitic plant genomes and transcriptomes reveal numerous genetic innovations, the functional-evolutionary relevance and roles of which open unprecedented research avenues.

Distribution of MITE family Monkey King in rapeseed (Brassica napus L) and its influence on gene expression

Miniature inverted-repeat transposable elements (MITEs) are a group of class II transposable elements. The MITE Monkey King (MK) was first discovered upstream of BnFLC.A10. In this study, genome resequencing of four selected B. napus accessions, revealed more than 4000 distributed copies of MKs constituting ~2.4 Mb of the B. napus genomic sequence and caused 677 polymorphisms among the four accessions. MK -polymorphism-related markers across 128 natural and 58 synthetic accessions revealed more polymorphic MKs in natural than synthetic accessions. Ten MK -induced indels significantly affected the expression levels of the nearest gene based on RNAseq analysis, six of these effects were subsequently confirmed using qRT-PCR. Decreased expression pattern of MK -derived miRNA-bna-miR6031 was also observed under various stress treatments. Further research focused on the MITE families should promote not only our understanding of gene regulatory networks but also inform crop improvement efforts.

A Practical Guide on Computational Tools and Databases for Transposable Elements in Plants

In the age of big data, obtaining precise information about the research topic of interesting is extremely important. Keeping this in mind, this chapter focuses on providing a practical knowledge guide about computational tools and databases of transposable elements (TE) in plants. For that, we organize and present this text in three sections: (1) a discussion about tools and databases on this theme; (2) hands-on of how to use a few of them; (3) an exploratory data analysis on public TE data. Finally, we are going deep to present the main challenges and possible solutions to improve resources and tools.

When junk DNA turns functional: transposon-derived non-coding RNAs in plants

Transposable elements (TEs) are major contributors to genome complexity in eukaryotes. TE mobilization may cause genome instability, although it can also drive genome diversity throughout evolution. TE transposition may influence the transcriptional activity of neighboring genes by modulating the epigenomic profile of the region or by altering the relative position of regulatory elements. Notably, TEs have emerged in the last few years as an important source of functional long and small non-coding RNAs. A plethora of small RNAs derived from TEs have been linked to the trans regulation of gene activity at the transcriptional and post-transcriptional levels. Furthermore, TE-derived long non-coding RNAs have been shown to modulate gene expression by interacting with protein partners, sequestering active small RNAs, and forming duplexes with DNA or other RNA molecules. In this review, we summarize our current knowledge of the functional and mechanistic paradigms of TE-derived long and small non-coding RNAs and discuss their role in plant development and evolution.

Role of non-coding RNAs in plant immunity

Crops are exposed to attacks by various pathogens that cause substantial yield losses and severely threaten food security. To cope with pathogenic infection, crops have elaborated strategies to enhance resistance against pathogens. In addition to the role of protein-coding genes as key regulators in plant immunity, accumulating evidence has demonstrated the importance of non-coding RNAs (ncRNAs) in the plant immune response. Here, we summarize the roles and molecular mechanisms of endogenous ncRNAs, especially microRNAs (miRNAs), long ncRNAs (lncRNAs), and circular RNAs (circRNAs), in plant immunity. We discuss the coordination between miRNAs and small interfering RNAs (siRNAs), between lncRNAs and miRNAs or siRNAs, and between circRNAs and miRNAs in the regulation of plant immune responses. We also address the role of cross-kingdom mobile small RNAs in plant-pathogen interactions. These insights improve our understanding of the mechanisms by which ncRNAs regulate plant immunity and can promote the development of better approaches for breeding disease-resistant crops.

Pan-genome: A promising resource for noncoding RNA discovery in plants

Plant genomes contain both protein-coding and noncoding sequences including transposable elements (TEs) and noncoding RNAs (ncRNAs). The ncRNAs are recognized as important elements that play fundamental roles in the structural organization and function of plant genomes. Despite various hypotheses, TEs are believed to be a major precursor of ncRNAs. Transposable elements are also prime factors that cause genomic variation among members of a species. Hence, TEs pose a major challenge in the discovery and analysis of ncRNAs. With the increase in the number of sequenced plant genomes, it is now accepted that a single reference genome is insufficient to represent the complete genomic diversity and contents of a species, and exploring the pan-genome of a species is critical. In this review, we summarize the recent progress in the field of plant pan-genomes. We also discuss TEs and their roles in ncRNA biogenesis and present our perspectives on the application of pan-genomes for the discovery of ncRNAs to fully explore and exploit their biological roles in plants.

Tissue-specific usage of transposable element-derived promoters in mouse development

BACKGROUND

Transposable elements (TEs) are a significant component of eukaryotic genomes and play essential roles in genome evolution. Mounting evidence indicates that TEs are highly transcribed in early embryo development and contribute to distinct biological functions and tissue morphology.

RESULTS

We examine the epigenetic dynamics of mouse TEs during the development of five tissues: intestine, liver, lung, stomach, and kidney. We found that TEs are associated with over 20% of open chromatin regions during development. Close to half of these accessible TEs are only activated in a single tissue and a specific developmental stage. Most accessible TEs are rodent-specific. Across these five tissues, 453 accessible TEs are found to create the transcription start sites of downstream genes in mouse, including 117 protein-coding genes and 144 lincRNA genes, 93.7% of which are mouse-specific. Species-specific TE-derived transcription start sites are found to drive the expression of tissue-specific genes and change their tissue-specific expression patterns during evolution.

CONCLUSION

Our results suggest that TE insertions increase the regulatory potential of the genome, and some TEs have been domesticated to become a crucial component of gene and regulate tissue-specific expression during mouse tissue development.

NAMS webserver: coding potential assessment and functional annotation of plant transcripts

Recent advances in transcriptomics have uncovered lots of novel transcripts in plants. To annotate such transcripts, dissecting their coding potential is a critical step. Computational approaches have been proven fruitful in this task; however, most current tools are designed/optimized for mammals and only a few of them have been tested on a limited number of plant species. In this work, we present NAMS webserver, which contains a novel coding potential classifier, NAMS, specifically optimized for plants. We have evaluated the performance of NAMS using a comprehensive dataset containing more than 3 million transcripts from various plant species, where NAMS demonstrates high accuracy and remarkable performance improvements over state-of-the-art software. Moreover, our webserver also furnishes functional annotations, aiming to provide users informative clues to the functions of their transcripts. Considering that most plant species are poorly characterized, our NAMS webserver could serve as a valuable resource to facilitate the transcriptomic studies. The webserver with testing dataset is freely available at http://sunlab.cpy.cuhk.edu.hk/NAMS/.

First Come, First Served: Sui Generis Features of the First Intron

Most of the transcribed genes in eukaryotic cells are interrupted by intervening sequences called introns that are co-transcriptionally removed from nascent messenger RNA through the process of splicing. In Arabidopsis, 79% of genes contain introns and more than 60% of intron-containing genes undergo alternative splicing (AS), which ostensibly is considered to increase protein diversity as one of the intrinsic mechanisms for fitness to the varying environment or the internal developmental program. In addition, recent findings have prevailed in terms of overlooked intron functions. Here, we review recent progress in the underlying mechanisms of intron function, in particular by focusing on unique features of the first intron that is located in close proximity to the transcription start site. The distinct deposition of epigenetic marks and nucleosome density on the first intronic DNA sequence, the impact of the first intron on determining the transcription start site and elongation of its own expression (called intron-mediated enhancement, IME), translation control in 5′-UTR, and the new mechanism of the trans-acting function of the first intron in regulating gene expression at the post-transcriptional level are summarized.

Identification and characterization of long non-coding RNA (lncRNA) in the developing seeds of Jatropha curcas

Long non-coding RNAs (lncRNAs) play critical roles in plant development. However, the information of lncRNAs in Jatropha curcas remains largely unexplored. Thus, an attempt has been made in J. curcas to identify 1,850 lncRNAs based on deep sequencing of developing seeds at three typical stages. About ten percent lncRNAs (196 lncRNAs) were differentially expressed lncRNAs during seed developing process. Together with reverse transcription quantitative real-time PCR, the lncRNA expression analyses revealed the stage-specific expression patterns of some novel lncRNAs in J. curcas. The target genes of lncRNAs were annotated for their roles in various biological processes such as gene expression, metabolism, and cell growth. Besides, 10 lncRNAs were identified as the precursors of microRNAs and 26 lncRNAs were predicted to be the targets of Jatropha miRNAs. A total of 31 key lncRNAs play critical roles in the seed developing process in the context of cell growth and development, lipid metabolism, and seed maturation. Our study provides the first systematic study of lncRNAs in the developing seeds of J. curcas and facilitates the functional research of plant lncRNAs and the regulation of seed development.

Plasmogamic Paternal Contributions to Early Zygotic Development in Flowering Plants

Flowering plant zygotes possess complete developmental potency, and the mixture of male and female genetic and cytosolic materials in the zygote is a trigger to initiate embryo development. Plasmogamy, the fusion of the gamete cytoplasms, facilitates the cellular dynamics of the zygote. In the last decade, mutant analyses, live cell imaging-based observations, and direct observations of fertilized egg cells by in vitro fusion of isolated gametes have accelerated our understanding of the post-plasmogamic events in flowering plants including cell wall formation, gamete nuclear migration and fusion, and zygotic cell elongation and asymmetric division. Especially, it has become more evident that paternal parent-of-origin effects, via sperm cytoplasm contents, not only control canonical early zygotic development, but also activate a biparental signaling pathway critical for cell fate determination after the first cell division. Here, we summarize the plasmogamic paternal contributions via the entry of sperm contents during/after fertilization in flowering plants.

Cobalt-induced retrotransposon polymorphism and humic acid protection on maize genome

Retrotransposon activity and genomic template stability (GTS) are one of the most significant rearranging mechanisms in environmental stress. Therefore, in this study, it is aimed to elucidate effecting of Cobalt (Co) on the instability of genomes and Long Terminal Repeat retrotransposon polymorphism in Zea mays and whether humic acid (HA) has any role on these parameters. For this purpose, Retrotransposon-microsatellite amplified polymorphism (REMAP) and Inter-Retrotransposon Amplified Polymorphism (IRAP) markers were applied to evaluate retrotransposon polymorphism and the GTS levels. It was found that IRAP and REMAP primers generate unique polymorphic band structures on maize plants treated with various doses of Co. Retrotransposon polymorphism increased and GTS decreased while increasing Co concentration. On the other hand, there was a reduction in negative effects of Co on retrotransposon GTS and polymorphism after treatment with HA. The results indicate that HA may be used effectively for the protection of maize seedlings from the destructive effects of Co toxicity.

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.

Comparative Analysis of miRNA Expression Profiles Between Heat-Tolerant and Heat-Sensitive Genotypes of Flowering Chinese Cabbage Under Heat Stress Using High-Throughput Sequencing

Heat stress disturbs cellular homeostasis, thus usually impairs yield of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). MicroRNAs (miRNAs) play a significant role in plant responses to different stresses by modulating gene expression at the post-transcriptional level. However, the roles that miRNAs and their target genes may play in heat tolerance of flowering Chinese cabbage remain poorly characterized. The current study sequenced six small RNA libraries generated from leaf tissues of flowering Chinese cabbage collected at 0, 6, and 12 h after 38 °C heat treatment, and identified 49 putative novel miRNAs and 43 known miRNAs that differentially expressed between heat-tolerant and heat-sensitive flowering Chinese cabbage. Among them, 14 novel and nine known miRNAs differentially expressed only in the heat-tolerant genotype under heat-stress, therefore, their target genes including disease resistance protein TAO1-like, RPS6, reticuline oxidase-like protein, etc. might play important roles in enhancing heat-tolerance. Gene Ontology (GO) analysis revealed that targets of these differentially expressed miRNAs may play key roles in responses to temperature stimulus, cell part, cellular process, cell, membrane, biological regulation, binding, and catalytic activities. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis identified their important functions in signal transduction, environmental adaptation, global and overview maps, as well as in stress adaptation and in MAPK signaling pathways such as cell death. These findings provide insight into the functions of the miRNAs in heat stress tolerance of flowering Chinese cabbage.

Engineering Multiple Abiotic Stress Tolerance in Canola, Brassica napus

Impacts of climate change like global warming, drought, flooding, and other extreme events are posing severe challenges to global crop production. Contribution of Brassica napus towards the oilseed industry makes it an essential component of international trade and agroeconomics. Consequences from increasing occurrences of multiple abiotic stresses on this crop are leading to agroeconomic losses making it vital to endow B. napus crop with an ability to survive and maintain yield when faced with simultaneous exposure to multiple abiotic stresses. For an improved understanding of the stress sensing machinery, there is a need for analyzing regulatory pathways of multiple stress-responsive genes and other regulatory elements such as non-coding RNAs. However, our understanding of these pathways and their interactions in B. napus is far from complete. This review outlines the current knowledge of stress-responsive genes and their role in imparting multiple stress tolerance in B. napus. Analysis of network cross-talk through omics data mining is now making it possible to unravel the underlying complexity required for stress sensing and signaling in plants. Novel biotechnological approaches such as transgene-free genome editing and utilization of nanoparticles as gene delivery tools are also discussed. These can contribute to providing solutions for developing climate change resilient B. napus varieties with reduced regulatory limitations. The potential ability of synthetic biology to engineer and modify networks through fine-tuning of stress regulatory elements for plant responses to stress adaption is also highlighted.

Insight into the Role of Epigenetic Processes in Abiotic and Biotic Stress Response in Wheat and Barley

Environmental stresses such as salinity, drought, heat, freezing, heavy metal and even pathogen infections seriously threaten the growth and yield of important cereal crops including wheat and barley. There is growing evidence indicating that plants employ sophisticated epigenetic mechanisms to fine-tune their responses to environmental stresses. Here, we provide an overview of recent developments in understanding the epigenetic processes and elements—such as DNA methylation, histone modification, chromatin remodeling, and non-coding RNAs—involved in plant responses to abiotic and biotic stresses in wheat and barley. Potentials of exploiting epigenetic variation for the improvement of wheat and barley are discussed.

Machine learning improves our knowledge about miRNA functions towards plant abiotic stresses

During the last two decades, human has increased his knowledge about the role of miRNAs and their target genes in plant stress response. Biotic and abiotic stresses result in simultaneous tissue-specific up/down-regulation of several miRNAs. In this study, for the first time, feature selection algorithms have been used to investigate the contribution of individual plant miRNAs in Arabidopsis thaliana response towards different levels of several abiotic stresses including drought, salinity, cold, and heat. Results of information theory-based feature selection revealed that miRNA-169, miRNA-159, miRNA-396, and miRNA-393 had the highest contributions to plant response towards drought, salinity, cold, and heat, respectively. Furthermore, regression models, i.e., decision tree (DT), support vector machines (SVMs), and Naïve Bayes (NB) were used to predict the plant stress by having the plant miRNAs' concentration. SVM with Gaussian kernel was capable of predicting plant stress (R2 = 0.96) considering miRNA concentrations as input features. Findings of this study prove the performance of machine learning as a promising tool to investigate some aspects of miRNAs' contribution to plant stress responses that have been undiscovered until today.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

Whole-transcriptome RNA sequencing reveals the global molecular responses and ceRNA regulatory network of mRNAs, lncRNAs, miRNAs and circRNAs in response to copper toxicity in Ziyang Xiangcheng (Citrus junos Sieb. Ex Tanaka)

BACKGROUND

Copper (Cu) toxicity has become a potential threat for citrus production, but little is known about related mechanisms. This study aims to uncover the global landscape of mRNAs, long non-coding RNAs (lncRNAs), circular RNAs (circRNAs) and microRNAs (miRNAs) in response to Cu toxicity so as to construct a regulatory network of competing endogenous RNAs (ceRNAs) and to provide valuable knowledge pertinent to Cu response in citrus.

RESULTS

Tolerance of four commonly used rootstocks to Cu toxicity was evaluated, and 'Ziyang Xiangcheng' (Citrus junos) was found to be the most tolerant genotype. Then the roots and leaves sampled from 'Ziyang Xiangcheng' with or without Cu treatment were used for whole-transcriptome sequencing. In total, 5734 and 222 mRNAs, 164 and 5 lncRNAs, 45 and 17 circRNAs, and 147 and 130 miRNAs were identified to be differentially expressed (DE) in Cu-treated roots and leaves, respectively, in comparison with the control. Gene ontology enrichment analysis showed that most of the DEmRNAs and targets of DElncRNAs and DEmiRNAs were annotated to the categories of 'oxidation-reduction', 'phosphorylation', 'membrane', and 'ion binding'. The ceRNA network was then constructed with the predicted pairs of DEmRNAs-DEmiRNAs and DElncRNAs-DEmiRNAs, which further revealed regulatory roles of these DERNAs in Cu toxicity.

CONCLUSIONS

A large number of mRNAs, lncRNAs, circRNAs, and miRNAs in 'Ziyang Xiangcheng' were altered in response to Cu toxicity, which may play crucial roles in mitigation of Cu toxicity through the ceRNA regulatory network in this Cu-tolerant rootstock.