Systematic Analysis and Functional Validation of Citrus XTH Genes Reveal the Role of Csxth04 in Citrus Bacterial Canker Resistance and Tolerance

In this study, we performed a comprehensive survey of xyloglucan endotransglucosylase/hydrolase (XTH) and a functional validation of Citrus sinensis (Cs) XTH genes to provide new insights into the involvement of XTHs in Xanthomonas citri subsp. citri (Xcc) infection. From the genome of sweet orange, 34 CsXTH genes with XTH characteristic domains were identified and clustered into groups I/II, IIIA, and IIIB. Except for chromosome 9, the CsXTH genes were unevenly distributed and duplicated among all chromosomes, identifying a CsXTH duplication hot spot on chromosome 4. With Xcc induction, a group of citrus canker-related CsXTHs were detected. CsXTH04 was identified as a putative candidate gene, which is up-regulated in citrus bacterial canker (CBC)-resistant varieties and induced by exogenous treatment with salicylic acid (SA) and methyl jasmonate (MeJA). CsXTH04 overexpression conferred CBC susceptibility to transgenic citrus, while CsXTH04 silencing conferred CBC resistance. Taken together, the annotation of the CsXTH family provides an initial basis for the functional and evolutionary study of this family as potential CBC-susceptible genes. CsXTH04, validated in this study, can be used in citrus breeding to improve CBC resistance.

Recently, XTHs have been shown to regulate the plant response to exogenous stress (Cho et al., 2006;Yang et al., 2011). For instance, AtXTH31 in roots is down-regulated by aluminum overexposure (Yang et al., 2011), and the overexpression of CaXTH3 improves salt and drought tolerance in transgenic Arabidopsis (Cho et al., 2006).
Along with the increased availability of the plant genome, the XTH family has been annotated in many plants and extensively researched (Yokoyama and Nishitani, 2001;Behar et al., 2018). In tomato, 56 XTH genes have been identified (Wang et al., 2018), 29 have been identified in rice, and 41 have been annotated in Populus (Yokoyama et al., 2004;Geisler-Lee et al., 2006;Saladié et al., 2006). A comprehensive cross-genome survey and phylogeny of the GH16 members revealed the evolutionary origin of EG16 and XTH proteins in a series of plant lineages including two citrus species (Behar et al., 2018). In the study, XTHs were extracted from genomic data (Wu et al., 2014) of Citrus clementina and Citrus sinensis based on Basic Local Alignment Search Tool (BLAST) strategy, which does not represent a specialized annotation process. The genomic data of C. sinensis used in studies from Behar and colleagues are yet to be assembled at the chromosome level, which brings uncertainty to subsequent research. In this study, we exhaustively annotated the CsXTHs from chromosome-level genomic data of C. sinensis (Xu et al., 2013) and performed an in-depth functional analysis.
Citrus bacterial canker (CBC) is caused by Xanthomonas citri subsp. citri (Xcc) (Schaad et al., 2005), which causes severe yield losses in citrus-producing regions globally. The cell wall is the initial barrier of the plant defense against bacterial infection and other stresses. In our previous long-term studies regarding the citrus transcriptomes induced by Xcc, XTH genes were identified as differentially expressed genes (DEGs) (Supplementary Table 1). In this study, we performed an in silico comprehensive detailed annotation and functional analysis of XTH genes in response to CBC by quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR) and reverse genetics.

