Abstract:HomeMolecular Plant-Microbe Interactions®Vol. 35, No. 9Sequence Data, Functional Annotation, and Relationship Analysis Between mRNAs and Long Noncoding RNAs from Tea Leaves During Infection by the Fungal Pathogen Epicoccum sorghinum Previous RESOURCE ANNOUNCEMENT OPENOpen Access licenseSequence Data, Functional Annotation, and Relationship Analysis Between mRNAs and Long Noncoding RNAs from Tea Leaves During Infection by the Fungal Pathogen Epicoccum sorghinumChen Huang, Hongke Huang, Zhongqiu Xia, Yuqin Yang, Xinyue Jiang, Yuanyou Yang, Delu Wang, Xiangyang Li, and Zhuo ChenChen HuangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, Hongke 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, ChinaSearch for more papers by this author, 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, ChinaSearch for more papers by this author, 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, ChinaSearch for more papers by this author, Xinyue JiangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, Yuanyou YangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, Delu WangCollege of Forestry, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, Xiangyang Li†Corresponding authors: X. Li; E-mail Address: xyli1@gzu.edu.cn, and Z. Chen; E-mail Address: gychenzhuo@aliyun.comKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, and Zhuo Chen†Corresponding authors: X. Li; E-mail Address: xyli1@gzu.edu.cn, and Z. Chen; E-mail Address: gychenzhuo@aliyun.comhttp://orcid.org/0000-0001-7130-8457Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this authorAffiliationsAuthors and Affiliations Chen Huang1 Hongke Huang1 2 Zhongqiu Xia1 2 Yuqin Yang1 2 Xinyue Jiang1 Yuanyou Yang1 Delu Wang3 Xiangyang Li1 † 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:15 Sep 2022https://doi.org/10.1094/MPMI-02-22-0048-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articleResource AnnouncementEpicoccum is one of the most species-rich genera in the family Didymellaceae (Q. Chen et al. 2017). Epicoccum sorghinum can cause leaf spot on wheat, tobacco, Oxalis debilis, lily, Bletilla striata, taro, Paris polyphylla, Brassica parachinensis, white chrysanthemum, rice, Atractylodes lancea, and Hemerocallis citrina (Q. H. Chen et al. 2021; X. L. Chen et al. 2017; Du et al. 2020; Fu et al. 2019; Liu et al. 2018; Liu et al. 2020; Ma et al. 2021; Perelló and Moreno 2005; Yu et al. 2019; Yuan et al. 2016; Zeng et al. 2018; Zhou et al. 2018). In addition, E. sorghinum was identified as the pathogen of leaf brown spot of Digitaria sanguinalis (Kang et al. 2019), leaf blight of Kadsura coccinea (Fan et al. 2020), soft rot of Dendrobium officinale (Xiao et al. 2021), lavender blackleg disease (Gu et al. 2021), and leaf sheath and leaf spot of maize (T.-M. Chen et al. 2021; Do Amaral et al. 2004). E. sorghinum and other fungi form a fungal complex, which can then cause low seed germination of Paspalum guenoarum (Gasparetto et al. 2017). E. sorghinum can cause sorghum grain mold, alone or as a member of another complex of fungi (Navi et al. 2005; Oliveira et al. 2018). Tea leaf spot caused by E. sorghinum results in huge losses of tea leaf productivity as well as decreases in tea leaf quality in tea plantations in Guizhou Province, China (Bao et al. 2019). Because of the complicated pathogenesis and the role of epidemiology in Epicoccum tea leaf spot, there are currently no effective control measures (Huang et al. 2021).Long noncoding RNAs (lncRNAs) are a number of non-protein-coding RNAs more than 200 nucleotides in length, which can function in transcriptional or posttranscriptional epigenetic regulation of gene expression (Chi et al. 2019). Analysis of lncRNAs can contribute to our understanding of the mechanisms of infection and resistance (Zaynab et al. 2018). A number of studies have indicated that lncRNAs are involved in plant–pathogen interactions. For instance, wheat lncRNAs participated in the response of wheat to infection by powdery mildew (Blumeria graminis f. sp. tritici) and stripe rust (Puccinia striiformis f. sp. tritici) (Xin et al. 2011). lncRNA ELF18-INDUCED LONG-NONCODING RNA1 (ELENA1) can interact with Mediator subunit 19a (MED19a), then affect the enrichment of MED19a on the PATHOGENESIS-RELATED PROTEIN 1 (PR1) gene promoter. So, the upregulation of the PR1 gene contributes to resistance to Pseudomonas syringae pv. tomato in Arabidopsis thaliana (Seo et al. 2017). The lncRNA GhlncLOX3 can modulate the expression of the cotton gene GhLOX3 involved in jasmonic acid biosynthesis and enhance the resistance of Gossypium hirsutum to Verticillium dahliae (Wang et al. 2021). In a previous study, we sequenced messenger RNAs (mRNAs) and lncRNAs originating from tea leaves during infection by the fungal pathogen Lasiodiplodia theobromae and predicted some cis regulation relationships between lncRNAs and mRNAs (Li et al. 2021). Nevertheless, no research has been published on the regulatory relationship between mRNAs and lncRNAs originating from tea leaves challenged by E. sorghinum, the pathogen of tea leaf spot. In this study, we sequenced and carried out functional annotation of mRNAs and lncRNAs originating from tea leaves during infection by E. sorghinum, and predicted the interaction relationship between putative target mRNAs and lncRNAs. These data will provide a resource for future research into host–pathogen interactions and disease resistance.Three-year-old potted tea plants (Camellia sinensis var. sinensis ‘Fuding-dabaicha’) were grown in the greenhouse under temperatures of 25°C (day) and 20°C (night), and a relative humidity of 70 to 80%. Single-spore E. sorghinum strain CGMCC3.20150 was inoculated onto potato dextrose agar (PDA) and incubated at 28°C in the dark to form colonies after 20 days. The third or fourth leaf on each detached, healthy tea twig was physically wounded, then inoculated with either a conidial suspension (106 conidia/ml, 20 to 30 μl of suspension per inoculation site) or mycelial plugs (6 mm in diameter). The adaxial leaf surface on either side of the main vein was punctured with four holes, using sterilized syringe needles, with the four holes being clustered in the range of 6 mm in diameter. For the uninfected treatment, sterile water or uninoculated PDA plugs were used to inoculate the tea leaves. For each of the infected and uninfected treatments, 15 healthy twigs were pooled to represent one replicate, and the experiment was conducted with three biological replicates for each treatment. At 48 h postinoculation, a lesion formed on each tea leaf site inoculated with a mycelial plug. The total RNA of the tea leaf was extracted at 48 h postinoculation, using a TRIzol Plus RNA Purification Kit (Invitrogen, Carlsbad, CA, U.S.A.), following the manufacturer's procedure. A cDNA library was prepared using a TruSeq Stranded mRNA Sample Preparation Kit (Illumina, San Diego, CA, U.S.A.), and sequenced on an Illumina NovaSeq 6000 platform 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). For mRNA analysis, valid data were aligned to the reference genome of C. sinensis var. sinensis (available from The Tea Plant Information Archive) using HISAT software (version 2.0.4) (Kim et al. 2015). The mapped reads of each sample were assembled, then expression was quantified in terms of the abundance of mRNAs by calculating fragments per kilobase of transcripts per million mapped read values at the level of log2 (fold change) > 1 or log2 (fold change) < −1 and with statistical significance (P value < 0.05), using StringTie software (version:stringtie-1.3.4d) (Frazee et al. 2015; Pertea et al. 2015) and Ballgown in R (Frazee et al. 2015). Functional analysis of the differentially expressed (DE) mRNAs was conducted by the gene ontology (GO) database (version 2016.04) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (version 2016.05), using the Blast2GO bioinformatics tool (Conesa et al. 2005).For lncRNA analysis, sequence quality was further verified using FastQC. Valid reads were used to map reads to the genome of C. sinensis var. sinensis, using Bowtie2 in R (Langmead and Salzberg 2012) and TopHat2 (Kim et al. 2013). The mapped reads of each sample were assembled using StringTie software (version stringtie-1.3.4d) (Pertea et al. 2015); then, the transcriptomes from the six samples consisting of the infected treatment and the uninfected treatment (three biological replicates of each) were merged to construct a comprehensive transcriptome using the Perl program. After the final transcriptome was generated, the abundance of all transcripts was estimated using StringTie software (Pertea et al. 2015) and Ballgown (Frazee et al. 2015). All transcripts were filtered as follows. First, transcripts that overlapped with known mRNAs or novel mRNAs with coding potential were filtered out; second, transcripts that were shorter than 200 bp and the coverage of reads < 3 were removed. Transcripts with coding potential were predicted using the Coding Potential Calculator (CPC) software (version 0.9-r2) (Kong et al. 2007) and Coding-Non-Coding Index (CNCI) software (version 2.0) (Sun et al. 2013). All transcripts with CPC score < −1 and CNCI score < 0 were removed. The remaining transcripts were considered to be lncRNAs. The expression level of lncRNAs was calculated using StringTie software (Pertea et al. 2015; Trapnell et al. 2010) and Ballgown in R (Frazee et al. 2015). The cis target genes of lncRNAs in the 100,000-bp upstream and downstream regions were predicted by the Python program. Then, functional analysis of the target genes of lncRNAs was conducted by using Blast2GO (Conesa et al. 2005). Significance was expressed as a P value < 0.05.Infection by E. sorghinum significantly up- or downregulated expression of 720 or 321 of the tea leaf mRNAs, respectively. When the GO database annotation of all mRNAs was carried out, the total number of genes (TB gene number) and the number of genes with significantly different expression between the infected treatment and the uninfected treatment (TS gene number) were 1,3585 and 523, respectively. At the level of biological process in GO, there was enrichment with particular respect to the number of genes related to defense response to fungal pathogens, with the number of genes with significantly different expression annotated in the specified database (S gene number) and the number of genes annotated in the specified GO term (B gene number) being 3 and 19, respectively. Expression of the genes Mcht-2 (TEA013240.1) and peroxidase 12 (TEA015011.1) was upregulated in the infected treatment relative to the uninfected treatment, whereas that of LRR receptor-like serine/threonine-protein kinase ERECTA (TEA016781.1) and long chain acyl-CoA synthetase 2 (TEA000885.1) was downregulated in the infected treatment relative to the uninfected treatment. KEGG enrichment analysis of DE mRNAs indicated that TB gene number/TS gene number was 16,706/605. Sixty S genes and 1,945 B genes were enriched in the pathway of plant–pathogen interaction (Table 1).Table 1. Messenger RNA (mRNA) and long noncoding RNA (lncRNA) sequences from tea (Camellia sinensis var. sinensis) leaves during infection by the fungal pathogen Epicoccum sorghinumTreatment sample sourceaSequenceUninfectedbInfectedcComparisonsMRNAs Raw reads52,155,31753,824,580− Valid reads49,393,45950,831,946− Valid bases7.41 G7.63 G− Infected vs. uninfected−−720 Upregulation, 321 downregulationlncRNAs Raw reads93,881,11691,514,902− Valid reads61,051,80565,195,740− Valid bases9.16 G9.78 G− Infected vs. uninfected−−242 Upregulation, 457 downregulationcis regulationd−−699 DE lncRNAs per 6 cis regulation−−1,165 DE mRNAs per 6 cis regulationaTea leaves were punctured with four holes using sterilized syringe needles, with the four holes being clustered within the range of 6 mm in diameter.bUninoculated potato dextrose agar (PDA) plugs were placed on the leaves.cMycelial PDA plugs were placed on the leaves.dcis regulation between differentially expressed (DE) lncRNAs and mRNAs (significant difference).Table 1. Messenger RNA (mRNA) and long noncoding RNA (lncRNA) sequences from tea (Camellia sinensis var. sinensis) leaves during infection by the fungal pathogen Epicoccum sorghinumView as image HTML The DE lncRNAs were analyzed in terms of gene expression. Infection by E. sorghinum significantly up- or downregulated expression of 242 or 457 DE lncRNAs, respectively. The DE mRNAs from the genes uncharacterized protein (TEA015043.1), cytochrome P450 CYP72A219-like (TEA018674.1), GLP1 (TEA005895.1), uncharacterized protein (TEA025952.1), hydroxymethylglutaryl-CoA synthase-like (TEA001661.1), and adenine phosphoribosyltransferase-like (TEA015523.1) occurred as a result of cis regulation by the DE lncRNAs MSTRG.63686.1 (TEA015043.1), MSTRG.23615.2 (TEA018674.1), MSTRG.38325.1 (TEA005895.1), MSTRG.38495.1 (TEA025952.1), MSTRG.26049.1 (TEA001661.1), and MSTRG.20281.2 (TEA015523.1), respectively. KEGG enrichment analysis of DE mRNAs targeted by lncRNAs indicated that the TB gene number and TS gene number were 693 and 23, respectively. At the level of plant–pathogen interaction, the S gene number and the B gene number of DE mRNAs were 4 and 86, respectively (Table 1).For tea leaf spot caused by E. sorghinum, mRNAs and lncRNAs originating from tea leaves during the infection of the pathogen of tea leaf spot (E. sorghinum) were sequenced. DE mRNAs and DE lncRNAs were analyzed and annotated using the GO and KEGG databases, and cis target genes of DE lncRNAs were also predicted. This dataset will provide an important resource for studying the mechanism underlying the interaction between E. sorghinum and the tea plant.Data AvailabilityThe sequences of mRNAs and lncRNAs of tea have been deposited in GenBank under Sequence Read Archive accession PRJNA799849. In addition, bioinformatics analysis of mRNAs and lncRNAs are presented in the Supplementary Figure and Supplementary Tables.Author-Recommended Internet ResourcesBlast2GO: http://www.ncbi.nlm.nih.gov/blastFastQC: http://www.bioinformatics.babraham.ac.uk/projects/fastqc/Gene ontology: http://geneontology.org HISAT software: http://www.ccb.jhu.edu/software/hisat/index.shtmlKyoto Encyclopedia of Genes and Genomes: http://www.kegg.jp/keggPerl program: https://www.perl.orgPython program: https://www.python.orgStringTie software: https://ccb.jhu.edu/software/stringtie/ The Tea Plant Information Archive: http://tpdb.shengxin.ren/index.html 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 Link, ISI, 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, Medline, ISI, Google ScholarChen, Q. H., Li, J. X., Qi, Y., Liu, D. H., and Miao, Y. H. 2021. First report of leaf spot on white chrysanthemum (Chrysanthemum morifolium) caused by Epicoccum sorghinum in Hubei Province, China. Plant Dis. 105:1212. https://doi.org/10.1094/PDIS-09-20-1896-PDN Link, ISI, Google ScholarChen, T.-M., Xie, Y. J., Sun, Q., Shi, X.-C., Wang, S.-Y., and Laborda, P. 2021. First report of Epicoccum sorghinum causing leaf sheath and leaf spot on maize in China. Plant Dis. 105:3741. https://doi.org/10.1094/PDIS-04-21-0746-PDN Abstract, ISI, Google ScholarChen, X. L., Wang, Y. H., and Luo, T. 2017. First report of leaf spot caused by Phoma sorghina on Oxalis debilis in China. Plant Dis. 101:1047. https://doi.org/10.1094/PDIS-11-16-1614-PDN Link, ISI, 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, ISI, 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, Medline, ISI, Google ScholarDo Amaral, A. L., De Carli, M. L., Barbosa Neto, J. F., and Dal Soglio, F. K. 2004. Phoma sorghina, a new pathogen associated with Phaeosphaeria leaf spot on maize in Brazil. Plant Pathol. 53:259. Crossref, ISI, Google ScholarDu, Y.-X., Xu, Y., Jiang, J.-Y., Wan, X.-F., Kang, C.-Z., Wang, T.-L., and Guo, L.-P. 2020. First report of Epicoccum sorghinum causing leaf spot on Atractylodes lancea in Hubei Province, China. Plant Dis. 104:3082. https://doi.org/10.1094/PDIS-03-20-0672-PDN Link, ISI, Google ScholarFan, H.-Y., Chen, W.-M., Huang, Y.-C., Zeng, Q.-Q., and Huang, Y. 2020. First report of Epicoccum sorghinum causing leaf blight on Kadsura coccinea in China. Plant Dis. 104:3258. https://doi.org/10.1094/PDIS-04-20-0884-PDN Link, ISI, 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, Medline, ISI, Google ScholarFu, R., Chen, C., Wang, J., Lu, D., and Gong, X. 2019. First report of Epicoccum sorghinum causing leaf spot on Paris polyphylla in China. Plant Dis. 103:1426. https://doi.org/10.1094/PDIS-11-18-2076-PDN Link, ISI, Google ScholarGasparetto, B. F., Franke, L. B., Andrade, C. C. L., Dalbosco, M., Duarte, V., Moreira, S. I., and Alves, E. 2017. First report of Bipolaris micropus, Curvularia geniculata, Epicoccum sorghinum, and Fusarium incarnatum on Paspalum guenoarum seeds in Rio Grande do Sul, Brazil. Plant Dis. 101:1679. https://doi.org/10.1094/PDIS-12-16-1746-PDN Link, ISI, Google ScholarGu, C.-Y., Pan, R., Abid, M., Zang, H.-Y., Yang, X., and Chen, Y. 2021. First report of blackleg disease caused by Epicoccum sorghinum on lavender (Lavandula stoechas) in China. Plant Dis. 105:2733. https://doi.org/10.1094/PDIS-10-20-2280-PDN Link, ISI, Google ScholarHuang, H.-L., Yin, Q.-X., Jiang, S.-L., Bao, X.-T., Wu, X., Wang, D.-L., and Chen, Z. 2021. Identification and biological characteristics of Epicoccum sorghinum, pathogen of tea leaf spot. Chin. J. Trop. Crops 42:3269-3277. Google ScholarKang, Y., Zhang, J. Wan, X.-Q., Xu, T.-T., Li, C.-X., and Cao, H.-Q. 2019. First report of leaf brown spot caused by Epicoccum sorghinum on Digitaria sanguinalis in China. Plant Dis. 103:1787. https://doi.org/10.1094/PDIS-02-19-0332-PDN Link, ISI, 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, Medline, ISI, 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, Medline, ISI, 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, Medline, ISI, 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, Medline, ISI, Google ScholarLi, D.-X., Jiang, S.-L., Wen, X.-D., Song, X.-C., Yang, Y.-Y., Wang, Y., Huang, H.-L., Wang, D.-L., and Chen, Z. 2021. Sequencing and functional annotation of mRNAs and lncRNAs from tea (Camellia sinensis) leaves during infection by the fungal pathogen, Lasiodiplodia theobromae. PhytoFrontiers 1:364-367. https://doi.org/10.1094/PHYTOFR-03-21-0020-A Link, Google ScholarLiu, L.-M., Zhao, Y., Zhang, Y.-L., Wang, L., Hou, Y.-X., and Huang, S.-W. 2020. First report of leaf spot disease on rice caused by Epicoccum sorghinum in China. Plant Dis. 104:2735. https://doi.org/10.1094/PDIS-03-20-0488-PDN Link, ISI, Google ScholarLiu, P.-Q., Wei, M.-Y., Zhu, L., Wang, R.-B., Li, B.-J., Weng, Q.-Y., and Chen, Q.-H. 2018. First report of leaf spot on taro caused by Epicoccum sorghinum in China. Plant Dis. 102:682. https://doi.org/10.1094/PDIS-11-16-1621-PDN Link, ISI, Google ScholarMa, Y.-M., Zhou, J.-L., Hu, Z., Zhong, J., and Zhu, J.-Z. 2021. First report of Epicoccum sorghinum causing leaf spot on Hemerocallis citrina in China. Plant Dis. 105:2251. https://doi.org/10.1094/PDIS-10-20-2144-PDN Link, ISI, 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 ScholarNavi, S. S., Bandyopadhyay, R., Reddy, R. K., Thakur, R. P., and Yang, X.-B. 2005. Effects of wetness duration and grain development stages on sorghum grain mold infection. Plant Dis. 89:872-878. https://doi.org/10.1094/PD-89-0872 Link, ISI, Google ScholarOliveira, R. C., Carnielli-Queiroz, L., and Correa, B. 2018. Epicoccum sorghinum in food: Occurrence, genetic aspects and tenuazonic acid production. Curr. Opin. Food Sci. 23:44-48. https://doi.org/10.1016/j.cofs.2018.05.011 Crossref, ISI, Google ScholarPerelló, A. E., and Moreno, M. V. 2005. First report of Phoma sorghina (Sacc.) Boerema Dorenbosch & van Kest on wheat leaves (Triticum aestivum L.) in Argentina. Mycopathologia 159:75-78. Crossref, Medline, ISI, 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, Medline, ISI, Google ScholarSeo, J. S., Sun, H. X., Park, B. S., Huang, C. H., Yeh, S. D., Jung, C., and Chua, N. H. 2017. ELF18-INDUCED LONG-NONCODING RNA associates with mediator to enhance expression of innate immune response genes in Arabidopsis. Plant Cell 29:1024-1038. https://doi.org/10.1105/tpc.16.00886 Crossref, Medline, ISI, 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, Medline, ISI, 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, Medline, ISI, 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, Medline, ISI, Google ScholarXiao, C.-Y., Li, R.-Y., Song, X.-C., Tian, X.-J., and Zhao, Q.-J. 2021. First report of soft rot on Dendrobium officinale caused by Epicoccum sorghinum in China. Plant Dis. 106:1063. https://doi.org/10.1094/PDIS-04-21-0866-PDN Link, ISI, Google ScholarXin, M.-M., Wang, Y, Yao, Y.-Y., Song, N, Hu, Z.-R., Qin, D.-D., Xie, C.-J., Peng, H.-R., Ni, Z.-F., and Sun, Q.-X. 2011. 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. 11:61. https://doi.org/10.1186/1471-2229-11-61 Crossref, Medline, ISI, Google ScholarYu, L., She, X.-M., Lan, G.-B., Tang, Y.-F., Li, Z.-G., Deng, M.-G., and He, Z.-F. 2019. First report of leaf spot caused by Epicoccum sorghinum on Chinese flowering cabbage (Brassica parachinensis) in China. Plant Dis. 103:2966. https://doi.org/10.1094/PDIS-04-19-0881-PDN Link, ISI, Google ScholarYuan, G.-Q., Liao, T., Tan, H.-W., Li, Q.-Q., and Lin, W. 2016. First report of leaf spot caused by Phoma sorghina on tobacco in China. Plant Dis. 100:1790. https://doi.org/10.1094/PDIS-11-15-1377-PDN Link, ISI, 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, Medline, ISI, Google ScholarZeng, H., Lu, Q., and Li, R. 2018. First report of leaf spot of lily caused by Epicoccum sorghinum in China. Plant Dis. 102:2648. https://doi.org/10.1094/PDIS-05-18-0787-PDN Link, ISI, Google ScholarZhou, H., Liu, P.-P., Qiu, S., Wei, S.-J., Xia, K., and Gao, Q. 2018. Identity of Epicoccum sorghinum causing leaf spot disease of Bletilla striata in China. Plant Dis. 102:1039. https://doi.org/10.1094/PDIS-11-17-1757-PDN Link, ISI, Google ScholarFunding: The project was supported by the Post-subsidy project of National Key Research Development Program of China (2018-5262), National Natural Science Foundation of China (numbers 31860515 and 21977023), Program of Introducing Talents to Chinese Universities (D20023) and China Agriculture Research System of MOF and MARA (D09).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. 35, No. 9 September 2022ISSN:0894-0282e-ISSN:1943-7706 Download Metrics Article History Issue Date: 22 Sep 2022Published: 15 Sep 2022Accepted: 9 May 2022 Pages: 875-879 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: 2018-5262 National Natural Science Foundation of ChinaGrant/Award Number: 31860515Grant/Award Number: 21977023 Program of Introducing Talents to Chinese UniversitiesGrant/Award Number: D20023 China Agriculture Research System of MOF and MARAGrant/Award Number: D09 Keywordsdifferential gene expressiondisease resistanceEpicoccum sorghinumlong noncoding RNAsteaThe author(s) declare no conflict of interest.PDF download

