Next Article in Journal
SGCD Missense Variant in a Lagotto Romagnolo Dog with Autosomal Recessively Inherited Limb-Girdle Muscular Dystrophy
Next Article in Special Issue
Genome-Wide Identification and Characterization of the WRKY Gene Family in Cucurbita maxima
Previous Article in Journal
Biological Embedding of Early-Life Adversity and a Scoping Review of the Evidence for Intergenerational Epigenetic Transmission of Stress and Trauma in Humans
Previous Article in Special Issue
Genome-Wide Identification and Functional Analysis of NAP1 in Triticum aestivum
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection

1
College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
State Key Laboratory for Managing Biotic and Chemical Threats to the Quality and Safety of Agro-Products, Key Laboratory of Biotechnology in Plant Protection of Ministry of Agriculture and Zhejiang Province, Institute of Plant Virology, Ningbo University, Ningbo 315211, China
*
Authors to whom correspondence should be addressed.
Genes 2023, 14(8), 1640; https://doi.org/10.3390/genes14081640
Submission received: 15 July 2023 / Revised: 9 August 2023 / Accepted: 17 August 2023 / Published: 18 August 2023
(This article belongs to the Special Issue Genetic Regulation of Biotic Stress Responses)

Abstract

:
The whitefly Bemisia tabaci is one of the most destructive pests worldwide, and causes tremendous economic losses. Tobacco Nicotiana tabacum serves as a model organism for studying fundamental biological processes and is severely damaged by whiteflies. Hitherto, our knowledge of how tobacco perceives and defends itself against whiteflies has been scare. In this study, we analyze the gene expression patterns of tobacco in response to whitefly infestation. A total of 244 and 2417 differentially expressed genes (DEGs) were identified at 12 h and 24 h post whitefly infestation, respectively. Enrichment analysis demonstrates that whitefly infestation activates plant defense at both time points, with genes involved in plant pattern recognition, transcription factors, and hormonal regulation significantly upregulated. Notably, defense genes are more intensely upregulated at 24 h post infestation than at 12 h, indicating an increased immunity induced by whitefly infestation. In contrast, genes associated with energy metabolism, carbohydrate metabolism, ribosomes, and photosynthesis are suppressed, suggesting impaired plant development. Taken together, our study provides comprehensive insights into how plants respond to phloem-feeding insects, and offers a theoretical basis for better research on plant–insect interactions.

1. Introduction

Plants and insects have engaged in a coevolutionary relationship spanning hundreds of millions of years, leading to the development of diverse strategies for combatting each other [1]. In the face of herbivorous insect attacks, plants have evolved various defense mechanisms. On the one hand, they produce specialized physical barriers, secondary metabolites, and toxic compounds to counteract herbivore assaults [2]. On the other hand, plants have developed inducible defenses, which can be categorized as direct or indirect defense strategies [3]. Indirect defense involves the production and emission of herbivore-induced plant volatiles (HIPVs), which serve as crucial cues for parasitoids and predators, ultimately reducing damage inflicted on the plant [4]. Conversely, direct defense mechanisms encompass physical barriers like leaf surface wax, thorns, and trichomes, constituting the first line of defense against herbivores [1]. Additionally, plants activate specific genes associated with the synthesis of secondary metabolites. Noteworthy compounds such as gossypol, glucosinolates, alkaloids, phenolics, and proteinase inhibitors have been identified to impede the development of herbivores [5,6].
The plant defense system is intricately regulated by sophisticated signaling networks, which involve phytohormone signaling pathways, reactive oxygen species (ROS) generation, transcription factors (TFs), kinase cascades, and mitogen-activated protein kinases (MAPKs) [7]. Pattern recognition receptors (PRRs) serve as the first line of defense, capable of recognizing conserved microbe-associated molecular patterns (MAMPs) and damage-associated molecular patterns (DAMPs), thereby activating primary innate immunity [8,9,10]. These receptors, located on the surface of plant cells, comprise receptor-like kinases (RLKs) and receptor-like proteins (RLPs), which can sense environmental changes, and transmit this information via activated signaling pathways to trigger adaptive responses [11,12]. In addition to PRRs, TFs play a crucial role in plant defense systems, including WRKY TFs, zinc finger TFs, MYB TFs, and ethylene-responsive TFs (ERFs). Previous research has demonstrated the involvement of TFs in plant–insect interactions, regulating herbivore-induced plant defense responses. For instance, in Arabidopsis thaliana, AtMYC2 was identified as a regulator of JA-mediated plant herbivore resistance [13], while OsWRKY53 activates rice defenses against brown planthoppers by activating an H2O2 burst and suppressing ethylene biosynthesis [14].
The whitefly B. tabaci (Gennadius) (Hemiptera: Aleyrodidae) is a prominent piercing–sucking herbivore known for causing significant crop losses through direct feeding and transmitting viruses [15,16]. As a highly polyphagous species, B. tabaci exhibits a wide host range, including ornamentals and greenhouse crops such as tomato, pepper, beans, eggplant, cucumber, and tobacco [17]. Among these crops, tobacco (N. tabacum), which serves as a major non-food species and a model organism, produces an array of plant defense metabolites that are toxic to most herbivorous insects. However, the mechanisms by which whiteflies manipulate plant defenses to adapt to tobacco plants, as well as the tobacco response to whitefly infestation, remain unclear. The advent of transcriptome technologies has provided an ideal method for studying plant physiological processes in response to B. tabaci. In this study, we conducted deep RNA sequencing of tobacco plants at different time points following whitefly infestation, unveiling the complex interaction between whiteflies and their tobacco hosts.