Comprehensive Annotation of XTHs in the Genome of Sweet Orange
For data mining, the proteome and genome of sweet orange were downloaded from C. sinensis annotation project (CAP) database (Xu et al., 2013;Wang et al., 2014) and Phytozome (Goodstein et al., 2012;Wu et al., 2014). A total of 33 XTH sequences of Arabidopsis thaliana were also collected from the Phytozome (Yokoyama and Nishitani, 2001;Goodstein et al., 2012). Exhaustive data mining and annotations were performed through a detailed three-step semi-automatic process that avoids the errors of automatic predictions (Fawal et al., 2014;Li et al., 2015). During the annotation process, Fgenesh++ (Solovyev et al., 2006), Pfam (El-Gebali et al., 2018 and Simple Modular Architecture Research Tool (SMART) (Letunic and Bork, 2018) were used for functional and structural re-annotation, while Scipio (Keller et al., 2008) was used to retrieve mis-annotated XTH genes in the automatic process and its corresponding chromosomal positions, gene structures, and sequences for all XTHs annotated. The records of XTHs in sweet orange were initially named "CsXTH" and followed by a number indicating its chromosomal order. EG16 (endoglucanases), a member of C. sinensis (CsEG16), was also annotated and analyzed with CsXTHs due to its relevance to CsXTHs (Behar et al., 2018). The EG16 member of C. clementina (Behar et al., 2018) was used as the query to survey the CsEG16s.

Plant, Bacterial Materials, and Treatments
All plant materials were sampled from the National Citrus Germplasm Repository, Chongqing, China. Varieties Wanjincheng (C. sinensis) and calamondin (Citrus madurensis) were used for CBC and exogenous hormone assays. Wanjincheng was used for gene transformation. All plants were planted in a greenhouse at 28°C. An orange cultivation site in the Yunnan province, China, was the source of a variant of Xcc termed XccYN1 derived from sweet orange leaves that displayed a natural infection. Peptone-yeast extract-malt extract (PYM) with d-glucose 1.5% (w/v) and antibiotics at 28°C was used for the culture of Xcc bacteria.

CBC and Exogenous Hormone Assay
For the analysis of the expression patterns of the CsXTHs, uniform and healthy fresh leaves excised from the plants were placed in culture plates containing sterile deionized water and maintained at 28°C with a 16-h light/8-h dark photoperiod. A 1,000-fold dilution of Xcc (OD600 = 0.5) was used to inoculate Wanjincheng and calamondin leaves in a 16-h light/8-h dark incubator at 28°C. Samples were collected at 0, 12, 24, 36, and 48 h postinoculation (hpi). To perform exogenous hormone assays, leaf discs were soaked in 10 μmol/L of salicylic acid (SA) or 100 μmol/L of methyl jasmonate (MeJA) and collected at 0, 6, 12, 24, 36, and 48 h posttreatment (hpt). Collected samples were assessed by qRT-PCR.

Subcellular Localization Analysis
Subcellular localizations were analyzed with software prediction and transient expression. SignalP 4.0 (Petersen et al., 2011) and CELLO (Yu et al., 2006) were used for signal peptide and subcellular localization predictions, respectively. The coding sequence (CDS) of CsXTH04 lacking a stop codon was amplified with primers F-SC (GGGGTACCATGGCTT CTTTTCTATGGACTCT) and R-SC (TCCCCCGGGAATGTC ACGGTCTCGTTTGCA) and inserted into the pLGNe-GFP vector to construct the transient expression vector. Agrobacterium tumefaciens EHA105 was used to infect onion epidermal cells, and green fluorescent protein (GFP) fluorescence was assessed after 48 h.

Overexpression and Gene Silencing of CsXTH04
For overexpression plasmid construction, fulllength CDS of CsXTH04 was amplified with primers (restriction sites included) F-OEc (CGGGATCCATGGC TTCTTTTCTATGGACT) and R-OEc (CGGAATTCTTAAAT GTCACGGTCTCGTTTGCA) and inserted into pGLNe. For silencing vector construction, a 216-bp fragment was amplified with primers (restriction sites included) F-RIc ( G C T C TAG AG G C G C G C C T T C ATAG C C T C C TAC A AGGGGT) and R-RIc (CGGGATCCATTTAAATGCATTC TGGTGAAGGCGTAGGGA) and integrated into the PUC-RANi vector. RNAi sequence including forward, intron, and reverse sequences were isolated and inserted into pLGNe to produce the final vector.

Citrus Transformation and the Transgenic Characterization of Plants
Overexpression and silencing plasmids were introduced into A. tumefaciens EHA105 by heat shock. The shoot segments of Wanjincheng (C. sinensis) were transformed with A. tumefaciens using the Peng method (Peng et al., 2017). PCR and GUS (beta-glucuronidase) assays were used to confirm the presence of transgenic genes. Primers F-OEd (CGACACGCTTGTCTA CTCCA) and R-OEd (TTAAATGTCACGGTCTCGTT TGCA) were used for overexpression plants; and F-RId (CT TCACCAGAATGCATTTAAATGTGTAA) and R-RId (GG TCTTACGGGATCCAAATACCTGCAAA) were used for silencing plants. The GUS activity in transgenic plants was detected using histochemical procedures (Sendín et al., 2017). Briefly, the leaf discs of transgenic citrus were infiltrated in X-GLUC solution (5-bromo-4-chloro-3-indolyl-betad-glucoronide) and incubated overnight at 37°C. Leaf discs were cleared by soaking in 70% ethanol and imaged by microscopy. After genomic PCR and GUS assays, qRT-PCRs were performed to analyze CsXTH04 expression in the positive plants. Wild-type (WT) plants were taken as controls for PCR verification, GUS assays, and qRT-PCR analysis.

Assessment of Transgenic Plant Resistance Against Xcc
In vitro assays were performed to detect the Xcc resistance of transgenic plants as previously described (Peng et al., 2017). Six fully mature healthy leaves were used per transgenic plant. Briefly, six punctures were made in each of the leaves using a 0.5-mm-diameter pin, after which 1 µL of XccYN1 bacterial suspension (1 × 10 5 CFU/mL) was inoculated at the sites. Images of the diseased areas were taken at 10 days postinoculation (dpi) and measured on ImageJ 2.0 software (National Institutes of Health, Bethesda, MD). The lesion sizes (LSs) of diseased spots and the disease index (DI) were used to evaluate Xcc resistance. The DI was calculated as previously described (Peng et al., 2017).

RNA Isolation, cDNA Synthesis, and qRT-PCR Assay
Frozen tissues were ground in liquid nitrogen, and total RNA was extracted from leaf samples using miniprep kits purchased from AidLab Ltd (Beijing, China) based on the manufacturer's protocols. RNA was reverse transcribed with TaKaRa kits (Dalian, China). The QuantStudio 7 (Applied Biosystems, USA) system and SYBR Premix ExTaq Green PCR kit (Bio-Rad, USA) were used for qRT-PCR with the actin gene of citrus (GenBank accession: GU911361.1) used as the internal control. Primer sequences were F-actin (GTTGCAGCAAT  GCCAGTGAA) and R-actin (GCGGCAGATGT GTTTTGTGT). The qRT-PCR parameters were as follows: 95°C for 5 min, followed by 40 cycles at 95°C for 10 s, 56°C for annealing, and extension for 30 s. A reaction mixture of 20 μL was composed of 100-ng cDNA, 0.5 μM of primers, and 10 μL of SYBR Green PCR mix. The 2 −∆∆CT method was used to evaluate the relative expression levels of each examined gene (Livak and Schmittgen, 2001). The primers for qRT-PCR were designed with Primer BLAST (Supplementary Table 2). For each gene, three biological and technical replicates were used.

Statistical Analysis
SPSS V22 (SPSS, Inc., Chicago, USA) was used for all analyses. Gene expression differences were compared through an analysis of variance (ANOVA) via Fisher's least significant difference (LSD) tests. Differences were considered statistically significant when *P < 0.05 and **P < 0.01. Each value represents the mean ± SD.

Identification and Bioinformatics Analysis of CsXTHs
With exhaustive mining and detailed annotation, 34 XTH members were annotated and characterized, including 33 complete genes and a single pseudogene (CsXTH14) ( Table 1). Among the XTH family, all genes were predicted by CAP, while 31 genes were predicted using the Phytozome protocol. Twelve XTH genes that were incorrectly predicted were manually annotated with Fgenesh++ re-annotation (Solovyev et al., 2006) and EST (expressed sequence tag) supports (Table 1 and Supplementary Table 3). The members were named CsXTH01 to CsXTH34 based on the chromosomal order. To validate the annotation of CsXTHs, the EST database was downloaded from National Center for Biotechnology Information (NCBI) using the best EST hits. Finally, 17 genes were found with a total of 103 ESTs (Table 1 and Supplementary Table 3). The ESTs helped the annotation of incorrectly predicted genes and predicted gene functions. Information on the CsXTH family, including gene name, ID, chromosome location, molecular weight (MW), amino acid length, isoelectric point (PI), and annotation modes are detailed in Table 1.
The complete amino acid sequences of CsXTHs (Supplementary Table 4) and AtXTHs (Supplementary Table 5, Supplementary Files 1 and 2) were used for phylogenetic analysis. The rooted ML phylogenetic tree indicated that the CsXTHs can be divided into I/ II, IIIA, and IIIB groups (Figure 1) according to the clade identifiers used in AtXTHs (Yokoyama et al., 2010).
Exon-intronic analysis was performed for the CsXTH genes, which revealed no introns in three genes (CsXTH25, CsXTH26, and CsXTH31). In other genes, the intron numbers ranged from 2 (15 genes) to 3 (15 genes). The introns and exons were well conserved between closely related genes. For example, closely related CsXTH, CsXTH25, and CsXTH26 genes contained no introns (Supplementary Figure 1A). With respect to functional annotations, all CsXTHs excluding CsXTH25 and CsXTH26 possessed an N-terminal glyco hydro 16 domain   Figure 1B).

