Temperature expression patterns of genes and their coexpression with LncRNAs revealed by RNA-Seq in non-heading Chinese cabbage

Emerging roles and mechanisms of lncRNAs in fruit and vegetables

With the development of genome sequencing technologies, many long non-coding RNAs (lncRNAs) have been identified in fruit and vegetables. lncRNAs are primarily transcribed and spliced by RNA polymerase II (Pol II) or plant-specific Pol IV/V, and exhibit limited evolutionary conservation. lncRNAs intricately regulate various aspects of fruit and vegetables, including pigment accumulation, reproductive tissue development, fruit ripening, and responses to biotic and abiotic stresses, through diverse mechanisms such as gene expression modulation, interaction with hormones and transcription factors, microRNA regulation, and involvement in alternative splicing. This review presents a comprehensive overview of lncRNA classification, basic characteristics, and, most importantly, recent advances in understanding their functions and regulatory mechanisms.

Transcriptome Analysis Reveals the Mechanism by Which Exogenous Melatonin Treatment Delays Leaf Senescence of Postharvest Chinese Kale (Brassica oleracea var. alboglabra)

Melatonin, a pleiotropic small molecule, is employed in horticultural crops to delay senescence and preserve postharvest quality. In this study, 100 µM melatonin treatment delayed a decline in the color difference index h* and a*, maintaining the content of chlorophyll and carotenoids, thereby delaying the yellowing and senescence of Chinese kale. Transcriptome analysis unequivocally validates melatonin's efficacy in delaying leaf senescence in postharvest Chinese kale stored at 20 °C. Following a three-day storage period, the melatonin treatment group exhibited 1637 differentially expressed genes (DEGs) compared to the control group. DEG analysis elucidated that melatonin-induced antisenescence primarily governs phenylpropanoid biosynthesis, lipid metabolism, plant signal transduction, and calcium signal transduction. Melatonin treatment up-regulated core enzyme genes associated with general phenylpropanoid biosynthesis, flavonoid biosynthesis, and the α-linolenic acid biosynthesis pathway. It influenced the redirection of lignin metabolic flux, suppressed jasmonic acid and abscisic acid signal transduction, and concurrently stimulated auxin signal transduction. Additionally, melatonin treatment down-regulated RBOH expression and up-regulated genes encoding CaM, thereby influencing calcium signal transduction. This study underscores melatonin as a promising approach for delaying leaf senescence and provides insights into the mechanism of melatonin-mediated antisenescence in postharvest Chinese kale.

Current perspectives of lncRNAs in abiotic and biotic stress tolerance in plants

Abiotic/biotic stresses pose a major threat to agriculture and food security by impacting plant growth, productivity and quality. The discovery of extensive transcription of large RNA transcripts that do not code for proteins, termed long non-coding RNAs (lncRNAs) with sizes larger than 200 nucleotides in length, provides an important new perspective on the centrality of RNA in gene regulation. In plants, lncRNAs are widespread and fulfill multiple biological functions in stress response. In this paper, the research advances on the biological function of lncRNA in plant stress response were summarized, like as Natural Antisense Transcripts (NATs), Competing Endogenous RNAs (ceRNAs) and Chromatin Modification etc. And in plants, lncRNAs act as a key regulatory hub of several phytohormone pathways, integrating abscisic acid (ABA), jasmonate (JA), salicylic acid (SA) and redox signaling in response to many abiotic/biotic stresses. Moreover, conserved sequence motifs and structural motifs enriched within stress-responsive lncRNAs may also be responsible for the stress-responsive functions of lncRNAs, it will provide a new focus and strategy for lncRNA research. Taken together, we highlight the unique role of lncRNAs in integrating plant response to adverse environmental conditions with different aspects of plant growth and development. We envisage that an improved understanding of the mechanisms by which lncRNAs regulate plant stress response may further promote the development of unconventional approaches for breeding stress-resistant crops.

Integrated Transcriptome and Proteome Analysis Revealed the Regulatory Mechanism of Hypocotyl Elongation in Pakchoi

Hypocotyl length is a critical determinant for the efficiency of mechanical harvesting in pakchoi production, but the knowledge on the molecular regulation of hypocotyl growth is very limited. Here, we report a spontaneous mutant of pakchoi, lhy7.1, and identified its characteristics. We found that it has an elongated hypocotyl phenotype compared to the wild type caused by the longitudinal growth of hypocotyl cells. Different light quality treatments, transcriptome, and proteomic analyses were performed to reveal the molecular mechanisms of hypocotyl elongation. The data showed that the hypocotyl length of lhy7.1 was significantly longer than that of WT under red, blue, and white lights but there was no significant difference under dark conditions. Furthermore, we used transcriptome and label-free proteome analyses to investigate differences in gene and protein expression levels between lhy7.1 and WT. At the transcript level, 4568 differentially expressed genes (DEGs) were identified, which were mainly enriched in "plant hormone signal transduction", "photosynthesis", "photosynthesis-antenna proteins", and "carbon fixation in photosynthetic organisms" pathways. At the protein level, 1007 differentially expressed proteins (DEPs) were identified and were mainly enriched in photosynthesis-related pathways. The comprehensive transcriptome and proteome analyses revealed a regulatory network of hypocotyl elongation involving plant hormone signal transduction and photosynthesis-related pathways. The findings of this study help elucidate the regulatory mechanisms of hypocotyl elongation in lhy7.1.

BcWRKY22 Activates BcCAT2 to Enhance Catalase (CAT) Activity and Reduce Hydrogen Peroxide (H2O2) Accumulation, Promoting Thermotolerance in Non-Heading Chinese Cabbage (Brassica campestris ssp. chinensis)

WRKY transcription factors (TFs) participate in plant defense mechanisms against biological and abiotic stresses. However, their regulatory role in heat resistance is still unclear in non-heading Chinese cabbage. Here, we identified the WRKY-IIe gene BcWRKY22(BraC09g001080.1), which is activated under high temperatures and plays an active role in regulating thermal stability, through transcriptome analysis. We further discovered that the BcWRKY22 protein is located in the nucleus and demonstrates transactivation activity in both the yeast and plant. Additionally, our studies showed that the transient overexpression of BcWRKY22 in non-heading Chinese cabbage activates the expression of catalase 2 (BcCAT2), enhances CAT enzyme activity, and reduces Hydrogen Peroxide (H2O2) accumulation under heat stress conditions. In addition, compared to its wild-type (WT) counterparts, Arabidopsis thaliana heterologously overexpresses BcWRKY22, improving thermotolerance. When the BcWRKY22 transgenic root was obtained, under heat stress, the accumulation of H2O2 was reduced, while the expression of catalase 2 (BcCAT2) was upregulated, thereby enhancing CAT enzyme activity. Further analysis revealed that BcWRKY22 directly activates the expression of BcCAT2 (BraC08g016240.1) by binding to the W-box element distributed within the promoter region of BcCAT2. Collectively, our findings suggest that BcWRKY22 may serve as a novel regulator of the heat stress response in non-heading Chinese cabbage, actively contributing to the establishment of thermal tolerance by upregulating catalase (CAT) activity and downregulating H2O2 accumulation via BcCAT2 expression.

Comprehensive transcriptome analysis reveals heat-responsive genes in flowering Chinese cabbage (Brassica campestris L. ssp. chinensis) using RNA sequencing

Flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee, 2n=20, AA) is a vegetable species in southern parts of China that faces high temperatures in the summer and winter seasons. While heat stress adversely impacts plant productivity and survival, the underlying molecular and biochemical causes are poorly understood. This study investigated the gene expression profiles of heat-sensitive (HS) '3T-6' and heat-tolerant (HT) 'Youlu-501' varieties of flowering Chinese cabbage in response to heat stress using RNA sequencing. Among the 37,958 genes expressed in leaves, 20,680 were differentially expressed genes (DEGs) at 1, 6, and 12 h, with 1,078 simultaneously expressed at all time points in both varieties. Hierarchical clustering analysis identified three clusters comprising 1,958, 556, and 591 down-regulated, up-regulated, and up- and/or down-regulated DEGs (3205 DEGs; 8.44%), which were significantly enriched in MAPK signaling, plant-pathogen interactions, plant hormone signal transduction, and brassinosteroid biosynthesis pathways and involved in stimulus, stress, growth, reproductive, and defense responses. Transcription factors, including MYB (12), NAC (13), WRKY (11), ERF (31), HSF (17), bHLH (16), and regulatory proteins such as PAL, CYP450, and photosystem II, played an essential role as effectors of homeostasis, kinases/phosphatases, and photosynthesis. Among 3205 DEGs, many previously reported genes underlying heat stress were also identified, e.g., BraWRKY25, BraHSP70, BraHSPB27, BraCYP71A23, BraPYL9, and BraA05g032350.3C. The genome-wide comparison of HS and HT provides a solid foundation for understanding the molecular mechanisms of heat tolerance in flowering Chinese cabbage.

A long noncoding RNA HILinc1 enhances pear thermotolerance by stabilizing PbHILT1 transcripts through complementary base pairing