2. Materials and Methods

2.1. Plants and Whitefly

We utilized tobacco plants of the N. tabacum cv. NC89 variety, which were cultivated in climate chambers at a temperature of 24 ± 2 °C, under a 14 h light and 10 h dark cycle, and with 50–70% humidity. The seeds for these plants were provided by Li-Rong Pan from the Institute of Biotechnology, Zhejiang University.
The Mediterranean (MED) species of the B. tabaci complex were used in this study. The insects were established from a field whitefly sample collected in Suzhou, China, in 2019. The whitefly colony was maintained on tobacco plants of the N. tabacum cv. NC89 variety. The whiteflies were reared in cages placed inside a climate chamber at a temperature of 27 ± 2 °C, under a 14 h light and 10 h dark cycle, and with 50–70% humidity.

2.2. Plant Treatments

To investigate the impact of whitefly feeding on tobacco plants, we subjected the plants to different treatments. Briefly, the tobacco plants were transferred to a climate chamber (27 ± 2 °C) one day before treatments. MED1 whitefly adults were released onto tobacco plants at the 4–5 leaf stage. The whiteflies were allowed to feed on the plants for 12 h and 24 h, respectively, under the same growth conditions. Non-infested tobacco plants were used as controls. Three biological replicates were performed.

2.3. RNA Isolation and Illumina Sequencing

