Candidate Resistant Genes of Sand Pear (Pyrus pyrifolia Nakai) to Alternaria alternata Revealed by Transcriptome Sequencing

Pear black spot (PBS) disease, which is caused by Alternaria alternata (Aa), is one of the most serious diseases affecting sand pear (Pyrus pyrifolia Nakai) cultivation worldwide. To investigate the defense mechanisms of sand pear in response to Aa, the transcriptome of a sand pear germplasm with differential resistance to Aa was analyzed using Illumina paired-end sequencing. Four libraries derived from PBS-resistant and PBS-susceptible sand pear leaves were characterized through inoculation or mock-inoculation. In total, 20.5 Gbp of sequence data and 101,632,565 reads were generated, representing 44717 genes. Approximately 66% of the genes or sequenced reads could be aligned to the pear reference genome. A large number (5213) of differentially expressed genes related to PBS resistance were obtained; 34 microsatellites were detected in these genes, and 28 genes were found to be closely related to PBS resistance. Using a transcriptome analysis in response to PBS inoculation and comparison analysis to the PHI database, 4 genes (Pbr039001, Pbr001627, Pbr025080 and Pbr023112) were considered to be promising candidates for sand pear resistance to PBS. This study provides insight into changes in the transcriptome of sand pear in response to PBS infection, and the findings have improved our understanding of the resistance mechanism of sand pear to PBS and will facilitate future gene discovery and functional genome studies of sand pear.