As global warming intensifies, heat stress has become a major environmental constraint threatening crop production and quality worldwide. Here, we characterize Heat-induced long intergenic noncoding RNA 1 (HILinc1), a cytoplasm-enriched lincRNA that plays a key role in thermotolerance regulation of pear (Pyrus spp.). HILinc1 Target 1 (PbHILT1) which is the target transcript of HILinc1, was stabilized via complementary base pairing to upregulate its expression. PbHILT1 could bind to Heat shock transcription factor A1b (PbHSFA1b) to enhance its transcriptional activity, leading to the upregulation of a major downstream transcriptional regulator, Multiprotein bridging factor 1c (PbMBF1c), during heat response. Transient overexpressing of either HILinc1 or PbHILT1 increases thermotolerance in pear, while transient silencing of HILinc1 or PbHILT1 makes pear plants more heat sensitive. These findings provide evidences for a new regulatory mechanism by which HILinc1 facilitates PbHSFA1b activity and enhances pear thermotolerance through stabilizing PbHILT1 transcripts.

Heat stress in pear cultivar results in upregulation of long non-coding RNA HILinc1, which binds to and stabilizes PbHILT1 mRNA, which codes for a protein that interacts with heat shock factor A1b, improving thermotolerance.

A natural antisense RNA improves chrysanthemum cold tolerance by regulating the transcription factor DgTCP1

Long noncoding RNAs (lncRNAs) are widely involved in the regulation of plant growth and development, but their mechanism of action in response to cold stress in plants remains unclear. Here, we found an lncRNA transcribed from the antisense strand of DgTCP1 (class I Teosinte branched1/Cycloidea/Proliferating [TCP] transcription factor) of chrysanthemum (Chrysanthemum morifolium Ramat.), named DglncTCP1. During the response of chrysanthemum to cold stress, overexpression of DgTCP1 improved the cold tolerance of chrysanthemum, while the DgTCP1 editing line (dgtcp1) showed decreased tolerance to cold stress. Overexpression of DglncTCP1 also increased the cold tolerance of chrysanthemum, while the DglncTCP1 amiRNA lines (DglncTCP1 amiR-18/38) also showed decreased tolerance to cold stress. Additionally, the overexpression of DglncTCP1 upregulated the expression of DgTCP1. This indicated that DglncTCP1 may play a cis-regulatory role in the regulatory process of DgTCP1 in cold tolerance. DglncTCP1 acts as a scaffold to recruit the histone modification protein DgATX (ARABIDOPSIS TRITHORAX from chrysanthemum) to DgTCP1 to enhance H3K4me3 levels, thereby activating DgTCP1 expression. Moreover, DgTCP1 can directly target DgPOD (peroxidase gene from chrysanthemum) to promote its expression and reduce reactive oxygen species accumulation, thereby improving the cold tolerance of chrysanthemum. In conclusion, these results suggest that natural antisense lncRNA plays a key role in improving the cold tolerance of chrysanthemum.

Transcriptomics and Genomics Analysis Uncover the Differentially Expressed Chlorophyll and Carotenoid-Related Genes in Celery

Celery (Apium graveolens L.), a plant from Apiaceae, is one of the most important vegetables and is grown worldwide. Carotenoids can capture light energy and transfer it to chlorophyll, which plays a central role in photosynthesis. Here, by performing transcriptomics and genomics analysis, we identified and conducted a comprehensive analysis of chlorophyll and carotenoid-related genes in celery and six representative species. Significantly, different contents and gene expression patterns were found among three celery varieties. In total, 237 and 290 chlorophyll and carotenoid-related genes were identified in seven species. No notable gene expansion of chlorophyll biosynthesis was detected in examined species. However, the gene encoding ζ-carotene desaturase (ZDS) enzyme in carotenoid was expanded in celery. Comparative genomics and RNA-seq analyses revealed 16 and 5 key genes, respectively, regulating chlorophyll and carotenoid. An intriguing finding is that chlorophyll and carotenoid-related genes were coordinately regulated by transcriptional factors, which could be distinctively classified into positive- and negative-regulation groups. Six CONSTANS (CO)-like transcription factors co-regulated chlorophyll and carotenoid-related genes were identified in celery. In conclusion, this study provides new insights into the regulation of chlorophyll and carotenoid by transcription factors.

The Functional Interplay between Ethylene, Hydrogen Sulfide, and Sulfur in Plant Heat Stress Tolerance

Plants encounter several abiotic stresses, among which heat stress is gaining paramount attention because of the changing climatic conditions. Severe heat stress conspicuously reduces crop productivity through changes in metabolic processes and in growth and development. Ethylene and hydrogen sulfide (H2S) are signaling molecules involved in defense against heat stress through modulation of biomolecule synthesis, the antioxidant system, and post-translational modifications. Other compounds containing the essential mineral nutrient sulfur (S) also play pivotal roles in these defense mechanisms. As biosynthesis of ethylene and H2S is connected to the S-assimilation pathway, it is logical to consider the existence of a functional interplay between ethylene, H2S, and S in relation to heat stress tolerance. The present review focuses on the crosstalk between ethylene, H2S, and S to highlight their joint involvement in heat stress tolerance.

Global identification and integrated analysis of heat-responsive long non-coding RNAs in contrasting rice cultivars

Identified 2743 rice lncRNAs LncRNAs in response to heat stress Function prediction of HRLs Network among HRLs, genes and miRNAs co-localization of HRLs with QTLs Significant motifs in HRL sequences Long non-coding RNAs (lncRNAs) play vital roles in plant responses to environmental challenges. A better understanding of the gene regulation mediated by lncRNAs and their systematic identification would provide great benefits for modern agriculture. In this study, we performed strand-specific RNA sequencing for two rice varieties, heat-tolerant ZS97B and heat-susceptible SYD2 under heat stress. In total, 2743 putative lncRNAs were identified, and their expression profiles in response to heat treatments were established. We identified 231 differentially expressed lncRNAs (DELs) under heat stress, including 31 DELs common to both varieties and 103 and 97 specific to ZS97B and SYD2, respectively, all defined as heat-responsive lncRNAs (HRLs). The target-coding genes of HRLs were predicted, and GO and KEGG annotations of HRL targets revealed functions in which HRLs might be involved. The interaction network between HRLs, target genes and relevant miRNAs was constructed. The HRLs and their targets were compared with publicly available QTLs for rice seedling growth under heat stimulus. Ten HRLs and twelve target genes were linked with five heat stress-relevant QTLs. Sequence analysis revealed several motifs significantly enriched within the 231 HRL sequences. Our findings provide a valuable resource for further characterization of lncRNAs in terms of heat response and plant heat tolerance improvement.

Large-scale analyses of heat shock transcription factors and database construction based on whole-genome genes in horticultural and representative plants