Chromosomal Distribution and Duplication Refer to a Duplication Hot Spot on Chromosome 4
The CsXTH genes were distributed among all chromosomes excluding chromosome 9 based on their positions in the genome of sweet orange (Figure 2A). Chromosome 4 contained the highest number of CsXTHs (17 members). By contrast, only a single CsXTH on chromosomes 1, 3, 5, and 7 was identified. Chromosome 4 showed the highest density of CsXTH genes (0.85 per Mb), while chromosomes 1 and 3 showed the lowest density. In addition, the distribution of CsXTH genes was uneven. A higher density of CsXTH was found in specific chromosomal regions, mainly on chromosome 4 (Figure 2A). To further understand how CsXTH genes evolved, gene duplication events were investigated in sweet orange. In this study, 10 pairs of duplications were arranged in blocks of duplicates, among which were a single pair of segmental duplications (SDs), six pairs of tandem duplications (TDs), and three pairs of whole-genome duplications (WGDs) (Figure 2A). These results strongly indicated that TD and WGD make major contributions to the expansion of the XTH family in sweet orange. According to the phylogenetic tree of CsXTHs (Figure 1) and the identified duplications, a large region on chromosome 4 was detected. In the 58,964-bp region (chromosome 04: 1564646 to 1623610), 11 CsXTHs were involved in the segment (CsXTH07 to CsXTH17) (Figure 2A). A phylogeny (ML) with 10 CsXTHs (pseudogene CsXTH14 was not used for phylogeny) showed five pairs of duplicated genes (CsXTH07-CsXTH12, CsXTH08-CsXTH13, CsXTH09-CsXTH15, CsXTH10-CsXTH16, and CsXTH11-CsXTH17) ( Figure 2B). The duplicated genes had the same directions on chromosomes, similar gene sizes, and the same number of introns except for CsXTH11-CsXTH17 ( Figure 2B). The duplicated genes in this region produced two duplicated segments ( Figure 2C). There was one unique case with an interval between CsXTH08 and CsXTH09 that differed from the interval between CsXTH13 and CsXTH15. We reasoned that these differences were due to the pseudogenization process of CsXTH14. In conclusion, a hypothesis for the evolution of the 10 CsXTH genes was constructed. One ancestral gene duplicated four times and formed the CsXTH07-to-CsXTH11 segment. A segment of CsXTH07 to CsXTH11 underwent TD and was inserted into a nearby loci on chromosome 4 to form the segment CsXTH12 to CsXTH17 (Figure 2D).

