Abstract
The Asian citrus psyllid Diaphorina citri is the transmission vector of Huanglongbing (HLB), a devastating disease of citrus plants. The bacterium “Candidatus Liberibacter asiaticus” (CLas) associated with HLB is transmitted between host plants by D. citri in a circulative manner. Understanding the interaction between CLas and its insect vector is key for protecting citrus cultivation from HLB damage. Here, we used RNA sequencing and liquid chromatography-mass spectrometry (LC-MS) to analyze the transcriptome and metabolome of D. citri interacting with CLas. We identified 662 upregulated and 532 downregulated genes in CLas-infected insects. These genes were enriched in pathways involving carbohydrate metabolism, the insects’ immune system, and metabolism of cofactors and vitamins. We also detected 105 differential metabolites between CLas-infected and non-infected insects, including multiple nucleosides and lipid-related molecules. The integrated analysis revealed nine pathways—including those of the glycine, serine, threonine, and purine metabolism—affected by the differentially expressed genes from both groups. The network for these pathways was subsequently constructed. Our results thus provide insights regarding the cross-talk between the transcriptomic and metabolomic changes in D. citri in response to CLas infection, as well as information on the pathways and genes/metabolites related to the CLas–D. citri interaction.
Introduction
The Asian citrus psyllid Diaphorina citri Kuwayama (Hemiptera: Psyllidae) is the main citrus pest, and it is the most widely distributed insect vector of huanglongbing (HLB, or citrus greening disease). HLB is a serious threat for the citrus industry worldwide (; ; ) as it shortens plant life, sharply reduces fruit yield, and may result in plant death. HLB is present in most citrus-producing countries of Asia, Africa, and America (, ). More than 20% of citrus trees in Brazil and 90% in Florida are currently destroyed by HLB, and HLB is currently the most serious infectious disease affecting global citrus yield (; ).
HLB is associated with the gram-negative bacterium, Candidatus Liberibacter asiaticus (CLas). Diaphorina citri typically feeds on plants infected with CLas for some time and then migrates to healthy plants; this accelerates the spread of CLas (). Feeding on diseased plants has been shown to increase the fecundity and shorten nymph development duration in D. citri (; ). A study using an electrical penetration graph showed that the salivary time of D. citri adults feeding on diseased plants was significantly prolonged compared to that of D. citri not feeding on diseased plants (). Moreover, infected D. citri have been shown to migrate to healthy plant earlier and have the ability to fly father than uninfected ones, and as bacterial contents increased, infected females became more attractive to males (; ). CLas infection downregulated the gene expression level of a protease and its transporter complex in D. citri nymphs, and the immune system’s ability to synthesize and release immune-related proteases was decreased, allowing the reproduction of CLas (). In addition, CLas can affect the behavior, gene expression, and substance metabolism of D. citri, thus promoting its transmission and spread. However, the genetic mechanisms regulating these phenomena are largely unknown.
The draft genome sequence of D. citri is available in the National Center for Biotechnology Information (NCBI). Accordingly, recent studies have continuously provided more information on the transcriptomic and proteomic changes regarding the species’ hemolymph and midgut (; ; , ; ), revealing the effects of CLas infection on D. citri’s nutrition metabolism, nymphal development, and immunity. Most aforementioned previous studies have focused on specific tissues of D. citri due to their implication in the insect’s immune system. In the present study, we compared the whole-body transcriptomes and metabolomes between CLas-infected and non-infected D. citri, and we analyzed the DEGs and differentially expressed metabolites to better understand the interaction between CLas and D. citri.
Materials and Methods
Insect Material
The D. citri specimens used in this study were collected from Sun Yat-sen University, Guangzhou, China, in 2014. The strain was reared on Citrus reticulata Blanco cv. Shatangju in a greenhouse under a 14 h: 10 h light: dark cycle at 26 ± 1°C and 65–70% relative humidity for more than 25 generations. The newly emerged female adults were fed for 14 days on the citrus plants infected with HLB, and they were confirmed to carry CLas through sampling detection. Meanwhile, the female adults that was continuously fed on healthy citrus plants for 14 days as control.
Quantitative PCR Detection of CLas
Total RNA (1 μg) was reverse-transcribed into first-strand cDNA using a PrimeScriptTM RT reagent kit (Takara Bio, Inc., Otsu, Shiga, Japan). Quantitative PCR (qPCR) was performed using a 10 μL reaction containing 1 μL cDNA, 0.3 μL each of 10 μmol⋅L–1 forward and reverse primers (Supplementary Table S1), and 5 μL SYBR FAST Universal qPCR mix (KAPA Biosystems, Woburn, MA, United States) in a LightCycler 480 system (Roche Diagnostics GmbH, Mannheim, Germany), with the following amplification conditions: 5 min at 95°C, followed by 45 cycles at 95°C for 10 s, 60°C for 20 s, and 72°C for 20 s. The qPCR experiments were performed for each sample using three biological and three technical replicates. The expression levels of selected genes were normalized to the expression levels of D. citri actin-1. The differential gene expression was calculated using the 2–△△Ct method (), and the results were expressed as mean ± SE. The differences in gene expression levels between infected and uninfected D. citri were analyzed using t-tests in the SPSS 18.0 statistical software (P < 0.01).
RNA-Seq and Differential Expression Analysis
All 20 insect individuals were pooled and ground in liquid nitrogen. Total RNA was extracted using a total RNA extract kit (Omega, Norcross, GA, United States) according to the manufacturer’s instructions. The cDNA libraries from triplicates of each insect line were generated using the Illumina TruSeq Stranded mRNA kit (Illumina, San Diego, CA, Untied States) and were size-selected to an average fragment size of 450 bp. The libraries were sequenced on the Illumina HiSeq Xten platform producing paired-end (PE) 150 bp reads.
Low quality reads from Illumina sequencing were filtered out using Cutadapt software v1.16 with parameters of—discard-trimmed-n3-e0.1 (). The obtained clean reads were mapped to the D. citri genome (NCBI accession No. GCA_000475195.1) with the Hisat2 software (). Read counts for all genes were extracted with HTSeq-count () and normalized using the R package DESeq2 (). The differentially expressed genes (DEGs) between the CLas-infected and non-infected lines were estimated by DESeq2 according to the threshold of | log2 ratio| > 1 and adjusted P < 0.05 (BH adjustment). The Blast2GO software was used to generate GO annotation of the DEGs. The DEGs were further annotated in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database by KEGG Automatic Annotation Server. The GO functional enrichment and KEGG pathway enrichment of the DEGs were analyzed using the Fisher’s exact test with Benjamini and Hochberg (BH) adjustment.
Validation of DEGs by qPCR
Twelve genes were selected for qPCR validation based on the differentially expressed level and specific functions (Supplementary Table S1). Following the user manual, the first cDNA strand fragments were synthesized from total RNA using the PrimeScriptTM RT Master Mix kit (Takara, Japan). Specific primers for selected genes are listed in Supplementary Table S1. qRT-PCR was performed on an ABI7500 real-time fluorescent quantitative PCR instrument (Applied Biosystems, Foster City, United States). The relative mRNA levels of the selected genes were estimated according to the aforementioned method.
Metabolite Extraction and LC-MS/MS Analysis
For metabolite extraction, 20 female adults from each insect line were pooled and surface washes of sterile water for 1 min were performed three times. Then the samples were ground to a powder form with an electronic grinding machine. A total of 20 mg of ground sample was placed in an EP tube, and a 1,000 μL extract solution (acetonitrile: methanol: water = 2: 2: 1, with isotopically labelled internal standard mixture) was added. After 30 s of vortex mixing, the samples were homogenized at 35 Hz for 4 min and sonicated for 5 min in ice-water bath. The homogenization and sonication cycle were repeated three times. Then, the samples were incubated for 1 h at −40°C and centrifuged at 12,000 rpm for 15 min at 4°C. The resulting supernatant was transferred to a fresh glass vial for analysis. The quality control sample was prepared by mixing an equal aliquot of supernatant from all the samples. Six biological replications were performed for each insect line.
LC-MS/MS analyses were performed using an UHPLC system (Vanquish, Thermo Fisher Scientific) with a UPLC BEH Amide column (2.1 mm × 100 mm, 1.7 μm) coupled to a Q Exactive HFX mass spectrometer (Orbitrap MS, Thermo Fisher Scientific). The mobile phase consisted of 25 mmol/L ammonium acetate and 25 ammonia hydroxide in water (pH = 9.75) (A) and acetonitrile (B). The analysis was carried with elution gradient as follows: 0–0.5 min, 95% B; 0.5–7.0 min, 95–65% B; 7.0–8.0 min, 65–40% B; 8.0–9.0 min, 40% B; 9.0–9.1 min, 40–95% B; 9.1–12.0 min, 95% B. Column temperature was 30°C. Auto-sampler temperature was 4°C, and the injection volume was 2 μL. The QE HFX mass spectrometer was used for its ability to acquire MS/MS spectra on information-dependent acquisition mode in the control of the acquisition software (Xcalibur, Thermo Fisher Scientific). In this mode, the acquisition software continuously evaluates the full scan MS spectrum. The electrospray ionization (ESI) source conditions were set as follows: sheath gas flow rate 50 Arb, Aux gas flow rate 10 Arb, capillary temperature 320°C, full MS resolution 60,000, MS/MS resolution 7,500, collision energy 10/30/60 in NCE mode, spray voltage 3.5 kV (positive) or −3.2 kV (negative).
Metabolome Data Processing and Annotation
For peak detection, extraction, alignment, and integration, the raw MS/MS spectra data were converted to the mzXML format using ProteoWizard and processed with an in-house program, developed using R and based on XCMS. Then, an in-house MS2 database (BiotreeDB) was used for metabolite annotation, with a cutoff value of 0.3.
The difference in metabolite contents between the two groups was calculated using analysis of variance (ANOVA) and the Mann-Whitney U-test in SPSS 13.0. Statistical significance was set at P < 0.05. MetaboAnalyst1 was further used for the hierarchical cluster analysis of all samples and overlapped metabolites. PCA and partial least squares-discriminant analysis (PLS-DA) were conducted to investigate the relationships among the test samples. For the pathway enrichment analysis, the differential metabolites between tested groups were assigned to metabolic pathways using the tool MetPA in MetaboAnalyst. The −log(P) and impact values of each metabolic pathway were calculated using the hypergeometric test.
Integrated Analysis
The transcriptome and metabolome pathways that overlapped were considered, and the corresponding KEGG xml files (KGML) were downloaded from the KEGG pathway database. The interactions among these pathways were estimated using a custom Perl script and the KGML files as input files (; ). Finally, the pathway–pathway network was visualized using the Cytoscape software (V.3.6.0).
Results
CLas-Infected D. citri
The CLas-specific fragments had 1,167 bp, and the gel imaging results showed that D. citri carried CLas after feeding on CLas-infected plants for 14 days (Figure 1A). The qPCR test also confirmed the CLas infection in the infected insect line, which showed a much higher level of CLas mRNA than the control non-infected line (Figure 1B). Thus, the target insect lines could be used for the subsequent transcriptomic and metabolomic analyses.
FIGURE 1
Differences in the Transcriptomic Profile Between CLas-Infected and Non-infected Insects
The RNA-seq experiments yielded an average 42,029,244 raw paired-end reads (range from 39,283,042 to 44,690,980) for each sample with an average Q30 value of 93.94% (Supplementary Table S2). After quality filtering, we obtained an average 39,538,101 high-quality clean reads for each sample. The average mapping rate of the clean reads to the D. citri genome was 81.91% (range from 81.10 to 82.66%), and the mapped reads were used for subsequent analysis.
The correlation test and principal component analysis (PCA) results were consistent regarding the triplicate RNA sequencing (RNA-seq) datasets of CLas-infected and non-infected insect lines (Supplementary Figure S1). The comparative analysis revealed that CLas infection resulted in 1,194 DEGs (4.97% of the total detected genes in the D. citri genome), including 662 up- and 532 downregulated genes (Figure 2A and Supplementary Table S3). Twelve randomly selected DEGs were validated by qPCR (Figure 2B), suggesting the reliability of the RNA-seq results. Among the upregulated genes, we found eight vitellogenin genes and two hemocyanin genes (Supplementary Table S2); the upregulation of these proteins has been observed in previous proteomic studies (; ), and our findings further confirmed this at a transcriptomic level.
FIGURE 2
DEGs were mostly enriched in gene ontology (GO) functional categories of carboxylic acid metabolic process, organic acid metabolic process, and cofactor binding (Figure 3A and Supplementary Table S3). However, the functional enrichment of upregulated and downregulated genes was largely different. For example, the downregulated genes were mostly related to functions of photoreceptor activity, protein-chromophore linkage, and hydroxymethyl-, formyl-, and related transferase activity, whereas the upregulated genes were preferentially enriched in the functions of lipid transporter activity and lipid transporter activity (Supplementary Table S4).
FIGURE 3
We further analyzed the pathways affected by the DEGs. The analysis showed the enrichment of DEGs in multiple carbohydrate metabolism pathways, including the citrate cycle (TCA cycle), ascorbate and aldarate metabolism, pentose, and glucuronate interconversions, and pyruvate metabolism (Figure 3B and Supplementary Table S5). This indicated that CLas infection may affect the carbon utilization of D. citri. Moreover, two immune system-related pathways, antigen processing, and presentation and NOD-like receptor signaling pathway, were found to be significantly enriched with DEGs, suggesting that CLas inevitably modulated D. citri immunity. Additionally, we noticed that several cofactor and vitamin metabolism pathways were significantly affected by the DEGs. This is in accordance with the GO enrichment in the functional category of cofactor binding. The affected cofactors and vitamins included folate, riboflavin, and retinol.
Metabolic Differences Conferred by CLas Infection
The metabolome of insects from CLas-infected and non-infected lines was measured using LC-MS/MS, which identified a total of 315 metabolites across all samples (Supplementary Table S6). The PCA showed that almost all samples were located within the 95% confidence interval, and these samples were distinctly separated into two groups according to their grouping situation (Figure 4A). The substantial difference between CLas-infected and non-infected insect lines was further confirmed by OPLS-DA score plots (Supplementary Figure S2). Of all detected metabolites, 105 showed significant difference in abundance between CLas-infected and non-infected lines (Mann-Whitney U-test, P < 0.05, Figure 4B). This result revealed that infection with CLas led to higher level of steroids (ergosterol), flavonoids (castavinol), and nucleosides such as cytosine, inosine, 5′-methylthioadenosine and 2-methylguanosine. In addition, we observed the upregulation of fatty acid molecules (linoleic acid, polyoxyethylene dioleate, 2-methylbutyroylcarnitine, L-acetylcarnitine, butyrylcarnitine, and polyoxyethylene dioleate) but the downregulation of fatty amides (linoleamide, palmitic amide, and oleamide) in CLas-infected insects (Figure 4B).
FIGURE 4
The metabolic pathway analysis revealed that the pathways of vitamin B6 metabolism, glycerophospholipid metabolism, and purine metabolism were largely affected by the differentially expressed metabolites under CLas infection (Figure 4C and Supplementary Table S7). Notably, most of the detected metabolites in these pathways were significantly upregulated, indicating the potentially positive effect of CLas in the metabolism of these compounds.
Regulatory Network Associated With CLas Infection
We further determined the gene regulatory network associated with CLas infection using an integrated analysis of the transcriptome and metabolome. Nine pathways were found to be affected by both differentially expressed genes and metabolites (Supplementary Table S8), and these pathways could be linked to each other through differential factors except the tryptophan metabolism pathway (Figure 5). The glycine, serine, and threonine metabolism contained the most differential factors, and this pathway interacted with those of the cysteine and methionine metabolism, glycerophospholipid metabolism, sphingolipid metabolism, and purine metabolism. The purine metabolism pathway also harbored multiple differential factors, most of which were significantly upregulated under CLas infection.
FIGURE 5
Among all related differential factors, the downregulated gene LOC103515380 (branched-chain-amino-acid aminotransferase) was involved in both pantothenate and CoA biosynthesis and cysteine and methionine metabolism pathways, indicating its reduced function in the regulatory network under CLas infection. In addition, we observed increased levels of the metabolite phosphatidylcholine, an intermediate product of the glycerophospholipid metabolism and the origin of the arachidonic acid metabolism (Supplementary Figure S3). Overall, the network analysis showed that CLas infection mainly impacted the metabolism of amino acids, lipids, and cofactors of D. citri.
Discussion
As a key vector of CLas transmission, D. citri interacts with CLas in a circulative and propagative manner; moreover, CLas can be distributed in various insect tissues, including the hemolymph and salivary glands (; ). Previous studies on the transcriptomic and proteomic changes in CLas-infected and non-infected D. citri have indicated the involvement of CLas infection in the nutrition metabolism, nymphal development, and immunity of D. citri (; ; , ; ). Here, we focused on the transcriptomic changes in the metabolism pathways conferred by CLas infection in D. citri by using a combined transcriptomic and metabolic approach. The transcriptomic analysis revealed 662 up- and 532 downregulated genes. The DEG pattern was similar to that previously obtained from the gut transcriptome (; ). The LC-MS/MS analysis identified 315 metabolites in D. citri, a number much larger than those obtained in previous GC-MS analyses (). These differential factors represent the body-wide transcriptomic and metabolomic responses to CLas infection, which should be compared to the values from specific tissues.
Previous transcriptomic and proteomic analyses of D. citri showed an upregulation of transcripts and proteins involved in defense and immunity (including phenoloxidase, vitellogenin, and hemocyanin) in response to CLas infection (, ; ). Our data further showed that genes in the pathways of antigen processing and presentation and NOD-like receptor signaling were significantly affected (Supplementary Table S5). The heat shock proteins (HSP) coding genes in these two pathways were upregulated, as reported by and . HSPs have been recognized as important immune response factors to biological and environmental stresses in multiple insect systems (; ; ). Thus, the upregulation of HSPs indicated the activation of D. citri’s immune system conferred by CLas infection. The draft genome sequence indicated that D. citri presented the Toll immune signaling pathway, as shown by previous studies (; ; ). In the transcriptome herein obtained, we also observed the upregulation of Ras-related C3 botulinum toxin substrate (LOC103512003), a gene induced by Toll-like receptor 2 in the Toll immune signaling pathway ().
Another significant change conferred by CLas infection was with respect to amino acid metabolism. At the transcriptomic level, the metabolism pathways of glycine, serine, threonine, tryptophan, phenylalanine, tyrosine, and glutathione were significantly affected (Supplementary Table S5). Genomic comparisons of CLas with cultivable L. crescens revealed a lack of genes related to the production of fatty acids and aromatic amino acids (phenylalanine, tryptophan, and tyrosine) (; ; ). This suggests that CLas must rely on either the host plant or the insect vector for these essential nutrients, and therefore they may use these amino acids by modulating the corresponding metabolism pathways in D. citri. The metabolomic analysis also revealed the significant change in the tryptophan metabolism pathway (Figure 4). Moreover, CLas infection elevated the abundance of several aliphatic amino acids (glycine, L-asparagine, and beta-alanine), the metabolism of which is closely associated with the fatty acid metabolism.
Our results showed that CLas-infected D. citri tended to be enriched in fatty acids (Supplementary Table S6), which is in accordance with the results of previous studies (; ). The increased abundance of linoleic acid was also observed in CLas-infected D. citri nymphs (). In contrast, the abundance of some fatty amides, derivatives of fatty acids, was significantly reduced. This means that CLas infection promotes the accumulation of fatty acids in D. citri through the modulation of fatty acid biosynthesis and the reduced transformation of fatty acids to the respective derivatives. The increasing abundance of fatty acids may also occur because of their important roles in the innate immunity of insects (). Remarkably, we observed the upregulation of lipid transporter activity and the significant overexpression of almost all vitellogenin genes (Supplementary Tables S3, S4). Vitellogenin is a lipoprotein with functions in insect fecundity and immunity, and it is synthesized in the fat body and secreted into the hemolymph (). Vitellogenin plays an important nutritional role during vitellogenesis () and has been shown to function in insect immunity by binding to pathogen lipopolysaccharides and peptidoglycans (). Previous proteomic data have shown that CLas infection causes the upregulation of vitellogenin proteins in the D. citri hemolymph, gut, or whole body (; , ). Our finding further supports the upregulation of vitellogenin in the CLas-infected D. citri, and this is in accordance with a previous report that showed the increasing reproductive fitness in CLas-infected D. citri ().
The transcriptomic results also suggested that CLas modulated the cofactor and vitamin metabolism in D. citri (Figures 3B, 4C). Pantothenate, the starting compound in CoA biosynthesis, was found to be significantly reduced by CLas infection. This may result from the interaction between CLas metabolism and host proteins (). Our study also indicated the decrease in riboflavin in CLas-infected insects. Riboflavin provisioning has been previously reported in bacteria–insect endosymbiotic relationships (; ). Therefore, the observed reduction in riboflavin may result from the impact of CLas on other symbionts that confer fitness in riboflavin provisioning to D. citri, as CLas infection not only modulates host immunity but is also related to other microbes in the host insect (; ). In fact, bacteria benefiting from the host insect’s cofactor and vitamin biosynthesis is a common phenomenon (; ; ), and this may also explain the potential downregulation of folate biosynthesis and retinol metabolism pathways inferred from the transcriptomic analysis herein shown (Supplementary Table S5).
Conclusion
In conclusion, by using transcriptomic and metabolomic analyses, we identified the genes and metabolites that differentially expressed between CLas-infected and non-infected D. citri. The integrated analysis of these two omics data revealed that CLas infection had a remarkable impact on pathways associated with the metabolism of amino acids, lipids, and cofactors. The present study provides a foundation for further investigation on the mechanisms of CLas–D. citri interaction.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
KL, RP, and QH conceived and designed the study and wrote the manuscript. JH, ZG, and MZ contributed to materials. KL and JH performed the experiments. RP performed the data analysis. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the Science and Technology Planning Project of Guangdong Province, China (2016A050502063), Characteristic Innovation Project of Ordinary University of Guangdong Province, China (KA1548807), and Ecological Control of Agricultural Pests Construction Fund (KA200540404).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2020.630037/full#supplementary-material
Supplementary Figure 1The correlation of the RNA-seq datasets from CLas-infected and non-infected D. citri. The correlation was calculated according to Spearman’s rank correlation coefficient (A) and principal component analysis (B).
Supplementary Figure 2Orthogonal projections tolatent structures- discriminant analysis (OPLS-DA) plots of metabolite composition of samples from CLas-infected and non-infected D. citri. Each symbol represents a sample, different symbol shapes denote different groups. The circle line indicates the 95% confidence interval (Hotelling’s T-squared ellipse).
Supplementary Figure 3Pathways of Glycerophospholipid metabolism (A) and Arachidonic acid metabolism (B) affected by CLas infection. Red indicates that the gene/compound is up-regulated, blue indicates that the gene/compound is down-regulated.
Supplementary Table 1The sequences of primers used in this study.
Supplementary Table 2The summary of RNA-seq datasets.
Supplementary Table 3The expression statistics of all genes obtained by RNA-seq.
Supplementary Table 4Results of the GO enrichment analysis for the DEGs.
Supplementary Table 5Results of the KEGG pathway enrichment analysis for the DEGs.
Supplementary Table 6The information of compounds identified by the metabolome analysis.
Supplementary Table 7The pathways affected by the differential metabolites.
Supplementary Table 8The overlapped pathways identified from transcriptome and metabolome analysis.
Footnotes
References
1
AguilarR.JedlickaA. E.MintzM.MahairakiV.ScottA. L.DimopoulosG. (2005). Global gene expression analysis of Anopheles gambiae responses to microbial challenge.Insect Biochem. Mol. Biol.35709–719. 10.1016/j.ibmb.2005.02.019
2
AmmarE.-D.ShattersR. G.Jr.HallD. G. (2011). Localization of Candidatus Liberibacter asiaticus, associated with Citrus Huanglongbing Disease, in its Psyllid vector using fluorescence in situ hybridization.J. Phytopathol.159726–734. 10.1111/j.1439-0434.2011.01836.x
3
AndersS.PylP. T.HuberW. (2015). HTSeq–a Python framework to work with high-throughput sequencing data.Bioinformatics31166–169. 10.1093/bioinformatics/btu638
4
ArbibeL.MiraJ. P.TeuschN.KlineL.GuhaM.MackmanN.et al (2000). Toll-like receptor 2-mediated NF-kappa B activation requires a Rac1-dependent pathway.Nat. Immunol.1533–540. 10.1038/82797
5
ArreseE. L.SoulagesJ. L. (2010). Insect fat body: energy, metabolism, and regulation.Annu. Rev. Entomol.55207–225. 10.1007/978-1-4615-9221-1_8
6
BauerE.SalemH.MarzM.VogelH.KaltenpothM. (2014). Transcriptomic immune response of the cotton stainer Dysdercus fasciatus to experimental elimination of vitamin-supplementing intestinal symbionts.PLoS One9:e114865. 10.1371/journal.pone.0114865
7
BovéJ. M. (2006). Huanglongbing: a destructive, newly-emerging, century-old disease of citrus.J. Plant Pathol.887–37.
8
CenY.YangC.HolfordP.BeattieG.a.C.Spooner-HartR. N.LiangG.et al (2012). Feeding behaviour of the Asiatic citrus psyllid, Diaphorina citri, on healthy and huanglongbing-infected citrus.Entomol. Exp. Appl.14313–22. 10.1111/j.1570-7458.2012.01222.x
9
ChenD.ZhangH.LuP.LiuX.CaoH. (2016). Synergy evaluation by a pathway-pathway interaction network: a new way to predict drug combination.Mol. Biosyst.12614–623. 10.1039/c5mb00599j
10
DuanY.ZhouL.HallD. G.LiW.DoddapaneniH.LinH.et al (2009). Complete genome sequence of citrus Huanglongbing bacterium, ‘Candidatus Liberibacter asiaticus’ obtained through metagenomics.Mol. Plant Microbe Interact.221011–1020. 10.1094/mpmi-22-8-1011
11
FerrareziR. S.VincentC. I.UrbanejaA.MachadoM. A. (2020). Editorial: unravelling citrus huanglongbing disease.Front. Plant. Sci.11:609655. 10.3389/fpls.2020.609655
12
GhanimM.AchorD.GhoshS.KontsedalovS.LebedevG.LevyA. (2017). ‘Candidatus Liberibacter asiaticus’ accumulates inside endoplasmic reticulum associated vacuoles in the gut cells of Diaphorina citri.Sci. Rep.7:16945.
13
GillT. A.ChuC.Pelz-StelinskiK. S. (2017). Comparative proteomic analysis of hemolymph from uninfected and Candidatus Liberibacter asiaticus-infected Diaphorina citri.Amino Acids49389–406. 10.1007/s00726-016-2373-2
14
GottwaldT. R. (2010). Current epidemiological understanding of citrus Huanglongbing.Annu. Rev. Phytopathol.48119–139. 10.1146/annurev-phyto-073009-114418
15
GotzM.PopovskiS.KollenbergM.GorovitsR.BrownJ. K.CiceroJ. M.et al (2012). Implication of Bemisia tabaci heat shock protein 70 in Begomovirus-whitefly interactions.J. Virol.8613241–13252. 10.1128/jvi.00880-12
16
GracaJ. V. (1991). Citrus greening disease.Annu. Rev. Phytopathol.29109–136.
17
GracaJ. V. (2008). “Biology, history and world status of huanglongbing,” in Proceedings of the I Taller Internacional Sobre Huanglongbing de los Cítricos (Candidatus Liberibacter spp.) y el Psílido Asiático de los Cítricos (Diaphorina citri), Hermosillo, 1–7. 10.1079/9781786394088.0001
18
JuJ. F.BingX. L.ZhaoD. S.GuoY.XiZ.HoffmannA. A.et al (2020). Wolbachia supplement biotin and riboflavin to enhance reproduction in planthoppers.ISME J.14676–687. 10.1038/s41396-019-0559-9
19
KillinyN.HijazF.El-SheshenyI.AlfaressS.JonesS. E.RogersM. E. (2017). Metabolomic analyses of the haemolymph of the Asian citrus psyllid Diaphorina citri, the vector of huanglongbing.Physiol. Entomol.42134–145. 10.1111/phen.12183
20
KillinyN.JonesS. E. (2018). Metabolic alterations in the nymphal instars of Diaphorina citri induced by Candidatus Liberibacter asiaticus, the putative pathogen of huanglongbing.PLoS One13:e0191871. 10.1371/journal.pone.0191871
21
KimD.LangmeadB.SalzbergS. L. (2015). HISAT: a fast spliced aligner with low memory requirements.Nat. Methods12357–360. 10.1038/nmeth.3317
22
KruseA.Fattah-HosseiniS.SahaS.JohnsonR.WarwickE.SturgeonK.et al (2017). Combining ‘omics and microscopy to visualize interactions between the Asian citrus psyllid vector and the Huanglongbing pathogen Candidatus Liberibacter asiaticus in the insect gut.PLoS One12:e0179531. 10.1371/journal.pone.0179531
23
KruseA.RamseyJ. S.JohnsonR.HallD. G.MacCossM. J.HeckM. (2018). Candidatus Liberibacter asiaticus minimally alters expression of immunity and metabolism proteins in hemolymph of Diaphorina citri, the insect vector of Huanglongbing.J. Proteome Res.172995–3011. 10.1021/acs.jproteome.8b00183
24
LaiK. K.Davis-RichardsonA. G.DiasR.TriplettE. W. (2016). Identification of the genes required for the culture of Liberibacter crescens, the closest cultured relative of the Liberibacter plant pathogens.Front. Microbiol.7:547. 10.3389/fmicb.2016.00547
25
LiZ.ZhangS.ZhangJ.LiuM.LiuZ. (2009). Vitellogenin is a cidal factor capable of killing bacteria via interaction with lipopolysaccharide and lipoteichoic acid.Mol. Immunol.463232–3239. 10.1016/j.molimm.2009.08.006
26
LivakK. J.SchmittgenT. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCT method.Methods25402–408. 10.1006/meth.2001.1262
27
LoveM. I.HuberW.AndersS. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2.Genome Biol.15550.
28
MannR. S.AliJ. G.HermannS. L.TiwariS.Pelz-StelinskiK. S.AlbornH. T.et al (2012). Induced release of a plant-defense volatile ‘deceptively’ attracts insect vectors to plants infected with a bacterial pathogen.PLoS Pathog.8:e1002610. 10.1371/journal.ppat.1002610
29
MartinM. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads.EMBnet J.1710–12. 10.14806/ej.17.1.200
30
MartínezY.LlaugerR.BatistaL.LuisM.IglesiaA.CollazoC. (2009). First report of ‘Candidatus Liberibacter asiaticus’ associated with Huanglongbing in Cuba.Plant Pathol.58389–389.
31
MartiniX.HoffmannM.CoyM. R.StelinskiL. L.Pelz-StelinskiK. S. (2015). Infection of an insect vector with a bacterial plant pathogen increases its propensity for dispersal.PLoS One10:e0129373. 10.1371/journal.pone.0129373
32
MoriyamaM.NikohN.HosokawaT.FukatsuT. (2015). Riboflavin provisioning underlies Wolbachia’s fitness contribution to its insect host.mBio6e01732-15.
33
PangR.ChenM.YueL.XingK.LiT.KangK.et al (2018). A distinct strain of Arsenophonus symbiont decreases insecticide resistance in its insect host.PLoS Genet.14:e1007725. 10.1371/journal.pgen.1007725
34
Pelz-StelinskiK. S.BrlanskyR. H.EbertT. A.RogersM. E. (2010). Transmission parameters for Candidatus liberibacter asiaticus by Asian citrus psyllid (Hemiptera: Psyllidae).J. Econ. Entomol.1031531–1541. 10.1603/ec10123
35
Pelz-StelinskiK. S.KillinyN. (2016). Better together: association with ‘Candidatus Liberibacter Asiaticus’ increases the reproductive fitness of its insect vector, Diaphorina citri (Hemiptera: Liviidae).Ann. Entomol. Soc. Am.109371–376. 10.1093/aesa/saw007
36
PockleyA. G. (2003). Heat shock proteins as regulators of the immune response.Lancet362469–476. 10.1016/s0140-6736(03)14075-5
37
RamseyJ. S.ChavezJ. D.JohnsonR.HosseinzadehS.MahoneyJ. E.MohrJ. P.et al (2017). Protein interaction networks at the host-microbe interface in Diaphorina citri, the insect vector of the citrus greening pathogen.R. Soc. Open Sci.4:160545. 10.1098/rsos.160545
38
RamseyJ. S.JohnsonR. S.HokiJ. S.KruseA.MahoneyJ.HilfM. E.et al (2015). Metabolic interplay between the Asian citrus psyllid and its Profftella symbiont: an Achilles’ heel of the citrus greening insect vector.PLoS One10:e0140826. 10.1371/journal.pone.0140826
39
RenS. L.LiY. H.ZhouY. T.XuW. M.CuthbertsonA. G. S.GuoY. J.et al (2016). Effects of Candidatus Liberibacter asiaticus on the fitness of the vector Diaphorina citri.J. Appl. Microbiol.1211718–1726.
40
SalemH.BauerE.StraussA. S.VogelH.MarzM.KaltenpothM. (2014). Vitamin supplementation by gut symbionts ensures metabolic homeostasis in an insect host.Proc. Biol. Sci.281:20141838. 10.1098/rspb.2014.1838
41
StanleyD.HaasE.MillerJ. (2012). Eicosanoids: exploiting insect immunity to improve biological control programs.Insects3492–510. 10.3390/insects3020492
42
TiwariS.GondhalekarA. D.MannR. S.ScharfM. E.StelinskiL. L. (2011). Characterization of five CYP4 genes from Asian citrus psyllid and their expression levels in Candidatus Liberibacter asiaticus-infected and uninfected psyllids.Insect Mol. Biol.20733–744. 10.1111/j.1365-2583.2011.01103.x
43
TufailM.TakedaM. (2008). Molecular characteristics of insect vitellogenins.J. Insect Physiol.541447–1458. 10.1016/j.jinsphys.2008.08.007
44
TylerH. L.RoeschL. F. W.GowdaS.DawsonW. O.TriplettE. W. (2009). Confirmation of the sequence of ‘Candidatus Liberibacter asiaticus’ and assessment of microbial diversity in Huanglongbing-infected citrus phloem using a metagenomic approach.Mol. Plant Microbe Interact.221624–1634. 10.1094/mpmi-22-12-1624
45
VyasM.FisherT. W.HeR.NelsonW.YinG.CiceroJ. M.et al (2015). Asian citrus psyllid expression profiles suggest Candidatus Liberibacter Asiaticus-mediated alteration of adult nutrition and metabolism, and of nymphal development and immunity.PLoS One10:e0130328. 10.1371/journal.pone.0130328
46
WangN.TrivediP. (2013). Citrus huanglongbing: a newly relevant disease presents unprecedented challenges.Phytopathology103652–665. 10.1094/phyto-12-12-0331-rvw
47
YuH. Z.LiN. Y.ZengX. D.SongJ. C.YuX. D.SuH. N.et al (2020). Transcriptome analyses of Diaphorina citri midgut responses to Candidatus Liberibacter Asiaticus Infection.Insects11:171. 10.3390/insects11030171
Summary
Keywords
Diaphorina citri, Candidatus Liberibacter asiaticus, interaction, transcriptome, metabolome
Citation
Liu K, He J, Guan Z, Zhong M, Pang R and Han Q (2021) Transcriptomic and Metabolomic Analyses of Diaphorina citri Kuwayama Infected and Non-infected With Candidatus Liberibacter Asiaticus. Front. Physiol. 11:630037. doi: 10.3389/fphys.2020.630037
Received
16 November 2020
Accepted
21 December 2020
Published
24 February 2021
Volume
11 - 2020
Edited by
Bin Tang, Hangzhou Normal University, China
Reviewed by
Yanjing Wang, Zhejiang Agriculture and Forestry University, China; Dandan Wei, Southwest University, China; Sheila Ons, National University of La Plata, Argentina
Updates
Copyright
© 2021 Liu, He, Guan, Zhong, Pang and Han.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rui Pang, pr839@163.com; pangr6@mail.sysu.edu.cnQunxin Han, hqx99@163.com
†These authors have contributed equally to this work
This article was submitted to Invertebrate Physiology, a section of the journal Frontiers in Physiology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.