Heat shock transcription factor (Hsf) plays a critical role in regulating heat resistance. Here, 2950 Hsf family genes were identified from 111 horticultural and representative plants. More Hsf genes were detected in higher plants than lower plants. Based on all Hsf genes, we constructed a phylogenetic tree, which indicated that Hsf genes of each branch evolved independently after species differentiation. Furthermore, we uncovered the evolutionary trajectories of Hsf genes by motif analysis. There were only 6 motifs (M1 to M6) in lower plants, and then 4 novel motifs (M7-M10) appeared in higher plants. However, the motifs of some Hsf genes were lost in higher plant, indicating that Hsf genes have undergone sequence variation during the evolution. The number of Hsf gene loss was more than duplication after whole-genome duplication in higher plants. The heat response network was constructed using 24 Hsf genes, 2421 downstream, and 222 upstream genes of Arabidopsis. Further enrichment analysis revealed that Hsf genes and other transcription factors interacted with each other to response heat resistance. The global expression maps were illustrated for Hsf genes under various abiotic, biotic stresses, and several developmental stages in Arabidopsis. The syntenic and phylogenetic analyses were conducted using Hsf genes of Arabidopsis and Pan-genome of 18 Brassica rapa accessions. We also performed the expression pattern analysis of Hsf and six Hsp family genes using expression values from different tissues and heat treatments in B. rapa. The interaction network between Hsf and Hsp gene families was constructed in B. rapa, and several core genes were detected in the network. Finally, we constructed a Hsf database (http://hsfdb.bio2db.com) for researchers to retrieve Hsf gene family information. Therefore, our study will provide rich resources for the evolution and functional study of Hsf genes.

Bidirectional lncRNA Transfer between Cuscuta Parasites and Their Host Plant

Dodder species (Cuscuta spp.) are holoparasites that have extensive material exchange with their host plants through vascular connections. Recent studies on cross-species transfer have provided breakthrough insights, but little is known about the interaction mechanisms of the inter-plant mobile substances in parasitic systems. We sequenced the transcriptomes of dodder growing on soybean hosts to characterize the long non-coding RNA (lncRNA) transfer between the two species, and found that lncRNAs can move in high numbers (365 dodder lncRNAs and 14 soybean lncRNAs) in a bidirectional manner. Reverse transcription-polymerase chain reaction further confirmed that individual lncRNAs were trafficked in the dodder-soybean parasitic system. To reveal the potential functions of mobile transcripts, the Gene Ontology terms of mobile lncRNA target genes were predicted, and mobile dodder target genes were found to be mainly enriched in "metabolic process", "catalytic activity", "signaling", and "response to stimulus" categories, whereas mobile soybean target genes were enriched in organelle-related categories, indicating that specific mobile lncRNAs may be important in regulating dodder parasitism. Our findings reveal that lncRNAs are transferred between dodder and its host soybean plants, which may act as critical regulators to coordinate the host-dodder interaction at the whole parasitic level.

Genetic and Physiological Responses to Heat Stress in Brassica napus

Given the current rise in global temperatures, heat stress has become a major abiotic challenge affecting the growth and development of various crops and reducing their productivity. Brassica napus, the second largest source of vegetable oil worldwide, experiences a drastic reduction in seed yield and quality in response to heat. This review outlines the latest research that explores the genetic and physiological impact of heat stress on different developmental stages of B. napus with a special attention to the reproductive stages of floral progression, organogenesis, and post flowering. Several studies have shown that extreme temperature fluctuations during these crucial periods have detrimental effects on the plant and often leading to impaired growth and reduced seed production. The underlying mechanisms of heat stress adaptations and associated key regulatory genes are discussed. Furthermore, an overview and the implications of the polyploidy nature of B. napus and the regulatory role of alternative splicing in forming a priming-induced heat-stress memory are presented. New insights into the dynamics of epigenetic modifications during heat stress are discussed. Interestingly, while such studies are scarce in B. napus, opposite trends in expression of key genetic and epigenetic components have been identified in different species and in cultivars within the same species under various abiotic stresses, suggesting a complex role of these genes and their regulation in heat stress tolerance mechanisms. Additionally, omics-based studies are discussed with emphasis on the transcriptome, proteome and metabolome of B. napus, to gain a systems level understanding of how heat stress alters its yield and quality traits. The combination of omics approaches has revealed crucial interactions and regulatory networks taking part in the complex machinery of heat stress tolerance. We identify key knowledge gaps regarding the impact of heat stress on B. napus during its yield determining reproductive stages, where in-depth analysis of this subject is still needed. A deeper knowledge of heat stress response components and mechanisms in tissue specific models would serve as a stepping-stone to gaining insights into the regulation of thermotolerance that takes place in this important crop species and support future breeding of heat tolerant crops.

The Dynamism of Transposon Methylation for Plant Development and Stress Adaptation

Plant development processes are regulated by epigenetic alterations that shape nuclear structure, gene expression, and phenotypic plasticity; these alterations can provide the plant with protection from environmental stresses. During plant growth and development, these processes play a significant role in regulating gene expression to remodel chromatin structure. These epigenetic alterations are mainly regulated by transposable elements (TEs) whose abundance in plant genomes results in their interaction with genomes. Thus, TEs are the main source of epigenetic changes and form a substantial part of the plant genome. Furthermore, TEs can be activated under stress conditions, and activated elements cause mutagenic effects and substantial genetic variability. This introduces novel gene functions and structural variation in the insertion sites and primarily contributes to epigenetic modifications. Altogether, these modifications indirectly or directly provide the ability to withstand environmental stresses. In recent years, many studies have shown that TE methylation plays a major role in the evolution of the plant genome through epigenetic process that regulate gene imprinting, thereby upholding genome stability. The induced genetic rearrangements and insertions of mobile genetic elements in regions of active euchromatin contribute to genome alteration, leading to genomic stress. These TE-mediated epigenetic modifications lead to phenotypic diversity, genetic variation, and environmental stress tolerance. Thus, TE methylation is essential for plant evolution and stress adaptation, and TEs hold a relevant military position in the plant genome. High-throughput techniques have greatly advanced the understanding of TE-mediated gene expression and its associations with genome methylation and suggest that controlled mobilization of TEs could be used for crop breeding. However, development application in this area has been limited, and an integrated view of TE function and subsequent processes is lacking. In this review, we explore the enormous diversity and likely functions of the TE repertoire in adaptive evolution and discuss some recent examples of how TEs impact gene expression in plant development and stress adaptation.

Transcriptional profiling reveals changes in gene regulation and signaling transduction pathways during temperature stress in wucai (Brassica campestris L.)

Background

Wucai (Brassica campestris L. ssp. chinensis var. rosularis Tsen) is a cold-tolerant plant that is vulnerable to high temperature. This study explored the response mechanism of wucai to low temperature. In this study, wucai seedlings were treated with different temperatures, including low temperature (LT), high temperature (HT), and a control.

Results

According to transcriptomics analysis, the number of differentially expressed genes (DEGs) in HT and LT was 10,702 and 7267, respectively, compared with the control. The key genes associated with the physiological response of wucai to the treatments were analyzed. The Kyoto Encyclopedia of Genes and Genomes and Gene Ontology annotations indicated the importance of the photosynthesis and photosynthetic-antenna protein pathways. We found that a high-temperature environment greatly inhibited the expression of important genes in the photosynthetic pathway (BrLhc superfamily members, PsaD, PsaE, PsaD, PsaD, PsbO, PsbP, PsbQ, PsbR, PsbS, PsbW, PsbY, Psb27, and Psb28), whereas low temperature resulted in the expression of certain key genes (BrLhc superfamily members, Psa F, Psa H, Psb S, Psb H, Psb 28). In addition, the wucai seedlings exhibited better photosynthetic performance under low-temperature conditions than high-temperature conditions.

Conclusions

Based on the above results, we speculate that upon exposure to low temperature, the plants developed higher cold tolerance by upregulating the expression of genes related to photosynthesis. Conversely, high-temperature stress inhibited the expression of pivotal genes and weakened the self-regulating ability of the plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07981-9.

Transcriptome and Small RNA Combined Sequencing Analysis of Cold Tolerance in Non-heading Chinese Cabbage

Background

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis) is an important leaf vegetable grown worldwide. However, there has currently been not enough transcriptome and small RNA combined sequencing analysis of cold tolerance, which hinders further functional genomics research.

Results

In this study, 63.43 Gb of clean data was obtained from the transcriptome analysis. The clean data of each sample reached 6.99 Gb, and the basic percentage of Q30 was 93.68% and above. The clean reads of each sample were sequence aligned with the designated reference genome (Brassica rapa, IVFCAASv1), and the efficiency of the alignment varied from 81.54 to 87.24%. According to the comparison results, 1,860 new genes were discovered in Pak-choi, of which 1,613 were functionally annotated. Among them, 13 common differentially expressed genes were detected in all materials, including seven upregulated and six downregulated. At the same time, we used quantitative real-time PCR to confirm the changes of these gene expression levels. In addition, we sequenced miRNA of the same material. Our findings revealed a total of 34,182,333 small RNA reads, 88,604,604 kinds of small RNAs, among which the most common size was 24 nt. In all materials, the number of common differential miRNAs is eight. According to the corresponding relationship between miRNA and its target genes, we carried out Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the set of target genes on each group of differentially expressed miRNAs. Through the analysis, it is found that the distributions of candidate target genes in different materials are different. We not only used transcriptome sequencing and small RNA sequencing but also used experiments to prove the expression levels of differentially expressed genes that were obtained by sequencing. Sequencing combined with experiments proved the mechanism of some differential gene expression levels after low-temperature treatment.

Conclusion

In all, this study provides a resource for genetic and genomic research under abiotic stress in Pak-choi.

The transcriptomic response to heat stress of a jujube (Ziziphus jujuba Mill.) cultivar is featured with changed expression of long noncoding RNAs

Long non-coding RNA (lncRNA) of plant species undergoes dynamic regulation and acts in developmental and stress regulation. Presently, there is little information regarding the identification of lncRNAs in jujube (Ziziphus jujuba Mill.), and it is uncertain whether the lncRNAs could respond to heat stress (HS) or not. In our previous study, a cultivar (Hqing1-HR) of Z. jujuba were treated by HS (45°C) for 0, 1, 3, 5 and 7 days, and it was found that HS globally changed the gene expression by RNA sequencing (RNA-seq) experiments and informatics analyses. In the current study, 8260 lncRNAs were identified successfully from the previous RNA-seq data, and it indicated that lncRNAs expression was also altered globally, suggesting that the lncRNAs might play vital roles in response to HS. Furthermore, bioinformatics analyses of potential target mRNAs of lncRNAs with cis-acting mechanism were performed, and it showed that multiple differentially expressed (DE) mRNAs co-located with DElncRNAs were highly enriched in pathways associated with response to stress and regulation of metabolic process. Taken together, these findings not only provide a comprehensive identification of lncRNAs but also useful clues for molecular mechanism response to HS in jujube.

Genome-wide analysis of long noncoding RNAs, 24-nt siRNAs, DNA methylation and H3K27me3 marks in Brassica rapa

Long noncoding RNAs (lncRNAs) are RNA fragments that generally do not code for a protein but are involved in epigenetic gene regulation. In this study, lncRNAs of Brassica rapa were classified into long intergenic noncoding RNAs, natural antisense RNAs, and intronic noncoding RNAs and their expression analyzed in relation to genome-wide 24-nt small interfering RNAs (siRNAs), DNA methylation, and histone H3 lysine 27 trimethylation marks (H3K27me3). More than 65% of the lncRNAs analyzed consisted of one exon, and more than 55% overlapped with inverted repeat regions (IRRs). Overlap of lncRNAs with IRRs or genomic regions encoding for 24-nt siRNAs resulted in increased DNA methylation levels when both were present. LncRNA did not overlap greatly with H3K27me3 marks, but the expression level of intronic noncoding RNAs that did coincide with H3K27me3 marks was higher than without H3K27me3 marks. The Brassica genus comprises important vegetables and oil seed crops grown across the world. B. rapa is a diploid (AA genome) thought to be one of the ancestral species of both B. juncea (AABB genome) and B. napus (AACC) through genome merging (allotetrapolyploidization). Complex genome restructuring and epigenetic alterations are thought to be involved in these allotetrapolyploidization events. Comparison of lncRNAs between B. rapa and B. nigra, B. oleracea, B. juncea, and B. napus showed the highest conservation with B. oleracea. This study presents a comprehensive analysis of the epigenome structure of B. rapa at multi-epigenetic levels (siRNAs, DNA methylation, H3K27me3, and lncRNAs) and identified a suite of candidate lncRNAs that may be epigenetically regulated in the Brassica genus.