Abstract:HomeMolecular Plant-Microbe Interactions®Vol. 35, No. 5Analysis of Competing Endogenous RNAs and MicroRNAs in Tea (Camellia sinensis) Leaves During Infection by the Leaf Spot Pathogen Pestalotiopsis trachicarpicola Previous RESOURCE ANNOUNCEMENT OPENOpen Access licenseAnalysis of Competing Endogenous RNAs and MicroRNAs in Tea (Camellia sinensis) Leaves During Infection by the Leaf Spot Pathogen Pestalotiopsis trachicarpicolaYuqin Yang, Qiaoxiu Yin, Changlong Qiu, Zhongqiu Xia, Hongke Huang, Chen Huang, Xinyue Jiang, Yuanyou Yang, Delu Wang, and Zhuo ChenYuqin 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, Qiaoxiu YinKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Changlong QiuKey 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, Hongke 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, Chen HuangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Xinyue JiangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Yuanyou YangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China, Delu Wang†Corresponding authors: D. L. Wang; E-mail Address: dlwang@gzu.edu.cn, and Z. Chen; E-mail Address: gychenzhuo@aliyun.comCollege of Forestry, Guizhou University, Guiyang, Guizhou 550025, China, and Zhuo Chen†Corresponding authors: D. L. Wang; E-mail Address: dlwang@gzu.edu.cn, and Z. Chen; E-mail Address: gychenzhuo@aliyun.comhttps://orcid.org/0000-0001-7130-8457Key Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, China AffiliationsAuthors and Affiliations Yuqin Yang1 2 Qiaoxiu Yin1 Changlong Qiu1 Zhongqiu Xia1 2 Hongke Huang1 2 Chen Huang1 Xinyue Jiang1 Yuanyou Yang1 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:23 Mar 2022https://doi.org/10.1094/MPMI-10-21-0262-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articleGenome AnnouncementPestalotiopsis trachicarpicola Y. M. Zhang & K. D. Hyde is a plant-pathogenic fungus that can cause leaf spot on tea (Camellia sinensis), seriously decreasing the production and quality of tea leaves. Competing endogenous RNAs (ceRNAs) include messenger RNAs (mRNAs), circular RNAs (circRNAs), and long noncoding RNAs (lncRNAs), which can regulate gene expression through binding to response elements in microRNAs (miRNAs), which, in turn, operate posttranscriptional regulation of target mRNAs. In the current study, ceRNAs and miRNAs originating from infected and uninfected leaves of tea (C. sinensis var. sinensis cv. Fuding-dabaicha) were sequenced and analyzed, and the interaction between individual ceRNAs and miRNAs was predicted, with target genes being identified. Of the ceRNA interaction pairs, 31, 24, six, and two were predicted in the form of miRNA-mRNA, miRNA-lncRNA, miRNA-circRNA, and mRNA-lncRNA interactions, respectively. Three types of interactions were predicted to occur in the regulation pairs of lncRNA/circRNA-miRNA-mRNA: i) circRNA-miRNA-mRNA (type 1 interactions), ii) lncRNA-miRNA-mRNA (type 2 interactions), and iii) circRNA-lncRNA-miRNA-mRNA (type 3 interactions). Some interaction pairs in lncRNA/circRNA-miRNA-mRNA in infected tea leaves ultimately changed the abundance of mRNAs, such as TEA015316.1:TEA015316.1, TEA020085.1:TEA020085.1, TEA027239.1:TEA027239.1, and TEA001447.1:TEA001447.1. Analysis of the types and expression abundance of ceRNAs and miRNAs originating from tea leaves during infection by P. trachicarpicola, and their interactive relationship will provide a valuable resource for future research into this host-pathogen interaction.Tea [Camellia sinensis (L.) O. Kuntze] is a perennial, evergreen, small tree or shrub and is an important cash crop grown in plantations in tropical and subtropical regions of Asia (Li et al. 2019b). Due to its perennial nature, widespread area of cultivation, and the complex, ecological environment of tea plantations, tea suffers from a number of diseases, largely caused by fungal phytopathogens (Pandey et al. 2021). These diseases mainly occur on the leaves of tea and reduce the production and quality of tea leaves (Pandey et al. 2021).In recent years, tea has been grown on a very large scale in Guizhou Province, China. The area of the tea plantations in Guizhou Province reached 350,000 square hectometers in 2017, making this province the top tea-growing region in China (Bao et al. 2019; Li et al. 2019a; Wang et al. 2020; Yin et al. 2021). Tea plantations in Guizhou Province combine the characters of high altitude (from 800 to 1,800 m above sea level), lower temperatures during the early spring, cold spells in late spring, and a rainy climate. As a consequence, tea suffers from a number of important diseases caused by pathogenic fungi, which are widespread and cause serious damage (Wang et al. 2020; Yin et al. 2021). In previous studies, the fungal pathogens causing tea leaf spot have been identified as Lasiodiplodia theobromae (Li et al. 2019a), Epicoccum sorghinum (Bao et al. 2019), Epicoccum nigrum (Yin et al. 2022), Didymella bellidis (Wang et al. 2020), and Alternaria longipes (Yin et al. 2021). In addition, we found that P. trachicarpicola can also cause an economically important leaf spot on tea plantations in Yuqing County, Guizhou Province, China. Early symptoms of infection of tea leaves by P. trachicarpicola are light brown spots that gradually increase in size and develop into brown, scattered, elliptical, or irregular lesions. The lesions later become grayish at the center and become surrounded by an obvious bright yellow halo. Adjacent lesions can coalesce into larger lesions, which cover the margin of the leaves. For the leaf spots with typical symptoms, the disease incidence was estimated to range between 75 and 83%. The disease severity was estimated to range from 75 to 86, depending on the field. In tea plantations with severe disease, the loss of yield can reach more than 68% (Li et al. 2021).P. trachicarpicola can also cause leaf spot of Trachycarpus fortunei (Maharachchikumbura et al. 2012; Zhang et al. 2012), fruit rot and trunk diseases of grapevine (Jayawardena et al. 2015), leaf spot of Eucommia ulmoides (Li et al. 2018), mango leaf spot (Shu et al. 2020), twig blight of Pinus bungeana (Qi et al. 2021), and gray blight disease of tea (Tsai et al. 2021).The components of ceRNAs, namely, circRNAs, lncRNAs, and mRNAs, can bind with specific sequence elements of miRNAs, resulting in posttranscriptional regulation of the expression of genes during growth, development, and environmental response (Hansen et al. 2013; Tay et al. 2014). ceRNAs play a role in the interaction between plant hosts and fungal pathogens (Cui et al. 2020; Fan et al. 2020; Jiang et al. 2020; Sun et al. 2020). For instance, infection of rice by Magnaporthe oryzae can induce differential expression of circRNAs originating from the host, and a large number of circRNAs have been shown to be involved in the disease defense response of the rice host (Fan et al. 2020). The present study identified tea leaf ceRNAs that were differentially expressed (DE) in response to infection and determined potential target genes that could be affected by interactions between these DE ceRNAs and the mRNAs of the target genes.In this study, a single-spore isolate of P. trachicarpicola ACCC39730 was cultured on potato dextrose agar (PDA) at 25°C in the dark, to form colonies for use as inoculum after incubation for 10 days. Three-year-old bushes of C. sinensis var. sinensis cv. Fuding-dabaicha were grown in the greenhouse under day and night conditions of 25°C and 20°C, respectively, and a relative humidity of 70 to 80%. Detached tea bush twigs were placed into 50-ml conical flasks (one twig per flask) containing 30 ml of sterile water, and the third leaf on each tea twig was inoculated. Four 6-mm diameter holes were punctured into the adaxial leaf surface, two on each side of the main vein of each leaf, using a sterile syringe needle. Mycelial plugs (6 mm diameter) were cut from the edge of a 10-day-old single-spore colony of P. trachicarpicola ACCC39730 on PDA, using a sterile cork borer, and were then placed, mycelial face downward, on the inoculation sites, one plug per site. For the uninoculated controls, plugs of sterile PDA were used instead. Damp cotton was placed over the plugs on the inoculation sites, and the tea twigs were then covered, using self-sealing plastic bags (20 cm × 30 cm) to retain high humidity, for 48 h. The tea twigs were then placed in the plant light incubator (NingboJiangnan Instrument Factory, Ningbo City, China) and were maintained under conditions of 25°C in the day and 20°C at night (200 μmol m−2⋅s−1) and a relative humidity of 70 to 80%. At 48 h postinoculation, the inoculated tea leaves showed symptoms but there were no symptoms on the control leaves. The experiment was set up with two different samples as follows. Control (uninfected) tea leaves were punctured by sterilized syringe needles, and tissue samples were harvested at 0 days postinoculation. Treatment (infected) tea leaves were punctured by sterilized syringe needles and were inoculated with the mycelial plugs and tissue samples were harvested at 1 day postinoculation. Each control and treatment sample contained 15 technical replicates, each replicate consisting of 15 tea twigs. The third leaf from each of the 15 twigs was pooled to carry out total RNA isolation, representing one biological replicate. The tea twigs were randomly assigned to the control or treatment groups, and each experiment was carried out with three independent biological replicates.Fifteen tea leaves from 15 tea twigs, representing one biological replicate of a sample, were pooled and were ground into a fine powder in the presence of liquid nitrogen. Then, 100 mg of the fine powder was used to extract total RNA, using TRIzol reagent (Invitrogen, Carlsbad, CA, U.S.A.), following manufacturer procedures. The RNA amount and purity of each sample were quantified and assessed using a NanoDrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, U.S.A.) and an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, U.S.A.), with the samples required to meet a threshold RNA integrity number >7.0.The cDNA libraries of all RNA species were prepared using the TruSeq Stranded Total RNA HT sample prep kit (Illumina, San Diego, CA, U.S.A.) and the TruSeq Small RNA sample prep kit (Illumina). Following the manufacturer-recommended protocol, ceRNAs were sequenced on an Illumina HiSeq 4000 at Lc-Bio Technologies Co., Ltd. (Hangzhou, China). Cutadapt software was used to remove the low-quality reads, including adapters, low-quality bases, and undetermined bases (Martin 2011). Then, the sequence quality was verified, using FastQC software to obtain clean reads for subsequent downstream analyses.The sequence analysis was divided into four parts. First, for the mRNA component, clean reads were aligned to the reference genome of C. sinensis var. sinensis using HISAT software (version 2.0.4) (Kim et al. 2015). The mapped reads of each sample were assembled, then, their abundance was quantified by calculating the value of fragments per kilobase of transcript per million mapped reads at the level of log2 (fold change) > 1 or < −1 and with statistical significance (P < 0.05), using StringTie software (version stringtie-1.3.4d) (Frazee et al. 2015; Pertea et al. 2015) and Ballgown software package in R (Frazee et al. 2015). The mRNAs targeted by the DE ceRNAs were annotated in the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, using BLAST (Buchfink et al. 2015; Love et al. 2014).Second, for the lncRNA component, sequence quality was verified using FastQC. Valid reads were used to map reads to the genome of C. sinensis var. sinensis, using the Bowtie 2 package in R (Langmead and Salzberg 2012) and TopHat2 (Kim et al. 2013). The mapped reads of each sample were assembled using StringTie software (version stringtie-1.3.4d) (Pertea et al. 2015), then, the transcriptomes from all replicates of both the infected and control samples were merged, to construct a comprehensive transcriptome, using the Perl program. After the final transcriptome was generated, the abundance of each transcript was estimated using StringTie software (Pertea et al. 2015) and the Ballgown package in R (Frazee et al. 2015). All unrelated transcripts, following identification according to the following criteria, were removed: transcripts that overlapped with known mRNAs or with novel mRNAs with coding potential, transcripts that were shorter than 200 bp, and transcripts with coverage of reads <3. The remaining transcripts with coding potential were used for prediction, using the Coding Potential Calculator (CPC) software (version 0.9-r2) (Kong et al. 2007) and Coding-Non-Coding Index (CNCI) software (version 2.0) (Sun et al. 2013), with transcripts with a CPC score <−1 and a CNCI score <0 being removed. The remaining transcripts were considered to be lncRNAs. The expression level of lncRNAs was calculated using StringTie software (Pertea et al. 2015; Trapnell et al. 2010) and the Ballgown package in R (Frazee et al. 2015).Third, for the circRNA component, Cutadapt was used to remove the reads that contained adapter contamination, low-quality bases, and undetermined bases (Martin 2011). Then, sequence quality was verified using FastQC. Bowtie2 (version 2.0.0-beta4) and TopHat2 software were used to map reads to the genome of C. sinensis var. sinensis (Kim et al. 2013; Langmead and Salzberg 2012). Any remaining unmapped reads were still mapped to the genome, using TopHat-Fusion software (Kim and Salzberg 2011). CIRCexplorer2 software was used to de novo assemble the mapped reads to circRNAs (Zhang et al. 2014, 2016). Then, back-splicing reads were identified in unmapped reads, using TopHat-Fusion software and CIRCexplorer2 software (Kim and Salzberg 2011; Zhang et al. 2014, 2016). All samples were analyzed to generate unique circRNAs, and the empirical analysis of digital gene expression data package in R (edgeR, version 3.13) was used to analyze the DE CircRNAs under the conditions of log2 (fold change) >1 or < −1 and with statistical significance (P < 0.05) (Robinson et al. 2010).Finally, for the miRNA component, raw reads were screened using ACGT101-miR software (LC Sciences, Houston, TX, U.S.A.), and adapter dimers and low-complexity sequences were removed, with the remaining reads being filtered by aligning with the databases of mRNA, RNA Families (ribosomal RNA, transfer RNA, small nuclear RNA, and small nucleolar RNA), and Repbase. Reads with lengths of 18 to 25 nt (miRNA-like) were retained, then, unique sequences with lengths in the range of 18 to 25 nt were mapped to tea species–specific precursors in miRBase 22.0 by BLAST search. Length variation at both the 3′ and 5′ ends and one mismatch within the sequence were allowed in the alignment. The unique sequences mapping to specific species mature miRNAs in hairpin arms were identified as known miRNAs. The unique sequences mapping to the other arm of known specific species precursor hairpin arms, opposite to the annotated mature miRNA-containing arm, were considered to be novel 5p- or 3p-derived miRNA candidates. The remaining sequences were mapped to other selected species precursors (with the exclusion of specific species) in miRBase 22.0 by BLAST search, and the mapped precursor-miRNAs were further BLAST searched against the specific species genomes to determine their genomic locations. Based on the above two steps, known miRNAs were identified. The unmapped sequences were searched for by BLAST against the specific genomes, and the hairpin RNA structures containing the sequences aligning in specific genomes were predicted from the flanking 120-nt sequences, using RNAfold software. The DE miRNAs, based on normalized deep-sequencing counts, were analyzed for significance using Student’s t test at the threshold of P value <0.05 (Li et al. 2016).The regulation mode of miRNA-mRNA interaction pairs was predicted using TargetFinder software (Allen et al. 2005; Manke et al. 2008; Mrowka et al. 2008; Wasserman and Sandelin 2004). The alignment score, as the screening criterion, was used to predict the interaction between mRNA and miRNA. When the sequences of miRNA and mRNA were aligned, the effect of the complementary sequences was evaluated using the alignment score. A mismatch or missing base was calculated as a “1” score, and a pair of “G: U” mismatches was calculated as a “0.5” score. The alignment score would double if the mismatch of “G: U” took place in the core region of miRNA; in general, the region from the second nucleotide to the thirteenth nucleotide was known as the core region. The alignment score indicated the degree of match between the target gene and miRNA. As the alignment score decreased, the degree of matching became more complete and the results became more reliable. The default maximum penalty in our research was set to a score of “4”.The regulation mode of miRNA-lncRNA/circRNA was predicted using the ssearch36 (36.3.6) software (Farrar 2007; Pearson 2016; Smith and Waterman 1981). In this study, the specific rules are as follows: i) the bump must be on the lncRNA, and the bump is in the middle of the miRNA. The middle region of the miRNA was defined as follows: #19-nt, 9 to 11 bits; #20-nt, 10 to 11 bits; #21-nt, 10 to 12 bits; #22-nt, 11 to 12 bits; #23-nt, 11 to 13 bits; #24-nt, 12 to 13 bits; #25-nt, 12 to 14 bits; #26-nt, 13 to 14 bits; ii) in addition to the middle region of the miRNA, the maximum values of mismatches were ≤4; iii) bumps are not allowed at non-middle positions of the miRNA. lncRNAs may play a cis role, acting on neighboring target genes. The cis target genes of lncRNAs were used to predict the mRNA-lncRNA interaction. In this study, coding genes in the 100,000-bp upstream and downstream base pairs were selected using Python script. Then, we carried out functional analysis of the target genes for lncRNAs by using a BLAST search (Conesa et al. 2005). Significance was expressed as P < 0.05.A total of 31, 24, six, and two of the ceRNA or miRNA interaction pairs were predicted in the miRNA-mRNA, miRNA-lncRNA, miRNA-circRNA, and mRNA-lncRNA interactions, respectively. For alignment scores of 2 and 2.5, miRNA-mRNA interactions were predicted for PC-3p-97705_24/TEA004183.1:TEA004183.1 and PC-3p-125521_18/TEA022449.1:TEA022449.1, respectively. A probable ribose-5-phosphate isomerase 2 (TEA004183.1) was annotated with the GO term “transmembrane transport” and the KEGG pathway “biosynthesis of amino acids”. The uncharacterized gene LOC100251837 (TEA022449.1) was annotated with the GO term “trehalose-phosphatase activity” and the KEGG pathway “starch and sucrose metabolism”. From the prediction of miRNA-lncRNA interactions, nta-MIR172e-p3_2ss15CG19TA/.:MSTRG.5527.1, nta-MIR172e-p3_2ss15CG19TA/.:MSTRG.30017.1, nta-MIR172e-p3_2ss15CG19TA/.:MSTRG.1706.1, nta-MIR172e-p3_2ss15CG19TA/TEA024896.1:MSTRG.5459.1, and PC-3p-103122_23/.:MSTRG.1183.1 exhibited significant interactions, and Smith-Waterman (SW) scores of these pairs were at the top of the list. Interestingly, nta-MIR172e-p3_2ss15CG19TA may target the four lncRNAs described above. In the prediction of miRNA-circRNA interactions, vvi-MIR396b-p3_2ss7AG20AG and circRNA292 may exhibit an interaction, with a SW score of 117 and an E-value of 1.10. The interaction of PC-5p-189883_10 and circRNA776 exhibited a higher SW score of 77 but a lower E-value of 0.00. For the prediction of mRNA-lncRNA interactions, TEA033765.1/MSTRG.24999.1, and TEA030398.1/MSTRG.12855.1 may exhibit interactions, with the Pearson’s correlation coefficients being 0.87 and 0.97, respectively. Integrated analysis of interactions among these four kinds of noncoding RNAs with mRNAs showed that the number of miRNA-mRNA and miRNA-lncRNA interactions was significantly higher than that of miRNA-circRNA interactions (Table 1).Table 1. The targeted prediction analysis of CeRNA from tea [Camellia sinensis (L.) O. Kuntze] leaves during challenge by the fungal pathogen Pestalotiopsis trachicarpicolaaTargeted prediction analysisNumber of correlation analysesmiRNA-mRNA31 pairs23 miRNAs/20 mRNAsmiRNA-lncRNA24 pairs11 miRNAs/16 lncRNAsmiRNA-circRNA6 pairs5 miRNAs/3 circRNAsmRNA-lncRNA2 pairs2 mRNA/2 lncRNAAll lncRNA/circRNA/miRNA/mRNA4/0/1/1 pairsType 1: circRNA-miRNA-mRNAcircRNA292/PC-5p-4673_314/ TEA020085.1:TEA020085.1circRNA791/PC-5p-4673_314/ TEA020085.1:TEA020085.1Type 2: lncRNA-miRNA-mRNA.:MSTRG.5527.1/nta-MIR172e-p3_2ss15CG19TA/ TEA015316.1:TEA015316.1.:MSTRG.1706.1/nta-MIR172e-p3_2ss15CG19TA/ TEA015316.1:TEA015316.1TEA024896.1:MSTRG.5459.1/nta-MIR172e-p3_2ss15CG19TA/TEA015316.1:TEA015316.1.:MSTRG.30017.1/nta-MIR172e-p3_2ss15CG19TA/ TEA015316.1:TEA015316.1Type 3: circRNA-lncRNA-miRNA-mRNAcircRNA292/vvi-MIR396b-p3_2ss7AG20AG/ TEA027239.1:TEA027239.1circRNA292/vvi-MIR396b-p3_2ss7AG20AG/ TEA001447.1:TEA001447.1amiRNA = microRNA, mRNA = messenger RNA, lncRNA = long noncoding RNA, and circRNA = circular RNA.Table 1. The targeted prediction analysis of CeRNA from tea [Camellia sinensis (L.) O. Kuntze] leaves during challenge by the fungal pathogen Pestalotiopsis trachicarpicolaaView as image HTML For analysis of the interaction pairs in lncRNA/circRNA-miRNA-mRNA, the regulation modes were divided into three types, namely circRNA-miRNA-mRNA (type 1), lncRNA-miRNA-mRNA (type 2), and circRNA-lncRNA-miRNA-mRNA (type 3). For type 1, two ceRNA pairs were predicted in the regulation modes and are listed as circRNA292/PC-5p-4673_314/TEA020085.1:TEA020085.1 and circRNA791/PC-5p-4673_314/TEA020085.1:TEA020085.1. For type 2, four pairs of ceRNAs, composed of four lncRNAs, one miRNA, and one mRNAs, were predicted and are listed as .:MSTRG.5527.1/nta-MIR172e-p3_2ss15CG19TA/TEA015316.1:TEA015316.1, .:MSTRG.1706.1/nta-MIR172e-p3_2ss15CG19TA/TEA015316.1:TEA015316.1, TEA024896.1:MSTRG.5459.1/nta-MIR172e-p3_2ss15CG19TA/TEA015316.1:TEA015316.1, and .:MSTRG.30017.1/nta-MIR172e-p3_2ss15CG19TA/TEA015316.1:TEA015316.1. GO annotation of TEA015316.1:TEA015316.1 indicated that it was involved with GO:0006810 (transport), GO:0005524 (ATP binding), GO:0016021 (integral component of membrane), and GO:0016887 (ATPase activity). For type 3, two ceRNA pairs were predicted in the regulation modes, listed as circRNA292/vvi-MIR396b-p3_2ss7AG20AG/TEA027239.1:TEA027239.1 and circRNA292/vvi-MIR396b-p3_2ss7AG20AG/TEA001447.1:TEA001447.1 (Table 1).Based on the analysis of interactions between ceRNAs and miRNAs, the abundance of mRNAs originating from tea leaves was predicted. For instance, lncRNAs .:MSTRG.5527.1, .:MSTRG.1706.1, TEA024896.1:MSTRG.5459.1, and .:MSTRG.30017.1 competitively bound miRNA nta-MIR172e-p3_2ss15CG19TA, and the abundance of mRNA TEA015316.1:TEA015316.1 was significantly upregulated through the comparison of the mRNA abundance of infected and uninfected samples. CircRNA292 and circRNA791 competitively bound miRNA PC-5p-4673_314, so that the abundance of mRNA TEA020085.1:TEA020085.1 was significantly upregulated. CircRNA292 competitively bound the miRNA vvi-MIR396b-p3_2ss7AG20AG, resulting in the abundance of mRNAs TEA027239.1:TEA027239.1 and TEA001447.1:TEA001447.1 being significantly down- and upregulated, respectively. These three mRNAs were not annotated in the GO database (Table 1).Data AvailabilityIn this study, the sequences and the interaction analysis of ceRNAs and miRNAs originating from infected tea leaves have provided a significant data resource for use in understanding the pathogenicity of P. trachicarpicola and the response mechanism of tea. The sequences of ceRNAs and miRNA of tea during infection by P. trachicarpicola have been deposited in GenBank under SRA accession number PRJNA756805 for lncRNAs and circRNAs, PRJNA756832 for mRNAs, and PRJNA742222 for small RNAs. In addition, the bioinformatics analysis of ceRNAs is presented in the Supplementary Tables.Author-Recommended Internet ResourcesCIRCexplorer2 software: https://github.com/YangLab/CIRCexplorer2Cutadapt: https://github.com/marcelm/cutadaptedgeR: http://bioconductor.org/packages/release/bioc/html/edgeR.htmlFastQC: https://www.bioinformatics.babraham.ac.uk/projects/fastqcHISAT software (version 2.0.4): http://www.ccb.jhu.edu/software/hisat/index.shtmlPerl program: https://www.perl.orgRepbase: https://www.girinst.org/repbaseRNA Families database: http://rfam.janelia.orgStringTie software: https://ccb.jhu.edu/software/stringtieTargetFinder software: https://github.com/carringtonlab/TargetFinderTopHat-Fusion software: http://tophat-fusion.sourceforge.netThe author(s) declare no conflict of interest.Literature CitedAllen, E., Xie, Z., Gustafson, A. M., and Carrington, J. C. 2005. microRNA-directed phasing during trans-acting siRNA biogenesis in plants. Cell 121:207-221. https://doi.org/10.1016/j.cell.2005.04.004 Crossref, Medline, ISI, Google ScholarBao, 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 Link, ISI, Google ScholarBuchfink, B., Xie, C., and Huson, D. H. 2015. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12:59-60. https://doi.org/10.1038/nmeth.3176 Crossref, Medline, ISI, 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, Medline, ISI, Google ScholarCui, J., Jiang, N., Hou, X., Wu, S., Zhang, Q., Meng, J., and Luan, Y. 2020. Genome-wide identification of lncRNAs and analysis of ceRNA networks during tomato resistance to Phytophthora infestans. Phytopathology 110:456-464. https://doi.org/10.1094/PHYTO-04-19-0137-R Link, ISI, Google ScholarFan, J., Quan, W., Li, G.-B., Hu, X.-H., Wang, Q., Wang, H., Li, X.-P., Luo, X., Feng, Q., Hu, Z.-J., Feng, H., Pu, M., Zhao, J.-Q., Huang, Y.-Y., Li, Y., Zhang, Y., and Wang, W.-M. 2020. circRNAs are involved in the rice–Magnaporthe oryzae interaction. Plant Physiol. 182:272-286. https://doi.org/10.1104/pp.19.00716 Crossref, Medline, ISI, Google ScholarFarrar, M. 2007. Striped Smith-Waterman speeds database searches six times over other SIMD implementations. Bioinformatics 23:156-161. https://doi.org/10.1093/bioinformatics/btl582 Crossref, Medline, ISI, 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, Medline, ISI, Google ScholarHansen, T. B., Jensen, T. I., Clausen, B. H., Bramsen, J. B., Finsen, B., Damgaard, C. K., and Kjems, J. 2013. Natural RNA circles function as efficient microRNA sponges. Nature 495:384-388. https://doi.org/10.1038/nature11993 Crossref, Medline, ISI, Google ScholarJayawardena, R. S., Zhang, W., Liu, M., Maharachchikumbura, S. S., Zhou, Y., Huang, J., Nilthong, S., Wang, Z., Li, X., Yan, J., and Hyde, K. D. 2015. Identification and characterization of Pestalotiopsis-like fungi related to grapevine diseases in China. Fungal Biol. 119:348-361. https://doi.org/10.1016/j.funbio.2014.11.001 Crossref, Medline, ISI, Google ScholarJiang, N., Cui, J., Hou, X., Yang, G., Xiao, Y., Han, L., Meng, J., and Luan, Y. 2020. Sl-lncRNA15492 interacts with Sl-miR482a and affects Solanum lycopersicum immunity against Phytophthora infestans. Plant J. 103:1561-1574. https://doi.org/10.1111/tpj.14847 Crossref, Medline, ISI, 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, Medline, ISI, 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-