Total RNA was isolated from the plant samples using a TRIzol Total RNA Isolation Kit (Takara, Dalian, China), according to the manufacturer’s instructions. Three biological replicates were performed for each treatment.
The extracted RNA samples were sent to Novogene (Tianjin, China) for RNA sequencing. RNA integrity was assessed using the Fragment Analyzer 5400 (Agilent Technologies, Palo Alto, CA, USA). After quantifying the total RNA, nine cDNA libraries were constructed according to the manufacturer’s recommendations of NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, Boston, MA, USA) and index codes were added to attribute sequences to each sample. Briefly, the mRNA was isolated from the total RNA using poly-T oligo-attached magnetic beads. The fragmentation of the mRNA was performed under elevated temperatures using divalent cations in NEBNext First Strand Synthesis Reaction Buffer (5 × X). The first strand cDNA synthesis was carried out using a random hexamer primer and M-MuLV Reverse Transcriptase. Subsequently, the second strand cDNA synthesis was performed using DNA Polymerase I and RNase H. The remaining overhangs were then converted into blunt ends using exonuclease/polymerase activities. Following the adenylation of the 3’ ends of DNA fragments, NEBNext Adaptors with hairpin loop structures were ligated to prepare the fragments for hybridization. To select cDNA fragments with a preferred length of 250–300 bp, the library fragments were purified using the AMPure XP system (Beckman Coulter, Beverly, MA, USA). Subsequently, 3 µL of USER Enzyme (NEB, Ipswich, MA, USA) was applied to the size-selected, adaptor-ligated cDNA and incubated at 37 °C for 15 min, followed by 5 min at 95 °C, prior to PCR amplification. PCR was performed using Phusion High-Fidelity DNA polymerase, Universal PCR primers, and Index (X) Primer. Finally, the PCR products were purified using the AMPure XP system, and the quality of the libraries was assessed using the Agilent Bioanalyzer 2100 system.
The index-coded samples were clustered using the TruSeq PE Cluster Kit v3-cBot-HS (Illumina, San Diego, CA, USA) on a cBot Cluster Generation System, following the manufacturer’s instructions. The sequencing was performed using the Illumina Novaseq 6000 platform with 150-bp paired-end reads. For each sample, cDNA library construction and Illumina sequencing were carried out with three technical replicates. The output data were submitted to the National Genomics Data Center under accession number PRJCA018862.
The internal software Fastp version 0.23.1 [18] was used to obtain clean reads by removing low-quality reads that contained adapters, empty reads or reads with unknown sequences “N” from the raw data. The quality control of reads was assessed using the FastQC program (version 0.11.3) [19]. This program executed a series of analysis modules on raw data and generated a report with statistics for the data analyzed, including the error rate of average base sequencing and the percentage of bases with Phred quality scores > 20 (Q20) or >30 (Q30). Subsequently, the clean reads from each cDNA library were aligned to the reference genome sequences of N. tabacum (https://solgenomics.net/ftp/genomes/Nicotiana_tabacum/edwards_et_al_2017/assembly/, accessed on 1 September 2022) using Hierarchical Indexing for Spliced Alignment of Transcripts (HISAT2) [20]. The low-quality alignments were filtered with Sequence Alignment/Map tools (SAMtools) [21].

2.4. DEGs Expression Analysis

Transcripts per million (TPM) expression values were calculated using cufflink [22]. The DESeq2 (http://bioconductor.org/packages/release/bioc/html/DESeq2) was used for identifying the differentially expressed genes with default parameters. Based on these statistical analyses, genes meeting the criteria of p-value < 0.01 and an absolute value of the log2 ratio > 1 were deemed significant differentially expressed genes (DEGs). To gain insights into the potential functions and metabolic pathways of the DEGs, further analysis was conducted for GO and KEGG enrichments using TBtools software v1.0697 [23]. In this software, enriched p-values were calculated according to one-sided hypergeometric test P = 1 i = 0 m 1 M i N M n i N n , with N representing the number of the gene with KEGG/GO annotation, n representing the number of DEGs in N, M representing the number of genes in each KEGG/GO term, and m representing the number of DEGs in each KEGG/GO term. The significant pathways were defined based on a corrected p-value ≤ 0.05. PCA was used to reveal overall differences in gene expression patterns among different transcriptomes, and R function plotPCA (github.com/franco-ye/TestRepository/blob/main/PCA_by_deseq2.R, accessed on 1 September 2022) was also used.

2.5. Quantitative Real-Time PCR (qRT-PCR)

Candidate transcriptomic genes were validated with qRT-PCR analysis using three biological replicates. Total RNA was extracted as described above. cDNAs were synthesized using the HiScript II Reverse Transcriptase (Vazyme, Nanjing, China). qRT-PCR was performed on a Roche Light Cycler® 480 Real-Time PCR System (Roche Diagnostics, Mannheim, Germany) using the Hieff qPCR SYBR Green Master Mix (Low Rox Plus) (Yeasen Biotechnology, Shanghai, China). The PCR procedure was as follows: denaturation for 5 min at 95 °C, followed by 40 cycles at 95 °C for 10 s and 60 °C for 30 s. The primers used in qRT-PCR were designed using Primer Premier v6.0 (Table S2). The reference gene Tubulin was used for normalization. The relative expression levels of the target genes were calculated using the 2−∆∆Ct method.

3. Results

3.1. Transcriptomic Profile of Tobacco Plants Infested with Whiteflies

To evaluate the global transcriptomic profile of tobacco in response to whitefly infestation, we performed deep RNA-Seq of the tobacco NC89 cultivars following whitefly infestation for 12 and 24 h (hereafter referred to as I-12 h and I-24 h). The non-infested tobacco plants were used as controls (hereafter referred to as NI), and three biological replicates were conducted for each treatment. In total, we generated approximately 2.3 hundred million paired-end (PE) reads from these nine libraries, with approximately 23–28 million reads generated per library. The raw reads underwent sequence duplication and GC content analysis using FastQC software (Table 1). The summaries of the total clean/raw reads, the total number of bases in the clean/raw data, and the error rate of average base sequencing, Q20, and Q30 are presented in Table 1. Additionally, principal component analysis (PCA) revealed significant changes in the tobacco transcriptomic profile induced by whitefly infestation. PC1 and PC2 explained 91% and 4% variance, respectively, and the samples from the three biological replicates exhibited good reproducibility (Figure S1).

3.2. Differently Expressed Genes (DEGs) between Whitefly-Infested and Non-Infested Tobacco

To evaluate alterations in the gene expression in tobacco plants in response to infestation by B. tabaci, we conducted a comprehensive analysis of DEGs. We performed a comparative analysis of RNA libraries from non-infested tobacco plants and those infested with whiteflies for 12 and 24 h, resulting in two distinct categories: (1) early response, a comparison between the control group and tobacco plants infested for 12 h, and (2) late response, a comparison between the control group and tobacco plants infested for 24 h.
For the early response, a total of 244 genes were differentially expressed, including 128 upregulated genes and 116 downregulated genes (Table S3, Figure 1a). For the late response, there were 2417 DEGs, with 983 upregulated and 1434 downregulated (Table S4, Figure 1a). To gain insights into the functions of these differentially expressed genes, we performed gene ontology (GO) enrichment analysis for both types of DEGs. We found that the upregulated genes in both the early and late stages were primarily associated with plant defense pathways, such as response to stimuli, response to biotic stimuli, response to hormones, response to reactive oxygen species, and innate immune responses (Figure 2a,b). These results suggest that whitefly infestation activates the plant innate immune system, particularly the pathways involving reactive oxygen species (ROS) and hormones, which are crucial components of the plant’s sophisticated signaling networks. In contrast, the downregulated genes in both early and late stages were mainly associated with cellular components of plant organelles, such as cytoplasm, chloroplast, and thylakoid (Figure 2c,d). Several genes related to photosynthesis were also downregulated. Previous studies have highlighted the critical role of plant organelles, such as chloroplasts and cytoplasm, in plant immunity [24,25]. The downregulated expression of organelle-related genes suggests that B. tabaci potentially inhibits organelle-mediated immunity to facilitate feeding. There were 89 DEGs specifically differentially expressed at 12 h, but not 24 h, when compared with the non-infested control. Notably, the majority of these genes (68, 76.4%) were specifically upregulated at the early stage, including ABC transporters, serine threonine tyrosine kinase, calmodulin-binding protein, BHLH transcription factor, and UDP-glucosyltransferase. These genes potentially act as immediate-early expressed genes [26] that efficiently respond to whitefly infestation.

3.3. Plant Signal Networks Response to Whitefly Attacks

To gain a better understanding of the plant immune response triggered by B. tabaci infestation and the key resistance genes involved in this process, we conducted a comparative analysis of DEGs between the two post-whitefly infestation time points (12 and 24 h). We found that a total of 54 genes were upregulated and 91 genes were downregulated (Figure 1b). Among the upregulated genes, the WRKY TF1 exhibited the highest upregulation among all upregulated genes. Additionally, other WRKY TFs were upregulated in both early and late stages. Specifically, in the late stage, WRKY TFs, including WRKY1, WRKY2, WRKY3, WRKY5, WRKY6, WRKY35, and WRKY78, were upregulated (Figure 3b). These findings indicate that WRKY genes play important roles in tobacco resistance against whitefly attacks. Interestingly, we also observed the upregulation of five ERFs in response to whitefly infestation, including ERF1, ERF1a, ERF2b, ERF4, and ERF5 (Figure 3c).
In addition to TFs, a significant number of genes associated with PRRs exhibited notable changes. PRRs consist of RLKs and RLPs and are involved in specifically recognizing herbivory-associated cues, such as oral secretions and eggs, thus activating plant immune responses [27,28]. We found that 49 PRR genes were altered, including 26 upregulated and 23 downregulated genes, upon whitefly infestation (Figure 3a). Overall, these findings suggest the activation of TFs and cell surface receptor protein kinases involved in regulating plant immunity in response to whitefly infestation.

3.4. Metabolic Processes Affected by Whitefly Attack

To investigate how B. tabaci regulates plant defenses, we examined DEGs in relation to metabolic pathways, particularly focusing on downregulated DEGs. The results of the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that the downregulated DEGs were primarily associated with metabolic pathways. The commonly downregulated DEGs at both early and late stages were enriched in pathways such as metabolism, energy metabolism, ribosomes, and carbohydrate metabolism (Figure 4a). Furthermore, whitefly infestation resulted in enrichment in other metabolism-related pathways, particularly at a late stage (Figure 4b,c). These results indicate that feeding by B. tabaci can impact fundamental metabolic and energy processes in plants. Metabolic pathways are essential for plant growth, development, and defense against biotic and abiotic stresses. The downregulation of genes associated with metabolism suggests that B. tabaci may inhibit basal metabolism and defense-related metabolic compounds as a strategy to overcome plant defenses.

3.5. Validation of DEGs Using Quantitative Real-Time PCR (qRT-PCR)

To confirm the expression of the DEGs, 10 genes that were associated with plant defenses or showed dramatic upregulation/downregulation in the RNA-Seq analysis were selected for qRT-PCR analysis (Table S1). The validation was conducted using three independent biological replicates. We observed the upregulation of zinc finger proteins and ATP-binding cassette (ABC) transporters. In contrast, two MYB TFs were downregulated. MYB TFs are a large protein family involved in regulating plant growth, development, and stress adaptation. Furthermore, peroxidase 40, L-ascorbate peroxidase and cryptochrome 1b, which are related to plant defense, were downregulated. Additionally, rubisco activate 1 was upregulated, while two late embryogenesis abundant protein 5s were downregulated. Overall, the expression profile of the genes selected for qRT-PCR validation was consistent with the DEG data (Figure 5), which indicates that the RNA-Seq data are reliable and reproducible.

4. Discussion

Whiteflies are among the most economically important pests of crops worldwide, causing significant losses in agricultural productivity each year. Despite previous studies on the whitefly–plant interaction, the ongoing evolutionary arms race between plants and whiteflies necessitates a deeper understanding of their interactions to develop novel strategies for controlling whitefly infestations. Next-generation sequencing technologies have revolutionized the study of insect–plant interactions. In this study, we employed transcriptome sequencing to investigate the gene expression changes in tobacco plants subjected to different treatments, specifically infestation by the whitefly B. tabaci. To gain insight into tobacco gene expression, we constructed cDNA libraries at different time points following whitefly infestation and analyzed the DEGs. Our statistical analysis revealed that the majority of DEGs were concentrated at 24 h after tobacco infestation by B. tabaci. Interestingly, we observed a higher number of downregulated genes compared to upregulated genes at this time point, indicating that whitefly infestation predominantly suppresses gene expression in tobacco.
Previous research has demonstrated that plants respond to pathogen exposure by PRRs, leading to downstream defense responses such as MAPK cascades, calcium flux, ROS bursts, transcriptional reprogramming, and phytohormone signaling pathways [29,30,31,32,33]. In our study, we found that the genes upregulated in response to whitefly infestation were associated with ROS and hormonal pathways. Additionally, GO enrichment pathway analysis highlighted their involvement in the response to stimulus, response to biotic stimulus, and innate immune responses. These findings suggest that the infestation of tobacco by B. tabaci triggers the activation of the plant innate immune system and responses to external stimuli. Notably, our study revealed significant changes in nearly 50 DEGs related to PRRs following whitefly infestation. Plants utilize PRRs to rapidly detect potential dangers posed by microbes and pests [11]. Based on this observation, we hypothesize that whitefly infestation in tobacco activates plant PRRs, thereby initiating downstream defense responses to counter the infestation.
TFs are key components of plant signaling networks that regulate gene expression in response to various factors, including plant growth, development, and biotic and abiotic stresses [34,35]. In our study, several TFs, including WRKY TFs and ERFs, were identified among the DEGs. Previous research has highlighted the crucial role of WRKY TFs in plant defenses against sap-sucking insects like aphids, planthoppers, and whitefly [36]. For example, WRKY89 enhanced resistance to the white-backed planthopper Sogatella furcifera in rice, and altered expression levels of six WRKY genes were observed in cotton following whitefly infestation [37,38]. Additionally, the upregulation of N. tabacum WRKY6 was reported when tobacco plants were infested with whiteflies [39]. In our study, with the exception of the downregulated WRKY 26 gene at 24 h, all other WRKY TFs, including WRKY6, were significantly upregulated. ERFs, another group of plant-specific TFs, are involved in regulating plant development and stress responses. Notably, ERF1 has been shown to play a key role in activating plant defenses against necrotrophic pathogens, while OsERF3 in rice regulates resistance to chewing or piercing/sucking insects through metabolic regulation [40,41,42]. Additionally, ERF5 has been implicated in plant innate immunity [43]. In our study, we observed a significant upregulation of ERF1, ERF1a, ERF2b, ERF4, and ERF5. Collectively, these results indicate that TFs, particularly WRKY TFs and ERFs, play a crucial role in regulating gene expression related to plant defense against whitefly infestation.

5. Conclusions

In conclusion, our study sheds light on the complex interplay between whiteflies and host plants, providing valuable insights into the gene expression changes in tobacco following B. tabaci infestation. The observed activation of the plant immune system, the involvement of PRRs, and the significant upregulation of TFs, especially WRKY TFs and ERFs, underscore their essential roles in plant defense against whitefly infestations. Understanding these molecular mechanisms will contribute to the development of effective strategies for whitefly control and crop protection. Further research should focus on unraveling the intricate details of these interactions and exploring the functional significance of specific genes involved in the plant defense response to whitefly infestations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14081640/s1. Figure S1: Principal component analysis (PCA) of gene expression patterns in tobacco plants with different treatments. The first two principal components (PC1 and PC2) based on transcriptomic results were shown; Table S1: Ten genes utilized for qRT-PCR validation; Table S2: Primers used in this study; Table S3: Differentially expressed genes (DEGs) at 12 h post B. tabaci infestation; Table S4: Differentially expressed genes (DEGs) at 24 h post B. tabaci infestation.

Author Contributions

H.-J.H. and J.-P.C. planned and designed the research. X.W., Z.-X.Y., Y.-Z.W. and X.-J.W. performed experiments and analyzed data. H.-J.H. and X.W. drafted the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 31801734.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequencing data generated in this study were submitted to the National Genomics Data Center under accession number PRJCA018862.

Acknowledgments

We would like to thank Novogene Company (Tianjin, China) for their technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Furstenberg-Hagg, J.; Zagrobelny, M.; Bak, S. Plant defense against insect herbivores. Int. J. Mol. Sci. 2013, 14, 10242–10297. [Google Scholar] [CrossRef] [PubMed]
  2. War, A.R.; Paulraj, M.G.; Ahmad, T.; Buhroo, A.A.; Hussain, B.; Ignacimuthu, S.; Sharma, H.C. Mechanisms of plant defense against insect herbivores. Plant Signal. Behav. 2012, 7, 1306–1320. [Google Scholar] [CrossRef] [PubMed]
  3. Chen, M.S. Inducible direct plant defense against insect herbivores: A review. Insect Sci. 2008, 15, 101–114. [Google Scholar] [CrossRef]
  4. War, A.R.; Sharma, H.C.; Paulraj, M.G.; War, M.Y.; Ignacimuthu, S. Herbivore induced plant volatiles: Their role in plant defense for pest management. Plant Signal. Behav. 2011, 6, 1973–1978. [Google Scholar] [CrossRef]
  5. Arimura, G.; Ozawa, R.; Maffei, M.E. Recent advances in plant early signaling in response to herbivory. Int. J. Mol. Sci. 2011, 12, 3723–3739. [Google Scholar] [CrossRef]
  6. Wu, J.Q.; Baldwin, I.T. New insights into plant responses to the attack from insect herbivores. Annu. Rev. Genet. 2010, 44, 1–24. [Google Scholar] [CrossRef]
  7. Verma, V.; Ravindran, P.; Kumar, P.P. Plant hormone-mediated regulation of stress responses. BMC Plant Biol. 2016, 16, 86. [Google Scholar] [CrossRef]
  8. Zipfel, C.; Felix, G. Plants and animals: A different taste for microbes? Curr. Opin. Plant Biol. 2005, 8, 353–360. [Google Scholar] [CrossRef]
  9. Nurnberger, T.; Brunner, F.; Kemmerling, B.; Piater, L. Innate immunity in plants and animals: Striking similarities and obvious differences. Immunol. Rev. 2004, 198, 249–266. [Google Scholar] [CrossRef]
  10. de Wit, P.J. How plants recognize pathogens and defend themselves. Cell. Mol. Life Sci. 2007, 64, 2726–2732. [Google Scholar] [CrossRef]
  11. Zipfel, C. Plant pattern-recognition receptors. Trends Immunol. 2014, 35, 345–351. [Google Scholar] [CrossRef] [PubMed]
  12. He, Y.; Zhou, J.; Shan, L.; Meng, X. Plant cell surface receptor-mediated signaling—A common theme amid diversity. J. Cell Sci. 2018, 131, jcs209353. [Google Scholar] [CrossRef] [PubMed]
  13. Dombrecht, B.; Xue, G.P.; Sprague, S.J.; Kirkegaard, J.A.; Ross, J.J.; Reid, J.B.; Fitt, G.P.; Sewelam, N.; Schenk, P.M.; Manners, J.M.; et al. MYC2 differentially modulates diverse jasmonate-dependent functions in Arabidopsis. Plant Cell 2007, 19, 2225–2245. [Google Scholar] [CrossRef] [PubMed]
  14. Hu, L.; Ye, M.; Li, R.; Lou, Y. OsWRKY53, a versatile switch in regulating herbivore-induced defense responses in rice. Plant Signal. Behav. 2016, 11, e1169357. [Google Scholar] [CrossRef]
  15. Oliveira, M.; Henneberry, T.; Anderson, P. History, current status, and collaborative research projects for Bemisia tabaci. Crop Prot. 2001, 20, 709–723. [Google Scholar] [CrossRef]
  16. Navas-Castillo, J.; Fiallo-Olive, E.; Sanchez-Campos, S. Emerging virus diseases transmitted by whiteflies. Annu. Rev. Phytopathol. 2011, 49, 219–248. [Google Scholar] [CrossRef]
  17. Cock, M.J.W. Bemisia tabaci—A Literature Survey on the Cotton Whitefly with an Annotated Bibliography; C.A.B International Institute of Biological Control: Ascot, UK, 1986. [Google Scholar]
  18. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  19. Andrews, S. FastQC a Quality Control Tool for High Throughput Sequence Data; Babraham Institute: Babraham, UK, 2014. [Google Scholar]
  20. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef]
  21. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics 2009, 25, 2078–2079. [Google Scholar] [CrossRef]
  22. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef]
  23. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  24. Nomura, H.; Komori, T.; Uemura, S.; Kanda, Y.; Shimotani, K.; Nakai, K.; Furuichi, T.; Takebayashi, K.; Sugimoto, T.; Sano, S.; et al. Chloroplast-mediated activation of plant immune signaling in Arabidopsis. Nat. Commun. 2012, 3, 926. [Google Scholar] [CrossRef] [PubMed]
  25. Nomura, H.; Shiina, T. Calcium signaling in plant endosymbiotic organelles: Mechanism and role in physiology. Mol. Plant 2014, 7, 1094–1104. [Google Scholar] [CrossRef] [PubMed]
  26. Muthamilarasan, M.; Prasad, M. Plant innate immunity: An updated insight into defense mechanism. J. Biosci. 2013, 38, 433–449. [Google Scholar] [CrossRef]
  27. Tang, D.Z.; Wang, G.X.; Zhou, J.M. Receptor kinases in plant-pathogen interactions: More than pattern recognition. Plant Cell 2017, 29, 618–637. [Google Scholar] [CrossRef]
  28. Reymond, P. Receptor kinases in plant responses to herbivory. Curr. Opin. Biotechnol. 2021, 70, 143–150. [Google Scholar] [CrossRef] [PubMed]
  29. Yuan, M.H.; Ngou, B.P.M.; Ding, P.T.; Xiu-Fan, X. PTI-ETI crosstalk: An integrative view of plant immunity. Curr. Opin. Plant Biol. 2021, 62, 102030. [Google Scholar] [CrossRef]
  30. Devendrakumar, K.T.; Li, X.; Zhang, Y.L. MAP kinase signalling: Interplays between plant PAMP- and effector-triggered immunity. Cell. Mol. Life Sci. 2018, 75, 2981–2989. [Google Scholar] [CrossRef]
  31. Rasmussen, S.B.; Reinert, L.S.; Paludan, S.R. Innate recognition of intracellular pathogens: Detection and activation of the first line of defense. Apmis 2009, 117, 323–337. [Google Scholar] [CrossRef]
  32. Peng, Y.J.; van Wersch, R.; Zhang, Y.L. Convergent and divergent signaling in pamp-triggered immunity and effector-triggered immunity. Mol. Plant Microbe Interact. 2018, 31, 403–409. [Google Scholar] [CrossRef]
  33. Tsuda, K.; Katagiri, F. Comparing signaling mechanisms engaged in pattern-triggered and effector-triggered immunity. Curr. Opin. Plant Biol. 2010, 13, 459–465. [Google Scholar] [CrossRef] [PubMed]
  34. Phukan, U.J.; Jeena, G.S.; Shukla, R.K. WRKY transcription factors: Molecular regulation and stress responses in plants. Front. Plant Sci. 2016, 7, 760. [Google Scholar] [CrossRef] [PubMed]
  35. Kundu, P.; Vadassery, J. Role of WRKY transcription factors in plant defense against lepidopteran insect herbivores: An overview. J. Plant Biochem. Biotechnol. 2021, 30, 698–707. [Google Scholar] [CrossRef]
  36. Kamphuis, L.G.; Zulak, K.; Gao, L.L.; Anderson, J.; Singh, K.B. Plant-aphid interactions with a focus on legumes. Funct. Plant Biol. 2013, 40, 1271–1284. [Google Scholar] [CrossRef]
  37. Li, J.; Zhu, L.; Hull, J.J.; Liang, S.; Daniell, H.; Jin, S.; Zhang, X. Transcriptome analysis reveals a comprehensive insect resistance response mechanism in cotton to infestation by the phloem feeding insect Bemisia tabaci (whitefly). Plant Biotechnol. J. 2016, 14, 1956–1975. [Google Scholar] [CrossRef]
  38. Wang, H.; Hao, J.; Chen, X.; Hao, Z.; Wang, X.; Lou, Y.; Peng, Y.; Guo, Z. Overexpression of rice WRKY89 enhances ultraviolet B tolerance and disease resistance in rice plants. Plant Mol. Biol. 2007, 65, 799–815. [Google Scholar] [CrossRef]
  39. Yao, D.M.; Zou, C.; Shu, Y.N.; Liu, S.S. WRKY transcription factors in nicotiana tabacum modulate plant immunity against whitefly via interacting with MAPK cascade pathways. Insects 2020, 12, 16. [Google Scholar] [CrossRef]
  40. Ohta, M.; Matsui, K.; Hiratsu, K.; Shinshi, H.; Ohme-Takagi, M. Repression domains of class II ERF transcriptional repressors share an essential motif for active repression. Plant Cell 2001, 13, 1959–1968. [Google Scholar] [CrossRef]
  41. Lorenzo, O.; Piqueras, R.; Sanchez-Serrano, J.J.; Solano, R. ETHYLENE RESPONSE FACTOR1 integrates signals from ethylene and jasmonate pathways in plant defense. Plant Cell 2003, 15, 165–178. [Google Scholar] [CrossRef]
  42. Lu, J.; Ju, H.P.; Zhou, G.X.; Zhu, C.S.; Erb, M.; Wang, X.P.; Wang, P.; Lou, Y.G. An EAR-motif-containing ERF transcription factor affects herbivore-induced signaling, defense and resistance in rice. Plant J. 2011, 68, 583–596. [Google Scholar] [CrossRef]
  43. Son, G.H.; Wan, J.R.; Kim, H.J.; Nguyen, X.C.; Chung, W.S.; Hong, J.C.; Stacey, G. Ethylene-responsive element-binding factor 5, ERF5, is involved in chitin-induced innate immunity response. Mol. Plant Microbe Interact. 2012, 25, 48–60. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Analysis of differentially expressed genes (DEG) upon B. tabaci infestation. (a) Number of DEGs in tobacco infested by B. tabaci at different time points. (b) The Venn diagram shows the number of specific and common DEGs that are upregulated and downregulated at different time points.