Comparative analysis of long noncoding RNAs in angiosperms and characterization of long noncoding RNAs in response to heat stress in Chinese cabbage

Long noncoding RNAs (lncRNAs) are widely present in different species and play critical roles in response to abiotic stresses. However, the functions of lncRNAs in Chinese cabbage under heat stress remain unknown. Here, we first conducted a global comparative analysis of 247,242 lncRNAs among 37 species. The results indicated that lncRNAs were poorly conserved among different species, and only 960 lncRNAs were homologous to 524 miRNA precursors. We then carried out lncRNA sequencing for a genome-wide analysis of lncRNAs and their target genes in Chinese cabbage at different stages of heat treatment. In total, 18,253 lncRNAs were identified, of which 1229 differentially expressed (DE) lncRNAs were characterized as being heat-responsive. The ceRNA network revealed that 38 lncRNAs, 16 miRNAs, and 167 mRNAs were involved in the heat response in Chinese cabbage. Combined analysis of the cis- and trans-regulated genes indicated that the targets of DE lncRNAs were significantly enriched in the "protein processing in endoplasmic reticulum" and "plant hormone signal transduction" pathways. Furthermore, the majority of HSP and PYL genes involved in these two pathways exhibited similar expression patterns and responded to heat stress rapidly. Based on the networks of DE lncRNA-mRNAs, 29 and 22 lncRNAs were found to interact with HSP and PYL genes, respectively. Finally, the expression of several critical lncRNAs and their targets was verified by qRT-PCR. Overall, we conducted a comparative analysis of lncRNAs among 37 species and performed a comprehensive analysis of lncRNAs in Chinese cabbage. Our findings expand the knowledge of lncRNAs involved in the heat stress response in Chinese cabbage, and the identified lncRNAs provide an abundance of resources for future comparative and functional studies.

Identification of long noncoding RNAs involved in resistance to downy mildew in Chinese cabbage

Brassica downy mildew, a severe disease caused by Hyaloperonospora brassicae, can cause enormous economic losses in Chinese cabbage (Brassica rapa L. ssp. pekinensis) production. Although some research has been reported recently concerning the underlying resistance to this disease, no studies have identified or characterized long noncoding RNAs involved in this defense response. In this study, using high-throughput RNA sequencing, we analyzed the disease-responding mRNAs and long noncoding RNAs in two resistant lines (T12-19 and 12-85) and one susceptible line (91-112). Clustering and Gene Ontology analysis of differentially expressed genes (DEGs) showed that more DEGs were involved in the defense response in the two resistant lines than in the susceptible line. Different expression patterns and proposed functions of differentially expressed long noncoding RNAs among T12-19, 12-85, and 91-112 indicated that each has a distinct disease response mechanism. There were significantly more cis- and trans-functional long noncoding RNAs in the resistant lines than in the susceptible line, and the genes regulated by these RNAs mostly participated in the disease defense response. Furthermore, we identified a candidate resistance-related long noncoding RNA, MSTRG.19915, which is a long noncoding natural antisense transcript of a MAPK gene, BrMAPK15. Via an agroinfiltration-mediated transient overexpression system and virus-induced gene silencing technology, BrMAPK15 was indicated to have a greater ability to defend against pathogens. MSTRG.19915-silenced seedlings showed enhanced resistance to downy mildew, probably because of the upregulated expression of BrMAPK15. This research identified and characterized long noncoding RNAs involved in resistance to downy mildew, laying a foundation for future in-depth studies of disease resistance mechanisms in Chinese cabbage.

Genome-Wide Identification and Comparative Analysis of ARF Family Genes in Three Apiaceae Species

The family Apiaceae includes many important vegetables and medicinal plants. Auxin response factors (ARFs) play critical roles in regulating plant growth and development. Here, we performed a comprehensive analysis of the ARF gene family in three Apiaceae species, celery, coriander, and carrot, and compared the results with the ARF gene family of lettuce, Arabidopsis, and grape. We identified 156 ARF genes in all six species and 89 genes in the three Apiaceae species, including 28, 34, and 27 in celery, coriander, and carrot, respectively. The paralogous gene number in coriander was far greater than that in carrot and celery. Our analysis revealed that ARF genes of the three Apiaceae species in 34 branches of the phylogenetic tree underwent significant positive selection. Additionally, our findings indicated that whole-genome duplication played an important role in ARF gene family expansion. Coriander contained a greater number of ARF genes than celery and carrot because of more gene duplications and less gene losses. We also analyzed the expression of ARF genes in three tissues by RNA-seq and verified the results by quantitative real-time PCR. Furthermore, we found that several paralogous genes exhibited divergent expression patterns. Overall, this study provides a valuable resource for exploring how ARF family genes regulate plant growth and development in other plants. Since this is the first report of the ARF gene family in Apiaceae, our results will serve as a guide for comparative and functional analyses of ARF and other gene families in Apiaceae.

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Screening and functional characterization of candidate resistance genes to powdery mildew from Dasypyrum villosum#4 in a wheat line Pm97033

KEY MESSAGE

Three genes designated DvLox, Pm21#4, and Pm21#4-H identified in a wheat-Dasypyrum villosum#4 T6V#4S·6DL translocation line Pm97033 conferred wheat for powdery mildew resistance. Powdery mildew (PM) caused by Blumeria graminis f. sp. tritici (Bgt) is one of the most devastating diseases in wheat. To date, only a few genes conferring resistance to wheat PM are cloned. Dasypyrum villosum is a wild relative of wheat, which provides Pm21 conferring wheat immunity to PM. In this study, we obtained many differentially expressed genes (DEGs) from a wheat-D. villosum#4 T6V#4S·6DL translocation line Pm97033 using RNA-sequencing. Among them, 7 DEGs associated with pathogen resistance were up-regulated in front of Bgt infection. Virus-induced gene silencing and transformation assays demonstrated that two of them, DvLox and Pm21#4 encoding a lipoxygenase and a encoding coiled-coil/nucleotide-binding site/leucine-rich repeat resistance protein, conferred wheat PM resistance. The transgenic wheat plants expressing DvLox enhanced PM resistance, and the transgenic wheat plants expressing Pm21#4 showed PM immunity. The Pm21#4-silenced Pm97033 plants by the cluster regularly interspaced short palindromic repeats-associated endonuclease (CRISPR/Cas9) system were susceptible to PM. Thus, Pm21#4 is a key gene contributing PM immune resistance in Pm97033. Constitutively expression of Pm21#4-H, which is silenced in Pm97033 and D. villosum#4, endowed a PM-susceptible wheat variety Fielder with PM immune resistance.

Long non-coding RNAs: emerging players regulating plant abiotic stress response and adaptation

BACKGROUND

The immobile nature of plants means that they can be frequently confronted by various biotic and abiotic stresses during their lifecycle. Among the various abiotic stresses, water stress, temperature extremities, salinity, and heavy metal toxicity are the major abiotic stresses challenging overall plant growth. Plants have evolved complex molecular mechanisms to adapt under the given abiotic stresses. Long non-coding RNAs (lncRNAs)-a diverse class of RNAs that contain > 200 nucleotides(nt)-play an essential role in plant adaptation to various abiotic stresses.

RESULTS

LncRNAs play a significant role as 'biological regulators' for various developmental processes and biotic and abiotic stress responses in animals and plants at the transcription, post-transcription, and epigenetic level, targeting various stress-responsive mRNAs, regulatory gene(s) encoding transcription factors, and numerous microRNAs (miRNAs) that regulate the expression of different genes. However, the mechanistic role of lncRNAs at the molecular level, and possible target gene(s) contributing to plant abiotic stress response and adaptation, remain largely unknown. Here, we review various types of lncRNAs found in different plant species, with a focus on understanding the complex molecular mechanisms that contribute to abiotic stress tolerance in plants. We start by discussing the biogenesis, type and function, phylogenetic relationships, and sequence conservation of lncRNAs. Next, we review the role of lncRNAs controlling various abiotic stresses, including drought, heat, cold, heavy metal toxicity, and nutrient deficiency, with relevant examples from various plant species. Lastly, we briefly discuss the various lncRNA databases and the role of bioinformatics for predicting the structural and functional annotation of novel lncRNAs.

CONCLUSIONS

Understanding the intricate molecular mechanisms of stress-responsive lncRNAs is in its infancy. The availability of a comprehensive atlas of lncRNAs across whole genomes in crop plants, coupled with a comprehensive understanding of the complex molecular mechanisms that regulate various abiotic stress responses, will enable us to use lncRNAs as potential biomarkers for tailoring abiotic stress-tolerant plants in the future.

Comparative transcriptome analysis reveals heat stress-responsive genes and their signalling pathways in lilies (Lilium longiflorum vs. Lilium distichum)

The lily, a famous bulbous flower, is seriously affected by high temperatures, which affect their growth and production. To date, the signalling pathways and the molecular mechanisms related to heat response in Lilium have not been elucidated. In this study, a comparative transcriptome analysis was performed in an important thermo-tolerant flower, L. longiflorum, and a thermo-sensitive flower, L. distichum. Lily seedlings were first exposed to heat stress at 42°C for different lengths of time, and the optimal time-points (2 h and 24 h) were selected for RNA sequencing (RNA-seq). Approximately 66.51, 66.21, and 65.36 Mb clean reads were identified from three libraries of L. longiflorum (LL_CK, LL_T2h and LL_T24h, respectively) and 66.18, 66.03, and 65.16 Mb clean reads were obtained from three libraries of L. distichum (LD_CK, LD_T2h and LD_T24h, respectively) after rRNA removing. A total of 34,301 unigenes showed similarity to known proteins in the database NCBI non-redundant protein (NR), Swiss-Prot proteins, InterPro proteins, Clusters of Orthologous Groups (COG) and Kyoto Encyclopedia of Genes and Genomes (KEGG). In addition, 1,621 genes were differentially expressed in the overlapping libraries between LL_DEGs and LD_DEGs; of these genes, 352 DEGs were obviously upregulated in L. longiflorum and downregulated in L. distichum during heat stress, including 4-coumarate, CoA ligase (4CL), caffeoyl-CoA O-methyltransferase (CCoAOMT), peroxidase, pathogenesis-related protein 10 family genes (PR10s), 14-3-3 protein, leucine-rich repeat receptor-like protein kinase, and glycine-rich cell wall structural protein-like. These genes were mainly involved in metabolic pathways, phenylpropanoid biosynthesis, plant-pathogen interactions, plant hormone signal transduction, and kinase signalling pathways. Quantitative RT-PCR was performed to validate the expression profiling of these DEGs in RNA-seq data. Taken together, the results obtained in the present study provide a comprehensive sequence resource for the discovery of heat-resistance genes and reveal potential key components that are responsive to heat stress in lilies, which may help to elucidate the heat signal transcription networks and facilitate heat-resistance breeding in lily.

Third-generation sequencing found LncRNA associated with heat shock protein response to heat stress in Populus qiongdaoensis seedlings

Background

As air temperatures increase globally, more and more plants are exposed to heat-stress conditions. Although many studies have explored regulation networks in plants with the aim of improving their heat-stress tolerance, only few have revealed them in trees. Here, individuals of Populus qiongdaoensis seedlings, which grows naturally in tropical areas, exposed to heat at 40 °C and the non-coding regulation networks were explored using the PacBio RSII and the Illumina sequencing platform.

Results

In total, we obtained 88,161 full-length transcripts representing 39,343 genes using 5,498,988 long reads and 350,026,252 clean reads, and also 216 microRNAs (miRNAs) via 95,794,107 reads. We then identified 928 putative long non-coding RNAs (lncRNAs), consisting of 828 sense lncRNAs (89.22%), 34 long intergenic non-coding RNAs (3.66%), 16 antisense (1.72%), and 50 sense intronic lncRNAs (5.39%). Under the dual criteria of |log₂fold-change| ≥ 1 and P-value < 0.05, 1690 genes, 25 lncRNAs, and 15 miRNAs were found differentially expressed under the heat stress treatment. Furthermore, 563 and 595 mRNAs were detected as target genes of 14 differently expressed miRNAs and 26 differentially expressed lncRNAs. Functional annotation analysis of these target genes demonstrated they were related to cell membrane stability, plant hormone signal transduction, antioxidation, and aldarate metabolism. Lastly, we uncovered a key interaction network of lncRNAs, miRNAs and mRNAs that consisted of miR1444d, miR482a.1, miR530a, lncHSP18.2, HSP18.1, and HSP18.2. Expression level analysis showed that miRNAs in the network were up-regulated, while mRNAs and lncRNA were down-regulated, and also found that lncHSP18.2 may cis-regulate HSP18.2.

Conclusions

Functional enrichment analysis of target genes of miRNAs and lncRNAs indicated that miRNAs and lncRNAs play an important role in the response to heat stress P. qiongdaoensis. Lastly, by investigating the miRNA–lncRNA–mRNA network of this species, we revealed that miRNAs may negatively regulate both lncRNAs and mRNAs in tree responses to heat stress, and found that lncHSP18.2 may cis-regulate HSP18.2.

The Arabidopsis lncRNA ASCO modulates the transcriptome through interaction with splicing factors

Alternative splicing (AS) is a major source of transcriptome diversity. Long noncoding RNAs (lncRNAs) have emerged as regulators of AS through different molecular mechanisms. In Arabidopsis thaliana, the AS regulators NSRs interact with the ALTERNATIVE SPLICING COMPETITOR (ASCO) lncRNA. Here, we analyze the effect of the knock-down and overexpression of ASCO at the genome-wide level and find a large number of deregulated and differentially spliced genes related to flagellin responses and biotic stress. In agreement, ASCO-silenced plants are more sensitive to flagellin. However, only a minor subset of deregulated genes overlaps with the AS defects of the nsra/b double mutant, suggesting an alternative way of action for ASCO. Using biotin-labeled oligonucleotides for RNA-mediated ribonucleoprotein purification, we show that ASCO binds to the highly conserved spliceosome component PRP8a. ASCO overaccumulation impairs the recognition of specific flagellin-related transcripts by PRP8a. We further show that ASCO also binds to another spliceosome component, SmD1b, indicating that it interacts with multiple splicing factors. Hence, lncRNAs may integrate a dynamic network including spliceosome core proteins, to modulate transcriptome reprogramming in eukaryotes.

Systematic identification and analysis of heat-stress-responsive lncRNAs, circRNAs and miRNAs with associated co-expression and ceRNA networks in cucumber (Cucumis sativus L.)

Researchers have shown that long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) act as competitive endogenous RNAs (ceRNAs) and are mutually regulated by competition for binding to common microRNA response elements (MREs). However, a comprehensive identification and analysis of lncRNAs and circRNAs as ceRNAs have not yet been completed in cucumber (Cucumis sativus L.) exposed to high-temperature stress. In our study, 32 663 coding transcripts, 2085 lncRNAs, 2477 circRNAs and 348 differentially expressed miRNAs were identified using RNA sequencing. In addition, six heat-stress-responsive miRNAs (five known and one novel miRNAs) and eight lncRNAs were selected for qPCR to confirm their expression profiles. By analyzing the cis effects of lncRNAs, we constructed a lncRNA-mRNA co-expression network. Based on the results, the corresponding lncRNAs play a regulatory role in the stress response in cucumber plants. In our study, the PatMatch software was used to predict the potential function of lncRNAs and circRNAs as ceRNAs. A total of 18 lncRNAs and seven circRNAs were predicted to bind to 114 differentially expressed miRNAs and compete with 359 mRNAs for miRNA binding sites. These mRNAs are predicted to be involved in various pathways, such as plant hormone signal transduction, plant-pathogen interaction and glutathione metabolism. Among them, TCONS_00031790, TCONS_00014332, TCONS_00014717, TCONS_00005674, novel_circ_001543 and novel_circ_000876 may interact with miR9748 by plant hormone signal transduction pathways in response to high-temperature stress. Moreover, indole-3-acetic acid (IAA) and 1-aminocyclopropane-l-carboxylic acid (ACC) levels decreased in the high-temperature treatment group, indicating that IAA and ethylene signaling might be involved in response to high-temperature stress. In this study, we conducted a full transcriptomic analysis in response to high-temperature stress in cucumber and, for the first time, integrated the potential ceRNA functions of lncRNAs/circRNAs. The results provide a basis for studying the potential functions of lncRNAs/circRNAs in response to high-temperature stress.

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.

Genome-wide analysis of long non-coding RNAs (lncRNAs) in two contrasting rapeseed (Brassica napus L.) genotypes subjected to drought stress and re-watering

BACKGROUND

Drought stress is a major abiotic factor that affects rapeseed (Brassica napus L.) productivity. Though previous studies indicated that long non-coding RNAs (lncRNAs) play a key role in response to drought stress, a scheme for genome-wide identification and characterization of lncRNAs' response to drought stress is still lacking, especially in the case of B. napus. In order to further understand the molecular mechanism of the response of B. napus to drought stress, we compared changes in the transcriptome between Q2 (a drought-tolerant genotype) and Qinyou8 (a drought-sensitive genotype) responding drought stress and rehydration treatment at the seedling stage.

RESULTS

A total of 5546 down-regulated and 6997 up-regulated mRNAs were detected in Q2 compared with 7824 and 10,251 in Qinyou8, respectively; 369 down-regulated and 108 up- regulated lncRNAs were detected in Q2 compared with 449 and 257 in Qinyou8, respectively. LncRNA-mRNA interaction network analysis indicated that the co-expression network of Q2 was composed of 145 network nodes and 5175 connections, while the co-expression network of Qinyou8 was composed of 305 network nodes and 22,327 connections. We further identified 34 transcription factors (TFs) corresponding to 126 differentially expressed lncRNAs in Q2, and 45 TFs corresponding to 359 differentially expressed lncRNAs in Qinyou8. Differential expression analysis of lncRNAs indicated that up- and down-regulated mRNAs co-expressed with lncRNAs participated in different metabolic pathways and were involved in different regulatory mechanisms in the two genotypes. Notably, some lncRNAs were co-expressed with BnaC07g44670D, which are associated with plant hormone signal transduction. Additionally, some mRNAs co-located with XLOC_052298, XLOC_094954 and XLOC_012868 were mainly categorized as signal transport and defense/stress response.

CONCLUSIONS

The results of this study increased our understanding of expression characterization of rapeseed lncRNAs in response to drought stress and re-watering, which would be useful to provide a reference for the further study of the function and action mechanisms of lncRNAs under drought stress and re-watering.

Global identification and analysis revealed differentially expressed lncRNAs associated with meiosis and low fertility in autotetraploid rice

BACKGROUND

Autotetraploid rice is a useful germplasm for polyploid rice breeding. Our previous research showed that non-coding RNAs might be associated with low fertility in autotetraploid rice. However, little information is available on long non-coding RNAs (lncRNAs) involved in the low fertility of autotetraploid rice. In the present study, RNA-seq was employed to detect the differentially expressed meiosis-related lncRNAs in autotetraploid rice, and gene overexpression and knock out experiments were used to validate the potential function of candidate lncRNA.

