Dynamic characteristics and functional analysis provide new insights into long non-coding RNA responsive to Verticillium dahliae infection in Gossypium hirsutum
BMC Plant Biology volume 21, Article number: 68 (2021)
Verticillium wilt is a widespread and destructive disease, which causes serious loss of cotton yield and quality. Long non-coding RNA (lncRNA) is involved in many biological processes, such as plant disease resistance response, through a variety of regulatory mechanisms, but their possible roles in cotton against Verticillium dahliae infection remain largely unclear.
Here, we measured the transcriptome of resistant G. hirsutum following infection by V. dahliae and 4277 differentially expressed lncRNAs (delncRNAs) were identified. Localization and abundance analysis revealed that delncRNAs were biased distribution on chromosomes. We explored the dynamic characteristics of disease resistance related lncRNAs in chromosome distribution, induced expression profiles, biological function, and these lncRNAs were divided into three categories according to their induced expression profiles. For the delncRNAs, 687 cis-acting pairs and 14,600 trans-acting pairs of lncRNA-mRNA were identified, which indicated that trans-acting was the main way of Verticillium wilt resistance-associated lncRNAs regulating target mRNAs in cotton. Analyzing the regulation pattern of delncRNAs revealed that cis-acting and trans-acting lncRNAs had different ways to influence target genes. Gene Ontology (GO) enrichment analysis revealed that the regulatory function of delncRNAs participated significantly in stimulus response process, kinase activity and plasma membrane components. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis indicated that delncRNAs participated in some important disease resistance pathways, such as plant-pathogen interaction, alpha-linolenic acid metabolism and plant hormone signal transduction. Additionally, 21 delncRNAs and 10 target genes were identified as being involved in alpha-linolenic acid metabolism associated with the biosynthesis of jasmonic acid (JA). Subsequently, we found that GhlncLOX3 might regulate resistance to V. dahliae through modulating the expression of GhLOX3 implicated in JA biosynthesis. Further functional analysis showed that GhlncLOX3-silenced seedlings displayed a reduced resistance to V. dahliae, with down-regulated expression of GhLOX3 and decreased content of JA.
This study shows the dynamic characteristics of delncRNAs in multiaspect, and suggests that GhlncLOX3-GhLOX3-JA network participates in response to V. dahliae invasion. Our results provide novel insights for genetic improvement of Verticillium wilt resistance in cotton using lncRNAs.
Cotton is one of the most important economic crops for its natural fiber and oil seed, and has been widely cultivated around the world. The quality and yield of cotton are frequently subjected to the serious threat from Verticillium wilt which is one of the most destructive diseases in cotton and deserves enormous researches to control efficiently . Verticillium wilt is caused by Verticillium dahliae, a soil-borne, xylem-invading, hemi-biotrophic fungal pathogen. More than 400 plant species were invaded by V. dahliae, causing an intractable vascular wilt disease and even death [1, 2]. Its microsclerotia can survive in soil for 15 years [3, 4], becoming a serious challenge to the control on Verticillium wilt.
Jasmonic acid (JA), as a vital plant hormone, plays an important regulatory role in plant responses to biotic and abiotic stresses. Increasing evidence has shown that the perception of pathogen attack promotes synthesis of JA from lipid, and the accumulation of JA activates downstream defensive genes expression in plant, thereby protects plant against the attacker [5, 6]. Using JA-insensitive mutants, it has been found that JA-dependent defense pathways in Arabidopsis contribute to resistance against the fungal pathogens Alternaria brassicicola and Botrytis cinerea . It has been reported that JA response can reduce tomato susceptibility to V. dahliae and Fusarium oxysporum . Importantly, it has been demonstrated that JA response pathway is vital to confer Verticillium wilt resistance in cotton [9,10,11].
In recent years, high-throughput technologies have been widely used to monitor expression profiles and identify a large number of differentially expressed genes/proteins in cotton inoculated with V. dahliae. Based on transcriptional analysis and proteomic analysis, it has been revealed that plant hormone signal transduction, phenylpropanoid biosynthesis pathway, lignin metabolism, cell wall relating enzymes/proteins, adenosine triphosphate (ATP)-binding cassette (ABC) proteins, reactive oxygen species and gossypol play important roles in the fight against V. dahliae in cotton [11,12,13,14,15,16,17,18]. In addition, a large number of Verticillium-resistance genes and proteins were identified. For example, previously we obtained 3027 Verticillium-resistance unigenes in G. barbadense cv Pima90–53 , and identified 1717 and 1476 differentially abundant proteins in resistant G. hirsutum cv. ND601 and susceptible G. hirsutum cv. CCRI8 after infection with V. dahliae, respectively . Moreover, there were some microRNAs (miRNAs) identified related to cotton defense against V. dahliae by high-throughput sequencing . These studies not only provide substantial information for research on molecular mechanism of cotton resistance to Verticillium wilt, but also lay a foundation for understanding the function of long non-coding RNA (lncRNA) in regulating plant response to V. dahliae infection.
LncRNAs are longer than 200 nucleotides (nt) RNAs without protein-coding ability . LncRNAs play important roles in flowering time regulation [21,22,23,24], photomorphogenesis , reproductive development , fruit development [27,28,29], biotic and abiotic stress responses [30,31,32]. Recently, the roles of lncRNAs in plant-pathogen interaction have attracted much attention. In wheat, lncRNAs participate in the response to powdery mildew and stripe rust infection [33, 34]. Several lncRNAs responding to F. oxysporum infection have been identified in Arabidopsis . A total of 110 lncRNAs responding to phytoplasma infection have been identified in Paulownia by high-throughput sequencing . It has been found that a number of lncRNAs are involved in response to Sclerotinia sclerotiorum infection in Brassica napus . LncRNA ELENA1 has been shown to regulate PATHOGENESIS-RELATED GENE1 (PR1) through interacting with Mediator subunits 19a (MED19a), resulting in enhancing Arabidopsis resistance to Pseudomonas syringe . More recently, it has been found that tomato lncRNA23468 can execute function as a competing endogenous RNA to modulate NBS-LRR genes by decoying miR482b in the interaction between tomato and Phytophthora infestans .
As was previously reported, a number of functional lncRNAs were detected in cotton. It was discovered that 1296 lncRNAs were related to fiber initiation, and 720 lncRNAs played a role in fiber elongation and secondary cell wall biosynthesis [39, 40]. Some lncRNAs were identified in G. hirsutum as being involved in regulating plant hormone against drought stress . A lncRNA (lncRNA973) was confirmed to play an important role in increasing salt tolerance of cotton . Furthermore, 1236 and 1907 lncRNAs were found in response to V. dahliae infection in resistant G. barbadense and susceptible G. hirsutum, respectively, and lineage-species specific (LS) lncRNAs were identified in response to the disease . However, the function and dynamic characteristics of lncRNAs in resistant G. hirsutum remain unknown.
In this study, we used RNA-seq technology to excavate Verticillium-resistance lncRNAs in resistant G. hirsutum and revealed the function of important lncRNAs by virus-induced gene silencing (VIGS) technology. We explored the dynamic characteristics of lncRNAs associated with disease response in some aspects, such as chromosome distribution, induced expression profiles, biological function, and linked the infection process by V. dahliae to lncRNAs dynamic expression profiles. Trans-acting was indicated as the main way of Verticillium wilt resistance-associated lncRNAs regulating the target mRNAs in cotton. The important function of GhlncLOX3 involved in cotton resistance to V. dahliae via GhlncLOX3-GhLOX3-JA network was revealed.
Up to 4277 lncRNAs were induced expression in cotton roots after inoculation with V. dahliae
In order to investigate the transcriptome of resistant G. hirsutum in response to V. dahliae infection, total RNAs isolated from V. dahliae-inoculated plants (cv. Nongda601) at 2, 6, 12, 24 and 48 h post inoculation (hpi) (VD2hpi, VD6hpi, VD12hpi, VD24hpi, and VD48hpi) were sequenced using a high-throughput RNA-seq approach. For comparison, the transcriptomes of mock-treated samples (MT2hpi, MT6hpi, MT12hpi, MT24hpi, and MT48hpi) were also sequenced. Through quality filters, about one billion clean reads were obtained. By transcript assembly, differential isoform and gene expression analysis using Cufflinks, approximately 77.1 to 83.3% of these clean reads were aligned to the reference genome of G. hirsutum . Approximately 89% of the aligned reads were uniquely mapped to a single genomic locus, attesting to the high quality of the RNA-seq reads and the reference cotton genome (Additional file 1: Table S1). By cufflinks assemblies, 380,559 transcripts were obtained ultimately, of which 16,876 were calculated lncRNAs.
LncRNAs with a greater than 2-fold expression change and q-value < 0.05 were recognized as differentially expressed. As a result, 4277 delncRNAs were identified (Fig. 1a). To validate the findings from sequencing data, the expression levels of six delncRNAs, randomly chosen, were assessed using qRT-PCR in all samples. The results from RNA sequencing and qRT-PCR showed a strong correlation (R = 0.917, p < 0.01) (Fig. 1b), suggesting the high quality of the RNA-seq.
DelncRNAs showed obvious dynamic characteristics in chromosome distribution and induced expression profiles
Through comparison of the number of delncRNAs at each of the five time points, we found that the abundance of delncRNAs differed evidently. A total of 1074, 2411, 1284, 1307 and 1802 lncRNAs were differentially expressed at 2 hpi, 6 hpi, 12 hpi, 24 hpi and 48 hpi, respectively. And these delncRNAs were mainly up-regulated, especially at 6 hpi, up to 85.9% (Table 1). These delncRNAs exhibited a high degree of temporal specificity: substantial number of delncRNAs were unique to 2 hpi, 6 hpi, 12 hpi, 24 hpi, or 48 hpi. Especially at 6 hpi, about 42.4% of delncRNAs (1022) were unique (Fig. 2a).
Localization analysis revealed that the distribution of lncRNAs induced by V. dahliae was not uniform on chromosomes. There were higher densities of delncRNAs on At05, At11 and Dt05, and lower on At06, Dt10 and At02. In addition, the delncRNAs tended to be distributed on chromosome ends. We also found that the distribution of delncRNAs on chromosomes varied with the infection process of V. dahliae (Fig. 2b).
All of the 4277 delncRNAs were divided into three categories according to their induced expression profiles, and each category had some distinctively induced expression patterns (Fig. 2c). Type I and type II clusters represented down- and up-regulated lncRNAs induced by V. dahliae, approximately 21.7 and 76.7%, respectively. Even though only 68 chimeric-regulated lncRNAs were classified as type III, they had the most complex expression patterns upon V. dahliae infection. For instance, there were 10 lncRNAs which were up-regulated at 6 hpi and then down-regulated at later one or two time points. Type I and type II clusters could be divided into two subclasses, respectively. We found that delncRNAs in the subtypes I-1 and II-1 clusters were down-regulated and up-regulated at only one time point, respectively. While delncRNAs in subtypes I-2 and II-2 clusters, more complicated, were down-regulated and up-regulated at least two time points and the induction rate varied with the infection times. These results suggested that delncRNAs had obvious dynamic characteristics in cotton after inoculation with V. dahliae.
Trans-acting was the main way of delncRNAs regulating the target mRNAs
It has been reported that lncRNAs can modulate gene activity through cis-acting and/or trans-acting [44, 45]. To analyze the potential regulation patterns of the delncRNAs, we selected the differentially expressed protein-coding genes that composed the cis- or trans- pairs with delncRNAs. In this study, 687 cis-acting pairs of lncRNA-mRNA were identified, and, more importantly, a large number of trans-acting pairs (14,600 pairs) were found. Obviously, trans-acting pairs were absolutely predominant in identified lncRNA-mRNA target pairs, accounting for 95.51% (Fig. 3a), which indicated that trans-acting was the main way of V. dahliae-induced lncRNAs regulating the target mRNAs in cotton.
As shown in the Fig. 3b, trans-acting pairs could be divided into three patterns. Pattern 1 contained 14,146 lncRNA-mRNA target pairs with 756 up-regulated lncRNAs and 2405 up-regulated mRNAs; Pattern 2 included 448 lncRNA-mRNA target pairs with 78 down-regulated lncRNAs and 290 down-regulated mRNAs; Pattern 3 only had 6 lncRNA-mRNA target pairs with 2 up-regulated lncRNAs and 6 down-regulated mRNAs. Obviously, the lncRNA-mRNA target pairs belonging to pattern 1 and 2 were positive regulatory action pairs, while those belonging to pattern 3 were negative regulatory action pairs, accounting for 99.96 and 0.04%, respectively (Fig. 3a and b). These suggested that the influence of trans-acting lncRNAs on its predicted target protein-coding genes was dominated by positive regulation.
The cis-acting pairs were put in four patterns, based on regulatory model of lncRNAs to its target mRNAs (Fig. 3c). Pattern 1 included 255 lncRNA-mRNA target pairs consisting of 214 up-regulated lncRNAs and 243 up-regulated mRNAs; Pattern 2 contained 129 lncRNA-mRNA target pairs composing of 92 down-regulated lncRNAs and 129 down-regulated mRNAs; Pattern 3 exhibited 223 lncRNA-mRNA target pairs comprising of 169 up-regulated lncRNAs and 209 down-regulated mRNAs; Pattern 4 displayed 80 lncRNA-mRNA target pairs with 66 down-regulated lncRNAs and 76 up-regulated mRNAs. Obviously, positive regulatory action pairs were clustered into pattern 1 and pattern 2, and negative regulatory action pairs were divided into pattern 3 and pattern 4, accounting for 55.90 and 44.10%, respectively (Fig. 3a and c). Comparing to trans-acting lncRNAs, negative regulation played an equally important role as positive regulation in cis-acting lncRNAs regulating their target mRNAs.
DelncRNAs participated in many disease resistance processes and exhibited a temporal specificity
In order to predict the potential function of delncRNAs, their target mRNAs were subjected to GO annotation and KEGG enrichment analysis. In the category of biological process, most of terms were related to stimulus response and protein modification (Fig. 4a). The top 10 terms of molecular function were related to kinase activity, protein binding and nucleotide binding (Fig. 4b). With regards to cellular components, biomembrane potentially played important role (Fig. 4c). In addition, the top 10 significant pathways were shown in Fig. 4d. It was interesting to note that glycolysis, plant-pathogen interaction, alpha-linolenic acid metabolism, fatty acid degradation, plant hormone signal transduction, terpenoid backbone biosynthesis, peroxisome and flavonoid biosynthesis were all associated with disease resistance, suggesting that these delncRNAs were implicated in responding to V. dahliae infection.
The GO enrichment analysis also showed that the function of unique delncRNAs varied with V. dahliae infection time (Fig. 2d). Concretely, it mainly related to stimulus response, kinase activity and phosphorylation at 2 hpi, stimulus response and hypoxia response at 6 hpi, stimulus response and single-organism process at 12 hpi, biomembrane at 24 hpi and cellular response to hypoxia at 48 hpi. These results indicated that lncRNAs could make a series of sequential responses rapidly to the infection of V. dahliae and further confirmed the dynamic characteristics of delncRNAs during cotton resistance to V. dahliae.
Regulatory roles and characteristics of functional candidate lncRNAs in JA biosynthesis
The result of KEGG enrichment analysis showed that alpha-linolenic acid metabolism was one of the significantly enriched pathways after infection with V. dahliae. Additionally, alpha-linolenic acid metabolism is associated with the biosynthesis of JA. Previous studies have showed that JA plays an important role in cotton resistance to V. dahliae . In this study, 22 lncRNA-mRNA pairs (21 lncRNAs and 10 mRNAs) were found involving in JA biosynthesis pathway, and they were assigned to four gene families: lipoxygenase (LOX), allene oxide synthase (AOS), allene oxide cyclase (AOC) and acyl-CoA oxidase (ACX) (Additional file 3: Table S2). Specifically, two lncRNAs were involved in regulating two LOX family members, five lncRNAs had regulatory relationships with three target genes of AOS family members, one lncRNA could be a potential regulatory factor of two AOC family members, 13 lncRNAs participated in the regulation of three ACX family members. The expression profiles of 21 delncRNAs involved in JA biosynthesis pathways were analyzed based on the normalized FPKM values. The results showed that all the delncRNAs were up-regulated in the cotton roots attacked by V. dahliae (Fig. 5). The correlation analysis indicated that all of the 22 lncRNA-mRNA pairs were positively regulated (Additional file 3: Table S2). In the regulatory relationships between lncRNAs and their target genes, only five lncRNA-mRNA pairs were in a one-to-one regulatory relationship, and the remaining 17 lncRNA-mRNA pairs were not in a one-to-one regulatory relationship, such as GhACX (Gh_A09G2434) regulated by 11 lncRNAs (Additional file 4: Fig. S2). These results suggested that the regulatory roles of these 21 lncRNAs on their target mRNAs were complex.
GhlncLOX3 was up-regulated in response to V. dahliae invasion and potentially related to JA biosynthesis
LOX is an important enzyme involved in JA biosynthesis via alpha-linolenic acid metabolism pathway in plant. According to the analysis of lncRNA-mRNA pairs, we found that GhLOX3 belonged to the target protein-coding gene of trans-acting lncRNA lnc_015092, designed as GhlncLOX3. Moreover, the expression of GhlncLOX3 was up-regulated in all time points and peaked at 6 hpi for about 30 folds. Thus, we inferred that GhlncLOX3 might be resistant to V. dahliae through regulating GhLOX3 expression to influence JA synthesis.
Further, quantitative real-time PCR experiments were carried out on eight cotton varieties with different resistance to V. dahliae. The results showed that GhlncLOX3 was up-regulated in all varieties at 6 hpi. The fold up-regulation of GhlncLOX3 was positively correlated with the disease resistance of cotton varieties (Fig. 6a). Moreover, the target gene GhLOX3 has a similar expression trend with GhlncLOX3 in all varieties (Additional file 5: Fig. S3). When GhlncLOX3 was silenced via VIGS (Fig. 6e), the plants showed enhanced susceptibility to V. dahliae with more serious wilting and chlorosis of leaf comparing to CK (Fig. 6c). Additionally, the disease index (DI) of GhlncLOX3-silenced plants was significantly higher than that of controls at 25 dpi (day post inoculation) (Fig. 6d). Compared with the control, GhlncLOX3-silenced plants had lower expression level of GhLOX3 and JA content (Fig. 6f and g). The correlation analysis indicated that the expression level of GhLOX3 was positively correlated with that of GhlncLOX3 (Fig. 6h). These results confirmed our hypothesis that GhlncLOX3 affected JA synthesis by regulating the expression of GhLOX3.
IntaRNA program was executed to calculate potential binding sites between GhlncLOX3 and GhLOX3 . As a result, a conceivable binding site with a free energy of − 17.54 k cal/mol was detected (Fig. 6b). LncRNA can exert its effects via direct mRNA binding and recruitment to some RNA binding proteins to regulate mRNA stability [47, 48]. Thus, we inferred that GhlncLOX3 might block GhLOX3 degradation via recruitment to some unknown RNA binding proteins. The result further suggested that GhlncLOX3 regulated the disease resistance of cotton to V. dahliae via GhlncLOX3-GhLOX3-JA model. Based on the results of this study and previous reports, we proposed a schematic diagram for the role of GhlncLOX3 in raising the resistance of cotton to V. dahliae (Fig. 6i).
LncRNA is a type of molecule with important functions in a wide range of biological processes across eukaryotes [44, 49]. It is a challenge to ascertain the biological function of lncRNAs due to the complex and diverse action mechanisms. Accumulating evidence suggests that lncRNAs can modulate some target genes via cis-acting and/or trans-acting [44, 45]. Nevertheless, most previous studies focused on cis-acting lncRNAs [1, 45, 50, 51]. For example, 92 disease resistance response cis-acting pairs of lncRNA-mRNA were identified in cotton . In our research, the function of trans-acting lncRNA to resist V. dahliae in cotton was firstly established; and 14,600 trans-acting pairs of lncRNA-mRNA were identified as related to disease resistance (Fig. 3b). The number of trans-acting pair was 21 times as much as that of cis-acting pair, approximately. Therefore, we suggested that more attention should be paid to plant trans-action lncRNAs in the future. Moreover, we studied the regulation patterns of cis-acting lncRNAs and trans-action lncRNAs. The results unveiled that the regulation patterns of cis-acting lncRNAs and trans-acting lncRNAs on their target mRNAs were not exactly the same. Trans-acting lncRNAs were mainly positive regulation, while the negative regulation played the same important role as positive regulation in cis-acting lncRNAs regulating their target mRNAs. These results represent the first to characterize the important status and the regulation pattern of trans-acting lncRNAs involved in plant responses to infection by V. dahliae. In addition, our study is mutually complementary with the previous research of lncRNAs to resist V. dahliae in cotton, for resistant G. hirsutum is complementary to susceptible G. hirsutum, and trans-action lncRNA is complementary to cis-action lncRNA.
Verticillium wilt responsive lncRNAs had obvious dynamic characteristics (Fig. 2), which probably closely related to the infection process of V. dahliae. Recently, Zhang and Zhao monitored the infection processes of V. dahliae in root tissue of Arabidopsis and cotton, respectively [18, 52]. At 6 hpi, the conidia covered the root surface of both the Arabidopsis and cotton; By 24 hpi, massive conidia began to germinate and extended to form neogenetic hyphae; At 48 hpi, a mass of hyphae were formed in cotton root [18, 52]. In this research, the unique delncRNAs were enriched to stimulus response and hypoxia response at 6 hpi, which might be due to the fact that the local hypoxia of the root perimeter caused by spores covering; At 24 hpi, it manly related to biomembrane, which might be due to the fact that the extension of neonatal hyphae affected the membrane structure; At 48 hpi, it enriched to cellular response to hypoxia, which might be due to the fact that oxygen was snatched by a large number of hyphae, resulting in intracellular hypoxia in root tissue.
LncRNAs have a variety of complex mechanisms to regulate genes expression at multiple levels [1, 49, 53]. Recent works in plants have identified a large number of lncRNAs and greatly enhanced our knowledge on lncRNAs biology. however, the detailed functions and action mechanisms of lncRNAs are still in its infancy [38, 39, 49]. In this study, lncRNAs have been categorized into several GO and pathway terms based on annotation of target protein-coding genes (Fig. 4), which will help us to understand the potential functions of lncRNAs in response to pathogen infection.
Plants have to encounter various biotic and abiotic stresses throughout their development process. To deal with these adversities, plants have evolved constitutive and inducible defense systems [7, 11]. It is now clear that several endogenous signal molecules, such as salicylic acid (SA), ethylene (ET) and JA, are synthesized and can activate a series of complex defense signaling networks [11, 54]. Previous studies have showed that JA plays an important role in cotton resistance to V. dahliae . Recently, 18 tea lncRNAs have been identified affecting JA biosynthesis through regulating the expression of JA-related genes . Moreover, lncRNAs show low conservation among species . In this study, 21 lncRNAs were found involving in JA biosynthesis pathway through complex regulatory networks.
LOX is a key enzyme of JA biosynthesis in plant and plays an important role in response to biotic stresses . For example, it was indicated that TomLoxD might be resistant to B. cinerea through regulating JA biosynthesis to influence the expression of plant defense genes ; CmLOX09 played a positive role in JA biosynthesis through the AOS pathway, and then increased resistance to Podosphaera xanthii . In this study, transcriptomes data analyses showed a linkage between GhlncLOX3 and GhLOX3, which prompted us to further investigate the function of GhlncLOX3. In order to figure out whether GhlncLOX3 is involved in cotton resistance to Verticillium wilt and has a positive regulatory relationship with GhLOX3, we tested the expression changes of GhlncLOX3 and GhLOX3 in eight different cotton varieties with V. dahliae infection. Moreover, GhlncLOX3 was silenced using the VIGS technique, and we detected that GhLOX3 and JA were significant down-regulated in GhlncLOX3-silenced plants. Thus, we preliminarily suggest that GhlncLOX3-GhLOX3-JA is an important network for cotton resistance to V. dahliae, although the further researches by RNA interference (RNAi) and over-expression need to be supplemented. In addition, a potential binding site between GhlncLOX3 and GhLOX3 was predicted (Fig. 6b). Recent studies have shown that lncRNAs can exert its effects via direct mRNA binding and recruitment to some RNA binding proteins to regulate mRNA stability [47, 48]. Therefore, we speculated that GhlncLOX3 might enhance the stability of GhLOX3 via nucleic acid binding and recruitment to some unknown RNA binding proteins. Although this conjecture needs further study to verify.
In this study, using a high-throughput RNA-seq approach, 4277 delncRNAs are identified. And these delncRNAs show obvious dynamic characteristics in temporal specificity, chromosome distribution, induced expression profiles and biological function. Although previous studies pay more attention to cis-acting, this research reveals that trans-acting is the main way of V. dahliae-induced lncRNAs regulating the target mRNAs in cotton. Silencing of GhlncLOX3 in cotton suppresses the accumulation of JA and reduces disease resistance, with down-regulation of GhLOX3. These results suggest that GhlncLOX3-GhLOX3-JA network plays an important role in cotton resistance to V. dahliae infection. Our results will extend the current view on lncRNA defence against V. dahliae infection and provide novel insights for genetic improvement of Verticillium wilt resistance in cotton.
Plant material and fungal pathogen
Resistant G. hirsutum cv. Nongda601 and Nongdamian8 were bred by our laboratory. Previously we also screened 419 accessions representing a core collection of cotton germplasm resources . The disease index (DI) of 419 accessions and Nongda601 have been identified at 25 days post inoculation (dpi) by our laboratory using at least 30 seedlings, which have not been published. The seedlings of Nongda601 were used to build strand-specific cDNA libraries for transcriptome analysis. We selected 8 varieties (Nongdamian8, Nianmian9, Shaansanyuan78–782, Zhong078, Lumian6, Dongting3, E82–6078 and Xin547) with different resistance to V. dahliae from 419 core germplasm resources for further functional validation.
The cotton seedlings were grown in 50% Hoagland’s solution under greenhouse conditions of 25 °C for 4 weeks. Changed the nutrient solution, once every 4 days, in order to ensure the healthy growth of seedlings. A highly aggressive V. dahliae strain Linxi2–1 , was cultured on Potato Dextrose Agar (PDA) for 2 weeks from storage at 4 °C, and then sub-cultured into Czapek’s medium for 10 days on a shaker (130 rpm) at 25 °C. The conidial suspension was adjusted to a density of 107 spores per milliliter with 50% Hoagland’s solution, as the final concentration for inoculation.
RNA isolation and sequencing
When the fourth true leaf spread, the cotton seedlings of Nongda601 were inoculated with the spores of V. dahliae. According to our previous research about the infection process of V. dahlia , roots were harvested at 2, 6, 12, 24 and 48 hpi. As mock treatments, the plants treated with 50% Hoagland’s solution were also collected at the corresponding time points. High-quality total RNAs of all samples were extracted using the RNA prep Pure Plant Kit (TIANGEN Biotech, Beijing, China). rRNA was removed by Ribo-zero rRNA Removal Kit (Epicentre, USA). The strand-specific cDNA libraries were generated using the rRNA-depleted RNA by NEBNext Ultra Directional RNA Library Prep Kit for Illumina (NEB, USA) according to the manufacturer’s recommendations. Sequencing was performed on the Illumina Hiseq 4000 platform and 150 bp paired-end reads were generated.
All sequence data were firstly processed by removing reads containing adapter, reads containing ploy-N and low-quality reads. The reference genome and the annotation files of G. hirsutum were downloaded from the CottonGen database (http://www.cottongen.org). Index of the reference genome was built using Bowtie2 v2.2.8 and paired-end clean reads were aligned to the reference genome using TopHat2 . The mapped reads of each sample were assembled by Cufflinks . According to the characteristics of lncRNA, we filtered the background noise based on Fragments Per Kilobase of exon per Million fragments mapped (FPKM), length, exon and coding potential (FPKM > 0.5; length > 200; exon ≥1; coding: NO).
Cufflinks provided computing program to calculate FPKM of both lncRNA and coding gene in each sample, and Cuffdiff provided statistical routines for determining differential expression data using a model based on the negative binomial distribution . The differential expression transcripts were obtained by the following criteria: adjusted P value < 0.05 and at least two-fold FPKM change .
All pathway enrichment and GO terms of listed genes were annotated using KOBAS 3.0 (http://kobas.cbi.pku.edu.cn/) by comparing to the reference genome background (P < 0.01).
Identification of trans-acting and cis-acting pairs
The correlation between lncRNA and mRNA was analyzed by Pearson correlation coefficient method to predict the target gene of lncRNA. The trans-acting pairs were obtained by the following criteria: absolute value of correlation > 0.95. Cis-acting is lncRNA acting on neighboring target genes. We set the threshold of cis-acting pairs to 100 kb, then lncRNA and mRNA with the same expression pattern were screened out by STME software .
As previously reported [1, 60], specific primers were designed to amplify GhlncLOX3 fragment and cloned into the TRV2:00 vector. Primer sequences are listed in Additional file 6: Table S3. Empty vector TRV2:00 and containing candidate fragment vector TRV2:GhlncLOX3 were, respectively, co-infiltrated with equal amount TRV1 via Agrobacterium tumefaciens GV3101 into cotton seedlings of Nongda601 when the cotyledons had spread. TRV2:CLA1 (CLOROPLASTOS ALTERADOS 1) and TRV2:00 were utilized as positive and negative control, respectively.
VIGS plants inoculation
When the bleaching phenotype of positive controls appeared, we started to perform the inoculation with at least 30 plants for each treatment using V. dahliae strain Linxi2–1 with at least three biological replicates. At 6 hpi, leaves were harvested from several plants to assess GhlncLOX3 silence degree and analyze changes in GhLOX3 expression level and JA content. At 25 dpi, the phenotype was observed and the DI for plant populations was calculated according to previously described .
Quantitative real-time (qRT) PCR analysis
Total plant RNA was extracted from cotton root or leaf using RNAprep Pure Plant Kit (TIANGEN Biotech, Beijing, China). The first stranded cDNA was synthesized from 0.5 μg RNA using the PrimeScript RT-PCR kit (TaKaRa, Dalian, China). Specific primers (Additional file 6: Table S3) were designed to implement the qRT-PCR on an ABI 7500 Real Time PCR system (Applied Biosystems, Waltham, Massachusetts, USA) with the SYBR Premix Ex Taq II system (TaKaRa, Dalian, China). The program of qRT-PCR was as follows: 95 °C for 2 min, followed by 40 cycles at 95 °C for 20 s, 55–60 °C for 20 s and 72 °C for 20 s. Gene expression levels were normalized to GhUBQ14 (ubiquitin 14) expression .
Availability of data and materials
The data generated or analyzed during the current study are included in this published article and its supplemental data files and available from the corresponding author on reasonable request.
Long non-coding RNA
Differentially expressed lncRNA
Kyoto Encyclopedia of Genes and Genomes
- PR1 :
Mediator subunits 19a
Virus-induced gene silencing
Hours post inoculation
Day post inoculation
Potato Dextrose Agar
- CLA1 :
CLOROPLASTOS ALTERADOS 1
Zhang L, Wang M, Li N, Wang H, Qiu P, Pei L, Xu Z, Wang T, Gao E, Liu J, et al. Long noncoding RNAs involve in resistance to Verticillium dahliae, a fungal disease in cotton. Plant Biotechnol J. 2018;16(6):1172–85.
Zhang T, Jin Y, Zhao JH, Gao F, Zhou BJ, Fang YY, Guo HS. Host-induced gene silencing of the target gene in fungal cells confers effective resistance to the cotton wilt disease pathogen Verticillium dahliae. Mol Plant. 2016;9(6):939–42.
Zhang DD, Wang J, Wang D, Kong ZQ, Zhou L, Zhang GY, Gui YJ, Li JJ, Huang JQ, Wang BL, et al. Population genomics demystifies the defoliation phenotype in the plant pathogen Verticillium dahliae. New Phytol. 2019;222(2):1012–29.
Fradin EF, Thomma BP. Physiology and molecular aspects of Verticillium wilt diseases caused by V. dahliae and V. albo-atrum. Mol Plant Pathol. 2006;7(2):71–86.
Wasternack C, Hause B. Jasmonates: biosynthesis, perception, signal transduction and action in plant stress response, growth and development. An update to the 2007 review in Annals of Botany. Ann Bot. 2013;111(6):1021–58.
Zhang L, Zhang F, Melotto M, Yao J, He SY. Jasmonate signaling and manipulation by pathogens and insects. J Exp Bot. 2017;68(6):1371–85.
Thomma BP, Eggermont K, Penninckx IA, Mauch-Mani B, Vogelsang R, Cammue BP, Broekaert WF. Separate jasmonate-dependent and salicylatedependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens. Proc Natl Acad Sci USA. 1998;95(25):15107–11.
Thaler JS, Owen B, Higgins VJ. The role of the jasmonate response in plant susceptibility to diverse pathogens with a range of lifestyles. Plant Physiol. 2004;135(1):530–8.
Li TG, Wang BL, Yin CM, Zhang DD, Wang D, Song J, Zhou L, Kong ZQ, Klosterman SJ, Li JJ, et al. The Gossypium hirsutum TIR-NBS-LRR gene GhDSC1 mediates resistance against Verticillium wilt. Mol Plant Pathol. 2019;20(6):857–76.
Li X, Pei Y, Sun Y, Liu N, Wang P, Liu D, Ge X, Li F, Hou Y. A cotton cyclin-dependent kinase E confers resistance to Verticillium dahliae mediated by jasmonate-responsive pathway. Front Plant Sci. 2018;9:642.
Gao W, Long L, Zhu LF, Xu L, Gao WH, Sun LQ, Liu LL, Zhang XL. Proteomic and virus-induced gene silencing (VIGS) analyses reveal that gossypol, brassinosteroids, and jasmonic acid contribute to the resistance of cotton to Verticillium dahliae. Mol Cell Proteomics. 2013;12(12):3690–703.
Xu L, Zhu L, Tu L, Liu L, Yuan D, Jin L, Long L, Zhang X. Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011;62(15):5607–21.
Zhang Y, Wang XF, Ding ZG, Ma Q, Zhang GR, Zhang SL, Li ZK, Wu LQ, Zhang GY, Ma ZY. Transcriptome profiling of Gossypium barbadense inoculated with Verticillium dahliae provides a resource for cotton improvement. BMC Genomics. 2013;14:637.
Fang W, Xie D, Zhu H, Li W, Xu Z, Yang L, Li Z, Sun L, Wang J, Nie L, et al. Comparative proteomic analysis of Gossypium thurberi in response to Verticillium dahliae inoculation. Int J Mol Sci. 2015;16(10):25121–40.
Li PT, MHO R, Chen TT, Lu QW, Ge Q, Gong WK, Liu AY, Gong JW, Shang HH, Deng XY, et al. Transcriptomic and biochemical analysis of upland cotton (Gossypium hirsutum) and a chromosome segment substitution line from G. hirsutum x G. barbadense in response to Verticillium dahliae infection. BMC Plant Biol. 2019;19(1):19.
Dong Q, Magwanga RO, Cai X, Lu P, Nyangasi Kirungu J, Zhou Z, Wang X, Wang X, Xu Y, Hou Y, et al. RNA-sequencing, physiological and RNAi analyses provide insights into the response mechanism of the ABC-mediated resistance to Verticillium dahliae infection in cotton. Genes (Basel). 2019;10(2):110.
Yang J, Wang X, Xie M, Wang G, Li Z, Zhang Y, Wu L, Zhang G, Ma Z. Proteomic analyses on xylem sap provides insights into the defense response of Gossypium hirsutum against Verticillium dahliae. J Proteome. 2020;213:103599.
Zhang Y, Wang X, Rong W, Yang J, Li Z, Wu L, Zhang G, Ma Z. Histochemical analyses reveal that stronger intrinsic defenses in Gossypium barbadense than in G. hirsutum are associated with resistance to Verticillium dahliae. Mol Plant-Microbe Interact. 2017;30(12):984–96.
Zhang Y, Wang W, Chen J, Liu J, Xia M, Shen F. Identification of miRNAs and their targets in cotton inoculated with Verticillium dahliae by high-throughput sequencing and degradome analysis. Int J Mol Sci. 2015;16(7):14749–68.
Mercer TR, Dinger ME, Mattick JS. Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009;10(3):155–9.
Heo JB, Sung S. Vernalization-mediated epigenetic silencing by a long intronic noncoding RNA. Science. 2011;331(6013):76–9.
Wu X, Shi T, Iqbal S, Zhang Y, Liu L, Gao Z. Genome-wide discovery and characterization of flower development related long non-coding RNAs in Prunus mume. BMC Plant Biol. 2019;19(1):64.
Swiezewski S, Liu F, Magusin A, Dean C. Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target. Nature. 2009;462(7274):799–802.
Liu F, Marquardt S, Lister C, Swiezewski S, Dean C. Targeted 3′ processing of antisense transcripts triggers Arabidopsis FLC chromatin silencing. Science. 2010;327(5961):94–7.
Wang H, Chung PJ, Liu J, Jang IC, Kean MJ, Xu J, Chua NH. Genome-wide identification of long noncoding natural antisense transcripts and their responses to light in Arabidopsis. Genome Res. 2014;24(3):444–53.
Wang M, Wu HJ, Fang J, Chu C, Wang XJ. A long noncoding RNA involved in rice reproductive development by negatively regulating Osa-miR160. Sci Bull. 2017;62(7):470–5.
Li R, Fu D, Zhu B, Luo Y, Zhu H. CRISPR/Cas9-mediated mutagenesis of lncRNA1459 alters tomato fruit ripening. Plant J. 2018;94(3):513–24.
Wang X, Ai G, Zhang C, Cui L, Wang J, Li H, Zhang J, Ye Z. Expression and diversification analysis reveals transposable elements play important roles in the origin of Lycopersicon-specific lncRNAs in tomato. New Phytol. 2016;209(4):1442–55.
Wang M, Zhao W, Gao L, Zhao L. Genome-wide profiling of long non-coding RNAs from tomato and a comparison with mRNAs associated with the regulation of fruit ripening. BMC Plant Biol. 2018;18(1):75.
Zhu QH, Stephen S, Taylor J, Helliwell CA, Wang MB. Long noncoding RNAs responsive to Fusarium oxysporum infection in Arabidopsis thaliana. New Phytol. 2014;201(2):574–84.
Qin T, Zhao H, Cui P, Albesher N, Xiong L. A nucleus-localized long non-coding RNA enhances drought and salt stress tolerance. Plant Physiol. 2017;175(3):1321–36.
Wang J, Yang Y, Jin L, Ling X, Liu T, Chen T, Ji Y, Yu W, Zhang B. Re-analysis of long non-coding RNAs and prediction of circRNAs reveal their novel roles in susceptible tomato following TYLCV infection. BMC Plant Biol. 2018;18(1):104.
Xin M, Wang Y, Yao Y, Song N, Hu Z, Qin D, Xie C, Peng H, Ni Z, Sun Q. Identification and characterization of wheat long non-protein coding RNAs responsive to powdery mildew infection and heat stress by using microarray analysis and SBS sequencing. BMC Plant Biol. 2011;11:61.
Zhang H, Chen X, Wang C, Xu Z, Wang Y, Liu X, Kang Z, Ji W. Long non-coding genes implicated in response to stripe rust pathogen stress in wheat (Triticum aestivum L.). Mol Biol Rep. 2013;40(11):6245–53.
Wang Z, Zhai X, Cao Y, Dong Y, Fan G. Long non-coding RNAs responsive to witches’ broom disease in Paulownia tomentosa. Forests. 2017;8(9):348.
Joshi RK, Megha S, Basu U, Rahman MH, Kav NN. Genome wide identification and functional prediction of long non-coding RNAs responsive to Sclerotinia sclerotiorum infection in Brassica napus. PLoS One. 2016;11(7):e0158784.
Seo JS, Sun HX, Park BS, Huang CH, Yeh SD, Jung C, Chua NH. ELF18-INDUCED LONG-NONCODING RNA associates with mediator to enhance expression of innate immune response genes in Arabidopsis. Plant Cell. 2017;29(5):1024–38.
Jiang N, Cui J, Shi Y, Yang G, Zhou X, Hou X, Meng J, Luan Y. Tomato lncRNA23468 functions as a competing endogenous RNA to modulate NBS-LRR genes by decoying miR482b in the tomato-Phytophthora infestans interaction. Hortic Res. 2019;6:28.
Wang M, Yuan D, Tu L, Gao W, He Y, Hu H, Wang P, Liu N, Lindsey K, Zhang X. Long noncoding RNAs and their proposed functions in fibre development of cotton (Gossypium spp.). New Phytol. 2015;207(4):1181–97.
Hu H, Wang M, Ding Y, Zhu S, Zhao G, Tu L, Zhang X. Transcriptomic repertoires depict the initiation of lint and fuzz fibres in cotton (Gossypium hirsutum L.). Plant Biotechnol J. 2018;16(5):1002–12.
Lu X, Chen X, Mu M, Wang J, Wang X, Wang D, Yin Z, Fan W, Wang S, Guo L, et al. Genome-wide analysis of long noncoding RNAs and their responses to drought stress in cotton (Gossypium hirsutum L.). PLoS One. 2016;11(6):e0156723.
Zhang X, Dong J, Deng F, Wang W, Cheng Y, Song L, Hu M, Shen J, Xu Q, Shen F. The long non-coding RNA lncRNA973 is involved in cotton response to salt stress. BMC Plant Biol. 2019;19(1):459.
Zhang T, Hu Y, Jiang W, Fang L, Guan X, Chen J, Zhang J, Saski CA, Scheffler BE, Stelly DM, et al. Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat Biotechnol. 2015;33(5):531–7.
Kim ED, Sung S. Long noncoding RNA: unveiling hidden layer of gene regulatory networks. Trends Plant Sci. 2012;17(1):16–21.
Deng F, Zhang X, Wang W, Yuan R, Shen F. Identification of Gossypium hirsutum long non-coding RNAs (lncRNAs) under salt stress. BMC Plant Biol. 2018;18(1):23.
Mann M, Wright PR, Backofen R. IntaRNA 2.0: enhanced and customizable prediction of RNA-RNA interactions. Nucleic Acids Res. 2017;45(W1):W435–9.
Schmitz SU, Grote P, Herrmann BG. Mechanisms of long noncoding RNA function in development and disease. Cell Mol Life Sci. 2016;73(13):2491–509.
Kretz M, Siprashvili Z, Chu C, Webster DE, Zehnder A, Qu K, Lee CS, Flockhart RJ, Groff AF, Chow J, et al. Control of somatic tissue differentiation by the long non-coding RNA TINCR. Nature. 2013;493(7431):231–5.
Chekanova JA. Long non-coding RNAs and their functions in plants. Curr Opin Plant Biol. 2015;27:207–16.
Pauli A, Valen E, Lin MF, Garber M, Vastenhouw NL, Levin JZ, Fan L, Sandelin A, Rinn JL, Regev A, et al. Systematic identification of long noncoding RNAs expressed during zebrafish embryogenesis. Genome Res. 2012;22(3):577–91.
Rinn JL, Kertesz M, Wang JK, Squazzo SL, Xu X, Brugmann SA, Goodnough LH, Helms JA, Farnham PJ, Segal E, et al. Functional demarcation of active and silent chromatin domains in human HOX loci by noncoding RNAs. Cell. 2007;129(7):1311–23.
Zhao P, Zhao YL, Jin Y, Zhang T, Guo HS. Colonization process of Arabidopsis thaliana roots by a green fluorescent protein-tagged isolate of Verticillium dahliae. Protein Cell. 2014;5(2):94–8.
Cabili MN, Trapnell C, Goff L, Koziol M, Tazon-Vega B, Regev A, Rinn JL. Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. Genes Dev. 2011;25(18):1915–27.
Pieterse CMJ, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SCM. Hormonal modulation of plant immunity. Annu Rev Cell Dev Biol. 2012;28(1):489–521.
Zhu C, Zhang S, Fu H, Zhou C, Chen L, Li X, Lin Y, Lai Z, Guo Y. Transcriptome and phytochemical analyses provide new insights into long non-coding RNAs modulating characteristic secondary metabolites of oolong tea (Camellia sinensis) in solar-withering. Front Plant Sci. 2019;10:1638.
Yan Q, Wu F, Yan Z, Li J, Ma T, Zhang Y, Zhao Y, Wang Y, Zhang J. Differential co-expression networks of long non-coding RNAs and mRNAs in Cleistogenes songorica under water stress and during recovery. BMC Plant Biol. 2019;19(1):23.
Zhu J, Wang X, Guo L, Xu Q, Zhao S, Li F, Yan X, Liu S, Wei C. Characterization and alternative splicing profiles of the lipoxygenase gene family in tea plant (Camellia sinensis). Plant Cell Physiol. 2018;59(9):1765–81.
Yan L, Zhai Q, Wei J, Li S, Wang B, Huang T, Du M, Sun J, Kang L, Li CB, et al. Role of tomato lipoxygenase D in wound-induced jasmonate biosynthesis and plant immunity to insect herbivores. PLoS Genet. 2013;9(12):e1003964.
Ju LJ, Zhang C, Liao JJ, Li YP, Qi HY. An oriental melon 9-lipoxygenase gene CmLOX09 response to stresses, hormones, and signal substances. J Zhejiang Univ Sci B. 2018;19(8):596–609.
Ma Z, He S, Wang X, Sun J, Zhang Y, Zhang G, Wu L, Li Z, Liu Z, Sun G, et al. Resequencing a core collection of upland cotton identifies genomic variation and loci influencing fiber quality and yield. Nat Genet. 2018;50(6):803–13.
Wang GN, Zhao GY, Yue XW, Li ZK, Zhang Y, Zhang GY, Wu LQ, Wang XF, Ma ZY. Pathogenicity and ISSR genetic differentiation of Verticillium dahliae isolates from cotton growing areas of Hebei province. Cotton Sci. 2012;24(4):348–57.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinformatics. 2006;7:191.
Artico S, Nardeli SM, Brilhante O, Grossi-de-Sa MF, Alves-Ferreira M. Identification and evaluation of new reference genes in Gossypium hirsutum for accurate normalization of real-time quantitative RT-PCR data. BMC Plant Biol. 2010;10:49.
All authors are grateful to the laboratory members for help, advice and discussion.
This work was supported by the National Natural Science Foundation of China (31871672), the China Agriculture Research System (CARS-15-03), the Science and Technology Support Program of Hebei Province (16226307D), the Outstanding Youth Fund of Hebei Province (C2019204365), and the Science and Technology Program of Baoding, China (1911ZN011). The funding bodies provided the financial support to this research, including experimental design and implementation, sampling and data analysis. No funder played the role in data collection and analysis and writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Summary of RNA -seq data.
The induced expression patterns of delncRNAs that were not included in Fig. 2c are shown here.
Correlation analysis of lncRNA-mRNA pairs related to JA biosynthesis pathway in cotton.
The regulatory relationships between lncRNAs and their target genes.
GhLOX3 gene expression fold change and disease index of infected plants in 8 cultivars. ‘MT’ and ‘VD’ mean mock treatment and seedling roots inoculated with V. dahliae.
List of PCR primers used in this study.
About this article
Cite this article
Wang, G., Wang, X., Zhang, Y. et al. Dynamic characteristics and functional analysis provide new insights into long non-coding RNA responsive to Verticillium dahliae infection in Gossypium hirsutum. BMC Plant Biol 21, 68 (2021). https://doi.org/10.1186/s12870-021-02835-8
- Verticillium dahliae