CsXTHs Are Involved in the Response to Xcc Infection
To further explore CsXTH functions, the expression of these XTHs together with CsEG16 during the induction of biotic stress was investigated. To detect the involvement of CsXTHs in Xcc infection, their expression at 0, 12, 24, 36, and 48 hpi was detected by qRT-PCR. From the qRT-PCR data, different expression profiles were detected for each gene (Figure 3). During Xcc infection, a group of CsXTHs with up-regulated expression profiles in the citrus canker sensitive variety Wanjincheng were identified including CsXTH03, CsXTH04, CsXTH12, CsXTH13, CsXTH16, CsXTH20, CsXTH21, CsXTH22, and CsXTH28. Among these genes, CsXTH04 and CsXTH21 showed a down-regulated profile in the citrus canker-resistant variety calamondin. This gene, therefore, represents an Xcc susceptibility gene. Based on our data mining and annotation together with Xcc-induced expression analyses, we thus identified CsXTH04 as a candidate gene for further CBC studies.

CsXTH04 Is an Extracellular Protein
CsXTH04 was identified by qRT-PCR as a candidate gene involved in Xcc infection (Figure 3). We, therefore, investigated the function of CsXTH04. The subcellular localization of CsXTH04 was determined by software predictions and exogenous transient expression systems. The protein contained a 23-aa signal peptide (MVVSYEGCFLLVFSLLAVVASGL) at the N-terminal, suggesting it to be secreted from cells. Potential extracellular loci were also identified (Supplementary Table 6). To validate these predictions, we examined the subcellular localization of CsXTH04 following transient expression of a recombinant pLGNe-CsXTH04-GFP plasmid ( Figure 4A). Both the cytoplasm and nucleus of the control exhibited green fluorescence ( Figure 4B), but the strongest GFP fluorescence was observed at the plasma membrane and cell wall of the epidermal onion cells with plasmid pLGNe-CsXTH04-GFP ( Figure 4C). Both predictions and transient expression analysis, therefore, indicated that CsXTH04 performs extracellular functions.