RESULTS

A total of 444 differentially expressed lncRNAs (DEL) were detected during anther and ovary meiosis in autotetraploid rice. Of these, 328 DEL were associated with the transposable elements, which displayed low expression levels during meiosis in autotetraploid rice. We used rapid amplification of cDNA ends (RACE) assay to validate 10 DEL and found that the lncRNAs were not assembly artifacts, and six of them were conserved in tetraploid rice. Moreover, 237 and 20 lncRNAs were associated with pollen mother cell (PMC) and embryo sac mother cell (EMC) meiosis in autotetraploid rice, respectively. The differential expressions of some meiosis-related targets and its DEL regulator, including MEL1 regulated by TCONS_00068868, LOC_Os12g41350 (meiotic asynaptic mutant 1) by TCONS_00057811 in PMC, and LOC_Os12g39420 by TCONS_00144592 in EMC, were confirmed by qRT-PCR. TCONS_00057811, TCONS_00055980 and TCONS_00130461 showed anther specific expression patterns and were found to be highly expressed during meiosis. CRISPR/Cas9 editing of lncRNA57811 displayed similar morphology compared to wild type. The overexpression of lncRNA57811 resulted in low pollen fertility (29.70%) and seed setting (33%) in rice.

CONCLUSION

The differential expression levels of lncRNAs, associated with transposable elements and meiosis-regulated targets, might be endogenous noncoding regulators of pollen/embryo sac development that cause low fertility in autotetraploid rice. The results enhance our understanding about rice lncRNAs, and facilitate functional research in autotetraploid rice.

Identification, functional prediction, and key lncRNA verification of cold stress-related lncRNAs in rats liver

Cold stimulation reduces the quality of animal products and increases animal mortality, causing huge losses to the livestock industry in cold regions. Long non-coding RNAs (lncRNAs) take part in many biological processes through transcriptional regulation, intracellular material transport, and chromosome remodeling. Although cold stress-related lncRNAs have been reported in plants, no research is available on the characteristic and functional analysis of lncRNAs after cold stress in rats. Here, we built a cold stress animal model firstly. Six SPF male Wistar rats were randomly divided to the acute cold stress group (4 °C, 12 h) and the normal group (24 °C, 12 h). lncRNA libraries were constructed by high-throughput sequencing (HTS) using rat livers. 2,120 new lncRNAs and 273 differentially expressed (DE) lncRNAs were identified in low temperature environments. The target genes of DElncRNA were predicted by cis and trans, and then functional and pathway analysis were performed to them. GO and KEGG analysis revealed that lncRNA targets were mainly participated in the regulation of nucleic acid binding, cold stimulation reaction, metabolic process, immune system processes, PI3K-Akt signaling pathway and pathways in cancer. Next, a interaction network between lncRNA and its targets was constructed. To further reveal the mechanism of cold stress, DElncRNA and DEmRNA were extracted to reconstruct a co-expression sub-network. We found the key lncRNA MSTRG.80946.2 in sub-network. Functional analysis of key lncRNA targets showed that targets were significantly enriched in fatty acid metabolism, the PI3K-Akt signaling pathway and pathways in cancer under cold stress. qRT-PCR confirmed the sequencing results. Finally, hub lncRNA MSTRG.80946.2 was characterized, and verified its relationship with related mRNAs by antisense oligonucleotide (ASO) interference and qRT-PCR. Results confirmed the accuracy of our analysis. To sum up, our work was the first to perform detailed characterization and functional analysis of cold stress-related lncRNAs in rats liver. lncRNAs played crucial roles in energy metabolism, growth and development, immunity and reproductive performance in cold stressed rats. The MSTRG.80946.2 was verified by network and experiments to be a key functional lncRNA under cold stress, regulating ACP1, TSPY1 and Tsn.

Integrating long noncoding RNAs and mRNAs expression profiles of response to Plasmodiophora brassicae infection in Pakchoi (Brassica campestris ssp. chinensis Makino)

The biotrophic protist Plasmodiophora brassicae causes serious damage to Brassicaceae crops grown worldwide. However, the molecular mechanism of the Brassica rapa response remains has not been determined. Long noncoding RNA and mRNA expression profiles in response to Plasmodiophora brassicae infection were investigated using RNA-seq on the Chinese cabbage inbred line C22 infected with P. brassicae. Approximately 5,193 mRNAs were significantly differentially expressed, among which 1,345 were upregulated and 3,848 were downregulated. The GO enrichment analysis shows that most of these mRNAs are related to the defense response. Meanwhile, 114 significantly differentially expressed lncRNAs were identified, including 31 upregulated and 83 downregulated. Furthermore, a total of 2,344 interaction relationships were detected between 1,725 mRNAs and 103 lncRNAs with a correlation coefficient greater than 0.8. We also found 15 P. brassicaerelated mRNAs and 16 lncRNA interactions within the correlation network. The functional annotation showed that 15 mRNAs belong to defense response proteins (66.67%), protein phosphorylation (13.33%), root hair cell differentiation (13.33%) and regulation of salicylic acid biosynthetic process (6.67%). KEGG annotation showed that the vast majority of these genes are involved in the biosynthesis of secondary metabolism pathways and plant-pathogen interactions. These results provide a new perspective on lncRNA-mRNA network function and help to elucidate the molecular mechanism of P. brassicae infection.

Integrated dual RNA-seq and dual iTRAQ of infected tissue reveals the functions of a diguanylate cyclase gene of Pseudomonas plecoglossicida in host-pathogen interactions with Epinephelus coioides

The interactions between host and pathogen is exceedingly complex, which involves alterations at multiple molecular layers. However, research to simultaneously monitor the alterations of transcriptome and proteome between a bacterial pathogen and aquatic animal host through integrated dual RNA-seq and dual iTRAQ of tissue during infection is currently lacking. The important role of a diguanylate cyclase gene (L321_RS15240) in pathogenicity of Pseudomonas plecoglossicida against Epinephelus coioides was suggested by previous dual RNA-seq of our lab. Then L321_RS15240-RNAi strains of P. plecoglossicida were constructed with pCM130/tac, and the mutant with the best silencing effect was selected for follow-up study. The RNAi of L321_RS15240 resulted in a significant decrease in bacterial virulence of P. plecoglossicida. The E. coioides spleens infected by wild type strain or L321_RS15240-RNAi strain of P. plecoglossicida were subjected to dual RNA-seq and dual iTRAQ, respectively. The results showed that: RNAi of L321_RS15240 led to 1)alterations of host transcriptome associated with complement and coagulation cascades, ribosome, arginine and proline metabolism, and oxidative phosphorylation; 2)high expression of host proteins which related to phagosome and metabolism responses (metabolism of glutathione, amino sugar and nucleotide sugar); 3)the highly differentially expression of host lncRNAs and miRNAs. The differentially expressed proteins and mRNAs of pathogen were different after infection, but the functions of these proteins and mRNAs were mainly related to metabolism and virulence. This study provides a new insight to comprehensively understand the gene functions of pathogens and hosts at multiple molecular layers during in vivo infection.

Non-coding RNAs: Functional roles in the regulation of stress response in Brassica crops

Brassica crops face a combination of different abiotic and biotic stresses in the field that can reduce plant growth and development by affecting biochemical and morpho-physiological processes. Emerging evidence suggests that non-coding RNAs (ncRNAs), especially microRNAs (miRNAs) and long ncRNAs (lncRNAs), play a significant role in the modulation of gene expression in response to plant stresses. Recent advances in computational and experimental approaches are of great interest for identifying and functionally characterizing ncRNAs. While progress in this field is limited, numerous ncRNAs involved in the regulation of gene expression in response to stress have been reported in Brassica. In this review, we summarize the modes of action and functions of stress-related miRNAs and lncRNAs in Brassica as well as the approaches used to identify ncRNAs.

Comparative Proteomics Indicates That Redox Homeostasis Is Involved in High- and Low-Temperature Stress Tolerance in a Novel Wucai (Brassica campestris L.) Genotype

The genotype WS-1, previously identified from novel wucai germplasm, is tolerant to both low-temperature (LT) and high-temperature (HT) stress. However, it is unclear which signal transduction pathway or acclimation mechanisms are involved in the temperature-stress response. In this study, we used the proteomic method of tandem mass tag (TMT) coupled with liquid chromatography-mass spectrometry (LC-MS/MS) to identify 1022 differentially expressed proteins (DEPs) common to WS-1, treated with either LT or HT. Among these 1022 DEPs, 172 were upregulated in response to both LT and HT, 324 were downregulated in response to both LT and HT, and 526 were upregulated in response to one temperature stress and downregulated in response to the other. To illustrate the common regulatory pathway in WS-1, 172 upregulated DEPs were further analyzed. The redox homeostasis, photosynthesis, carbohydrate metabolism, heat-shockprotein, and chaperones and signal transduction pathways were identified to be associated with temperature stress tolerance in wucai. In addition, 35S:BcccrGLU1 overexpressed in Arabidopsis, exhibited higher reduced glutathione (GSH) content and reduced glutathione/oxidized glutathione (GSH/GSSG) ratio and less oxidative damage under temperature stress. This result is consistent with the dynamic regulation of the relevant proteins involved in redox homeostasis. These data demonstrate that maintaining redox homeostasis is an important common regulatory pathway for tolerance to temperature stress in novel wucai germplasm.

