Isolation of two new genes encoding heat shock protein 70 in Bemisia tabaci and analysis during thermal stress

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Changes in Gene Expression of Whiteflies, Bemisia tabaci MED Feeding on Tomato Plants Infected by One of the Criniviruses, Tomato Chlorosis Virus through Transcriptome Analysis

Tomato chlorosis virus (ToCV), transmitted by the whitefly, Bemisia tabaci (Gennadius; Hemiptera: Aleyrodidae) has been continuously emerging on tomato plants and causing a significant economic loss throughout China. In the current study, RNA-Seq analysis was used to explore the gene expression profiles of B. tabaci Mediterranean (MED) that fed on both ToCV-infected and -uninfected tomato plants for 6, 12, 24, and 48 hours, respectively. The results revealed that dynamic changes occurred in the gene expressions of whiteflies at different time intervals after they acquired the virus. A total of 1709, 461, 4548, and 1748 differentially expressed genes (DEGs) were identified after a 6, 12, 24, and 48 hours feeding interval for the viral acquisition, respectively. The least number of expressed genes appeared in whiteflies with the 12 hours feeding treatment, and the largest numbers of those found in those with 24 hours feeding treatment. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that B. tabaci MED responded to ToCV acquisition through altering its nerve system development, fertility, detoxification, glucose metabolism, and immune function before it lost its ability to transmit the virus. The number of DEGs, degree of differential gene expressions, expression level of the same gene, involved biological processes, and metabolic functions in whiteflies post the 12 hours feeding, and viral acquisition were different from those from other three feeding treatments, which could be a significant finding suggesting an effective control of B. tabaci MED should be done less than 12 hours after whiteflies started feeding on ToCV-infected tomatoes. Our results further provided a clarified understanding in how B. tabaci was protected from viral acquisitions through comparison of the differential profile of gene expressions in whiteflies feeding on plants that were infected by semipersistent viruses.

Comparative transcriptome analysis of Callosobruchus chinensis (L.) (Coleoptera: Chrysomelidae-Bruchinae) after heat and cold stress exposure

Callosobruchus chinensis is regarded as one of the cosmopolitan pests of legume crops and can cause tremendous losses to a variety of beans. This study focused on comparative transcriptome analyses of C. chinensis exposed to 45 °C (heat stress), 27 °C (ambient temperature) and -3 °C (cold stress) for 3 h to investigate the gene differences and underlying molecular mechanisms. There were 402 and 111 differentially expressed genes (DEGs) identified in the heat and cold stress treatments, respectively. "cell process", "cell" and "binding" were the main enriched functions and biological processes revealed by gene ontology (GO) analysis. The clusters of orthologous genes (COG) showed that DEGs were assigned to the categories: "posttranslational modification, protein turnover, chaperones", "lipid transport and metabolism", and "general function prediction only". With respect to the Kyoto Encyclopedia of Genes and Genomes (KEGG), the "longevity regulating pathway-multiple species", "carbon metabolism", "peroxisome", "protein processing in endoplasmic", "glyoxylate and dicarboxylate metabolism" pathways were significantly enriched. The annotation and enrichment analysis revealed that genes encoding heat shock proteins (Hsps) and cuticular proteins were significantly upregulated under high and low-temperature stresses, respectively. In addition, some DEGs encoding "Protein lethal essential for life", "Reverse transcriptase", "DnaJ domain", "Cytochrome" and "Zinc finger protein" were also upregulated to varying degrees. Transcriptomic data were validated using qRT‒PCR, which confirmed that they were consistent. In this paper, the temperature tolerance of C. chinensis adults was evaluated and the results showed that female adults were more sensitive to heat and cold stress than males, and the upregulation of heat shock protein and epidermal protein was the largest in DEGs after heat and cold stress, respectively. These findings provide a reference for further understanding the biological characteristics of C. chinensis adults and the molecular mechanisms underlying the response to low and high temperatures.

Molecular characterization of six heat shock protein 70 genes from Arma chinensis and their expression patterns in response to temperature stress

Arma chinensis is an important predatory enemy of many agricultural and forest pests. Heat shock protein 70 (Hsp70) plays an essential role in insect adaptation to various stress factors. To explore the functions of Hsp70s in relation to thermal tolerance of A. chinensis, full-length cDNAs of six Hsp70 genes (AcHsp70Ba, AcHsp70-4, AcHsp68a, AcHsp68b, AcHsp70-2, and AcHsc70-4) were cloned. Their open reading frames (ORFs) were 1902, 2454, 1884, 1905, 1872, and 1947 bp, respectively. Developmental expression profiles showed that AcHsp70Ba, AcHsp70-4, and AcHsc70-4 were extremely highly expressed in adult stages. AcHsp68a and AcHsp70-2 showed the highest level of expression in nymph stages, and AcHsp68b was mainly expressed in male adults. Tissue distribution analysis demonstrated that the AcHsp70s were ubiquitously expressed but showing gene-specific and sex-driven patterns of expression. High temperature induced the expression of the six AcHsp70s. Among them, AcHsp70Ba, AcHsp70-4, AcHsp68a, and AcHsc70-4 were significantly induced at 38 °C for 6 h, while all six AcHsp70s were significantly induced at 38 °C for 24 h. There were differences in responses of the six AcHsp70s to low-temperature stress. The expressions of AcHsp70-4, AcHsp68a, and AcHsp68b in male adults were significantly repressed at 4 °C for 6 h, whereas AcHsp70Ba and AcHsp70-2 were significantly induced. The levels of AcHsp70Ba, AcHsp68b, and AcHsp70-2 in female adults were significantly repressed at 4 °C for 24 h, whereas AcHsc70-4 was significantly induced. These results suggested that AcHsp70s play important roles in various developmental stages and tissue function, and contribute to the tolerance of A. chinensis to extreme temperatures.

