Next Article in Journal
Mapping Potential Distribution of Spodoptera frugiperda (Lepidoptera: Noctuidae) in Central Asia
Previous Article in Journal
Long-Tailed Silverfish (Ctenolepisma longicaudata) Control; Bait Choice Based on Primary and Secondary Poisoning
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 11
Issue 3
10.3390/insects11030171
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptome Analyses of Diaphorina citri Midgut Responses to Candidatus Liberibacter Asiaticus Infection

by
Hai-Zhong Yu
1,2,†,
Ning-Yan Li
1,†,
Xiang-Dong Zeng
1,
Jian-Chun Song
1,
Xiu-Dao Yu
1,2,
Hua-Nan Su
1,2,
Ci-Xiang Chen
3,
Long Yi
1,2 and
Zhan-Jun Lu
1,2,*
1
College of Life Sciences, Gannan Normal University, Ganzhou 341000, China
2
National Navel Orange Engineering Research Center, Gannan Normal University, Ganzhou 341000, China
3
Fruit Bureau of Ganzhou, Ganzhou 341000, China
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Insects 2020, 11(3), 171; https://doi.org/10.3390/insects11030171
Submission received: 4 February 2020 / Revised: 29 February 2020 / Accepted: 5 March 2020 / Published: 7 March 2020
Browse Figures
Review Reports Versions Notes

Abstract

:
The Asian citrus psyllid (ACP), Diaphorina citri Kuwayama (Hemiptera: Liviidae), is an important transmission vector of the citrus greening disease Candidatus Liberibacter asiaticus (CLas). The D. citri midgut exhibits an important tissue barrier against CLas infection. However, the molecular mechanism of the midgut response to CLas infection has not been comprehensively elucidated. In this study, we identified 778 differentially expressed genes (DEGs) in the midgut upon CLas infection, by comparative transcriptome analyses, including 499 upregulated DEGs and 279 downregulated DEGs. Functional annotation analysis showed that these DEGs were associated with ubiquitination, the immune response, the ribosome, endocytosis, the cytoskeleton and insecticide resistance. KEGG enrichment analysis revealed that most of the DEGs were primarily involved in endocytosis and the ribosome. A total of fourteen DEG functions were further validated by reverse transcription quantitative PCR (RT-qPCR). This study will contribute to our understanding of the molecular interaction between CLas and D. citri.

1. Introduction

Huanglongbing (HLB) is a destructive disease of citrus that represents a major threat to the world’s citrus industry. HLB is caused by the phloem-limited, Gram-negative bacterium “Candidatus Liberibacter spp” (CLas). HLB nearly destroyed the citrus industry in Florida (USA) and most citrus-producing regions of the world [1,2]. CLas can decrease plant growth vigor, and ultimately result in the death of the infected citrus tree [3]. There are no effective methods to control HLB once established. Recently, many antibiotics have been adopted to control CLas bacteria; examples include the trunk injection of penicillin, streptomycin and oxytetracycline hydrochloride [1,4]. However, these approaches do not work well in the field. The killing of vectors is an effective method to reduce HLB spread [5].
The Asian citrus psyllid (ACP), Diaphorina citri Kuwayama, is the principal transmission vector of HLB [6]. Vector control of D. citri is recognized as a key approach to preventing the spread of HLB. The application of insecticides is the most widely employed option for reducing D. citri populations in some citrus-growing regions [7]. However, the improper use of such chemicals has caused the poisoning of farmers, environmental pollution and insect resistance [8,9]. Some researchers have focused on the olfactory system of D. citri, revealing a suite of odorants that can be used to develop affordable and safe odor-based surveillance for D. citri control [10]. Yu et al. revealed that the silencing of D. citri muscle protein 20 (DcMP20) resulted in significant mortality and reduced the body weight of D. citri [11]. By RNA interference (RNAi) technology, the knockdown of D. citri tropomyosin1-X1 (DcTm1-X1) significantly increased the mortality rate of nymphs [12]. Although some success regarding the control of D. citri has been achieved, the threats caused by HLB are going to continue. It is essential to investigate the molecular mechanisms of HLB transmission in D. citri.
Both D. citri adults and nymphs can transmit CLas. Recent studies suggested that the pathogen multiplies in D. citri after acquisition and becomes distributed among various internal tissues, including the alimentary canal, salivary glands, hemolymph, filter chamber, midgut, fat body, muscle tissues and ovaries [13,14,15]. Among these tissues, the midgut is the first barrier that the bacterium must breach before entry into the hemolymph, indicating that the midgut plays an important role in CLas infection [16]. Pathogens can manipulate the host’s essential biological processes during the host-pathogen interaction, e.g., by changing protein translation, vesicular transport and protein metabolism [17,18]. In the midgut, many genes and proteins are involved in the defense against pathogens. Wang et al. identified 869 differentially expressed genes (DEGs) in the Bombyx mori larval midgut via comparative transcriptome analysis. The results showed that many of the DEGs were associated with protein metabolism, the cytoskeleton and apoptosis [19]. Bao et al. performed a transcriptome-wide analysis on the Nilaparvata lugens intestine, and a total of 33 digestion-related genes, 25 immune response genes and 27 detoxification-related genes were identified [20]. The interactions between D. citri and CLas have become a heavily researched topic regarding the molecular mechanisms of HLB transmission. In recent years, high-throughput omics technologies have been widely used to study the interactions of vector-pathogens. Combining transcriptomics with proteomics methods, Kruse et al. identified many DEGs involved in the TCA cycle, iron metabolism, insecticide resistance and the insect immune system in the D. citri gut after exposure to CLas [21]. Lu et al. also identified 62 DEPs from uninfected and CLas-infected adult D. citri, using two-dimensional gel electrophoresis. These DEPs were associated with energy metabolism, host detoxification processes and the cytoskeleton [22]. The hemolymph is a critical component of the D. citri immune system, where it coordinates the insect’s immune system activity. Kruse et al. investigated the effect of CLas exposure on the D. citri hemolymph, using high-resolution quantitative mass spectrometry, and the results revealed that proteins associated with fatty acid synthesis and energy metabolism were upregulated [23]. However, many confusing issues remain to be addressed, such as the mechanism underlying CLas entry, multiplication and dissemination in the midgut of D. citri.
In this study, we performed a comparative transcriptome analysis to determine the specific responses of D. citri to CLas infection in the midgut at the mRNA level. A total of 778 DEGs were identified, and these genes were mainly involved in ubiquitination, the immune response, the ribosome, endocytosis, the cytoskeleton and insecticide resistance. This study established a foundation for further understanding the interaction mechanisms between D. citri and CLas.

2. Materials and Methods

2.1. Insect Rearing, Tissue Collection and Total RNA Extraction

D. citri was reared in mesh cages on Murraya exotica at 27 ± 1 °C, 70% ± 5% relative humidity, and a 14:10 (light–dark) photocycle. Healthy D. citri adults were transferred to CLas-infected Newhall Navel oranges and maintained under the same rearing conditions. One hundred D. citri from CLas-infected Newhall Navel oranges were randomly collected for calculating the infection rate. In general, the head of D. citri was cut off using a razor blade, and the genome DNA was extracted using the TIANamp Micro DNA kit (TianGen Biotech, Beijing, China). Afterwards, a PCR reaction was performed to amplify the 16S ribosomal DNA fragment using the HLB pathogen specific primer set OI1/OI2c (forward primer 5′-GCGCGTATGCAATACGAGCGGCA-3′ and reverse primer 5′-GCCTCGCGACTTCGCAACCCAT-3′) based on a previous protocol [24]. Then, agarose gel electrophoresis was performed to investigate the CLas infection rate. When the CLas-infected percentage reached above 80%, the infected D. citri were selected. The uninfected and CLas-infected D. citri were collected and dissected to obtain the midgut, which was then washed with precooled DEPC-water to remove the remaining tissues. Genomic DNA was isolated from the D. citri midgut, and a PCR was performed to confirm the CLas infection rate.
For transcriptome sequencing, three groups, each consisting of three hundred CLas-infected or uninfected D. citri, were dissected for RNA extraction. In order to avoid RNA degradation, the collected midgut samples were stored in RNAlater and frozen at −80 °C. Total RNA was extracted from the midgut of uninfected and CLas-infected D. citri using an animal tissue total RNA kit (Simgen, Hangzhou, China). RNA concentration and purity were assayed using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, New York, NY, USA) at absorbance ratios of A260/280 and A260/230. The integrity of total RNA was confirmed using agarose gel electrophoresis.

