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-

Abstract:HomePhytoFrontiers™Vol. 2, No. 3Sequencing and Functional Annotation of Competing Endogenous RNAs and MicroRNAs in Tea Leaves During Infection by Lasiodiplodia theobromae PreviousNext RESOURCE ANNOUNCEMENT OPENOpen Access licenseSequencing and Functional Annotation of Competing Endogenous RNAs and MicroRNAs in Tea Leaves During Infection by Lasiodiplodia theobromaeDi Guo, Zhongqiu Xia, Xinyue Jiang, Hongke Huang, Yuqin Yang, Shilong Jiang, Delu Wang, and Zhuo ChenDi GuoKey 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, Xinyue JiangKey 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, 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, Shilong JiangKey Laboratory of Green Pesticide and Agricultural Bioengineering, Ministry of Education, Guizhou University, Guiyang, Guizhou 550025, ChinaAgricultural College, Guizhou University, Guiyang, Guizhou 550025, ChinaSearch for more papers by this author, 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, ChinaSearch for more papers by this author, and Zhuo Chen†Corresponding authors: D. L. Wang; E-mail Address: dlwang@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 Di Guo1 2 Zhongqiu Xia1 2 Xinyue Jiang1 Hongke Huang1 2 Yuqin Yang1 2 Shilong Jiang1 3 Delu Wang4 † 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 3Agricultural College, Guizhou University, Guiyang, Guizhou 550025, China 4College of Forestry, Guizhou University, Guiyang, Guizhou 550025, China Published Online:31 Mar 2022https://doi.org/10.1094/PHYTOFR-11-21-0075-AAboutSectionsView articlePDFSupplemental ToolsAdd to favoritesDownload CitationsTrack Citations ShareShare onFacebookTwitterLinked InRedditEmailWechat View articleAnnouncement of Sequences of Competing Endogenous RNAs and miRNAsLasiodiplodia theobromae is a plant-pathogenic fungus which can cause leaf spot in tea (Camellia sinensis (L.) O. Kuntze), seriously decreasing the productivity and quality of the crop. 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). In total, 40,197 candidate regulatory pairs isolated from L. theobromae-infected leaves were predicted in the circRNA-lncRNA-miRNA-mRNA network, with 39,092 of the candidate regulatory pairs being at the circRNA-lncRNA-miRNA-mRNA level (type 3), 1,103 being at the lncRNA-miRNA-mRNA level (type 2), and 2 at the circRNA-miRNA-mRNA level (type 1). Differentially expressed candidate genes targeted by ceRNAs originating from inoculated tea leaves were predicted, including pathogenesis-related gene transcriptional activator PTI6 (putative), probable receptor-like protein kinase At5g15080, disease resistance protein RPM1-like, and WRKY transcription factor 22-like, and were enriched in the plant–pathogen interaction pathway. Analysis of the types and expression (in terms of abundance) of ceRNAs and miRNAs originating from tea leaves during infection by L. theobromae and their interaction will provide a valuable resource for future research into the host–pathogen interactions.The Botryosphaeriaceae family (phylum Ascomycota) consists of saprobes, endophytes, and phytopathogens, the latter having largely woody plant hosts (Chethana et al. 2016; Morales-Cruz et al. 2015; Paolinelli-Alfonso et al. 2016; Phillips et al. 2013; Yan et al. 2013). In this family, Lasiodiplodia Ellis & Everh. was clustered in a clade containing 25 species, which can cause plant disease (Phillips et al. 2013; Yan et al. 2013). Within this genus, Lasiodiplodia theobromae (Pat.) Griffon & Maubl., whose anamorph is Botryosphaeria rhodina (Berk. & M. A. Curtis) Arx, can cause diseases in many plants, including canker, dieback, fruit spot, stem-end rot, postharvest fruit rot, gummosis, inflorescence blight, and leaf spot (Cardoso et al. 2002; Garibaldi et al. 2012; Li et al. 2013; Luo et al. 2011; Pérez et al. 2017; Phipps and Porter 1998; Rodríguez-Gálvez et al. 2015; Serrato-Diaz et al. 2013; Serrato-Diaz et al. 2014; Woodward et al. 2005; Yildiz et al. 2014; Zheng et al. 2021). Because of the high incidence and severity of the disease, tea leaf spot infected by L. theobromae results in serious decreases in the quality and yield of tea (C. sinensis (L.) O. Kuntze) leaves, causing substantial economic losses in some tea plantations in Guizhou Province, China (Li et al. 2019).In a previous study, the mRNAs and miRNAs originating from L. theobromae and tea (C. sinensis) leaves during infection by the pathogen were sequenced, and the differentially expressed (DE) mRNAs and DEmiRNAs originating from the pathogen and tea leaves were predicted (Jiang et al. 2021). The results indicated that miRNAs stu-MIR5303j-p5_2ss10GT17CT, cpa-MIR171a-p3_2ss12TG18TG, aly-miR157a-5p, gma-miR166a-3p, and gma-MIR167c-p3_2ss14TA20AG could potentially interact with the greatest number of target mRNAs of tea leaves. These target genes were significantly enriched in the pathways of amino sugar and nucleotide sugar metabolism, plant hormone signal transduction, and zeatin biosynthesis (Jiang et al. 2021). Additionally, the miRNAs gma-miR6300_R+1, PC-5p-130614_34, and cpa-MIR171a-p3_2ss12TG18TG could interact with the greatest number of target mRNAs of L. theobromae. These target genes were significantly enriched in the pathways of ubiquitin-mediated proteolysis, RNA degradation, and starch and sucrose metabolism (Jiang et al. 2021).Based on these studies, the sequences of the DEmRNAs and DElncRNAs in the comparison between infected leaves versus uninfected leaves were predicted and, in total, 2,503 and 5,278 of the DEmRNAs and 3,403 and 1,477 of the DEIncRNAs were up- and downregulated, respectively, by L. theobromae during the infection process. The effect of cis-regulation between the DElncRNAs and DEmRNAs was predicted (Li et al. 2021). The miRNAs can bind to specific sequences in target mRNAs, then block their translation. circRNAs are a type of lncRNAs with covalently closed circles, lacking both a 5′ cap and a 3′ tail, which include circular exonic RNAs (circRNA), circular intronic RNAs (ciRNA), and exon-intron circRNAs (EIciRNAs) (Jeck and Sharpless 2014; Li et al. 2015, 2016). lncRNAs are transcribed by RNA polymerase II, III, IV, and V; are longer than 200 nucleotides (nt); and are divided into intronic lncRNA, bidirectional lncRNA, sense lncRNA, and antisense lncRNA (Mercer et al. 2009; Ponting et al. 2009). miRNAs are small noncoding RNAs of 18 to 25 nt, which function posttranscriptionally by usually base-pairing to the mRNA to repress protein synthesis (Fabian et al. 2010). CircRNAs, lncRNAs, and mRNAs are known to be important components of ceRNAs, which can bind with specific elements of miRNAs, inhibiting the effect of the latter and, in turn, regulating the expression of genes during growth, development, and environmental response (Hansen et al. 2013; Tay et al. 2014). Thus far, there have been few reports on the role of ceRNAs and miRNAs in the interaction between plant hosts and fungal pathogens (Cui et al. 2020; Fan et al. 2020; Sun et al. 2020). A recent study indicated that infection by Magnaporthe oryzae can induce formation of DEcircRNAs originating from the host, with a large number of circRNAs being involved with defense by the rice host against the pathogen (Fan et al. 2020). Because tea foliar disease can seriously influence the production and quality of tea leaves, the sequence and the functional annotation of ceRNAs and miRNAs would contribute to learning about the pathogenic mechanism and how to cultivate disease-resistant varieties of tea. Here, we conducted sequence and functional annotation of ceRNAs and miRNAs in tea leaves during infection by L. theobromae, and targeted prediction analysis between ceRNAs and miRNAs.In the current study, a single conidium of L. theobromae CGMCC3.20151 was inoculated onto potato dextrose agar (PDA) and incubated at 28°C in the dark to form a single-spore colony after 10 days. Three-year-old bushes of C. sinensis var. sinensis ‘Fuding-dabaicha’ were grown in the greenhouse under conditions of 25 and 20°C (day and night, respectively), and a relative humidity of 70 to 80%. Detached tea twigs were placed into a 50-ml conical flask, each containing 30 ml of sterile water, and the third and fourth leaves on each tea twig were inoculated. Four holes (6 mm in diameter) were punctured into the adaxial surface of the leaf, two on either side of the main vein of each leaf, using a sterile syringe needle. Mycelial plugs (diameter = 6 mm) were cut from the edge of 10-day-old single-spore colonies of L. theobromae, using a sterile hole punch, which were placed, mycelial face downward, on the inoculation sites on the leaf adaxial surface. For the uninfected treatment, the inoculation sites on the tea leaves were covered with plugs of sterile PDA. Damp cotton wool was placed over the agar plugs on the inoculation sites; then, each tea twig was enclosed in a self-sealing Ziplock plastic bag to raise the relative humidity for 48 h. The flasks containing the tea twigs were placed in an illuminated plant growth chamber incubator and maintained under day and night temperatures of 25 and 20°C, respectively, with a cycle of 14 h/day and 8 h/night and a relative humidity of 70 to 80%. At 48 h postinoculation, the pathogen-inoculated tea leaves showed symptoms but there were no symptoms on the leaves inoculated with sterile PDA plugs. Each treatment contained 15 tea twigs. Each experiment was carried out with three independent biological replicates.In all, 30 tea leaves from 15 tea twigs were pooled to represent one treatment replicate, then ground into a fine powder under the protection 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 the manufacturer's procedure. The amount and purity of the RNA in 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, with samples required to meet the threshold RNA integrity number > 7.0. Approximately 5 μg of total RNA was used to deplete ribosomal RNA (rRNA) according to the protocol of the Ribo-Zero rRNA Removal Kit (Illumina, San Diego, CA, U.S.A.). The cDNA libraries of all RNA species were prepared using a TruSeq Stranded Total RNA HT Sample Prep Kit (Illumina) and TruSeq Small RNA Sample Prep Kit (Illumina).Following the manufacturer's recommended protocol, ceRNAs were sequenced on an Illumina HiSeq 4000 platform at Lc-bio Technologies Co., Ltd. (Hangzhou, China). Cutadapt (https://github.com/marcelm/cutadapt) 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 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) 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 (http://tpdb.shengxin.ren/index.html), using HISAT (hierarchical indexing for spliced alignment of transcripts) software (version 2.0.4) (http://www.ccb.jhu.edu/software/hisat/index.shtml) (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 log2 (fold change) < −1 and with statistical significance (P value < 0.05), using StringTie software (version stringtie-1.3.4d) (https://ccb.jhu.edu/software/stringtie/) (Frazee et al. 2015; Pertea et al. 2015) and the Ballgown R package (Frazee et al. 2015). The mRNAs targeted by DEceRNAs were annotated in the databases of the Gene Ontology (GO) (http://geneontology.org) (version 2016.04) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) (https://www.kegg.jp/kegg) (version 2016.05) using BLAST2GO software (Conesa et al. 2005).Second, for the lncRNA component, sequence quality was 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 HISAT2 software (https://daehwankimlab.github.io/hisat2/) (Kim et al. 2015, 2019; Langmead and Salzberg 2012) and TopHat2 (http://ccb.jhu.edu/software/tophat) (Kim et al. 2013). The mapped reads of each sample were assembled using StringTie software (version stringtie-1.3.4d) (https://ccb.jhu.edu/software/stringtie/) (Pertea et al. 2015); then, the transcriptomes from the replicates of the infected or the uninfected treatment samples were merged to construct a comprehensive transcriptome, using the Perl program (https://www.perl.org). 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 were removed based on the following criteria: (i) transcripts which overlapped with known mRNAs or with novel mRNAs with coding potential, (ii) transcripts which were shorter than 200 bp, and (iii) transcripts with coverage of reads < 3. Based on the above steps, 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 all 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 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Bowtie2 (version 2.0.0-beta4) and TopHat2 software (http://ccb.jhu.edu/software/tophat) were used to map reads to the genome of C. sinensis var. sinensis (http://tpdb.shengxin.ren/index.html) (Kim et al. 2013; Langmead and Salzberg 2012). As unmapped reads, remaining reads were still mapped to the genome using TopHat-Fusion software (http://tophat-fusion.sourceforge.net/) (Kim and Salzberg 2011). CIRCexplorer2 software (https://github.com/YangLab/CIRCexplorer2) was used to assemble the mapped reads to circRNAs de novo (Zhang et al. 2014, 2016). Then, back-splicing reads were identified in unmapped reads by TopHat-Fusion software (http://tophat-fusion.sourceforge.net/) and CIRCexplorer2 software (https://github.com/YangLab/CIRCexplorer2) (Kim and Salzberg 2011; Zhang et al. 2014, 2016). All samples were analyzed to generate unique circRNAs, and the R package Empirical Analysis of Digital Gene Expression data in R (edgeR) (version 3.13) (http://bioconductor.org/packages/release/bioc/html/edgeR.html) was used to analyze the DEcircRNAs under the conditions of log2 (fold change) > 1 or log2 (fold change) < −1 and with statistical significance (P value < 0.05) (Robinson et al. 2010).Finally, for the miRNA component, raw reads were treated using ACGT101-miR software (LC Sciences, Houston, TX, U.S.A.), adapter dimers and low-complexity sequences were removed, and the remaining reads were filtered by aligning the databases of mRNA, RFam (rRNA, transfer RNA, small nuclear RNA, and small nucleolar RNA) (http://rfam.janelia.org), and Repbase (http://www.girinst.org/repbase). Reads with lengths of 18 to 25 nt (miRNA-like) were retained; then, unique sequences with lengths in the range 18 to 25 nt were mapped to tea-species-specific precursors in miRBase 22.0 (ftp://mirbase.org/pub/mirbase/CURRENT/) by BLAST search. Length variation at both the 3’ and 5’ ends and one mismatch inside of the sequence were allowed in the alignment. Unique sequences mapping to species-specific mature miRNAs in hairpin arms were identified as known miRNAs. Unique sequences mapping to the other arm of known species-specific precursor hairpin sequences opposite to the annotated mature miRNA-containing arm were considered to be novel 5’p- or 3’p-derived miRNA candidates. The remaining unmapped sequences were mapped to the genomic sequences of C. sinensis in search of candidate novel miRNAs by BLAST search.Based on the above two steps, known miRNAs were determined. The unmapped sequences were searched by BLAST against the genomes of C. sinensis, and the hairpin RNA structures containing the sequences aligning to specific genomes were predicted from the flanking 120-nt sequences using RNAfold software (http://rna.tbi.univie.ac.at/cgi-bin/RNAWebSuite/RNAfold.cgi). The noncoding sequences that could form a stem-loop structure and meet the standard of miRNA prediction (Axtell et al. 2018) were regarded as the true miRNAs of C. sinensis. The DEmiRNAs, based on normalized deep-sequencing counts, were analyzed using Student's t test and the significance threshold was set to be 0.05.The regulation mode of miRNA-mRNA was predicted using TargetFinder software (https://github.com/carringtonlab/TargetFinder) (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” was calculated as a “0.5” score. The alignment score would increase twofold when the mismatch of G:U took place in the core region of the miRNA; in general, the distance from nucleotides 2 to 13 was known as the core region. The alignment score indicated the degree of complementarity between the target gene and the miRNA. As the alignment score decreased, the degree of complementarity became more complete, and the results became more reliable. The default maximum penalty in our research was set to a score of “4”.In addition, the regulation mode of miRNA-lncRNA/circRNA was predicted using the ssearch36 software (36.3.6) (Farrar 2007; Pearson 2016; Smith and Waterman 1981). In the current study, the specific rules were listed as follows. (i) The bump must be on the ncRNA, and the bump is in the middle of the miRNA. The middle region of the miRNA was defined: nucleotide 19: 9 to 11 bits, nucleotide 20: 10 to 11 bits, nucleotide 21: 10 to 12 bits, nucleotide 22: 11 to 12 bits, nucleotide 23: 11 to 13 bits, nucleotide 24: 12 to 13 bits, nucleotide 25: 12 to 14 bits, and nucleotide 26: 13 to 14 bits. (ii) In addition to the middle region of the miRNA, the maximum values of mismatches was ≤ 4. (iii) Bumps were not allowed at nonmiddle positions of the miRNA.In total, 40,197 candidate regulatory pairs were predicted in the regulation mode of circRNA-lncRNA-miRNA-mRNA of the tea leaf infected by L. theobromae. Among them, 39,092, 1,103, and 2 candidate regulatory pairs were predicted in circRNA-lncRNA-miRNA-mRNA, lncRNA-miRNA-mRNA, and circRNA-miRNA-mRNA, respectively. For instance, lncRNA MSTRG.28003.1 competitively bound miRNA PC-3p-14126704_1, and lncRNA MSTRG.52547.1 competitively bound miRNA PC-5p-6832841_2. These two lncRNA-miRNA pairs could regulate the expression (in terms of mRNA abundance) of WRKY transcription factor 22-like (TEA028913.1). When expression of lncRNA MSTRG.28003.1 was upregulated, the expression level of the WRKY transcription factor 22-like gene, targeted by miRNA PC-3p-14126704_1, was significantly upregulated. When lncRNA MSTRG.52547.1 bound to miRNA PC-5p-6832841_2, the expression level of the WRKY transcription factor 22-like gene was significantly upregulated (Table 1).Table 1. Competing endogenous RNA sequences from tea (Camellia sinensis var. sinensis (L.) O. Kuntze) leaves during challenge by the fungal pathogen Lasiodiplodia theobromaeSample informationDataaSequenceUninfectedbInfectedcCorrelation analysesRaw datacircRNA-lncRNARaw reads141,647,262145,772,279−Valid reads139,330,285143,076,549−Valid bases20.90 G21.46 G−mRNARaw reads46,471,88949,257,845−Valid reads45,582,99247,825,308−Valid bases6.84 G7.17 G−miRNARaw reads16,673,70016,134,186−Valid reads15,140,4109,506,213−Valid bases−−−Targeted prediction analysismiRNA-mRNA−−−362 pairs (151 miRNAs/222 mRNAs)miRNA-lncRNA−−−2,489 pairs (314 miRNAs/1,705 lncRNAs)miRNA-circRNA−−−115 pairs (55 miRNAs/87 circRNAs)All lncRNA-circRNA-miRNA-mRNA−−−40,197 pairsType 1: circRNA-miRNA-mRNA−−−2 pairsType 2: lncRNA-miRNA-mRNA−−−1,103 pairsType 3: circRNA-lncRNA-miRNA-mRNA−−−39,092 pairsaAbbreviations: circRNA = circular RNA, lncRNA = long noncoding RNA, mRNA = messenger RNA, and miRNA = microRNA.bUninfected treatments represent tea leaves wounded by puncturing four holes on the adaxial side of each leaf, using sterilized syringe needles, with four holes being clustered on the range of 6-mm diameter; then, sterile potato dextrose agar plugs (diameter = 6 mm) were placed on the puncture sites on the tea leaves.cInfected treatments represent tea leaves wounded by puncturing four holes on the adaxial side of each leaf, using sterilized syringe needles, with four holes clustered within the range of 6 mm in diameter; then, the mycelial plugs (diameter = 6 mm) were placed on the puncture sites on the tea leaves.Table 1. Competing endogenous RNA sequences from tea (Camellia sinensis var. sinensis (L.) O. Kuntze) leaves during challenge by the fungal pathogen Lasiodiplodia theobromaeView as image HTML When the KEGG database annotation of genes targeted by the DEceRNAs was carried out, the total number of genes (TB) and the number of genes with significantly different expression between the infected and uninfected treatments (TS) were 17,235 and 137, respectively. The genes targeted by the DEceRNAs were functionally enriched in the pathway of plant–pathogen interaction, including WRKY transcription factor 22-like, kinase protein with adenine nucleotide α hydrolase-like domain, protein NSP-INTERACTING KINASE1, probable serine/threonine-protein kinase At1g18390, chaperone protein ClpB3 (chloroplastic-like), probable receptor-like protein kinase At5g15080, probable leucine-rich repeat (LRR) receptor-like serine/threonine-protein kinase At3g47570, receptor-like protein kinase FERONIA, LysM domain receptor-like kinase 4, probable LRR receptor-like serine/threonine-protein kinase At1g06840, pathogenesis-related gene transcriptional activator PTI6 (putative), disease resistance protein RPM1-like, and putative disease resistance protein At1g50180. When the GO database annotation of the genes targeted by the DEceRNAs was carried out, the TB and TS gene numbers were 13,989 and 111, respectively. At the level of molecular function, the genes targeted by the DEceRNAs included those enriched with respect to ATP binding, including kinase protein with adenine nucleotide α hydrolases-like domain, protein NSP-INTERACTING KINASE1, and probable LRR receptor-like serine/threonine-protein kinase At1g06840.In this study, the sequences and the subsequent analysis of bioinformatics of ceRNAs and miRNAs provide a significant resource for improving our understanding of the pathogenicity of L. theobromae and the response mechanism of tea. The sequences of ceRNAs and miRNAs in tea during the infection process of L. theobromae have been deposited in GenBank under Sequence Read Archive accession numbers PRJNA689677 for circRNAs, PRJNA689677 for lncRNAs, PRJNA637492 for mRNAs, and PRJNA558429 for miRNAs. In addition, the bioinformatics analysis of ceRNAs is presented in Supplementary Tables S1 through S39.The 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, Google ScholarAxtell, M. J., and Meyers, B. C. 2018. Revisiting criteria for plant microRNA annotation in the era of big data. Plant Cell 30:272-284. https://doi.org/10.1105/tpc.17.00851 Crossref, Google ScholarCardoso, J. E., Vidal, J. C., Santos, A. A., Freire, F. C. O., and Viana, F. M. P. 2002. First report of black branch dieback of cashew caused by Lasiodiplodia theobromae in Brazil. Plant Dis. 86:558. https://doi.org/10.1094/PDIS.2002.86.5.558B Link, Google ScholarChethana, K. W. T., Li, X.-H., Zhang, W., Hyde, K. D., and Yan, J.-Y. 2016. Trail decryption of molecular research on Botryosphaeriaceae in woody plants. Phytopathol. Mediterr. 55:147-171. 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 ScholarCui, J., Jiang, N., Hou, X.-X., Wu, S.-H., Zhang, Q., Meng, J., and Luan, Y.-S. 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, Google ScholarFabian, M. R., Sonenberg, N., and Filipowicz, W. 2010. Regulation of mRNA translation and stability by microRNAs. Annu. Rev. Biochem. 79:351-379. https://doi.org/10.1146/annurev-biochem-060308-103103 Crossref, 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, 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, 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 ScholarGaribaldi, A., Bertetti, D., Amatulli, M. T., Cardinale, J., and Gullino, M. L. 2012. First report of postharvest fruit rot in avocado (Persea americana) caused by Lasiodiplodia theobromae in Italy. Plant Dis. 96:460. https://doi.org/10.1094/PDIS-10-11-0886 Link, 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, Google ScholarJeck, W. R., and Sharpless, N. E. 2014. Detecting and characterizing circular RNAs. Nat. Biotechnol.. 32:453-461. https://doi.org/10.1038/nbt.2890 Crossref, Google ScholarJiang, S.-L., Yin, Q.-X., Li, D.-X., Wu, X., Wang, Y., Wang, D.-L., and Chen, Z. 2021. Integrated mRNA and small RNA sequencing for analyzing tea leaf spot pathogen Lasiodiplodia theobromae, under in vitro conditions and the course of infection. Phytopathology 111:882-885. https://doi.org/10.1094/PHYTO-07-20-0297-A Link, 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., Paggi, J. M., Park, C., Bennett, C., and Salzberg, S. L. 2019. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 37:907-915. https://doi.org/10.1038/s41587-01

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

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

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 Leaf Spot Pathogen Didymella Segeticola and Its Host, Tea (camellia Sinensis), During Infection

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

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

Functional Annotation of Circrnas of Tea Leaves During Infection by the Tea Leaf Spot Pathogen Didymella Segeticola

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

Sequencing and Functional Annotation of Competing Endogenous RNAs and MicroRNAs in Tea Leaves During Infection by Lasiodiplodia Theobromae

Identification of Cerna-Vsirna-mrna Network for Exploring the Mechanism Underlying Pathogenesis of Sugarcane Mosaic Virus in Resistant and Susceptible Maize Inbred Lines

MicroRNAs and Their Targeted Genes Associated with Phase Changes of Stem Explants During Tissue Culture of Tea Plant

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

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

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

Whole Genome Sequences of the Tea Leaf Spot Pathogen Didymella Segeticola.

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

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

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

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

Whole Genome Sequence and Gene Annotation Resource for Didymella Bellidis, Associated with Tea Leaf Spot.

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