Sequences and Genome-Wide Analysis of Mrna and Microrna Expression in Tea (camellia Sinensis) Leaves in Response to Epicoccum Sorghinum Infection

Sequencing and Functional Annotation of Mrnas and Lncrnas from Tea (camellia Sinensis) Leaves During Infection by the Fungal Pathogen Lasiodiplodia Theobromae

Integrated Mrna and Small RNA Sequencing for Analyzing Leaf Spot Pathogen Didymella Segeticola and Its Host, Tea (camellia Sinensis), During Infection

Integrated Sequencing Data, Annotation, and Targeting Analysis of Mrnas and MicroRNAs from Tea Leaf During Infection by Tea Leaf Spot Pathogen, Epicoccum Nigrum

Sequencing, Functional Annotation, and Interaction Prediction of Mrnas and Candidate Long Noncoding RNAs Originating from Tea Leaves During Infection by the Fungal Pathogen Causing Tea Leaf Spot, Didymella Bellidis

Transcriptome Profiling Data of Tea Leaves During Infection by the Dark Brown Tea Leaf Spot Pathogen, Epicoccum mackenziei

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

The Sequence and Integrated Analysis of Competing Endogenous RNAs Originating from Tea Leaves Infected by the Pathogen of Tea Leaf Spot, Didymella Segeticola.