2.2. Library Preparation and Illumina Sequencing

Transcriptome sequencing was conducted on the Illumina HiSeq platform (Novogene Bioinformatics Technology Co., Ltd. Tianjin, China). Sequencing libraries were generated using the NEBNext1Ultra™ RNA Library Prep Kit for Illumina1 (NEB, Ipswich, MA, USA). The library quality was assessed on the Agilent Bioanalyzer 2100 system. The clustering of the index-coded samples was performed on a cBot Cluster Generation System using TruSeq PE Cluster Kitv3-cBot-HS (Illumina, San Diego, CA, USA). After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform, and 125 bp paired-end reads were generated. The sequencing fragments were translated into fastq format raw reads using the CASAVA software. The raw reads of fastq format were firstly processed through in-house Perl scripts. The clean reads were obtained by removing reads containing adapter, reads containing poly-N, and low quality reads in which a base number of Qphred less than or equal to twenty accounted for more than fifty percent of the entire read length. At the same time, Q20, Q30 and the GC content of the clean data were calculated.

2.3. Read Mapping and Identification of DEGs

The reference genome and gene model annotation files were downloaded from the genome website (ftp://ftp.citrusgreening.org/annotation/OGSv2.0/) [25]. Hisat2 v2.0.5 (https://anaconda.org/biobuilds/hisat2) was used to build the index of the reference genome and align the paired-end clean reads with the reference genome [26]. This generated a database of splice junctions based on the gene model annotation file. The expression levels of the genes were calculated by fragments per kilobase of transcript per million fragments mapped (FPKM). Differential expression analysis between CLas-free groups and CLas-infected groups were performed using the DESeq2 R package (1.16.1). DESeq2 provides statistical routines for determining differential expression according to the negative binomial distribution. The resulting P-values were adjusted using the Benjamini–Hochberg method. A Corrected Genes P-value of 0.05 and an absolute |log2 (fold change)| of 0 were set as the threshold for significantly differential expression. The hierarchical cluster analysis of differentially expressed genes (DEGs) was conducted using Genesis software (http://genome.tugraz.at/genesisclient_download.shtml).

2.4. Gene Ontology (GO) and KEGG Enrichment Analysis of DEGs

Gene Ontology (GO) is a tool used for gene annotation by collecting defined, structured, controlled vocabulary [27]. KEGG (Kyoto Encyclopedia of Genes and Genomes) is a database used to categorize associated gene sets into appropriate pathways [28]. The clusterProfiler R package, which implements the GO terms, was used for the enrichment analysis of DEGs, in which the gene lengths of these DEGs were corrected [29]. KEGG pathway enrichment analysis for DEGs was performed using the clusterProfiler R package [30]. A P-value of < 0.01 was set as the threshold.

2.5. Reverse Transcription Quantitative PCR (RT-qPCR) Analysis

In order to verify the reliability of the transcriptome data, the expression levels of 14 randomly selected DEGs were examined by RT-qPCR. All of the primers used are listed in Table S1. PCR reactions were prepared containing 10 μL of SYBR Ⅱ, 8 μL of ddH2O, 0.5 μL of forward primer, 0.5 μL of reverse primer and 1.0 μL of cDNA template. The reactions were performed on the LightCycle®96 PCR Detection System (Roche, Basel, Switzerland). The relative expression level of each gene was calculated using the 2−∆∆Ct method. Three biological replicates were conducted for each sample. D. citri glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as a reference gene. The expression patterns were compared by ANOVAs, followed by Fisher’s protected least significant difference (LSD) tests [31].

3. Results

3.1. Detection of CLas Infection in the Midgut of D. citri

D. citri midguts were dissected in DEPC-water using a dissecting needle. The results showed that the detached midguts were relatively intact (Figure 1A). To ensure the integrity of the midgut samples for transcriptome sequencing, genomic DNA was isolated, and PCR was performed using the OI1 primer. The primer sequences are shown in Table 1. Analysis of the PCR products on agarose gels showed that the CLas-infected midguts exhibited a clear OI1 band, of 1160 bp in length, that was not present in the control group (Figure 1B) [24]. These results indicated that the selected midgut samples could be used for further transcriptome sequencing.

3.2. Illumina Sequencing and Read Assembly

According to the transcriptome analysis, we generated 62,582,728, 58,399,847, 56,575,636, 57,180,938, 63,335,414 and 54,887,010 raw reads from the CLas-infected and uninfected groups, respectively. The raw sequencing data were deposited in NCBI SRA under the accession numbers SRX7026384, SRX7026385, SRX7026386, SRX7026387, SRX7026388 and SRX7026389. After stringent quality assessment and data filtering, 61,521,780, 57,557,156, 55,768,652, 55,707,048, 62,309,800 and 53,907,044 clean reads were obtained. The Q30 (sequencing error rate < 0.1%) and Q20 (sequencing error rate < 1%) were more than 92% and 97%, respectively. The GC contents were 40.57%, 40.36%, 40.91%, 37.74%, 39.94% and 39.39%, respectively (Table 1). Therefore, the accuracy of the sequencing data was sufficient for further analysis.

3.3. Identification of DEGs in Response to CLas Infection

In total, 798 DEGs were identified in the midgut of CLas-infected D. citri compared with the uninfected D. citri. Among these DEGs, 499 were upregulated and 279 were downregulated (Figure 2A). Basing on the log10(RPKM+1) values of the two groups, we performed hierarchical clustering of all the DEGs to determine the expression patterns of the identified genes (Figure 2B). These results indicated that CLas infection altered the transcriptional profiles of the DEGs.

3.4. GO and KEGG Enrichment Analysis

GO and KEGG enrichment analyses were conducted to further investigate the functions of the DEGs [32]. In this study, a total of 106 upregulated DEGs and 73 downregulated DEGs were performed for GO enrichment analysis. The GO enrichment analysis revealed that a total of 18 upregulated DEGs were mainly associated with the cytoplasm, 20 upregulated DEGs were associated with the organonitrogen compound metabolic process and 14 upregulated DEGs were related to the carbohydrate metabolic process (Figure 3A, Table S2); a total of 20 downregulated DEGs were mainly involved in transmembrane transport and 14 downregulated DEGs were associated with transmembrane transporter activity (Figure 3B, Table S2). KEGG pathway analysis is useful for researching the complex biological functions of genes [33]. According to the KEGG pathway enrichment analysis, a total of seven, three and four upregulated DEGs were significantly enriched in ribosome, proteasome and oxidative phosphorylation, respectively (Figure 4A, Table S3). A total of six downregulated DEGs were significantly enriched in endocytosis (Figure 4B, Table S3).

3.5. Validation of DEGs at the Transcriptional Level

To validate the reliability of the transcriptome sequencing data, the relative expression levels of 14 DEGs involved in different functions were analyzed by RT-qPCR (Figure 5). These results were consistent with the transcriptome data. For example, the gene E3 ubiquitin-protein ligase (DcitrP055490.1) was upregulated in both the transcriptome data and the RT-qPCR analysis, with a similar fold change. However, 40S ribosomal protein S9 (DcitrP013835.1) and heat shock protein 70 (DcitrP017415.1) were upregulated in the midgut after CLas infection, while they had no obvious change in the RT-qPCR analysis (Figure 5A). The linear regression analysis of the correlation between RT-qPCR and transcriptome showed an R2 value of 0.936 and a corresponding slope of 1.0208 (Figure 5B). These results suggested a significant correlation between RT-qPCR and the transcriptome data.

3.6. Analysis of DEGs Associated with Ubiquitination, the Immune Response and the Ribosome

Based on the transcriptome analysis, many DEGs associated with ubiquitination, the immune response and ribosomes were altered in the uninfected groups and the CLas-infected groups (Table 2, Figure 6). For the ubiquitination analysis, a total of 13 DEGs were identified. Among them, eight DEGs (61.5%) were upregulated in the midgut after CLas infection. In addition, five DEGs were downregulated, including ubiquitin-associated domain-containing protein 1, ubiquitin-like modifier-activating enzyme ATG7, E3 ubiquitin-protein ligase RNF8, E3 ubiquitin-protein ligase MARCH2 and ubiquitin-conjugating enzyme E2J2. For the immune response analysis, 11 DEGs were obtained after CLas infection in the midgut; nine genes (81.8%) were upregulated except for Toll-like 8B and Beat protein. For the ribosome analysis, a total of 15 DEGs were identified, and most of the genes were upregulated, except for 28S ribosomal protein S27.

3.7. Analysis of DEGs Involved in Endocytosis, the Cytoskeleton and Insecticide Resistance

Based on the GO and KEGG enrichment analyses, the DEGs involved in endocytosis, the cytoskeleton and insecticide resistance were identified (Table 3, Figure 7). For endocytosis, six (85.7%) DEGs were downregulated, and only one gene was upregulated. For the cytoskeleton analysis, a total of 15 genes were screened, and most of the proteins were upregulated after CLas infection. However, tyrosine protein kinase receptor torso and dynein heavy chain were downregulated. For insecticide resistance, a total of 19 DEGs were found; of these DEGs, 13 (68.4%) were upregulated and six (31.6%) were upregulated, in the midgut after CLas infection.

4. Discussion

HLB is devastating citrus production worldwide and belongs to a phloem-limited α-proteobacterium. However, the study of this bacterium is highly difficult due to the difficulties in culturing it in vitro [34]. The management of HLB depends on the removal of infected trees and the control of the D. citri. Chemical insecticides are currently employed as the primary management strategy, but the widespread use of insecticides leads to serious resistance in D. citri [35]. The study of the interaction between CLas and D. citri is promising for the control of HLB. The transmission of CLas in D. citri is a long-term process that is composed of an acquisition access period (APP), a period of latency and an inoculation access period (IAP) [13]. During CLas entry into the D. citri, the midgut is an important immune barrier for defense against bacterial infection. In a previous report, draft genome sequencing revealed that D. citri contained the Toll and JAK/STAT pathways, but lacked genes for the IMD pathway response to Gram-negative bacterial infection [36]. D. citri transmits CLas, beginning from a few days to a week after acquisition, and for the lifetime of the vector, suggesting that CLas have developed mechanisms to avoid psyllid cellular and humoral immune defenses. In Galleria mellonella, Bacillus thuringiensis infection significantly influenced antioxidant activity and the level of lipid peroxidation in the larval midgut [37]. Therefore, we considered that the D. citri midgut would play critical roles in the interaction between CLas and D. citri. In this study, we performed RNA-pooling sequencing, and a total of 62,582,728, 58,399,847 and 56,575,636 raw reads were obtained from CLas-free groups, and 57,180,938, 63,335,414 and 54,887,010 raw reads were obtained from CLas-infected groups. These results showed that the number of raw reads for CLas-free groups and CLas-infected groups were similar. Previous studies revealed that pooling sequencing did not represent the population variations in gene expression levels, but it could save costs and limit starting material. However, stringent false discovery rates (FDRs) and the high-throughput validation of DEGs should be considered [38]. After assembly, a total of 778 DEGs were identified in the midgut after CLas infection according to comparative transcriptome analysis. Many pathogens can induce the upregulation of host gene expression. Xiong et al. identified many immunity-related genes encoding pattern recognition receptors, signal modulators and immune effectors after the injection of the fungal pathogen Beauveria bassiana and the Gram-negative bacterium Enterobacter cloacae [39]. Tang et al. performed transcriptome analysis in Musca domestica larvae inoculated with a mixture of Escherichia coli and Staphylococcus aureus. The results showed that many genes involved in innate immunity were induced following infection [40]. In this study, KEGG enrichment analysis showed that upregulated genes were significantly enriched in ribosomes, oxidative phosphorylation and proteasomes. Downregulated genes were significantly enriched in endocytosis.

4.1. Ubiquitination, the Immune Response and Ribosomes May Play Important Roles in the Midgut Response to CLas Infection

Ubiquitination occurs through a series of reactions catalyzed by different enzymes, including Ub-activating E1, Ub-conjugating E2 and Ub-ligase E3 [41]. Ubiquitination plays an important role in the recognition and clearance of some invading bacteria. Many bacterial effectors enable them to interfere with the host’s ubiquitination system, and thus to achieve successful infection [42]. Wang et al. confirmed that ubiquitin modification had multiple effects on the host immune system against Salmonella infection [43]. In this study, a total of 13 ubiquitination-related genes were identified; among them, eight genes were upregulated and five genes were downregulated after CLas infection. E3 ubiquitin-protein ligase was upregulated in the D. citri midgut following CLas infection. Zhang et al. revealed that zebrafish bloodthirsty member 20 with E3 ubiquitin ligase activity was involved in the immune response against bacterial infection [44]. Therefore, we speculated that E3 ubiquitin-protein ligase might play an important role in CLas infection. Many ubiquitinated proteins are recognized by the proteasome and are then further degraded [45]. The ubiquitin proteasome system is a key signaling pathway in the host response to bacterial or viral infection. Isaacson et al. revealed that host cells could utilize the ubiquitin-proteasome system to counteract viral infections through the generation of target structures recognized by T cells [46]. Yu et al. also identified eight differentially expressed proteins associated with ubiquitination in the B. mori midgut after B. mori nucleopolyhedrovirus (BmNPV) infection [47]. In this study, four proteasome-related proteins were upregulated after CLas infection, including proteasome activator complex subunit 3-like, proteasome subunit beta type, proteasome subunit alpha type and 26S proteasome non-ATPase regulatory subunit 12. We speculated that CLas bacteria invading the midgut might be recognized by the D. citri immune system. Meanwhile, CLas could activate host ubiquitination to eliminate immune-related proteins.
CLas, as Gram-negative bacteria, activate the host immune system after invading the midgut [48]. Thus, the midgut immune response may play an important role in the defense against pathogen infection [49]. In total, 11 genes related to immune responses were differentially expressed in the CLas-infected groups. Serine proteases (SPs) are involved in both the prophenoloxidase (PPO) activation cascade and the Toll immune signaling pathway, especially those with a clip domain. The proteases are secreted into the hemolymph as inactive precursors and transform into the active protein involved in the melanization reaction requiring specific proteolytic cleavage [50]. In this study, we found that three genes belonging to the serine protease family were upregulated after CLas infection, including CLIP domain-containing serine protease 2-like, CLIBP-serine protease 5 and CLIBP-serine protease 4, partial. In insects, PPO is activated through a serine protease cascade upon the recognition of pathogen-associated molecular patterns. In previous research, two PPOs were identified from the genome database of D. citri. We considered that CLIP domain-containing serine protease 2-like might regulate D. citri PPOs to activate melanization to defend against CLas infection [27].
In the process of host invasion by pathogens, many proteins are utilized to achieve replication. Ribosomal proteins, in conjunction with rRNA, make up the ribosomal subunits involved in the cellular process of translation [51]. In this study, a total of 15 DEGs related to ribosomes were screened, comparing the CLas-infected group to the CLas-free group. Among these genes, 14 were upregulated after CLas infection, except for 28S ribosomal protein S27 (mitochondrial). Yu et al. also found that some ribosome-associated proteins were upregulated after BmNPV infection [47]. These results indicate that ribosomal proteins are likely to play crucial roles in the response to CLas infection.

4.2. Endocytosis, the Cytoskeleton and Insecticide Resistance May Play Crucial Roles in the Midgut Response to CLas Infection

Clathrin-mediated endocytosis is mainly involved in the selective and facilitated internalization of cell surface receptors [52]. Many pathogens can take advantage of the endocytosis machinery in the cytosol for replication [53,54,55]. However, the specific mechanism of endocytosis between viruses and bacteria is different. Bacteria can secrete some proteins or components that allow the modification of the internalization vacuole to permit an intravacuolar lifestyle with concomitant replication [56]. For some viruses, endocytes can help transport incoming particles deep into the cytoplasm unobstructed by cytoplasmic crowding and obstacles such as the cytoskeleton [57]. Interestingly, seven genes associated with endocytosis were found, and most of them were downregulated after CLas infection, except for heat shock protein 70. During CLso infection of the potato psyllid, CLso cells were observed between the basal lamina and the basal epithelial cell membranes [58]. For example, ADP ribosylation factor proteins comprise a group of five Ras-related GTPases that are thought to function as regulators of membrane traffic. Interestingly, CLas are Gram-negative bacteria, but they are not pathogenic to D. citri. In addition, many studies revealed that CLas could be detected in the hemolymph from CLas-infected D. citri [23]. Therefore, we speculated that after CLas invading, D. citri might inhibit the expression of endocytosis-related genes in the midgut to prevent the further transmission of CLas. Moreover, CLas might avoid the host immune system by endocytosis.
The cytoskeleton is critical for the maintenance of cell shape, cell motility and intracellular transport, and bacterial and viral infections require the cytoskeleton [59]. Bacteria in the process of invasion and proliferation, at all stages of the intracellular bacterial life cycle, have the same three-dimensional cytosolic space containing the cytoskeleton [60]. In this instance, a total of 15 DEGs associated with the cytoskeleton were identified, and most were upregulated in the CLas-infected groups compared to the CLas-free groups. As a major protein constituent of cytoskeletal filaments, tubulin is involved in many vital cellular processes, including cell motility, cellular division and cytokinesis [61]. These results indicated that cytoskeleton-related genes played an important role in the process of CLas infection in D. citri.

5. Conclusions

By using transcriptome sequencing, we identified 778 genes differentially expressed in the midgut between the CLas-infected groups and CLas-free groups. KEGG and GO enrichment analyses revealed that 80 DEGs were associated with ubiquitination, the immune response, ribosomes, endocytosis, the cytoskeleton and insecticide resistance. This study provides a foundation for further research to investigate the mechanisms of CLas invasion of the D. citri midgut.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/3/171/s1, Table S1: All primers used for RT-qPCR in this study, Table S2: GO enrichment analysis of DEGs, Table S3: KEGG enrichment analysis of DEGs.

Author Contributions

Conceptualization, H.-Z.Y. and Z.-J.L.; methodology, N.-Y.L.; software, H.-Z.Y., H.-N.S. and L.Y.; investigation, X.-D.Z., C.-X.C. and J.-C.S.; writing—original draft preparation, H.-Z.Y. and N.-Y.L.; writing—review and editing, H.-Z.Y., X.-D.Y. and Z.-J.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Natural Science Foundation of China (31960116), Ministry of Science and Technology of China (2018YFD0201504), the Educational Commission of Jiangxi Province of China (GJJ180747), and Science and Technology Department of Jiangxi Province (20171BCB23074), the earmarked fund for Jiangxi Agriculture Research System (JXARS-07) and Key Research and Development Plan of Jiangxi Province (20192BBF60057).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bove, J.M. Huanglongbing: A destructive, newly-emerging, century-old disease of citrus. J. Plant Pathol. 2006, 88, 7–37. [Google Scholar]
  2. Wang, N.; Pierson, E.A.; Setubal, J.C.; Xu, J.; Levy, J.G.; Zhang, Y.; Li, J.; Rangel, L.T.; Martins, J.J. The Candidatus Liberibacter-host interface: Insights into pathogenesis mechanisms and disease control. Annu. Rev. Phytopathol. 2017, 55, 451–482. [Google Scholar] [CrossRef] [PubMed]
  3. Slisz, A.M.; Breksa, A.P.; Mishchuk, D.O.; McCollum, G.; Slupsky, C.M. Metabolomic analysis of citrus infection by ‘Candidatus Liberibacter’ reveals insights into pathogenicity. J. Proteome Res. 2012, 22, 4223–4230. [Google Scholar] [CrossRef]
  4. Hu, J.H.; Jiang, J.X.; Wang, N. Control of citrus Huanglongbing (HLB) via trunk injection of plant activators and antibiotics. Phytopathology 2018, 108, 186–195. [Google Scholar] [CrossRef] [Green Version]
  5. Qureshi, J.A.; Kostyk, B.C.; Stansly, P.A. Insecticidal suppression of Asian citrus psyllid Diaphorina citri (Hemiptera: Liviidae) vector of huanglongbing pathogens. PLoS ONE 2014, 9, e112331. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kolora, L.D.; Powell, C.M.; Hunter, W.; Bextine, B.; Lauzon, C.R. Internal extracellular bacteria of Diaphorina citri Kuwayama (Hemiptera: Psyllidae), the Asian citrus psyllid. Curr. Microbiol. 2015, 70, 710–715. [Google Scholar] [CrossRef] [PubMed]
  7. Boina, D.R.; Bloomquist, J.R. Chemical control of the Asian citrus psyllid and of huanglongbing disease in citrus. Pest Manag. Sci. 2015, 71, 808–823. [Google Scholar] [CrossRef] [PubMed]
  8. Boina, D.R.; Onagbola, E.O.; Salyani, M.; Stelinski, L.L. Antifeedant and sublethal effects of imidacloprid on Asian citrus psyllid, Diaphorina citri. Pest Manag. Sci. 2010, 65, 870–877. [Google Scholar] [CrossRef] [PubMed]
  9. Huang, J.K.; Hu, R.F.; Pray, C.; Qiao, F.B.; Rozelle, S. Biotechnology as an alternative to chemical pesticides: A case study of Bt cotton in China. Agric. Econ. 2015, 29, 55–67. [Google Scholar] [CrossRef]
  10. Coutinho-Abreu, I.V.; Forster, L.; Guda, T.; Ray, A. Odorants for surveillance and control of the Asian citrus psyllid (Diaphorina citri). PLoS ONE 2014, 9, e109236. [Google Scholar] [CrossRef]
  11. Yu, X.D.; Gowda, S.; Killiny, N. Double-stranded RNA delivery through soaking mediates silencing of the muscle protein 20 and increases mortality to the Asian citrus psyllid. Diaphorina citri. Pest Manag. Sci. 2017, 73, 1846–1853. [Google Scholar] [CrossRef] [PubMed]
  12. Lu, Z.J.; Zhou, C.H.; Yu, H.Z.; Huang, Y.L.; Liu, Y.X.; Xie, Y.X.; Wang, J.; Hu, W.; Huang, A.J.; Su, H.N.; et al. Potential roles of insect Tropomyosin1-X1 isoform in the process of Candidatus Liberibacter asiaticus infection of Diaphorina citri. J. Insect Physiol. 2019, 114, 125–135. [Google Scholar] [CrossRef] [PubMed]
  13. Inoue, H.; Ohnishi, J.; Ito, T.; Tomimura, K.; Miyata, S.; Iwanami, T.; Ashihara, W. Enhanced proliferation and efficient transmission of Candidatus Liberibacter asiaticus by adult Diaphorina citri after acquisition feeding in the nymphal stage. Ann. Appl. Biol. 2010, 155, 29–36. [Google Scholar] [CrossRef]
  14. Coletta-Filho, H.D.; Daugherty, M.P.; Ferreira, C.; Lopes, J.R. Temporal progression of ‘Candidatus Liberibacter asiaticus’ infection in citrus and acquisition efficiency by Diaphorina citri. Phytopathology 2014, 104, 416–421. [Google Scholar] [CrossRef] [Green Version]
  15. Ammar, E.D.; Shatters, R.G., Jr.; Hall, D.G. Localization of Candidatus Liberibacter asiaticus, associated with citrus huanglongbing disease, in its psyllid vector using fluorescence in situ hybridization. J. Phytopathol. 2011, 159, 726–734. [Google Scholar] [CrossRef]
  16. Gakhar, S.K.; Shandilya, H.K. Midgut specific immune response of vector mosquito Anopheles stephensi to malaria parasite Plasmodium. Indian J. Exp. Biol. 2011, 39, 287–290. [Google Scholar]
  17. Vincendeau, P.; Gobert, A.P.; Daulouede, S.; Moynet, D.; Mossalayi, M.D. Arginases in parasitic diseases. Trends Parasitol. 2003, 19, 9–12. [Google Scholar] [CrossRef]
  18. Rider, M.A.; Zou, J.; Vanlandingham, D.; Nuckols, J.T.; Higgs, S.; Zhang, Q.; Lacey, M.; Kim, J.; Wang, G.; Hong, Y.S. Quantitative proteomic analysis of the Anopheles gambiae (Diptera: Culicidae) midgut infected with o’nyong-nyong virus. J. Med. Entomol. 2013, 50, 1077–1088. [Google Scholar] [CrossRef] [Green Version]
  19. Wang, X.Y.; Yu, H.Z.; Geng, L.; Xu, J.P.; Yu, D.; Zhang, S.Z.; Ma, Y.; Fei, D.Q. Comparative transcriptome analysis of Bombyx mori (Lepidoptera) larval midgut response to BmNPV in susceptible and near-isogenic resistant strains. PLoS ONE 2016, 11, e0155341. [Google Scholar] [CrossRef]
  20. Bao, Y.Y.; Wang, Y.; Wu, W.J.; Zhao, D.; Xue, J.; Zhang, B.Q.; Shen, Z.C.; Zhang, C.X. De novo intestine-specific transcriptome of the brown planthopper Nilaparvata lugens revealed potential functions in digestion, detoxification and immune response. Genomics 2012, 99, 256–264. [Google Scholar] [CrossRef] [Green Version]
  21. Kruse, A.; Fattah-Hosseini, S.; Saha, S.; Johnson, R.; Warwick, E.; Sturgeon, K.; Mueller, L.; MacCoss, M.J.; Shatters, R.G.J.; Cilia Heck, M. 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 ONE 2017, 12, e0179531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Lu, Z.J.; Killiny, N. Huanglongbing pathogen Candidatus Liberibacter asiaticus exploits the energy metabolism and host defence responses of its vector Diaphorina citri. Physiol. Entomol. 2017, 42, 319–335. [Google Scholar] [CrossRef]
  23. Kruse, A.; Ramsey, J.S.; Johnson, R.; Hall, D.G.; MacCoss, M.J.; Heck, M. Candidatus Liberibacter asiaticus’ minimally alters expression of immunity and metabolism proteins in the hemolymph of Diaphorina citri, the insect vector of Huanglongbing. J. Proteome Res. 2018, 17, 2995–3011. [Google Scholar] [CrossRef] [PubMed]
  24. Shokrollah, H.; Abdullah, T.L.; Sijam, K.; Abdullah, S.N.A. Ultrastructure of Candidatus Liberibacter asiaticus and its damage in huanglongbing (HLB) infected citrus. Afr. J. Biotechnol. 2010, 9, 5897–5901. [Google Scholar]
  25. Saha, S.; Hosmani, P.S.; Villalobos-Ayala, K.; Miller, S.; Shippy, T.; Flores, M.; Rosendale, A.; Cordola, C.; Bell, T.; Mann, H.; et al. Improved annotation of the insect vector of citrus greening disease: Biocuration by a diverse genomics community. Database 2017, 2017, bax032. [Google Scholar] [CrossRef]
  26. Trapnell, C.; Pachter, L.; Salzberg, S.L. TopHatL discovering splice junctions with RNA-Seq. Bioinformatics 2009, 25, 1105–1111. [Google Scholar] [CrossRef]
  27. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwiqht, S.S.; Eppiq, J.T.; et al. Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [Green Version]
  28. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  29. Zheng, Q.; Wang, X.J. GOEAST: A web-based software toolkit for Gene Ontology enrichment analysis. Nucleic Acids Res. 2008, 36, 358–363. [Google Scholar] [CrossRef] [Green Version]
  30. Yu, G.C.; Wang, L.G.; Han, Y.Y.; He, Q.Y. Clusterprofiler: An R Package for comparing biological themes among gene clusters. Omics 2012, 16, 284–287. [Google Scholar] [CrossRef]
  31. Cheraghi, E.; Mehranjani, M.S.; Shariatzadeh, S.M.A.; Esfahani, M.H.N.; Alani, B. N-acetylchsteine compared to metformin, improves the expression profiles of growth differentiation factor-9 and receptor tyrosine kinase c-kit in the oocytes of patients with polycystic ovarian syndrome. Int. J. Fertil. Steril. 2018, 11, 270–278. [Google Scholar] [PubMed]
  32. Camon, E.; Magrane, M.; Barrell, D.; Lee, V.; Dimmer, E.; Maslen, J.; Binns, D.; Harte, N.; Lopez, R.; Apweiler, R. The gene ontology annotation (GOA) database: Sharing knowledge in uniprot with gene ontology. Nucleic Acids Res. 2004, 32, D262–D266. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Kanehisa, M.; Goto, S.; Kawashima, S.; Nakaya, A. The KEGG databases at GenomeNet. Nucleic Acids Res. 2002, 30, 42–46. [Google Scholar] [CrossRef] [PubMed]
  34. Zheng, Z.; Xu, M.R.; Bao, M.L.; Wu, F.N.; Chen, J.C.; Deng, X.L. Unusual five copies and dual Forms of nrdB in “Candidatus Liberibacter asiaticus”: Biological implications and PCR detection application. Sci. Rep. 2016, 6, 39020. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Grafton-Cardwell, E.E.; Stelinski, L.L.; Stansly, P.A. Biology and management of Asian citrus psyllid, vector of the huanglongbing pathogens. Annu. Rev. Entomol. 2013, 58, 413–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Arp, A.P.; Hunter, W.B.; Pelz-Stelinski, K.S. Annotation of the Asian citrus psyllid genome reveals a reduced innate immune system. Front. Physiol. 2016, 7, 570. [Google Scholar] [CrossRef] [Green Version]
  37. Dubovskiy, I.M.; Martemyanov, V.V.; Vorontsova, Y.L.; Rantala, M.J.; Gryzanova, E.V.; Glupov, V.V. Effect of bacterial infection on antioxidant activity and lipid peroxidation in the midgut of Galleria mellonella L. larvae (Lepidoptera, Pyralidae). Comp. Biochem. Physiol. C Toxicol. Pharmacol. 2008, 48, 1–5. [Google Scholar] [CrossRef]
  38. Rajkumar, A.P.; Qvist, P.; Lazarus, R.; Lescai, F.; Ju, J.; Nyegaard, M.; Mors, O.; Børglum, A.D.; Li, Q.B.; Christensen, J.H. Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq. BMC Genom. 2015, 16, 548. [Google Scholar] [CrossRef] [Green Version]
  39. Xiong, G.H.; Xing, L.S.; Lin, Z.; Saha, T.T.; Wang, C.; Jiang, H.B.; Zou, Z. High throughput profiling of the cotton bollworm Helicoverpa armigera immunotranscriptome during the fungal and bacterial infections. BMC Genom. 2015, 16, 321. [Google Scholar] [CrossRef] [Green Version]
  40. Tang, T.; Li, X.; Yang, X.; Yu, X.; Wang, J.H.; Liu, F.S.; Huang, D.W. Transcriptional response of Musca domestica larvae to bacterial infection. PLoS ONE 2014, 9, e104867. [Google Scholar] [CrossRef] [Green Version]
  41. Koegl, M.; Hoppe, T.; Schlenker, S.; Ulrich, H.D.; Mayer, T.U.; Jentsch, S. A novel ubiquitination factor, E4, is involved in multiubiquitin chain assembly. Cell 1999, 96, 635–644. [Google Scholar] [CrossRef] [Green Version]
  42. Dupont, N.; Temime-Smaali, N.; Lafont, F. How ubiquitination and autophagy participate in the regulation of the cell response to bacterial infection. Biol. Cell 2010, 102, 621–634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Wang, L.D.; Yan, J.; Niu, H.; Huang, R.; Wu, S.Y. Autophagy and ubiquitination in Salmonella infection and the related inflammatory responses. Front. Cell Infect Microbiol. 2018, 8, 78. [Google Scholar] [CrossRef] [PubMed]
  44. Zhang, X.S.; Zhao, H.; Chen, Y.Y.; Luo, H.Y.; Yang, P.L.; Yao, B. A zebrafish (Danio rerio) bloodthirsty member 20 with E3 ubiquitin ligase activity involved in immune response against bacterial infection. Biochem. Biophys. Res. Commun. 2015, 457, 83–89. [Google Scholar] [CrossRef] [PubMed]
  45. Grice, G.L.; Nathan, J.A. The recognition of ubiquitinated proteins by the proteasome. Cell Mol. Life Sci. 2016, 73, 3497–3506. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Isaacson, M.K.; Ploegh, H.L. Ubiquitination, ubiquitin-like modifiers, and deubiquitination in viral infection. Cell Host Microbe 2009, 5, 559–570. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Yu, H.Z.; Wang, X.Y.; Xu, J.P.; Ma, Y.; Zhang, S.Z.; Yu, D.; Fei, D.Q.A. Muhammad, iTRAQ-based quantitative proteomics analysis of molecular mechanisms associated with Bombyx mori (Lepodoptera) larval midgut response to BmNPV in susceptible and near-isogenic strains. J. Proteom. 2017, 8, 35–50. [Google Scholar] [CrossRef]
  48. Wang, N.; Trivedi, P. Citrus huanglongbing: A newly relevant disease presents unprecedented challenges. Phytopathology 2013, 103, 652–665. [Google Scholar] [CrossRef] [Green Version]
  49. Crava, C.M.; Jacubowska, A.K.; Escriche, B.; Herrero, S.; Bel, Y. Dissimilar regulation of antimicrobial proteins in the midgut of Spodoptera exigua larvae challenged with Bacillus thuringiensis toxins or baculovirus. PLoS ONE 2015, 10, e0125991. [Google Scholar] [CrossRef]
  50. Chu, Y.; Liu, Y.; Shen, D.X.; Hong, F.; Wang, G.R.; An, C.J. Serine proteases SP1 and SP13 mediate the melanization response of Asian corn borer, Ostrinia furnacalis, against entomopathogenic fungus Beauveria bassiana. J. Invertebr. Pathol. 2015, 128, 64–72. [Google Scholar] [CrossRef]
  51. Rodnina, M.V.; Wolfgang, W. The ribosome as a molecular machine: The mechanism of tRNA-mRNA movement in translocation. Biochem. Soc. Trans. 2011, 39, 658–662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Brodsky, F.M.; Guagliardi, L.E.; Galloway, C.J.; Koppelman, B.; Acton, S.L.; Wong, D.H.; Nathke, I. Endocytosis: Molecular mechanism and role in the immune response. Endocytosis 1992, 62, 343–353. [Google Scholar]
  53. Nonnenmacher, M.; Weber, T. Adeno-associated virus 2 infection requires endocytosis through the CLIC/GEEC pathway. Cell Host Microbe 2011, 10, 563–576. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Pauza, C.D.; Price, T.M. Human immunodeficiency virus infection of T cells and monocytes proceeds via receptor-mediated endocytosis. J. Cell Biol. 1988, 107, 959–968. [Google Scholar] [CrossRef] [PubMed]
  55. Branza-Nichita, N.; Macovei, A.; Lazar, C. Caveolae-Dependent Endocytosis in Viral Infection; Ceresa, B., Ed.; Intech: London, UK, 2012. [Google Scholar]
  56. Kumar, Y.; Valdivia, R.H. Leading a sheltered life: Intracellular pathogens and maintenance of vacuolar compartments. Cell Host Microbe 2009, 5, 593–601. [Google Scholar] [CrossRef] [Green Version]
  57. Vazquez-Calvo, A.; Sobrino, F.; Martin-Acebes, M.A. Plasma membrane phosphatidylinositol 4, 5 bisphosphate is required for internalization of foot-and mouth disease virus and vesicular stomatitis virus. PLoS ONE 2012, 7, e45172. [Google Scholar] [CrossRef] [Green Version]
  58. Cicero, J.M.; Fisher, T.W.; Qureshi, J.A.; Stansly, P.A.; Brown, J.K. Colonization and intrusive invasion of potato psyllid by ‘Candidatus Liberibacter solanacearum’. Phytopathology 2017, 107, 36–47. [Google Scholar] [CrossRef] [Green Version]
  59. Harries, P.A.; Schoelz, J.E.; Nelson, R.S. Intracellular transport of viruses and their components: Utilizing the cytoskeleton and membrane highways. Mol. Plant Microbe Interact. 2010, 23, 1381–1393. [Google Scholar] [CrossRef]
  60. Colonne, P.M.; Winchell, C.G.; Voth, D.E. Hijacking host cell highways: Manipulation of the host actin cytoskeleton by obligate intracellular bacterial pathogens. Front. Cell. Infect. Microbiol. 2016, 6, 107. [Google Scholar] [CrossRef] [Green Version]
  61. Lu, Z.J.; Hu, H.; Killiny, N. Proteomic maps of subcellular protein fractions of the Asian citrus psyllid Diaphorina citri, the vector of citrus huanglongbing. Physiol. Entomol. 2017, 42, 36–64. [Google Scholar] [CrossRef] [Green Version]
Insects 11 00171 g001 550
Figure 1. The structure of the D. citri midgut, and the PCR detection of Candidatus Liberibacter asiaticus (CLas). (A) The structure of the whole D. citri midgut was observed under a Lecia S8AP0 stereomicroscope. (B) Agarose gel electrophoresis analysis in CLas-free groups and CLas-infected groups. M, 2000 DNA marker; PC, positive control; NC, negative control.
Figure 1. The structure of the D. citri midgut, and the PCR detection of Candidatus Liberibacter asiaticus (CLas). (A) The structure of the whole D. citri midgut was observed under a Lecia S8AP0 stereomicroscope. (B) Agarose gel electrophoresis analysis in CLas-free groups and CLas-infected groups. M, 2000 DNA marker; PC, positive control; NC, negative control.
Insects 11 00171 g001
Insects 11 00171 g002 550
Figure 2. The identification and hierarchical cluster analysis of differentially expressed genes. (A) A scatter diagram for each gene. The blue, red and green points represent no difference in expression, upregulated genes and downregulated genes, respectively. (B) Hierarchical clustering of differentially expressed genes (DEGs) between the CLas-free groups and CLas-infected groups. Columns indicate different samples. Rows represent different DEGs. Blue bands indicate a low expression level, and red bands indicate a high gene expression level.
Figure 2. The identification and hierarchical cluster analysis of differentially expressed genes. (A) A scatter diagram for each gene. The blue, red and green points represent no difference in expression, upregulated genes and downregulated genes, respectively. (B) Hierarchical clustering of differentially expressed genes (DEGs) between the CLas-free groups and CLas-infected groups. Columns indicate different samples. Rows represent different DEGs. Blue bands indicate a low expression level, and red bands indicate a high gene expression level.
Insects 11 00171 g002
Insects 11 00171 g003 550
Figure 3. Gene Ontology (GO) enrichment analysis of DEGs. A scatter diagram of GO categories. The X-axis indicates the gene ratio. The Y-axis indicates different categories. (A) Upregulated DEGs. (B) Downregulated DEGs.
Figure 3. Gene Ontology (GO) enrichment analysis of DEGs. A scatter diagram of GO categories. The X-axis indicates the gene ratio. The Y-axis indicates different categories. (A) Upregulated DEGs. (B) Downregulated DEGs.
Insects 11 00171 g003
Insects 11 00171 g004 550
Figure 4. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of DEGs. A scatter diagram of KEGG pathways. The X-axis indicates the gene ratio. The Y-axis indicates different pathways. (A) Upregulated DEGs. (B) Downregulated DEGs.
Figure 4. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis of DEGs. A scatter diagram of KEGG pathways. The X-axis indicates the gene ratio. The Y-axis indicates different pathways. (A) Upregulated DEGs. (B) Downregulated DEGs.
Insects 11 00171 g004
Insects 11 00171 g005 550
Figure 5. The correlation between the gene expression ratios obtained from the transcriptome data and RT-qPCR data. (A) The differential expression levels of 14 differentially expressed genes in the CLas-free and CLas-infected D. citri midgut. The relative expression levels were calculated using the 2−∆∆Ct method. Statistical analysis was performed using the SPSS software. The significant differences are indicated by * (p < 0.05) or ** (p < 0.01); (B) Lineage analysis between the transcriptome and RT-qPCR data. The ratios obtained by RT-qPCR (Y-axis) were plotted against the ratios obtained by the transcriptome (X-axis).
Figure 5. The correlation between the gene expression ratios obtained from the transcriptome data and RT-qPCR data. (A) The differential expression levels of 14 differentially expressed genes in the CLas-free and CLas-infected D. citri midgut. The relative expression levels were calculated using the 2−∆∆Ct method. Statistical analysis was performed using the SPSS software. The significant differences are indicated by * (p < 0.05) or ** (p < 0.01); (B) Lineage analysis between the transcriptome and RT-qPCR data. The ratios obtained by RT-qPCR (Y-axis) were plotted against the ratios obtained by the transcriptome (X-axis).
Insects 11 00171 g005
Insects 11 00171 g006 550
Figure 6. Hierarchical analysis for DEGs related to ubiquitination, the immune response and ribosomes between CLas-free and CLas-infected D. citri midgut. DEG expression is shown with a pseudocolor scale (from −3 to 3), with a red color indicating high expression levels and a green color indicating low expression. Each group represents three biological replicates.
Figure 6. Hierarchical analysis for DEGs related to ubiquitination, the immune response and ribosomes between CLas-free and CLas-infected D. citri midgut. DEG expression is shown with a pseudocolor scale (from −3 to 3), with a red color indicating high expression levels and a green color indicating low expression. Each group represents three biological replicates.
Insects 11 00171 g006
Insects 11 00171 g007 550
Figure 7. Hierarchical analysis for DEGs associated with endocytosis, the cytoskeleton and insecticide resistance between CLas-free and CLas-infected D. citri midgut. A hierarchical cluster analysis was performed using the Genesis software. DEG expression is shown with a pseudocolour scale (from −3 to 3), with the red colour indicating high expression levels and the green color indicating low expression. Each group represents three biological replicates.
Figure 7. Hierarchical analysis for DEGs associated with endocytosis, the cytoskeleton and insecticide resistance between CLas-free and CLas-infected D. citri midgut. A hierarchical cluster analysis was performed using the Genesis software. DEG expression is shown with a pseudocolour scale (from −3 to 3), with the red colour indicating high expression levels and the green color indicating low expression. Each group represents three biological replicates.
Insects 11 00171 g007
Table
Table 1. A summary of the transcriptomes in the different treatments in D. citri.
Table 1. A summary of the transcriptomes in the different treatments in D. citri.
SampleCLas-Free-1CLas-Free-2CLas-Free-3CLas-Infected-1CLas-Infected-2CLas-Infected-3
Total raw reads62,582,72858,399,84756,575,63657,180,93863,335,41454,887,010
Total clean reads61,521,78057,557,15655,768,65255,707,04862,309,80053,907,044
Q2097.25%97.14%97.37%97.59%97.16%97.48%
Q3092.45%92.20%92.72%93.13%92.24%92.90%
GC percent40.57%40.46%40.91%37.74%39.94%39.39%
Total map43,131,403 (70.11%)40,754,962 (70.81%)40,200,207 (72.08%)32,534,301 (58.4%)41,738,828 (66.99%)35,207,811 (65.31%)
Unique map34,876,548 (56.69%)32,834,545 (57.05%)32,347,622 (58.0%)26,881,942 (48.26%)33,556,647 (53.85%)28,376,185 (52.64%)
Multiple map8,254,855 (13.42%)7,920,417 (13.76%)7,852,585 (14.08%)5,652,359 (10.15%)8,182,181 (13.13%)6,831,626 (12.67%)
Read 1 map17,503,161 (28.45%)16,493,012 (28.66%)16,205,676 (29.06%)13,455,363 (24.15%)16,865,739 (27.07%)14,200,948 (26.34%)
Read 2 map17,373,387 (28.24%)16,341,533 (28.39%)16,141,946 (28.94%)13,426,579 (24.1%)16,690,908 (26.79%)14,175,237 (26.3%)
Positive map17,505,192 (28.45%)16,467,596 (28.61%)16,228,967 (29.1%)13,459,583 (24.16%)16,829,779 (27.01%)14,230,468 (26.4%)
Table
Table 2. Differentially expressed genes upon CLas infection involved in ubiquitination, the immune response and the ribosomes comparing CLas-infected groups and CLas-free groups.
Table 2. Differentially expressed genes upon CLas infection involved in ubiquitination, the immune response and the ribosomes comparing CLas-infected groups and CLas-free groups.
Gene IDGene DescriptionCLas-Free FPKMCLas-Infected FPKMLog2Fold Change
Ubiquitination
DcitrP024565.1Proteasome subunit alpha type8.59062418.910981.142678
DcitrP057820.1S-phase kinase-associated protein 1 5.0743512.833141.338850457
DcitrP057935.1Ubiquitin-fold modifier-conjugating enzyme 17.19537114.553071.015273471
DcitrP026550.1Ubiquitin-associated domain-containing protein 110.655573.702037−1.531855619
DcitrP066965.1Ubiquitin-like modifier-activating enzyme ATG729.4318414.80442−0.993594896
DcitrP055490.1E3 ubiquitin-protein ligase0.3330952.286882.782229875
DcitrP062105.1Proteasome activator complex subunit 3-like protein2.29000213.834692.595512129
DcitrP059695.1Proteasome subunit beta type1.7535263.5879192.414010482
DcitrP081240.1Ubiquitin domain-containing protein UBFD18.09510315.501530.93331621
DcitrP077930.126S proteasome non-ATPase regulatory subunit 1211.590520.462840.818448568
DcitrP069975.1E3 ubiquitin-protein ligase MARCH278.948138.47541−1.034157481
DcitrP069805.1E3 ubiquitin-protein ligase RNF8203.47397.33118−1.064416343
DcitrP080605.1Ubiquitin-conjugating enzyme E2 J224.04910.95165−1.132484873
Immune response
DcitrP022770.1CLIP domain-containing serine protease 2-like0.150492.3549213.964562534
DcitrP022775.1CLIPB-serine protease 40.5699215.4859223.275711152
DcitrP046955.1Aminopeptidase0.2232171.0707382.276052354
DcitrP016435.1Lysozyme i-type4.80354619.426172.011353707
DcitrP097315.1CLIPB-serine protease 52.7448376.3969241.223892409
DcitrP029315.1CLIPB-serine protease 5-RA1.5519073.4888691.163879263
DcitrP022220.1Apolipophorin-Ⅲ82.6109184.93531.162289806
DcitrP036745.1Toll-like 8B66.7779239.1133−0.771504993
DcitrP077945.1Beat protein1.7320960.063739−4.609706023
DcitrP077305.1Cathespin F-like protein 5-RA0.0145750.05788435.280394143
DcitrP036095.1Spondin-2-like0.0360880.6121924.060891064
Ribosome
DcitrP073005.160S ribosomal protein L2328.24808126.53572.162205397
DcitrP093430.1Ribosomal protein L27A98.15065435.65472.152761801
DcitrP013835.140S ribosomal protein S94.20966213.430441.666951923
DcitrP008150.160S ribosomal protein L3417.6527651.353971.536118663
DcitrP059870.140S ribosomal protein S24155.3869427.22291.458686099
DcitrP062275.1Ribosomal protein L15101.5892271.75711.418811677
DcitrP059845.140S ribosomal protein S2421.1147355.262061.385121791
DcitrP068410.128S ribosomal protein S18a, mitochondrial43.2703102.0781.235994788
DcitrP007765.1Ribosomal RNA processing protein 36 like protein4.52550910.31741.186528712
DcitrP019375.160S ribosomal protein L21206.8503463.48051.16374505
DcitrP027780.140S ribosomal protein S9419.7915883.73741.073779521
DcitrP025515.160S ribosomal protein L32533.56681100.2741.044015195
DcitrP002745.140S ribosomal protein S11102.8649207.19521.009873361
DcitrP094130.160S ribosomal protein L1348.4580296.809030.997246375
DcitrP037015.128S ribosomal protein S27, mitochondrial6.2450213.043301−1.04281238
Table
Table 3. Differentially expressed genes upon CLas infection involved in endocytosis, the cytoskeleton and insecticide resistance, comparing CLas-infected groups and CLas-free groups.
Table 3. Differentially expressed genes upon CLas infection involved in endocytosis, the cytoskeleton and insecticide resistance, comparing CLas-infected groups and CLas-free groups.
Gene IDGene DescriptionCLas-Free FPKMCLas-Infected FPKMLog2Fold Change
Endocytosis
DcitrP017415.1Heat shock protein 704.97065834.092952.779406588
DcitrP097875.1Multivesicular body subunit 12B24.385668.356045−1.541534978
DcitrP031105.1WASH complex subunit 7-like protein17.458947.028736−1.315473724
DcitrP034560.1Partioning defective 6 like protein gamma26.12179.444002−1.46961334
DcitrP007925.1ADP-ribosylation factor16.665526.056527−1.459442252
DcitrP079685.1Arf-GAP with SH3 domain6.1497812.706943−1.176262825
novel.2039Uncharacterized LOC11346621948.3890123.34858−1.053542718
Cytoskeleton
DcitrP027020.1Paramyosin-like protein6.93305615.753051.185025088
DcitrP028800.1Troponin I29.919867.876141.18133546
DcitrP031990.1Tyrosine protein kinase receptor torso2.7711041.295568−1.104829478
DcitrP084355.1Troponin I0.4333393.1635642.881484544
DcitrP056675.1Alpha tubulin 1 chain-like9.99305247.739732.256687099
DcitrP062610.1Myosin light chain kinase0.720032.9303622.021891929
DcitrP046925.1Transgelin9.71144822.083031.186590264
DcitrP015360.1Kinesin-like protein1.971324.459141.180135096
DcitrP016785.1Alpha tubulin 1 chain-like 6.09683153.220633.126504327
DcitrP034805.1Gelsolin 17.55555120.65492.780374779
DcitrP011480.1Troponin C0.9288416.2300012.755828695
DcitrP027235.1Dynein light chain Tctex-type 10.844665.5036712.705884305
DcitrP081065.1Alpha tubulin 1 chain-like 1.5226677.6398562.330954816
DcitrP062790.1Keratin, type I cytoskeletal 93.79760610.281381.432892752
DcitrP002545.1Dynein heavy chain1.3178320.439357−1.557423656
Insecticide resistance
DcitrP079275.1Cytochrome P450 CYP4C40.3318856.3045644.238246725
DcitrP045800.1Probable cytochrome P450 6a140.1728212.7344634.010391446
DcitrP079270.1CYP301A1 partial1.01965911.57463.500092719
DcitrP079255.1CYP301A1 partial2.04563213.910842.763107035
DcitrP039115.1ATP-binding cassette sub-family G member 50.3495811.9146872.457939089
DcitrP019570.1Glutathione S-transferase108.29310.52191.519243483
DcitrP022100.1CYP3167A2 partial3.103118.0452371.366858366
DcitrP071995.1CYP3172A5 partial1.8477494.734521.355174187
DcitrP019590.1Glutathione S-transferase14.9414337.157581.311635196
DcitrP069025.1ATP-binding cassette sub-family G member 43.4126528.3487521.294131073
DcitrP042025.1ATP-binding cassette sub-family A member21.795379.547165−1.187067336
DcitrP090960.1Multidrug resistance-associated protein 156.7160422.4627−1.335398633
DcitrP013660.1CYP6KB117.943546.599974−1.448082102
DcitrP078955.1ATP-binding cassette sub-family A member 381.791969.226527−1.549025159
DcitrP051650.1Carboxylic ester hydrolase8.5384772.508855−1.773370392

Share and Cite

MDPI and ACS Style

Yu, H.-Z.; Li, N.-Y.; Zeng, X.-D.; Song, J.-C.; Yu, X.-D.; Su, H.-N.; Chen, C.-X.; Yi, L.; Lu, Z.-J. Transcriptome Analyses of Diaphorina citri Midgut Responses to Candidatus Liberibacter Asiaticus Infection. Insects 2020, 11, 171. https://doi.org/10.3390/insects11030171

AMA Style

Yu H-Z, Li N-Y, Zeng X-D, Song J-C, Yu X-D, Su H-N, Chen C-X, Yi L, Lu Z-J. Transcriptome Analyses of Diaphorina citri Midgut Responses to Candidatus Liberibacter Asiaticus Infection. Insects. 2020; 11(3):171. https://doi.org/10.3390/insects11030171

Chicago/Turabian Style

Yu, Hai-Zhong, Ning-Yan Li, Xiang-Dong Zeng, Jian-Chun Song, Xiu-Dao Yu, Hua-Nan Su, Ci-Xiang Chen, Long Yi, and Zhan-Jun Lu. 2020. "Transcriptome Analyses of Diaphorina citri Midgut Responses to Candidatus Liberibacter Asiaticus Infection" Insects 11, no. 3: 171. https://doi.org/10.3390/insects11030171

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop