Functional Annotation, Prediction of Cis Target Gene for the Sequence of Mrnas, and Candidate Long Noncoding RNAs from Tea (camellia Sinensis Var. Sinensis) Leaves During Infection by the Fungal Pathogen Epicoccum Nigrum
Hongke Huang,Yuanyou Yang,Chen Huang,Zhongqiu Xia,Yuqin Yang,Xinyue Jiang,Delu Wang,Zhuo Chen
DOI: https://doi.org/10.1094/phytofr-03-22-0015-a
2022-01-01
PhytoFrontiers™
Abstract:HomePhytoFrontiers™Vol. 2, No. 4Functional Annotation, Prediction of Cis Target Gene for the Sequence of mRNAs, and Candidate Long Noncoding RNAs from Tea (Camellia sinensis var. sinensis) Leaves During Infection by the Fungal Pathogen Epicoccum nigrum Previous RESOURCE ANNOUNCEMENT OPENOpen Access licenseFunctional Annotation, Prediction of Cis Target Gene for the Sequence of mRNAs, and Candidate Long Noncoding RNAs from Tea (Camellia sinensis var. sinensis) Leaves During Infection by the Fungal Pathogen Epicoccum nigrumHongke Huang, Yuanyou Yang, Chen Huang, Zhongqiu Xia, Yuqin Yang, Xinyue Jiang, Delu Wang, and Zhuo ChenHongke HuangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaCollege of Tea Science, Guizhou University, Guiyang, Guizhou 550025, China, Yuanyou YangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Chen HuangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Zhongqiu XiaKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaCollege of Tea Science, Guizhou University, Guiyang, Guizhou 550025, China, Yuqin YangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaCollege of Tea Science, Guizhou University, Guiyang, Guizhou 550025, China, Xinyue JiangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Delu Wang†Corresponding authors: D. Wang; E-mail Address: [email protected], and Z. Chen; E-mail Address: [email protected]College of Forestry, Guizhou University, Guiyang, Guizhou 550025, China, and Zhuo Chen†Corresponding authors: D. Wang; E-mail Address: [email protected], and Z. Chen; E-mail Address: [email protected]http://orcid.org/0000-0001-7130-8457Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaAffiliationsAuthors and Affiliations Hongke Huang1 2 Yuanyou Yang1 Chen Huang1 Zhongqiu Xia1 2 Yuqin Yang1 2 Xinyue Jiang1 Delu Wang3 † Zhuo Chen1 † 1Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China 2College of Tea Science, Guizhou University, Guiyang, Guizhou 550025, China 3College of Forestry, Guizhou University, Guiyang, Guizhou 550025, China Published Online:26 Jul 2022https://doi.org/10.1094/PHYTOFR-03-22-0015-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articlemRNA and lncRNA AnnouncementThe Didymellaceae family consists of 19 monophyletic generic clades, each representing an individual genus such as Epicoccum, in which some species are phytopathogens (Chen et al. 2017). Epicoccum nigrum was found in India, China, the United States, and The Netherlands, and can cause leaf spot on Lablab purpureus, loquat, cowpea, and tea (Chen et al. 2017; Deepika et al. 2021; Mahadevakumar et al. 2014; Wu et al. 2017). In addition, as a minor pathogen, E. nigrum contributes to a fungal complex with the primary pathogens Alternaria alternata and Diaporthe ilicicola sp. nov. to cause fruit rot of deciduous holly (Lin et al. 2018).As an emerging disease, tea leaf spot caused by E. nigrum can reduce the productivity and quality of tea (Camellia sinensis var. sinensis) in Guizhou Province, southwestern China. In this study, high-throughput sequencing of messenger RNAs (mRNAs) and long noncoding RNAs (lncRNAs) was conducted using the Illumina Novaseq 6000 platform.lncRNAs are a heterogeneous group of transcripts which are longer than 200 nucleotides and have no protein-coding ability. lncRNAs resemble mRNA, with short, open reading frames. Most annotated lncRNAs are transcribed by RNA polymerase II (Pol II) and the plant-specific RNA polymerases Pol IV and V (Batista and Chang 2013; Chi et al. 2019; Matzke and Mosher 2014; Wierzbicki et al. 2008). lncRNAs can interact with mRNA, DNA, protein, and microRNAs and, consequently, regulate gene expression, being involved with epigenetic, transcriptional, post-transcriptional, and translational regulation as well as post-translational modification (Dykes and Emanueli 2017; Heo et al. 2013; Vance and Ponting 2014). Recently, lncRNAs in plants were found to be involved with plant disease resistance during the process of infection, with the interaction between lncRNAs and mRNAs from plants contributing to defense against infection (Nie et al. 2021; Song et al. 2021; Wang et al. 2021; Zaynab et al. 2018; Zhang et al. 2020).Tea brown leaf spot caused by E. nigrum was first found on tea (C. sinensis (L.) O. Kuntze) in Sinan County (27.74°N, 108.35°E) and Kaiyang County (27.96°N, 107.34°E), Guizhou Province, southwest China. The disease results in huge tea leaf production losses and decreases of tea leaf quality during cold spells in late spring and at high-altitude in Guizhou Province (Yin et al. 2022). Interestingly, previous studies indicated that E. sorghinum, E. nigrum, E. mackenziei, Didymella segeticola, and D. bellidis, as members in the Didymellaceae family, can cause tea leaf spot in Guizhou Province, and E. sorghinum, E. nigrum, D. segeticola, and D. bellidis represented more virulence and induced lesions with a relatively larger size compared with E. mackenziei under laboratory conditions (Bao et al. 2019; Wang et al. 2020; Yin et al. 2022; Zhao et al. 2018). To date, there have been no published reports on the sequence, annotation, and the interactions between lncRNAs and mRNAs from tea leaves infected by E. nigrum. Thus, analysis of the expression abundance of lncRNAs and mRNAs originating from tea leaves during infection by E. nigrum and their interactive relationship will provide a valuable resource for future research into this host–pathogen interaction.In this study, 3-year-old potted tea plants (C. sinensis var. sinensis 'Fuding-dabaicha') were grown in a greenhouse under day and night temperatures of 25 and 20°C, respectively, with a cycle of 14 h/day and 8 h/night and relative humidity of 70 to 80%. A single-conidium isolate of E. nigrum strain CGMCC3.20484 was cultured on potato dextrose agar at 28°C in the dark for 15 days to obtain a vigorously growing colony. The third or fourth leaf on each detached, healthy tea twig was used in the inoculation assay. Either side of the main vein on each tea leaf was punctured twice using sterilized syringe needles, making four clustered holes within a region 6 mm in diameter. For the infected treatment, a conidial suspension of E. nigrum strain CGMCC3.20484 (106 conidia/ml, 20 to 30 μl/inoculation site) or mycelial plugs (6 mm in diameter) were used to inoculate the physically wounded sites on the adaxial surface of tea leaves. For the uninfected treatment, sterile water was used to inoculate the wounded tea leaves. For both the infected and uninfected treatments, 15 tea leaves (1 from each twig) were pooled to represent one replicate and three biological replicates were carried out for each treatment. At 48 h postinoculation, a lesion formed rapidly on each leaf site inoculated with the conidial suspension or mycelial plugs but not at the sites inoculated with sterile water. The tissues neighboring the inoculation sites of each leaf were harvested using a sterile hole punch at 48 h after inoculation. The tissues were ground into a fine powder under the protection of liquid nitrogen. The total RNA was extracted at 48 h postinoculation using a TRIzol Plus RNA Purification Kit (Invitrogen, Carlsbad, CA, U.S.A.), following the manufacturer's procedure. Approximately 5 μg of total RNA was used to deplete ribosomal RNA according to the protocol of the Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, U.S.A.). A cDNA library was prepared using a TruSeq Stranded mRNA Sample Preparation Kit (Illumina), and sequenced on the Illumina Novaseq 6000 by Lc-Bio Technologies Co., Ltd. (Hangzhou, China), following the manufacturer's recommended protocol. Low-quality reads were removed using the Cutadapt program (Martin 2011). In order to analyze mRNA from tea leaves, valid reads were aligned to the reference genome of C. sinensis var. sinensis (http://tpdb.shengxin.ren/index.html) using HISAT2 software (https://daehwankimlab.github.io/hisat2/) (Kim et al. 2015). The mapped reads of each treatment were assembled using StringTie (http://ccb.jhu.edu/software/stringtie/) (version: stringtie-1.3.4d. Linux_x86_64). Then, the expression abundance of mRNAs was calculated as the fragment per kilobase of transcripts per million mapped reads (FPKM) value, using StringTie software and Ballgown in R (Frazee et al. 2015; Pertea et al. 2015). The differentially expressed mRNAs (DEmRNAs) were identified by the criteria of log2 (fold change) > 1 or log2 (fold change) < −1, and statistical significance (P value < 0.05), using the R package Empirical Analysis of Digital Gene Expression Data in R (edgeR) (version 3.13) (https://bioconductor.org/packages/release/bioc/html/edgeR.html). Functional analysis of the DEmRNAs was conducted by reference to the databases of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG), using the Blast2GO software (Conesa et al. 2005).In order to analyze lncRNA from tea leaves, sequence quality was further verified using FastQC (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Valid reads were used to map reads to the genome of C. sinensis var. sinensis (http://tpdb.shengxin.ren/index.html) using Bowtie2 in R and TopHat2 (Kim et al. 2013; Langmead and Salzberg 2012). The mapped reads of each sample were assembled using StringTie software (Pertea et al. 2015); then, the transcriptomes from the six samples from the infected and uninfected treatments were merged to construct a comprehensive transcriptome using an in-house Perl program (https://www.perl.org). After the final transcriptome was generated, the abundance of each of the transcripts was estimated using StringTie software and Ballgown (Frazee et al. 2015; Pertea et al. 2015). All transcripts were filtered to screen out the lncRNAs. First, transcripts that overlapped with known mRNAs and novel mRNAs with coding potential were filtered out. Second, to obtain a wider range of candidate lncRNAs, we only adopted the Coding Potential Calculator (CPC) software (version 0.9-r2) and Coding-Non-Coding Index (CNCI) software (version 2.0) for prediction in this study (Kong et al. 2007; Sun et al. 2013). Any transcript with a CPC score ≤ −1 and a CNCI score < 0 were removed (Kong et al. 2007; Sun et al. 2013). Therefore, based on such prediction methods, only candidate lncRNAs were obtained, and the predicted results are not completely lncRNAs. Finally, transcripts that were shorter than 200 bp and where the coverage of reads was < 3 were removed. The remaining transcripts were considered to be candidate lncRNAs. The abundance level of lncRNAs was calculated using StringTie software by calculating FPKM (Pertea et al. 2015; Trapnell et al. 2010) and the differentially expressed (DE) lncRNAs (DElncRNAs) were selected on the criteria of log2 (fold change) > 1 or log2 (fold change) < −1, and with statistical significance (P value < 0.05) by the R package Empirical Analysis of Digital Gene Expression Data in R (edgeR) (version 3.13) (https://bioconductor.org/packages/release/bioc/html/edgeR.html) (Frazee et al. 2015; Pertea et al. 2015). Cis target genes of lncRNAs in the region within 100,000 bp upstream and downstream were predicted by the Python program (https://www.python.org). Then, the functional annotation and enrichment analysis of the target genes identified for individual DElncRNAs was conducted using Blast2GO (Conesa et al. 2005).Infection by E. nigrum significantly up- or downregulated expression of 430 or 128, respectively, of the DEmRNAs originating from the infected tea leaf (Table 1). When the GO database annotation of all DEmRNAs was carried out, the total number of genes (TB gene number) and the number of genes with significantly differential expression between the infected treatment and the uninfected treatment (TS gene number) were 13,594 and 256, respectively. At the biological process (BP) aspect, there was enrichment with particular respect to the genes related to response to biotic stimulus, with the the number of genes with significantly differential expression annotated in the specified database (S gene number) and the number of genes annotated in the specified term (B gene number) being 3 and 9, respectively. The abundance of the DEmRNAs for zeamatin (TEA032349.1, TEA005276.1, TEA014254.1, TEA026507.1, TEA014256.1, and TEA014245.1) encoding an antifungal protein was upregulated in infected leaves relative to uninfected leaves.Table 1. Characteristics of messenger RNA (mRNA) and long noncoding RNA (lncRNA) sequences from tea (Camellia sinensis var. sinensis) leaves during infection by Epicoccum nigrumSample sourceaSequenceUninfected treatmentInfected treatmentlncRNAs and mRNAsRaw reads89,868,99193,224,757Valid reads58,544,88662,754,655Valid bases8.78 G9.42 GmRNAsInfected treatment versus−430 (upregulation)Uninfected treatment−128 (downregulation)lncRNAsInfected treatment versus−128 (upregulation)Uninfected treatment−931 (downregulation)Cis regulation between DElncRNAs andb−1,059 DElncRNAs/3 Cis regulation (diff)DEmRNAs−986 DEmRNAs/3 Cis regulation (diff)aUninfected and infected treatments represent a tea leaf that was punctured with four holes using sterilized syringe needles, with the four holes being clustered on the adaxial surface in the range of 6 mm in diameter, after which sterile water (uninfected) or the mycelial plug (6 mm in diameter per inoculation site) (infected) was applied to each inoculation site on the tea leaves.bCis regulation between the differentially expressed lncRNAs (DElncRNAs) and the differentially expressed mRNAs (DEmRNAs).Table 1. Characteristics of messenger RNA (mRNA) and long noncoding RNA (lncRNA) sequences from tea (Camellia sinensis var. sinensis) leaves during infection by Epicoccum nigrumView as image HTML Infection by E. nigrum significantly up- and downregulated expression of 128 and 931 DElncRNAs, respectively. For response to infection by E. nigrum, 1,059 DElncRNAs exhibited cis regulation of expression of 986 DEmRNAs (Table 1). GO enrichment analysis of DEmRNAs targeted by lncRNAs indicated that the TB and TS gene numbers were 545 and 25, respectively. At the level of BP, the DEmRNAs were mainly enriched with respect to response to biotic stimulus, with the S gene number being 1. KEGG analysis indicated that the DEmRNA for golgin subfamily B member 1-like (TEA003533.1), enriched in the pathway of nucleotide excision repair, was cis regulated by the DElncRNA MSTRG.37064.1. The DEmRNA histone acetyltransferase type B catalytic subunit (TEA005332.1), unenriched in the database of the GO and KEGG, was cis regulated by the DElncRNA MSTRG.50640.1. The DEmRNA uncharacterized LOC103646289 (TEA019367.1), enriched in the pathway of ribosome, was cis regulated by the DElncRNA MSTRG.32067.2.For tea leaf spot caused by E. nigrum, mRNAs and lncRNAs originating from tea leaves during the infection were sequenced, and the DEmRNAs and DElncRNAs were annotated and predicted the interaction. This dataset will provide an important resource for studying the mechanism underlying the interaction between E. nigrum and tea, and association studies will provide a valuable resource for future research into the host–pathogen interaction.Data AvailabilityThe sequences of mRNAs and lncRNAs of tea have been deposited in GenBank under Sequence Read Archive accession PRJNA806359. In addition, bioinformatics analysis of mRNAs and lncRNAs are presented in Supplementary Tables S1 through S34 and Supplementary Figure S1. E. nigrum strain CGMCC3.20484 has been deposited in the China General Microbiological Culture Collection.The author(s) declare no conflict of interest.Literature CitedBao, X.-T., Dharmasena, D. S. P., Li, D.-X., Wang, X., Jiang, S.-L., Ren, Y.-F., Wang, D.-L., Song, B.-A., and Chen, Z. 2019. First report of Epicoccum sorghinum causing leaf spot on tea in China. Plant Dis. 103:3282. https://doi.org/10.1094/PDIS-06-19-1296-PDN Google ScholarBatista, P. J., and Chang, H. Y. 2013. Long noncoding RNAs: Cellular address codes in development and disease. Cell 152:1298-1307. https://doi.org/10.1016/j.cell.2013.02.012 Crossref, Google ScholarChen, Q., Hou, L.-W., Duan, W.-J., Crous, P.-W., and Cai, L. 2017. Didymellaceae revisited. Stud. Mycol. 87:105-159. https://doi.org/10.1016/j.simyco.2017.06.002 Crossref, Google ScholarChi, Y.-J., Wang, D., Wang, J.-P., Yu, W.-D., and Yang, J.-C. 2019. Long non-coding RNA in the pathogenesis of cancers. Cells 8:1015. https://doi.org/10.3390/cells8091015 Crossref, Google ScholarConesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., and Robles, M. 2005. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21:3674-3676. https://doi.org/10.1093/bioinformatics/bti610 Crossref, Google ScholarDeepika, Y. S., Mahadevakumar, S., Amruthesh, K. N., and Lakshmidevi, N. 2021. First report of Epicoccum nigrum associated with leaf spot disease of cowpea (Vigna unguiculata) from India. J. Plant Pathol. 103:391-392. https://doi.org/10.1007/s42161-020-00718-9 Crossref, Google ScholarDykes, I. M., and Emanueli, C. 2017. Transcriptional and post-transcriptional gene regulation by long non-coding RNA. Genomics Proteomics Bioinf. 15:177-186. https://doi.org/10.1016/j.gpb.2016.12.005 Crossref, Google ScholarFrazee, A. C., Pertea, G., Jaffe, A. E., Langmead, B., Salzberg, S. L., and Leek, J. T. 2015. Ballgown bridges the gap between transcriptome assembly and expression analysis. Nat. Biotechnol. 33:243-246. https://doi.org/10.1038/nbt.3172 Crossref, Google ScholarHeo, J. B., Lee, Y. S., and Sung, S. 2013. Epigenetic regulation by long noncoding RNAs in plants. Chromosome Res. 21:685-693. https://doi.org/10.1007/s10577-013-9392-6 Crossref, Google ScholarKim, D., Langmead, B., and Salzberg, S. L. 2015. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 12:357-360. https://doi.org/10.1038/nmeth.3317 Crossref, Google ScholarKim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. 2013. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14:R36. https://doi.org/10.1186/gb-2013-14-4-r36 Crossref, Google ScholarKong, L., Zhang, Y., Ye, Z.-Q., Liu, X.-Q., Zhao, S.-Q., Wei, L., and Gao, G. 2007. CPC: Assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 35:W345-W349. https://doi.org/10.1093/nar/gkm391 Crossref, Google ScholarLangmead, B., and Salzberg, S. L. 2012. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9:357-359. https://doi.org/10.1038/nmeth.1923 Crossref, Google ScholarLin, S., Taylor, N. J., and Peduto Hand, F. 2018. Identification and characterization of fungal pathogens causing fruit rot of deciduous holly. Plant Dis. 102:2430-2445. https://doi.org/10.1094/PDIS-02-18-0372-RE Link, Google ScholarMahadevakumar, S., Jayaramaiah, K. M., and Janardhana, G. R. 2014. First report of leaf spot disease caused by Epicoccum nigrum on lablab purpureus in India. Plant Dis. 98:284. https://doi.org/10.1094/PDIS-07-13-0798-PDN Link, Google ScholarMartin, M. 2011. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 17:10-12. https://doi.org/10.14806/ej.17.1.200 Crossref, Google ScholarMatzke, M. A., and Mosher, R. A. 2014. RNA-directed DNA methylation: An epigenetic pathway of increasing complexity. Nat. Rev. Genet. 15:394-408. https://doi.org/10.1038/nrg3683 Crossref, Google ScholarNie, J.-T., Wang, H.-S., Zhang, W.-L., Teng, X., Yu, C., Cai, R., and Wu, G. 2021. Characterization of lncRNAs and mRNAs involved in powdery mildew resistance in Cucumber. Phytopathology 111:1613-1624. https://doi.org/10.1094/PHYTO-11-20-0521-R Link, Google ScholarPertea, M., Pertea, G. M., Antonescu, C. M., Chang, T. C., Mendell, J. T., and Salzberg, S. L. 2015. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33:290-295. https://doi.org/10.1038/nbt.3122 Crossref, Google ScholarSong, L., Fang, Y., Chen, L., Wang, J., and Chen, X.-W. 2021. Role of non-coding RNAs in plant immunity. Plant Commun. 2:100180. https://doi.org/10.1016/j.xplc.2021.100180 Crossref, Google ScholarSun, L., Luo, H.-T., Bu, D.-C., Zhao, G.-G., Yu, K.-T., Zhang, C.-H., Liu, Y.-N., Chen, R.-S., and Zhao, Y. 2013. Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts. Nucleic Acids Res. 41:e166. https://doi.org/10.1093/nar/gkt646 Crossref, Google ScholarTrapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Van Baren, M. J., Salzberg, S. L., Wold, B. J., and Pachter, L. 2010. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28:511-515. https://doi.org/10.1038/nbt.1621 Crossref, Google ScholarVance, K. W., and Ponting, C. P. 2014. Transcriptional regulatory functions of nuclear long noncoding RNAs. Trends Genet. 30:348-355. https://doi.org/10.1016/j.tig.2014.06.001 Crossref, Google ScholarWang, G.-N., Wang, X.-F., Zhang, Y., Yang, J., Li, Z.-K., Wu, L.-Z., Wu, J.-H., Wu, N., Liu, L.-X., Liu, Z.-W., Zhang, M., Wu, L.-Q., Zhang, G.-Y., and Ma, Z.-Y. 2021. 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. https://doi.org/10.1186/s12870-021-02835-8 Crossref, Google ScholarWang, X., Yin, Q.-X., Jiang, S.-L., Wu, X., Wang, D.-L., Song, B.-A., and Chen, Z. 2020. First report of Didymella bellidis causing tea leaf spot in China. Plant Dis. 104:1254. https://doi.org/10.1094/PDIS-08-19-1782-PDN Link, Google ScholarWierzbicki, A. T., Haag, J. R., and Pikaard, C. S. 2008. Noncoding transcription by RNA polymerase Pol IVb/Pol V mediates transcriptional silencing of overlapping and adjacent genes. Cell 135:635-648. https://doi.org/10.1016/j.cell.2008.09.035 Crossref, Google ScholarWu, D., Zhang, D.-H., Timko, M. P., Li, M.-Y., and Liang, G.-L. 2017. First report of Epicoccum nigrum causing brown leaf spot of loquat in Southwestern China. Plant Dis. 101:1553. https://doi.org/10.1094/PDIS-12-16-1840-PDN Link, Google ScholarYin, Q.-X., Jiang, S.-L., Li, D.-X., Huang, H.-L., Wang, Y., Wang, D.-L., and Chen, Z. 2022. First report of Epicoccum nigrum causing brown leaf spot in tea in Guizhou Province. China. Plant Dis. 106:321. https://doi.org/10.1094/PDIS-04-21-0815-PDN Link, Google ScholarZaynab, M., Fatima, M., Abbas, S., Umair, M., Sharif, Y., and Raza, M. A. 2018. Long non-coding RNAs as molecular players in plant defense against pathogens. Microb. Pathog. 121:277-282. https://doi.org/10.1016/j.micpath.2018.05.050 Crossref, Google ScholarZhang, H., Guo, H., Hu, W.-G., and Ji, W.-Q. 2020. The emerging role of long non-coding RNAs in plant defense against fungal stress. Int. J. Mol. Sci. 21:2659. https://doi.org/10.3390/ijms21082659 Crossref, Google ScholarZhao, X.-Z., Wang, Y., Li, D.-X., Ren, Y.-F., and Chen, Z. 2018. Morphological characterization and phylogenetic analysis of the pathogen Phoma segeticola var. camelliae causing a new tea disease. Acta Phytopathol. Sin. 48:556-559. https://doi.org/10.13926/j.cnki.apps.000184 Google ScholarFunding: This work was supported by the National Key Research Development Program of China (2017YFD0200308) and its post-subsidy project (2018-5262), National Natural Science Foundation of China (numbers 21977023 and 31860515), China Agriculture Research System of MOF and MARA (D09), Program of Introducing Talents to Chinese Universities (111 Program, D20023), and Guizhou Provincial Science and Technology Program (2020-5001).The author(s) declare no conflict of interest. Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.DetailsFiguresLiterature CitedRelated Vol. 2, No. 4 2022ISSN:2690-5442 Download Metrics Article History Issue Date: 23 Dec 2022Published: 26 Jul 2022Accepted: 19 Apr 2022 Pages: 404-407 InformationCopyright © 2022 The Author(s).This is an open access article distributed under the CC BY-NC-ND 4.0 International license.Funding National Key Research Development Program of ChinaGrant/Award Number: 2017YFD0200308Grant/Award Number: 2018-5262 National Natural Science Foundation of ChinaGrant/Award Number: 21977023Grant/Award Number: 31860515 China Agriculture Research System of MOF and MARAGrant/Award Number: D09 Program of Introducing Talents to Chinese UniversitiesGrant/Award Number: D20023 Guizhou Provincial Science and Technology ProgramGrant/Award Number: 2020-5001 Keywordsassociation studies of mRNA and lncRNAcis regulation of gene expressionleaf spotlncRNAmRNAThe author(s) declare no conflict of interest.PDF download