Sequence Data, Functional Annotation, and Relationship Analysis Between Mrnas and Long Noncoding RNAs from Tea Leaves During Infection by the Fungal Pathogen Epicoccum Sorghinum.

Identification of Regulatory Networks of MicroRNAs and Their Targets in Response to Colletotrichum gloeosporioides in Tea Plant (Camellia sinensis L.)

Analysis of Competing Endogenous RNAs and MicroRNAs in Tea (camellia Sinensis) Leaves During Infection by the Leaf Spot Pathogen Pestalotiopsis Trachicarpicola

Involvement Of Salicylic Acid In Anthracnose Infection In Tea Plants Revealed By Transcriptome Profiling

Functional Annotation of Circrnas in Tea Leaves after Infection by the Tea Leaf Spot Pathogen, Lasiodiplodia Theobromae.

Integrated mRNA and Small RNA Sequencing for Analyzing Tea Leaf Spot Pathogen <i>Lasiodiplodia theobromae</i>, Under In Vitro Conditions and the Course of Infection

Lifestyle Characteristics and Gene Expression Analysis of Colletotrichum Camelliae Isolated from Tea Plant [camellia Sinensis (L.) O. Kuntze] Based on Transcriptome.

First Report of Lasiodiplodia Theobromae Causing Leaf Spot on Tea Plant in Guizhou Province of China

Impact of Exogenous Caffeine on Regulatory Networks of Micrornas in Response to Colletotrichum Gloeosporioides in Tea Plant

Transcriptomic Analysis Reveals the Molecular Adaptation of Three Major Secondary Metabolic Pathways to Multiple Macronutrient Starvation in Tea (<i>Camellia sinensis</i>)

Quantification of the fungal pathogen Didymella segeticola in Camellia sinensis using a DNA-based qRT-PCR assay

HS-SPEM-GC-MS Analyses of Volatiles and Transcriptome Analysis of Camellia Tetracocca in Response to Feeding by Empoasca Onukii Matsuda

Multi-omics analysis to visualize the dynamic roles of defense genes in the response of tea plants to gray blight