Long noncoding RNAs in Brassica rapa L. following vernalization

Brassica rapa L. is an important agricultural crop that requires a period of prolonged cold for flowering. This process is known as vernalization. Studies have shown that long noncoding RNAs (lncRNAs) play important roles in abiotic stress responses and several cold-responsive noncoding RNAs have been suggested to be involved in vernalization. We examined the transcriptome of the Chinese cabbage inbred line (B. rapa L. var. pekinensis) RJKB-T24, and identified 1,444 long intergenic noncoding RNAs (lincRNAs), 551 natural antisense transcripts (NATs), and 93 intronic noncoding RNAs (incRNAs); 549 of the 2,088 lncRNAs significantly altered their expression in response to four weeks of cold treatment. Most differentially expressed lncRNAs did not lead to a change of expression levels in mRNAs covering or near lncRNAs, suggesting that the transcriptional responses to four weeks of cold treatment in lncRNA and mRNA are independent. However, some differentially expressed mRNAs had NATs with expression altered in the same direction. These genes were categorized as having an abiotic stress response, suggesting that the paired-expression may play a role in the transcriptional response to vernalization or cold treatment. We also identified short-term cold treatment induced NATs in BrFLC and BrMAF genes, which are involved in vernalization. The lncRNAs we identified differed from those reported in Arabidopsis thaliana, suggesting the role of lncRNAs in vernalization differ between these two species.

Genome-wide analysis of long non-coding RNAs unveils the regulatory roles in the heat tolerance of Chinese cabbage (Brassica rapa ssp.chinensis)

Long non-coding RNAs (lncRNAs) mediate important epigenetic regulation in various biological processes related to the stress response in plants. However, the systematic analysis of the lncRNAs expressed in Brassica rapa under heat stress has been elusive. In this study, we performed a genome-wide analysis of the lncRNA expression profiles in non-heading Chinese cabbage leaves using strand-specific RNA-sequencing. A total of 4594 putative lncRNAs were identified with a comprehensive landscape of dynamic lncRNA expression networks under heat stress. Co-expression networks of the interactions among the differentially expressed lncRNAs, mRNAs and microRNAs revealed that several phytohormones were associated with heat tolerance, including salicylic acid (SA) and brassinosteroid (BR) pathways. Of particular importance is the discovery of 25 lncRNAs that were highly co-expressed with 10 heat responsive genes. Thirty-nine lncRNAs were predicted as endogenous target mimics (eTMs) for 35 miRNAs, and five of them were validated to be involved in the heat tolerance of Chinese cabbage. Heat responsive lncRNA (TCONS_00048391) is an eTM for bra-miR164a, that could be a sponge for miRNA binding and may be a competing endogenous RNA (ceRNA) for the target gene NAC1 (Bra030820), affecting the expression of bra-miR164a in Chinese cabbage. Thus, these findings provide new insights into the functions of lncRNAs in heat tolerance and highlight a set of candidate lncRNAs for further studies in non-heading Chinese cabbage.

Genome-wide identification of oil biosynthesis-related long non-coding RNAs in allopolyploid Brassica napus

Background

Long noncoding RNAs (lncRNAs) are transcripts longer than 200 bp that do not encode proteins but nonetheless have been shown to play important roles in various biological processes in plants. Brassica napus is an important seed oil crop worldwide and the target of many genetic improvement activities. To understand better the function of lncRNAs in regulating plant metabolic activities, we carried out a genome-wide lncRNA identification of lncRNAs in Brassica napus with a focus on lncRNAs involved in lipid metabolism. Twenty ribosomal RNA depleted strand specific RNA-seq (ssRNA-seq) datasets were generatred using RNAs isolated from B. napus seeds at four developmental stages. For comparison we also included 30 publically available RNA-seq datasets generated from poly(A) enriched mRNAs isolated from from various Brassica napus tissues in our analysis.

Results

A total of 8905 lncRNA loci were identified, including 7100 long intergenic noncoding RNA (lincRNA) loci and 1805 loci generating long noncoding natural antisense transcript (lncNAT). Many lncRNAs were identified only in the ssRNA-seq and poly(A) RNA-seq dataset, suggesting that B. napus has a large lncRNA repertoire and it is necessary to use libraries prepared from different tissues and developmental stages as well as different library preparation approaches to capture the whole spectrum of lncRNAs. Analysis of coexpression networks revealed that among the regulatory modules are networks containing lncRNAs and protein-coding genes related to oil biosynthesis indicating a possible role of lncRNAs in the control of lipid metabolism. One such example is that several lncRNAs are potential regulators of BnaC08g11970D that encodes oleosin1, a protein found in oil bodies and involved in seed lipid accumulation. We also observed that the expression levels of B. napus lncRNAs is positively correlated with their conservation levels.

Conclusions

We demonstrated that the B. napus genome has a large number of lncRNA and that these lncRNAs are expressed broadly across many developmental times and in different tissue types. We also provide evidence indicating that specific lncRNAs appear to be important regulators of lipid biosynthesis forming regulatory networks with transcripts involved in lipid biosynthesis. We also provide evidence that these lncRNAs are conserved in other species of the Brassicaceae family.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5117-8) contains supplementary material, which is available to authorized users.

Integration of RNAi and RNA-seq uncovers the immune responses of Epinephelus coioides to L321_RS19110 gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is a threatening and important pathogen in aquaculture and causes heavy losses. Expression of L321_RS19110 was found significant up-regulated at 18 °C than at 28 °C, which was confirmed by quantitative real-time PCR. RNAi significantly reduced the content of L321_RS19110 mRNA of P. plecoglossicida, and exhibited the best efficiency of gene silencing with a reduction of 84.9%. Compared with the wild type strain, the infection of L321_RS19110-RNAi-1 strain resulted in the onset time delay, and 30% reduction in mortality of Epinephelus coioides, as well as alleviates in the symptoms of E. coioides spleen. Moreover, compared with wild type strain, the gene silence of L321_RS19110 in P. plecoglossicida resulted in a significant change in transcriptome of infected E. coioides. The results of KEGG analysis showed that genes of chemokine signaling pathway and cytokine-cytokine receptor interaction, digestion and absorption of proteins (non-specific immune) and antigen processing and presentation pathways were most affected by L321_RS19110-RNAi of P. plecoglossicida. In these immune pathways, the most affected immune genes were associated with different number of non-coding RNAs. Among them, most lncRNAs and miRNAs exerted negative regulation on immune genes. The results indicated that L321_RS19110 was an important virulent gene of P. plecoglossicida, the up-regulation of the immune pathways made E. coioides more likely to remove L321_RS19110-RNAi strain than the wild type strain of P. plecoglossicida, the expression of immune genes were regulated by miRNA and lncRNA in a complex way.

Integration of RNAi and RNA-seq Reveals the Immune Responses of Epinephelus coioides to sigX Gene of Pseudomonas plecoglossicida

Pseudomonas plecoglossicida is an important pathogen for aquaculture and causes high mortality in various marine fishes. Expression of sigX was found significantly up-regulated at 18°C than at 28°C, which was verified by quantitative real-time PCR (qRT-PCR). RNAi significantly reduced the content of sigX mRNA of P. plecoglossicida, whether in in vitro or in the spleen at all sampling time points. Compared with the wild-type strain, the infection of sigX-RNAi strain resulted in the onset time delay, and 20% reduction in mortality of Epinephelus coioides, as well as alleviates in the symptoms of E. coioides spleen. Compared with wild-type strain, the gene silence of sigX in P. plecoglossicida resulted in a significant change in transcriptome of infected E. coioides. The result of gene ontology and KEGG analysis on E. coioides showed that genes of serine-type endopeptidase and chemokine signaling pathway, coagulation and complement system, and intestinal immune network for IgA production pathway were mostly affected by sigX of P. plecoglossicida. Meanwhile, the immune genes were associated with different number of miRNA and lncRNA, and some miRNAs were associated with more than one gene at the same time. The results indicated that sigX was a virulent gene of P. plecoglossicida. The up-regulation of the immune pathways made E. coioides more likely to kill sigX-RNAi strain than the wild-type strain of P. plecoglossicida, while the immune genes were regulated by miRNA and lncRNA by a complex mode.

Re-analysis of long non-coding RNAs and prediction of circRNAs reveal their novel roles in susceptible tomato following TYLCV infection

BACKGROUND

Long Noncoding-RNAs (LncRNAs) are known to be involved in some biological processes, but their roles in plant-virus interactions remain largely unexplored. While circular RNAs (circRNAs) have been studied in animals, there has yet to be extensive research on them in a plant system, especially in tomato-tomato yellow leaf curl virus (TYLCV) interaction.

RESULTS

In this study, RNA transcripts from the susceptible tomato line JS-CT-9210 either infected with TYLCV or untreated, were sequenced in a pair-end strand-specific manner using ribo-zero rRNA removal library method. A total of 2056 lncRNAs including 1767 long intergenic non-coding RNA (lincRNAs) and 289 long non-coding natural antisense transcripts (lncNATs) were obtained. The expression patterns in lncRNAs were similar in susceptible tomato plants between control check (CK) and TYLCV infected samples. Our analysis suggested that lncRNAs likely played a role in a variety of functions, including plant hormone signaling, protein processing in the endoplasmic reticulum, RNA transport, ribosome function, photosynthesis, glulathione metabolism, and plant-pathogen interactions. Using virus-induced gene silencing (VIGS) analysis, we found that reduced expression of the lncRNA S-slylnc0957 resulted in enhanced resistance to TYLCV in susceptible tomato plants. Moreover, we identified 184 circRNAs candidates using the CircRNA Identifier (CIRI) software, of which 32 circRNAs were specifically expressed in untreated samples and 83 circRNAs in TYLCV samples. Approximately 62% of these circRNAs were derived from exons. We validated the circRNAs by both PCR and Sanger sequencing using divergent primers, and found that most of circRNAs were derived from the exons of protein coding genes. The silencing of these circRNAs parent genes resulted in decreased TYLCV virus accumulation.