Introduction
Pear is one of the most important fruit crops in the world and has a long history of commercial cultivation in approximately 50 temperate countries [1][2][3][4]. However, pear diseases have become a major constraint on the pear industry worldwide [5]. Pear black spot (PBS) disease, which is caused by Alternaria alternata (Aa), is one of the most serious diseases in sand pear. PBS disease causes both necrosis on pear leaves, twigs, and fruits and early leaf drop, and it reduces productivity and the quality of fruit. PBS disease decreases pear yields in areas with high temperatures and humidity. The control of PBS disease is mainly achieved through the use of fungicides; however, large quantities of fungicides can produce potential toxic effects on humans and wildlife and lead to environment pollution [6,7]. Moreover, the practice of applying large amounts of fungicides has increased the difficulty and cost of producing high-quality fruit [8]. Breeding disease-resistant varieties is the most effective and economical method of controlling the disease, although traditional breeding methods are hampered by the long growth period of pear [9]. However, the drawbacks of conventional breeding can be overcome by using plant resistance genes (R-gene) and molecular breeding approaches [10].
R-genes play an important role in plant genetic resistance mechanisms. A large number of R-genes have been isolated from several plant species, including rice, wheat, soybean, apple and grape [11][12][13][14][15]. Previous studies have attempted to isolate R-genes in pear, and four putative QTLs for fire-blight resistance were identified via the European pear genetic linkage map [16]. In addition, two genetic linkage maps of the pear scab resistance gene Vnk region of the cultivar Kinchaku have been constructed [17]. A leucine-rich repeat receptor-like protein kinase (LRPK) gene with possible involvement in scab resistance was identified from Pyrus pyrifolia cv. Kousui using reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends PCR [18]. Approximately 100 disease-resistance gene candidates were identified from six species of pear based on conserved nucleotide-binding sites and leucine-rich repeat domains in resistance genes using a homology-based cloning method [5]. However, R-genes specific for PBS disease have yet to be cloned in pears. Thus, new approaches in pears, such as transcriptome sequencing, are required to identify an R-gene for PBS disease.
Transcriptome sequencing has proven to be an effective approach for analyzing changes in developmentally and environmentally induced gene expression at the transcriptome level [19]. Transcriptomic information has been used in a wide range of biological studies and provides fundamental insights into biological processes and applications such as gene expression levels [20], gene expression functional annotation during development or after experimental treatments [21,22], gene discovery [23], SSR mining and SNP discovery [24]. This information can be used to help predict the roles and interactions of individual genes, discover more complex signal pathways activated in response to external stimuli and uncover potential cross-talk between these pathways. The latest paired-end tag sequencing strategy has further improved DNA sequencing efficiency and expands short-read lengths, thus providing a better depiction of the transcriptome [19].
In this study, differentially expressed genes (DEGs) of the pear in response to PBS infection were analyzed by paired-end transcriptome sequencing and available pear genome sequences (http://peargenome.njau.edu.cn:8004/default.asp?d=1&m=1) [25]. The objectives of this study were to screen R-genes that are differentially expressed in PBS-resistant and PBS-susceptible sand pear germplasms in response to Aa infection. The new dataset generated in this study will be a useful resource for future genetic and genomic studies of sand pear.

Plant Material and PBS Inoculation
Two sand pear germplasms that differ greatly in their resistance to Aa were used in this study. Pyrus pyrifolia Nakai. cv. 'Jinjing' pear (J), a cultivar with high resistance to Aa that was identified at the Wu Chang Sand Pear Garden National Fruit-tree Germplasm Resource (WCSPGNFGR) through three consecutive years (2009, 2010 and 2011) of field artificial inoculation experiments, was bred at the Research Institute of Fruit and Tea, Hubei Academy of Agricultural Science. Pyrus pyrifolia Nakai. cv. 'Hongfen' pear (H), a landrace germplasm resource with a high susceptibility to Aa, was also identified in the same field artificial inoculation experiments.
For the inoculations, inoculum production was conducted using a PBS strain (H) obtained from leaves of the Xiangnan pear from the germplasm bank of WCSPGNFGR. The H strain could infect the 'Hongfen' pear but not the 'Jinjing' pear. Stock cultures of the strain were stored at 4°C and sub-cultured in potato-sucrose-agar (PSA) medium at 28°C. After three weeks, 5 mL of sterilized distilled water was added to the surface of the colonies, which were scraped off using a small brush to remove the mycelium, and the conidial suspensions were filtered through four layers of cheesecloth. The conidial concentrate was centrifuged once for 20 min at 6000 g to remove mycelial fragments.
The conidial concentration was adjusted to 1 × 10 6 spores mL -1 . The spore suspension was sprayed onto detached young leaves of 'Jinjing' pear and 'Hongfen' pear with a glass atomizer. Control leaves of 'Jinjing' pear and 'Hongfen' pear were similarly sprayed with distilled water. The inoculated leaves were incubated in a moist chamber at 28°C. After inoculation with PBS, the leaves of 'Jinjing' pear and 'Hongfen' pear were sampled at 0 day, 1 day, 2 days, 3 days and 4 days. 'Jinjing' and 'Hongfen' pear leaves that were inoculated with distilled water were also sampled on these days.

Sample Preparation, Read Alignment, and Sequence Analysis
Total RNA was extracted from each sample using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. RNA purity was determined using a Nanodrop spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, DE, USA), 1% formaldehyde gel electrophoresis, and a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The total RNA of each sample was diluted to 750 ng μl -1 . The 200 μl mixture of RNA ('Jinjing' pear sample, J-P; and 'Hongfen' pear sample, H-P) used for sequencing included 40 μl of sample RNA of pear leaves inoculated with PBS at 0 day, 1 day, 2 days, 3 days and 4 days. Additionally, the 200 μl mixture of RNA ('Jinjing' pear sample, J-CK; and 'Hongfen' pear sample, H-CK) used for sequencing included 40 μl of sample RNA of pear leaves inoculated with distilled water at 0 day, 1 day, 2 days, 3 days and 4 days. The sample numbers H-CK, H-P, J-CK and J-P were used in the experimental design and data analysis.
The purity of the four RNA mixtures was determined using a Nanodrop spectrophotometer, 1% formaldehyde gel electrophoresis, and a 2100 Bioanalyzer. Oligo (dT)-coated beads were used to isolate poly (A) mRNA after the four RNA mixtures were collected. A fragmentation buffer (Ambion, Austin, TX, USA) was added to digest the mRNA to produce short fragments. First-strand cDNA synthesis was performed using random hexamer primers, and secondstrand cDNA synthesis was performed using buffer, dNTPs, RNaseH and DNA polymerase I. Short fragments were purified with the QIAquick PCR Purification kit (Qiagen, Valencia, CA) for the end repair and poly (A) addition reaction. Subsequently, sequencing adapters were connected to the short fragments; the cDNA library was then constructed and purified. Finally, the sequencing library was constructed by PCR amplification and sequenced by Biomarker Company using an Illumina HiSeq 2000 sequencing platform.
Sequence reads were aligned using Bowtie (http://bowtie-bio.sourceforge.net), an ultrafast short-read mapping program, using the pear genome as a reference [25]. TopHat (http://tophat. cbcb.umd.edu) was used to identify the splice junctions, and raw digital gene expression data were normalized as reads per kilobase pair of transcript per million mapped reads (RPKM), with genes that had P-values <0.001 selected for further analysis. Expression data were log 2 transformed and filtered at differences of 2-fold or greater in expression for each mixture sample. Differential patterns of gene expression from the various samples were represented by scatter diagrams.

Sequence Annotation
The optimal assembly results were selected according to the assembly evaluation, and clustering analysis was performed to produce a DEG database composed of potential alternative splicing transcripts. The DEGs were annotated by alignment with those deposited in diverse protein databases, including the National Center for Biotechnology Information (NCBI) nonredundant protein (Nr) database, the NCBI non-redundant nucleotide sequence (Nt) database, the Swiss-Prot database, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, the Gene Ontology (GO) database, the Cluster of Orthologous Groups of proteins (COG) database and the TrEMBL database.

SSR Mining and Primer Design
MIcroSAtellite (MISA) (http://pgrc.ipk-gatersleben.de/misa/) was employed for microsatellite mining. In this study, SSRs were considered to contain motifs with at least 10 contiguous repeats of one nucleotide, 6 repeats of two nucleotides, and 5 repeats of 3 to 6 nucleotides. Based on the MISA results, Primer3 software (http://primer3.sourceforge.net) was used to design the SSR primers, with the PCR product size ranging from 250 to 400 bp. In total, sixtyeight pairs of random SSR primers were used to analyze polymorphisms between PBS-resistant cultivar 'Jinjing' pear and PBS-susceptible 'Hongfen' pear cultivar.

Validation of RNA-Seq Data by Quantitative Real-Time RT-PCR
Total RNA was reverse transcribed into first-strand cDNA using an M-MLV First Strand Kit (Invitrogen) according to the manufacturer's instructions. Twenty-six genes were selected for confirmation by qRT-PCR with SYBR Premix Ex Taq (Takara, Japan). Primers for the chosen genes were designed using the Primer 6.0 program (PREMIER Biosoft International, Canada). The qRT-PCR gene expression analysis was performed on a StepOne real-time PCR System (Applied Biosystems, USA) using the glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene as an endogenous control [26]. Briefly, each reaction (final volume, 20 μl) contained 10 μl of 2 × TransStart Top Green qPCR SuperMix (TransGen, Beijing), 0.5 μl of passive reference dye II (50×), 0.4 μl (each) of forward and reverse primers (10 mM), 2 μl of cDNA template (corresponding to 50 ng total RNA), and 7 μl of RNase-free water. The reaction mixtures were heated to 94°C for 30 s, followed by 40 cycles at 94°C for 5 s, 55°C for 15 s, and 72°C for 34 s. A melting curve was generated for each sample at the end of each run to monitor the purity of the amplified products. Gene expression levels were analyzed using StepOne Software v2.0. The relative expression levels of target genes were calculated using the 2 −ΔΔCT method [27]. All of the samples were tested in triplicate, and the experiments were performed on three biological replicates.

Symptoms in Leaves after PBS Inoculation
The 'Jinjing' and 'Hongfen' pear cultivars showed different symptoms in response to PBS strain (H) after inoculation. 'Jinjing' pear leaves inoculated with strain H showed negligible necrosis on days 0, 1, 2, 3 and 4. However, 'Hongfen' pear leaves inoculated with strain H presented small brown necrotic spots on day 1. Leaf necrosis became more serious over time and developed into large necrotic lesions on day 4. 'Jinjing' and 'Hongfen' pear leaves inoculated with distilled water did not present necrotic spots on any day (Fig 1).

Sequencing and Assembly of Illumina Short Reads
Four libraries (H-CK, H-P, J-CK and J-P) were analyzed by RNA sequencing (RNA-seq). H-CK and H-P RNA were derived from the leaves of 'Hongfen' pear inoculated with distilled water and the H strain, respectively. J-CK RNA and J-P RNA were derived from the leaves of 'Jinjing' pear inoculated with distilled water and the H strain, respectively. More than 25 million 100-bp paired-end reads, encompassing over 5 Gb of sequence data, were generated for each library ( Table 1). The data sets supporting the results of this study are available at the NCBI SRA repository under accession number SRP051914.
Approximately 66% of the sequenced reads (33 million reads) were successfully mapped to the pear genome reference sequence using Bowtie and TopHat software. Of these, 87.23% were uniquely mapped reads and 12.77% were multiple mapped reads in the H-CK sample ( Table 2).

Transcriptome Analysis in Response to PBS Inoculation
The DEGs were analyzed using EBSeq software (Version 1.1.3) [28]. Using transcriptome sequencing, 5213 DEGs related to PBS resistance were obtained. A total of 909 DEGs were obtained between H-P and H-CK, and 501 DEGs were obtained between J-P and J-CK. A greater number of genes showed down-regulation following inoculation (Table 3), whereas a greater number of DEGs were screened between different varieties with either the water or PBS strain inoculation treatment. For example, 3460 DEGs were detected between H-CK and J-CK, whereas 3305 DEGs were obtained between H-P and J-P (Table 3). Scatter diagrams of gene expression are shown in S1 Fig.  The DEGs were grouped into distinct clusters based on expression patterns, and the results indicated that 3071 (88.8%) of the 3460 DEGs between H-CK and J-CK had significant matches in the Nr database, 3015 (87.1%) had significant matches in the Nt database, 2359 (68.2%) had significant matches in the Swiss-Prot database, 3071 (88.8%) had significant matches in the TrEMBL database, 2612 (75.5%) had significant matches in the GO database, 617 (17.8%) had significant matches in the KEGG database, and 1267 (36.6%) had significant matches in the COG database ( Table 4).
The COG database classifies orthologous gene products. Every protein in COG is assumed to be evolved from an ancestor protein, and the entire database is built on coding proteins with complete genomes. The DEGs were aligned to the COG database to predict and classify possible functions. Among the 25 COG categories, the cluster for 'General function prediction only' represented the largest group, followed by 'Transcription', 'Replication, recombination and repair' and 'Signal transduction mechanisms'. The cluster for 'Cell motility' represents the smallest group (S2, S3, S4 and S5 Figs).
Furthermore, GO functional annotation suggested that the DEGs can be categorized into 56 functional terms. A number of DEGs between the inoculation treatment and the water treatment in the same species were categorized into the 'cell killing' group. For example, 6 genes were categorized into the 'cell killing' group between H-P and H-CK, and 2 genes were categorized into the 'cell killing' group between J-CK and J-P (S6, S7, S8 and S9 Figs).
The KEGG PATHWAY database records networks of molecular interactions in these cells as well as their variants specific to particular organisms. Pathway-based analyses help to further elucidate the biological functions of genes. KEGG PATHWAY enrichment analysis showed that the most abundant categories in our analysis included 'Glycolysis/Gluconeogenesis,' 'Phenylalanine metabolism,' 'Plant hormone signal transduction,' 'Plant-pathogen interaction,' 'Ribosome,' and 'Spliceosome' (S10, S11, S12 and S13 Figs). These pathways might be significant in the resistance to Aa processes, especially in the 'Plant-pathogen interaction' metabolic pathway. Of the DEGs that participated in the 'Plant-pathogen interaction' metabolic pathway, 19 were between samples H-CK and J-CK, 23 were between samples H-P and J-P, 4 were between samples H-CK and H-P, and 1 was between samples J-CK and J-P.
We also compared the DEGs derived from 'Hongfen' (H-CK/H-P) with the DEGs from 'Jinjing' (J-CK/J-P), and 152 DEGs were observed for both varieties with PBS inoculation. These DEGs could represent the common genes induced by the PBS H strain. The DEGs derived from different varieties within the same treatment were compared, and 1987 DEGs were detected from the water treatment (H-CK/J-CK) and PBS inoculation (H-P/J-P) (Fig 2).
Resistance induced by the pathogen and the resistance of sand pear cultivars were the focus of the data analysis. The 152 DEGs derived from 'Jinjing' (J-CK/J-P) compared with the DEGs from 'Hongfen' (H-CK/H-P) and the 1987 DEGs derived from the water treatment (H-CK/ J-CK) and PBS inoculation (H-P/J-P) were categorized into 33 functional terms that were probably related to plant resistance to disease. For example, 11 genes were categorized into defense response to fungus (GO: 0050832), 13 were categorized into intracellular membranebounded organelle (GO: 0043231) and 4 were categorized into hydrolase activity, hydrolyzing O-glycosyl compounds (GO: 0004553) among the 152 DEGs (Fig 3). Furthermore, among the 1987 DEGs, 81 genes were categorized into defense response to fungus (GO: 0050832), 115 were categorized into intracellular membrane-bounded organelle (GO: 0043231) 14 genes were categorized into hydrolase activity, hydrolyzing O-glycosyl compounds (GO: 0004553) (Fig 4).

SSR Mining from the Sand Pear Transcriptome
In this study, 34 microsatellites were detected in the 152 DEGs derived from 'Jinjing' (J-CK/ J-P) compared with the DEGs from 'Hongfen' (H-CK/H-P), and 21 (61.8%) gene sequences contained 1 motif. The microsatellites included 9 (26.5%) dinucleotide motifs, 1 (2.9%) trinucleotide motif, and 3 (8.8%) compound motifs. A total of 68 primer pairs were successfully designed based on these 34 SSRs, 11 of which were able to distinguish between 'Jinjing' and  'Hongfen' pear (S1 Table). As the markers were developed from the 152 DEGs, which may represent genes related to resistance or susceptibility traits, these markers could be related to candidate genes involved in PBS resistance or susceptibility and could be useful for future efforts to breed resistance.  Table).

Analysis of the Sand Pear Resistance Gene to Aa
In addition, the 28 DEGs were annotated based on the GO and COG databases. The GO categories of the 28 genes showed that 'differences in response to stimulus', 'signal transduction' and 'immune process and death process' were the top 4 significantly enriched categories ( Fig  5). The classification was distinguished by orthologous gene products, and the functionality of each category was annotated using the COG database. Seven genes (Pbr023112, Pbr023278, Pbr025080, Pbr033741, Pbr039001, Pbr040608 and Pbr008283) were annotated with COG functions. Pbr025080 was associated with transcription, and the protein product of the other genes belonged to the unknown functional annotation in the COG database (S14 Fig). To further demonstrate the relationship between the 28 DEGs and PBS resistance, we identified metabolic pathways represented by the DEG collection. Annotations of 28 DEGs were fed into the KEGG Pathway Tools, which is an alternative approach to categorizing genes functions that emphasizes metabolic pathways. This process predicted the metabolic pathways represented by the 28 DEGs. An enrichment analysis of the 28 genes showed that the Pbr039001 gene encoded the disease resistance protein RPM1, which was significantly enriched in the plant pathogen interaction pathway (P-value = 0.0328) (Fig 6).   The Pathogen Host Interaction database (PHI-base) contains molecular and biological information on genes shown to affect the outcome of pathogen-host interactions (http://www. phi-base.org/). Twenty-two genes were successfully aligned to PHI-base R-gene sequences in the 28 genes related to PBS resistance, and the scores of 4 genes (Pbr001627 is 211, Pbr025080 is 193, Pbr023112 is 112, and Pbr025376 is 102) were higher than 100 (Table 5).

RNA-Seq Expression Validation by qPCR
To quantitatively determine the reliability of our transcriptome data, we monitored the expression of 26 DEGs by qRT-PCR. These genes included 22 DEGs related to PBS resistance and 4 regulatory proteins (Pbr020071 represents a RING finger protein, Pbr037418 represents a jasmonate-zim-domain protein, Pbr040066 represents a small heat shock protein, and Pbr042781 represents an NADP-dependent D-sorbitol-6-phosphate dehydrogenase) (S3 Table). The results of the qRT-PCR analysis revealed expression patterns that were consistent with the Illumina sequencing, although quantitative differences in expression levels were observed. The correlation between the RNA-seq and qRT-PCR results was measured by constructing a scatter plot of the log2 fold changes, and a positive correlation coefficient was observed (Pearson coefficient R 2 = 0.732) (Fig 7).

Discussion
The pear (P. bretschneideri Rehd. cv. 'Dangshansuli') genome, which encodes 42812 proteincoding genes, was released in 2013 [25]. In the present study, of the 44717 genes identified by transcriptome sequences, 41497 (92.8%) were successfully aligned to the pear genome A total of 396 nucleotide-binding site (NBS)-containing R genes were identified in the pear genome [25], and 4 new NBS-containing R genes (pear_newGene_2235, pear_newGene_1053, pear_newGene_1262 and pear_newGene_2229) were identified by transcriptome sequencing in the present study.
Pathogens proliferate within plant tissues through the secretion of effector virulence-associated proteins, RNAs and metabolites that down-regulate basal defenses. In turn, plants employ myriad defense mechanisms to prevent pathogen proliferation [29]. Plants have evolved sophisticated surveillance systems to recognize pathogen effectors delivered into host cells. In this study, a KEGG enrichment analysis showed that the Pbr039001 gene encoded disease resistance protein RPM1. RPM1 participates in plant immune responses, is an NBS-LRR immune receptor that recognizes the Pseudomonas syringae effectors AvrB and AvrRpm1, and may guard against pathogens by using AvrRpm1 and AvrB to manipulate RIN4 activity [30,31]. Pathogens inject pathogenic proteins directly into plant cells through the secretion system of the plant cell wall. However, plants are capable of activating a hypersensitivity response through the RPM1 protein stress response, and programmed cell death in the plant is employed to inhibit pathogen growth.
Based on the RNA-seq and qRT-PCR analyses, 16 genes were divided into 4 groups (Fig 8). The first group contained 4 genes (Pbr023112, Pbr007974, Pbr012560 and Pbr022874) that were highly expressed in the resistant variety ('Jinjing' pear) but were minimally expressed in the susceptible variety ('Hongfen' pear). The fold difference in expression was more than 4, especially for Pbr023112, where the fold difference in expression was more than 20. The second group contained 4 genes (Pbr025080, Pbr000681, Pbr12606 and Pbr038352) that were more highly expressed in the resistant variety ('Jinjing' pear) than in the susceptible variety ('Hongfen' pear), and the fold difference in expression was between 2 and 3. The third group contained 4 genes (Pbr001627, Pbr008283, Pbr033741 and Pbr034022) that were more highly expressed in the resistant variety ('Jinjing' pear) than in the susceptible variety ('Hongfen' pear), and the fold difference in expression was less than 2. The fourth group contained 4 genes (Pbr039001, Pbr001247, pear_newGene_1262 and Pbr022889) that were differentially expressed in response to PBS inoculation, especially Pbr039001, whose fold difference in expression was more than 10.
Based on the GO functional annotation, COG function classification, KEGG analysis, and R gene comparison in the PHI-base of the 28 DEGs, 4 genes (Pbr039001, Pbr001627, Pbr025080 and Pbr023112) were screened as candidates with resistance to PBS. Pbr039001 mainly participated in biological plant-type hypersensitive response processes (GO: 0009626); Pbr001627 mainly participated in biological regulation of cellular processes (GO: 0050794), defense responses to bacterium (GO: 0042742), primary metabolic processes (GO: 0044238), innate immune responses (GO: 0045087), aromatic compound biosynthetic processes (GO: 0019438), and small molecule biosynthetic processes (GO: 0044283); Pbr025080 mainly participated in biological defense responses (GO: 0006952), cellular processes (GO: 0009987), and multiorganism processes (GO: 0051704); and Pbr023112 mainly participated in biological defense responses (GO: 0006952). The methods by which these 4 candidate genes function in the sand pear can be verified through transgenic experiments in future research.

Conclusions
In conclusion, this work presents the first transcriptome sequencing analysis of resistant sand pear cultivars and susceptible germplasm leaves inoculated with PBS using the Illumina platform. In total, 20.5 Gbp of sequence data and 101,632,565 reads were generated, and approximately 66% of the sequenced reads were successfully aligned to the pear genome reference sequence. In summary, 44717 genes were identified by transcriptome sequences in this study, 5213 DEGs related to PBS resistance were obtained, 34 microsatellites were detected in the DEGs, and 107525 reliable SNPs were identified through transcriptome sequencing. From the transcriptome analysis of PBS inoculation, 28 DEGs related to PBS resistance were identified, and 4 genes (Pbr039001, Pbr001627, Pbr025080 and Pbr023112) were screened as candidates for sand pear resistance to Aa. These data will facilitate gene discovery and functional genomic studies in sand pear, and the findings will improve our understanding of the resistance mechanisms of sand pear to PBS. Therefore, this study provides insight into the complex transcriptome of the sand pear and establishes a biotechnological platform for future research.