Csxth04 Expression Was Induced by SA and MeJA
Proteins associated with plant disease are often regulated by phytohormones (Zuo et al., 2015;He et al., 2019). To assess the involvement of CsXTH04 in disease resistance-related signaling pathways, SA and MeJA inductions were performed, and CsXTH04 expression was detected by qRT-PCR. The results showed that the expression of CsXTH04 in calamondin in response to exogenous SA increased ~14-fold during treatment, while in Wanjincheng, the expression was unaffected within the 48-h treatment period ( Figure  4D). Regarding MeJA induction, the expression of CsXTH04 also increased ~3-fold in calamondin. After 36 hpt, the expression levels decreased to below starting levels. Within the 48-h MeJA treatment period, CsXTH04 showed no significant differences in Wanjincheng ( Figure 4E). From these results, we concluded that in calamondin, CsXTH04 expression is be induced by both SA and MeJA, while in Wanjincheng, CsXTH04 expression could not be induced.

Overexpression of CsXTH04 Confers CBC Susceptibility to Transgenic Plants
Transgenic citrus plants were produced to aid our understanding of the role of CsXTH04. The CsXTH04 overexpression plasmid contained a GUS coding sequence and CaMV 35S promoter ( Figure 5A). Three plants (OE1 to OE3) with the transgenic integration of CsXTH04 were confirmed by genomic PCR and GUS assays. Compared with WT, OE1 to OE3 had a 1.8-kb signature fragment that was detected by PCR ( Figure 5B). The OE1 to OE3 showed a blue color on the edge of the leaf discs, while WT leaves showed no color ( Figure 5C). Following September 2019 | Volume 10 | Article 1109 Frontiers in Plant Science | www.frontiersin.org qRT-PCR detection, the three overexpressed plants possessed high levels of CsXTH04 (300-, 100-, and 55-fold higher than WT, respectively) ( Figure 5D). Regarding phenotypes, the transgenic plants exhibited normal growth rates than did WT plants ( Figure 5E). To study the CBC resistance of the overexpressed transgenic plants OE1 to OE3, in vitro assays were performed through the inoculation with pinpricks. The sizes of the lesions on OE-transgenic leaves were larger than those of WT plants ( Figure 5F). The pustules of infection by Xcc were exacerbated by the overexpression of CsXTH04, while the highest levels of susceptibility were displayed by OE1, followed by OE2 and OE3. Regarding the LS, OE1 possessed the largest lesions, which were ~2.36-fold higher than those of WT, while OE3 possessed the lowest, which were 1.39-fold higher than those of WT ( Figure 5G). DI, as a measure of disease severity, was enhanced by 36.7-87.8% in the transgenic lines compared with WT ( Figure 5H). We, therefore, concluded that the overexpression of CsXTH04 conferred Xcc suppressibility to transgenic citrus.

Silencing of CsXTH04 Confers CBC Resistance to Transgenic Plants
To further elucidate the role of Wanjincheng, CsXTH04 was silenced using RNAi approaches. RNAi sequences, including forward, intron, and reverse sequences, were digested and inserted into pLGNe (Figure 6A). To verify the transgenic plants, genomic PCRs were performed, and three silenced plants were obtained ( Figure 6B). The plants were subjected to verification by GUS stain ( Figure 6C). The three plants (R1, R2, and R3) exhibited relatively low levels of expression of CsXTH04 than did WT plants (41%, 40%, and 23% of WT), which was validated by qRT-PCR ( Figure 6D). Compared with WT, the three silenced plants exhibited normal growth states (Figure 6E). The use of inoculation by pinprick was applied to study the CBC resistance of the three silenced mutant plants. Transgenic plants R1 to R3 exhibited smaller-sized pustules than did WT ( Figure 6F). We, therefore, conclude that silencing CsXTH04 significantly enhances its resistance to Xcc infection. Statistically diseased lesions in the three silenced plants (R1, R2, and R3) showed smaller dimensions when examined against WT (58%, 57%, and 42%, respectively) ( Figure 6G). The analysis of CBC severity revealed that the three silenced mutant plants had significantly higher DI than WT ( Figure 6H). Consequently, disease severity decreased by 29% (R1) to 65% (R3). As such, CsXTH04 mutant's exhibit increased CBC resistance, suggesting that CsXTH04 is an Xcc susceptibility gene.