CONCLUSION

In this study, we identified novel lncRNAs and circRNAs using bioinformatic approaches and showed that these RNAs function as negative regulators of TYLCV infection. Moreover, the expression patterns of lncRNAs in susceptible tomato plants were different from that of resistant tomato plants, while exonic circRNAs expression positively associated with their respective protein coding genes. This work provides a foundation for elaborating the novel roles of lncRNAs and circRNAs in susceptible tomatoes following TYLCV infection.

The DNA Methylome and Association of Differentially Methylated Regions with Differential Gene Expression during Heat Stress in Brassica rapa

Cytosine DNA methylation is a critical epigenetic mechanism in the silencing of transposable elements, imprinting and regulating gene expression. However, little is known about the potential role of mC in response to heat stress. To determine and explore the functions of the dynamic DNA methylome during heat stress, we characterized single-base resolution methylome maps of Brassica rapa and assessed the dynamic changes of mC under heat stress using whole genome bisulfite sequencing. On average, the DNA methylation levels of CG, CHG and CHH are 39.3%, 15.38% and 5.24% in non-heading Chinese cabbage (NHCC), respectively. We found that the patterns of methylation are similar to other eudicot plants, but with higher CHH methylation levels. Further comparative analysis revealed varying patterns for three sequence contexts (mCG, mCHG and mCHH) under heat stress indicating context- and position-dependent methylation regulation. DNA methylation near the TSS and TES may be closely associated with methylation-dependent transcriptional silencing. Association analysis of differential methylation and differential gene expression revealed a different set of methDEGs involved at early and late stages under heat stress. The systemic characterization of the dynamic DNA methylome during heat stress will improve our understanding of the mechanism of epigenetic regulation under heat stress.

Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells

Long non-coding RNAs (lncRNAs) are increasingly regarded as a key role in regulating diverse biological processes in various tissues and species. Although the cold responsive lncRNAs have been reported in plants, no data is available on screening and functional prediction of lncRNAs in cold acclimation in fish so far. Here we compared the expression profile of lncRNAs in cold acclimated zebrafish embryonic fibroblast cells (ZF4) cultured at 18°C for 30 days with that of cells cultured at 28°C as control by high-throughput sequencing. Totally 8,363 novel lncRNAs were identified. Including known and novel lncRNAs, there are 347 lncRNAs up-regulated and 342 lncRNAs down-regulated in cold acclimated cells. Among the differentially expressed lncRNAs, 74 and 61 were detected only in control cells or cold-acclimated cells, respectively. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses of adjacent genes to the differentially expressed lncRNAs showed that the enriched genes are involved in electron transport, cell adhesion, oxidation-reduction process, and so on. We also predicted the target genes of the differentially expressed lncRNAs by looking for interactions between lncRNAs and mRNAs, and constructed an interaction network. In summary, our genome-wide systematic identification and functional prediction of cold responsive lncRNAs in zebrafish cells suggests a crucial role of lincRNAs in cold acclimation in fish.

RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed

Perilla frutescen is used as traditional food and medicine in East Asia. Its seeds contain high levels of α-linolenic acid (ALA), which is important for health, but is scarce in our daily meals. Previous reports on RNA-seq of perilla seed had identified fatty acid (FA) and triacylglycerol (TAG) synthesis genes, but the underlying mechanism of ALA biosynthesis and its regulation still need to be further explored. So we conducted Illumina RNA-sequencing in seven temporal developmental stages of perilla seeds. Sequencing generated a total of 127 million clean reads, containing 15.88 Gb of valid data. The de novo assembly of sequence reads yielded 64,156 unigenes with an average length of 777 bp. A total of 39,760 unigenes were annotated and 11,693 unigenes were found to be differentially expressed in all samples. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 486 unigenes were annotated in the “lipid metabolism” pathway. Of these, 150 unigenes were found to be involved in fatty acid (FA) biosynthesis and triacylglycerol (TAG) assembly in perilla seeds. A coexpression analysis showed that a total of 104 genes were highly coexpressed (r > 0.95). The coexpression network could be divided into two main subnetworks showing over expression in the medium or earlier and late phases, respectively. In order to identify the putative regulatory genes, a transcription factor (TF) analysis was performed. This led to the identification of 45 gene families, mainly including the AP2-EREBP, bHLH, MYB, and NAC families, etc. After coexpression analysis of TFs with highly expression of FAD2 and FAD3 genes, 162 TFs were found to be significantly associated with two FAD genes (r > 0.95). Those TFs were predicted to be the key regulatory factors in ALA biosynthesis in perilla seed. The qRT-PCR analysis also verified the relevance of expression pattern between two FAD genes and partial candidate TFs. Although it has been reported that some TFs are involved in seed development, more direct evidence is still needed to verify their function. However, these findings can provide clues to reveal the possible molecular mechanisms of ALA biosynthesis and its regulation in perilla seed.

Development and comparative genomic mapping of Dasypyrum villosum 6V#4S-specific PCR markers using transcriptome data

Twenty-five Dasypyrum villosum 6V#4S-specific PCR markers were developed using transcriptome data and further assigned to comparative genomic maps of wheat chromosome 6A, 6B, and 6D and barley chromosome 6H contrasting their homologous genes in these genomes. Two Dasypyrum villosum accessions, D.v#2 and No. 1026 from England and Russia, respectively, contain Pm21 on chromosome 6V#2S and PmV on chromosome 6V#4S. Both genes confer high resistance to powdery mildew (PM) in wheat. Even though several molecular markers have been developed to detect Pm21 and PmV, only the MBH1 marker can simultaneously detect both Pm21 and PmV. In this study, we first used a high-throughput sequencing technique to obtain the transcriptome sequences of a wheat-D. villosum translocation line, Pm97033-which contains chromosome 6V#4S carrying the PmV locus, under wheat PM pathogen induction. Twenty-five 6V#4S chromosome-specific markers were developed. Three of them were able to clearly distinguish chromosomes 6V#4S and 6V#2S by product size, four amplified the product specific for chromosome 6V#4S only, and the remaining 18 markers identified chromosome 6VS in wheat backgrounds. Two different D. villosum accessions, their derived translocation lines and wheat varieties carrying different chromosome 6VS were identified using these specific markers. The 25 newly developed markers together with the known PM resistance gene Stpk-V were used to construct comparative genomic maps with the homoeologous chromosome arms of wheat and barley. The colinearity of the identified gene sequences amplified by the 25 markers among wheat chromosomes 6A, 6B, and 6D and barley chromosome 6H was not very conserved and interrupted frequently by inversion and insertion. Our markers have potential in marker assisted selection for PM resistance breeding, and for locating other potential important genes and cloning the PmV gene on chromosome 6V#4S.

Comparative transcriptome analysis between resistant and susceptible tomato allows the identification of lncRNA16397 conferring resistance to Phytophthora infestans by co-expressing glutaredoxin

The rapid development of omics sequencing technology has facilitated the identification of thousands of long non-coding (lnc)RNAs in plant species, but the role of lncRNAs in plant-pathogen interactions remains largely unexplored. We used comparative transcriptome analysis of Phytophthora infestans-resistant and -susceptible tomatoes to identify differentially expressed genes (DEGs) and lncRNAs (DELs), and examine lncRNA-mRNA networks. A total of 1037 DEGs and 688 DELs were identified between P. infestans-resistant and -susceptible tomatoes. The co-localization networks, including 128 DEGs and 127 DELs, were performed. We found that lncRNA16397 acted as an antisense transcript of SlGRX22 to regulate its expression, and also induced SlGRX21 expression when lncRNA16397 was overexpressed. In addition, disease symptoms and reactive oxygen species (ROS) accumulation in tomatoes overexpressing lncRNA16397 and SpGRX were fewer and lower than those in wild-type after P. infestans infection. This result suggests that tomato lncRNA16397 induces SlGRX expression to reduce ROS accumulation and alleviate cell membrane injury, resulting in enhanced resistance to P. infestans. Our results provide insight into lncRNAs involved in the response of tomato to P. infestans infection, demonstrate that the lncRNA16397-GRXs network is an important component of the P. infestans network in tomato, and provide candidates for breeding to enhance biotic stress-resistance in tomato.

The Role of Long Noncoding RNAs in Plant Stress Tolerance

Plants must adapt to multiple biotic and abiotic stresses ; thus, sensing and responding to environmental signals is imperative for their survival. Moreover, understanding these responses is imperative for efforts to improve plant yield and consistency. Regulation of transcript levels is a key aspect of the plant response to environmental signals. Long noncoding RNAs (lncRNAs) have gained widespread attention in recent years with the advance of high-throughput sequencing technologies. As important biological regulators, lncRNAs have been implicated in a wide range of developmental processes and diseases in animals. However, knowledge of the role that lncRNAs play in plant stress tolerance remains limited. Here, we review recent studies on the identification, characteristics, classification, and biological functions of lncRNAs in response to various stresses, including bacterial pathogens, excess light, drought, salinity, hypoxia, extreme temperatures, and nitrogen/phosphate deficiency. We also discuss possible directions for future research.