Transcriptome analysis and screening of putative sex-determining genes in the invasive pest, Frankliniella occidentalis (Thysanoptera: Thripidae)

The invasive insect pest, Frankliniella occidentalis, is a well-known vector that transmits a variety of ornamental and vegetable viruses. The mechanistic basis of sex determination in F. occidentalis is not well understood, and this hinders our ability to deploy sterile insect technology as an integrated pest management strategy. In this study, six cDNA libraries from female and male adults of F. occidentalis (three biological replicates each) were constructed and transcriptomes were sequenced. A total of 6000 differentially-expressed genes were identified in the two sexes including 2355 up- and 3645 down-regulated genes. A total of 149 sex-related genes were identified based on GO enrichment data and included transformer-2 (tra2), fruitless (fru), male-specific lethal (msl) and sex lethal (sxl); several of these exhibited sex-specific and/or sex-biased expression in F. occidentalis. This study contributes to our understanding of the sex-determined cascade in F. occidentalis and other members of the Thysanoptera.

Induced heat shock protein 70 confers biological tolerance in UV-B stress-adapted Myzus persicae (Hemiptera)

As an environmental stress factor, ultraviolet-B (UV-B) radiation directly affects insect growth, development, and reproduction. Heat shock protein 70s kDa (Hsp70s) plays an important role in the environmental adaptation of insects. To determine the role of MpHsp70s in the UV-B tolerance of Myzus persicae (Sulzer), we identified the complete complementary DNA sequences of seven MpHsp70s. They were found to be ubiquitously expressed during different developmental stages and were highly expressed in second-instar nymphs and wingless adults. The expression levels of the MpHsp70s were significantly upregulated when exposed to different durations of UV-B stress. Nanocarrier-mediated dsMpHsp70 suppressed the expression of the MpHsp70s and reduced the body length, weight, survival rate, and fecundity of M. persicae under UV-B exposure. When the combinational RNAi approach was adopted, the effects on the survival rate and fecundity were greater under UV-B stress, except for MpHsc70-4. These results suggest that MpHsp70s are essential for the resistance of M. persicae to UV-B stress.

Genome-Wide Identification and Analysis of the Heat-Shock Protein Gene Superfamily in Bemisia tabaci and Expression Pattern Analysis under Heat Shock

Simple Summary

Bemisia tabaci MED is an invasive pest that had caused considerable economic damage in the past decades. Its successful colonization is closely related to heat-shock proteins (HSPs), which are related to heat resistance. In this study, 33 BtaHsps were identified based on the sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 BtaHsps were newly identified. Analysis of the secondary structure and evolutionary relationship showed that they were all closely related. In addition, BtaHsp90A3 of the HSP90 family was screened by analyzing the expression level changes of these genes under 42 °C heat shock and RNAi was performed on the BtaHsp90A3. The results showed that the silencing of BtaHsp90A3 is closely related to the heat resistance of B. tabaci MED. Taken together, this study conducted an in-depth identification of BtaHsps that clarifies their evolutionary relationships and their response to thermal stress in B. tabaci MED.

Abstract

The thermal tolerance of Bemisia tabaci MED, an invasive whitefly species with worldwide distribution, plays an important role in its ecological adaptation during the invasion process. Heat-shock proteins (HSPs) are closely related to heat resistance. In this study, 33 Hsps (BtaHsps) were identified based on sequenced genome of B. tabaci MED belonging to six HSP families, among which 22 Hsps were newly identified. The secondary structures of a further 22 BtaHsps were also predicted. The results of RT-qPCR showed that heat shock could affect the expression of 14 of the 22 Hsps newly identified in this study. Among them, the expression level of six Hsps increased under 42 °C treatment. As the unstudied gene, BtaHsp90A3 had the highest increase rate. Therefore, BtaHsp90A3 was chosen for the RNAi test, and silencing BtaHsp90A3 by RNAi decreased the survival rate of adult B. tabaci at 42 °C. The results indicated that only a few Hsps were involved in the thermal tolerance of host whitefly although many Hsps would response under heat stress. This study conducted a more in-depth and comprehensive identification that demonstrates the evolutionary relationship of BtaHsps and illustrates the response of BtaHsps under the influence of thermal stress in B. tabaci MED.