Figure 1. Analysis of differentially expressed genes (DEG) upon B. tabaci infestation. (a) Number of DEGs in tobacco infested by B. tabaci at different time points. (b) The Venn diagram shows the number of specific and common DEGs that are upregulated and downregulated at different time points.
Genes 14 01640 g001
Figure 2. Scatterplot of Gene Ontology (GO) pathway enrichment analysis for differentially expressed genes (DEG). (a) Analysis of upregulated DEGs at 12 h post infestation. (b) Analysis of upregulated DEGs at 24 h post infestation. (c) Analysis of downregulated DEGs at 12 h post infestation. (d) Analysis of downregulated DEGs at 24 h post infestation. The rich factor is the ratio of DEGs annotated in a given pathway term to the total number of genes annotated in that pathway term. The gene ratio is represented on the horizontal axis, while the enriched pathway names are shown on the vertical axis. The color scale indicates various p-value thresholds, and the size of the dots corresponds to the number of genes associated with each pathway.
Figure 2. Scatterplot of Gene Ontology (GO) pathway enrichment analysis for differentially expressed genes (DEG). (a) Analysis of upregulated DEGs at 12 h post infestation. (b) Analysis of upregulated DEGs at 24 h post infestation. (c) Analysis of downregulated DEGs at 12 h post infestation. (d) Analysis of downregulated DEGs at 24 h post infestation. The rich factor is the ratio of DEGs annotated in a given pathway term to the total number of genes annotated in that pathway term. The gene ratio is represented on the horizontal axis, while the enriched pathway names are shown on the vertical axis. The color scale indicates various p-value thresholds, and the size of the dots corresponds to the number of genes associated with each pathway.
Genes 14 01640 g002
Figure 3. Analysis of differentially expressed genes (DEG) associated with plant resistance. The expression patterns of receptor-like kinases (RLKs) and receptor-like proteins (RLPs) (a), WRKY transcription factors (TFs) (b), and ethylene-responsive transcription factors (ERFs) (c) were displayed. The heatmap was drawn using TBtools software, and the function “scale” was used for the data by rows. The scale method was “zero to one”. A red color indicates high expression, while a blue color indicates low expression.
Figure 3. Analysis of differentially expressed genes (DEG) associated with plant resistance. The expression patterns of receptor-like kinases (RLKs) and receptor-like proteins (RLPs) (a), WRKY transcription factors (TFs) (b), and ethylene-responsive transcription factors (ERFs) (c) were displayed. The heatmap was drawn using TBtools software, and the function “scale” was used for the data by rows. The scale method was “zero to one”. A red color indicates high expression, while a blue color indicates low expression.
Genes 14 01640 g003
Figure 4. Enrichment analysis of downregulated, differentially expressed genes (DEG) upon B. tabaci infestation. (a) Analysis of DEGs commonly downregulated at 12 h and 24 h post whitefly infestation. (b,c) Analysis of DEGs downregulated at 12 h (b) and 24 h (c) post whitefly infestation.
Figure 4. Enrichment analysis of downregulated, differentially expressed genes (DEG) upon B. tabaci infestation. (a) Analysis of DEGs commonly downregulated at 12 h and 24 h post whitefly infestation. (b,c) Analysis of DEGs downregulated at 12 h (b) and 24 h (c) post whitefly infestation.
Genes 14 01640 g004
Figure 5. Correlation between transcriptomic data (green) and qRT-PCR results (blue). The relative expression level of each gene was determined using qRT-PCR and was compared with the expression of the transcriptomic data. p-values in transcriptomic data were determined with DESeq software, while p-values in qPCR were determined with a two-tailed unpaired Student’s t test. ** p < 0.01; *** p < 0.001; ns, not significant. Data are presented as mean values ± SD. (a) Zinc finger CCCH domain-containing protein; (b) Ribulose-1,5-bisphosphate carboxylase; (c) ABC transporter; (d) MYB transcription factor REVEILLE 8; (e) MYB transcription factor REVEILLE 1; (f) Protein kinase At1g79600; (g) Late embryogenesis abundant protein 5; (h) L-ascorbate peroxidase; (i) Cryptochrome b; (j) Late embryogenesis abundant protein 5.
Figure 5. Correlation between transcriptomic data (green) and qRT-PCR results (blue). The relative expression level of each gene was determined using qRT-PCR and was compared with the expression of the transcriptomic data. p-values in transcriptomic data were determined with DESeq software, while p-values in qPCR were determined with a two-tailed unpaired Student’s t test. ** p < 0.01; *** p < 0.001; ns, not significant. Data are presented as mean values ± SD. (a) Zinc finger CCCH domain-containing protein; (b) Ribulose-1,5-bisphosphate carboxylase; (c) ABC transporter; (d) MYB transcription factor REVEILLE 8; (e) MYB transcription factor REVEILLE 1; (f) Protein kinase At1g79600; (g) Late embryogenesis abundant protein 5; (h) L-ascorbate peroxidase; (i) Cryptochrome b; (j) Late embryogenesis abundant protein 5.
Genes 14 01640 g005
Table 1. Summary of statistics from Illumina sequencing.
Table 1. Summary of statistics from Illumina sequencing.
SampleRaw ReadsClean ReadsRaw Base (G) 1Clean Base (G) 2Error (%) 3Q20 (%) 4Q30 (%) 5GC Content (%)
NI-122825325223898776.856.720.0397.3292.5341.69
NI-226091282244683307.837.340.0397.8993.6943.02
NI-328556897266261768.577.990.0397.8993.7643.18
I-12 h-127446453258175328.237.750.0397.9793.8843.32
I-12 h-226080109243278437.827.30.0397.9593.8943.31
I-12 h-326357914244812417.917.340.0397.9893.9843.36
I-24 h-126781167249715088.037.490.0397.8693.6942.38
I-24 h-223722409233215137.1270.0397.2392.2841.91
I-24 h-324041972223845677.216.720.0397.9193.7941.82
1 Raw Base: the total number of bases in the raw data. 2 Clean Base: the total number of bases in the clean data. 3 Error: error rate of average base sequencing. 4 Q2: percentage of bases with Phred quality scores > 20. 5 Q30: percentage of bases with Phred quality scores > 30.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, X.; Ye, Z.-X.; Wang, Y.-Z.; Wang, X.-J.; Chen, J.-P.; Huang, H.-J. Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection. Genes 2023, 14, 1640. https://doi.org/10.3390/genes14081640

AMA Style

Wang X, Ye Z-X, Wang Y-Z, Wang X-J, Chen J-P, Huang H-J. Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection. Genes. 2023; 14(8):1640. https://doi.org/10.3390/genes14081640

Chicago/Turabian Style

Wang, Xin, Zhuang-Xin Ye, Yi-Zhe Wang, Xiao-Jing Wang, Jian-Ping Chen, and Hai-Jian Huang. 2023. "Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection" Genes 14, no. 8: 1640. https://doi.org/10.3390/genes14081640

APA Style

Wang, X., Ye, Z. -X., Wang, Y. -Z., Wang, X. -J., Chen, J. -P., & Huang, H. -J. (2023). Transcriptomic Analysis of Tobacco Plants in Response to Whitefly Infection. Genes, 14(8), 1640. https://doi.org/10.3390/genes14081640

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