DISCUSSION
Xyloglucan is a hemicellulose polysaccharide that is a major constituent of the primary cell walls of dicots and some monocots (Vogel, 2008;Scheller and Ulvskov, 2010;Pauly and Keegstra, 2016). XTHs can cut and rejoin xyloglucan to regulate the organization and composition of the cell wall during plant development and stress responses, including abiotic and biotic stresses (Yokoyama and Nishitani, 2001). To date, numerous XTHs genes have been identified from plants (Behar et al., 2018) including 33 genes in Arabidopsis (Rose et al., 2002), 29 in rice (Yokoyama et al., 2004), 56 in tomato (Wang et al., 2018), and 44 in Medicago  (Xuan et al., 2016). Despite the XTH family of two Citrus species (C. clementina and C. sinensis) being surveyed from genomic data, the XTHs have not been annotated in a detailed manner. The genes extracted from genomics using automated BLAST-based protocols typically contain errors and/or redundancy (Fawal et al., 2014;Li et al., 2015). In this study, 12/34 XTH genes were incorrectly predicted and manually corrected ( Table 1). At the chromosome level, assembled genomic datasets (Xu et al., 2013) are in contrast to those described in previous studies (Behar et al., 2018).
The number of XTH genes in sweet orange was higher than that of Arabidopsis and rice but lower than that of tomato and Medicago, suggesting that this was a result of plant-specific gains and losses (Song et al., 2015). On chromosome 4 of sweet orange, there is a CsXTH hot spot with 10 duplicated genes. This duplication hot spot contains TDs combined with segmental duplications (Figure 2A). The up-regulation of CsXTH was subject to protein concentrations for specific functional requirements or functional differentiation . This duplication hot spot of CsXTHs represents research interests on future functional and evolutionary studies.
Plant cell walls play an important role in the regulation of environmental stresses. Plant XTHs are important cell wallmodifying enzymes that participate in cell wall extension and degradation, maintaining the integrity and strength of the cell wall under normal and stressful environments (Zhu et al., 2012). Compared with abiotic stresses, the relationship between XTHs and plant disease, particularly CBCs, is poorly defined. In our experience of studying the transcriptomes of species infected with Xcc, XTHs were consistently represented (Supplementary Table 1), indicating the involvement of XTHs in Xcc infection. To investigate the XTH genes that potentially respond to CBC, we performed qRT-PCR to detect the induced expression of CsXTHs. We found that CsXTH04 was involved in Xcc infection and SA/MeJA signaling pathways, which preliminarily validated their links to disease. The differential induction of Wanjincheng and calamondin provide insight into the different cis-elements between these two species (Song et al., 2016). Our present work highlighted the involvement of three CsXTH genes in CBC resistance, further extending the list of CBC-related genes, confirming a role for XTHs in pathogen infection. CsXTH04 was investigated in depth using overexpression and silencing strategies. We found that the overexpression of CsXTH04 conferred CBC suppressibility to transgenic citrus, while the silencing of CsXTH04 conferred CBC resistance (Figures 5 and 6). However, many questions remain unsolved and require further investigation. These include the functional mechanisms of CsXTH04 in CBC infection and how bacteria induce the expression of CsXTH04.
Further physiological and molecular/genetic studies of CsXTH04 must now be performed to clarify its role in the plant response to CBC stress and tolerance.

CONCLUSIONS
To our knowledge, this is the first focused annotation and expression profiling of XTHs in citrus in response to Xcc induction. In the present study, 34 XTHs from the genome of sweet orange were annotated, providing an initial basis for the functional investigation and evolutionary study of this family. The expression profiling of CsXTH induced by Xcc provides insight into the involvement and development of CBC. The CBC identified genes, including CsXTH04, hold value for citrus breeding to improve CBC resistance.

DATA AVAILABILITY
All the datasets generated for this study are included in the Supplementary Files.

AUTHOR CONTRIBUTIONS
QLi, SC and YH designed the experiments; QLi, AH, JQ, WD, YH, QLo, and XZ performed the experiments; QLi, TL and LY performed data analysis; SC and YH contributed materials and tools for analysis and reagents; QLi wrote this manuscript; All authors read and approved the final manuscript.