Analysis of a native whitefly transcriptome and its sequence divergence with two invasive whitefly species

Comparative transcriptome analysis of Liriomyza trifolii (Burgess) and Liriomyza sativae (Blanchard) (Diptera: Agromyzidae) in response to rapid cold hardening

The ability of insets to react efficiently to fluctuation in temperature is crucial for them to survive in variable surroundings. Rapid cold hardening (RCH) is a process that increase cold tolerance in most insect species. The molecular mechanisms of RCH remain largely unknown, and whether it is associated with transcriptional changes is unclear. In this study, we compared the transcriptomes of Liriomyza trifolii and L. sativae exposed to RCH to investigate the transcript abundance due to RCH in both species. RNA-seq revealed 93,166 assembled unigenes, and 34,303 of these were annotated in the L. trifolii and L. sativae transcriptome libraries. After a 4-h treatment at 1°C (RCH) compared with control, 268 and 606 unigenes were differentially expressed in L. trifolii and L. sativae, respectively. When comparing pupae exposed to 2h cold shock directly with pupae went through 4h acclimation prior to 2h cold shock, 60 and 399 unigenes were differentially expressed in L trifolii and L sativae, respectively. Genes that were commonly expressed in both L. trifolii and L. sativae, included cytochrome P450, cuticular protein, glucose dehydrogenase, solute carrier family 22 and cationic amino acid transporter. Additionally, several pathways including galactose metabolism and peroxisome were significantly enriched during RCH. Our results show that the transcriptional response is correlated with RCH in the pupal stage of the two Liriomyza species, but more transcriptional changes were identified in L sativae than in L. trifolii.

Whole-genome single nucleotide polymorphism and mating compatibility studies reveal the presence of distinct species in sub-Saharan Africa Bemisia tabaci whiteflies

In sub-Saharan Africa cassava growing areas, two members of the Bemisia tabaci species complex termed sub-Saharan Africa 1 (SSA1) and SSA2 have been reported as the prevalent whiteflies associated with the spread of viruses that cause cassava mosaic disease (CMD) and cassava brown streak disease (CBSD) pandemics. At the peak of CMD pandemic in the late 1990s, SSA2 was the prevalent whitefly, although its numbers have diminished over the last two decades with the resurgence of SSA1 whiteflies. Three SSA1 subgroups (SG1 to SG3) are the predominant whiteflies in East Africa and vary in distribution and biological properties. Mating compatibility between SSA1 subgroups and SSA2 whiteflies was reported as the possible driver for the resurgence of SSA1 whiteflies. In this study, a combination of both phylogenomic methods and reciprocal crossing experiments were applied to determine species status of SSA1 subgroups and SSA2 whitefly populations. Phylogenomic analyses conducted with 26 548 205 bp whole genome single nucleotide polymorphisms (SNPs) and the full mitogenomes clustered SSA1 subgroups together and separate from SSA2 species. Mating incompatibility between SSA1 subgroups and SSA2 further demonstrated their distinctiveness from each other. Phylogenomic analyses conducted with SNPs and mitogenomes also revealed different genetic relationships among SSA1 subgroups. The former clustered SSA1-SG1 and SSA1-SG2 together but separate from SSA1-SG3, while the latter clustered SSA1-SG2 and SSA1-SG3 together but separate from SSA1-SG1. Mating compatibility was observed between SSA1-SG1 and SSA1-SG2, while incompatibility occurred between SSA1-SG1 and SSA1-SG3, and SSA1-SG2 and SSA1-SG3. Mating results among SSA1 subgroups were coherent with phylogenomics results based on SNPs but not the full mitogenomes. Furthermore, this study revealed that the secondary endosymbiont-Wolbachia-did not mediate reproductive success in the crossing assays carried out. Overall, using genome wide SNPs together with reciprocal crossings assays, this study established accurate genetic relationships among cassava-colonizing populations, illustrating that SSA1 and SSA2 are distinct species while at least two species occur within SSA1 species.

Lysine provisioning by horizontally acquired genes promotes mutual dependence between whitefly and two intracellular symbionts

Horizontal gene transfer is widespread in insects bearing intracellular symbionts. Horizontally transferred genes (HTGs) are presumably involved in amino acid synthesis in sternorrhynchan insects. However, their role in insect-symbiont interactions remains largely unknown. We found symbionts Portiera, Hamiltonella and Rickettsia possess most genes involved in lysine synthesis in the whitefly Bemisia tabaci MEAM1 although their genomes are reduced. Hamiltonella maintains a nearly complete lysine synthesis pathway. In contrast, Portiera and Rickettsia require the complementation of whitefly HTGs for lysine synthesis and have lysE, encoding a lysine exporter. Furthermore, each horizontally transferred lysine gene of ten B. tabaci cryptic species shares an evolutionary origin. We demonstrated that Hamiltonella did not alter the titers of Portiera and Rickettsia or lysine gene expression of Portiera, Rickettsia and whiteflies. Hamiltonella also did not impact on lysine levels or protein localization in bacteriocytes harboring Portiera and ovaries infected with Rickettsia. Complementation with whitefly lysine synthesis HTGs rescued E. coli lysine gene knockout mutants. Silencing whitefly lysA in whiteflies harboring Hamiltonella reduced lysine levels, adult fecundity and titers of Portiera and Rickettsia without influencing the expression of Hamiltonella lysA. Furthermore, silencing whitefly lysA in whiteflies lacking Hamiltonella reduced lysine levels, adult fecundity and titers of Portiera and Rickettsia in ovarioles. Therefore, we, for the first time, demonstrated an essential amino acid lysine synthesized through HTGs is important for whitefly reproduction and fitness of both obligate and facultative symbionts, and it illustrates the mutual dependence between whitefly and its two symbionts. Collectively, this study reveals that acquisition of horizontally transferred lysine genes contributes to coadaptation and coevolution between B. tabaci and its symbionts.

Molecular Characterization of Donacia provosti (Coleoptera: Chrysomelidae) Larval Transcriptome by De Novo Assembly to Discover Genes Associated with Underwater Environmental Adaptations

Simple Summary

Donacia provosti is one of the major pests of aquatic crops. It has been widely distributed in the world causing extensive damage to lotus and rice plants. The larvae generally live-in water; however, little is known about the evolution and molecular mechanisms underlying the adaptation. Here, we generated the first larval transcriptome of D. provosti in order to identify potential genetic mechanisms of aquatic adaptation. About 5036 orthologous clusters were identified among four species and 494 unique clusters were identified from D. provosti larvae including the visual perception. Moreover, 93 orthologous gene pairs were found evolving under positive selection. Our results also showed that 4 gene pairs out of the 93 gene pairs were associated with the “mTOR signaling pathway”, which are predicted to be involved in the molecular mechanism of D. provosti adaptation to the underwater environment. In the light of the increasing availability of transcriptomic information for beetle underwater habitat and evolutionary analyses, it is expected that this paper will provide us with some novel insights into aquatic adaptation in beetles and serves as a foundation for future studies aiming to identify candidate genes underlying the genetic basis of aquatic adaptation in beetles.

Abstract

Donacia provosti (Fairmaire, 1885) is a major pest of aquatic crops. It has been widely distributed in the world causing extensive damage to lotus and rice plants. Changes in gene regulation may play an important role in adaptive evolution, particularly during adaptation to feeding and living habits. However, little is known about the evolution and molecular mechanisms underlying the adaptation of D. provosti to its lifestyle and living habits. To address this question, we generated the first larval transcriptome of D. provosti. A total of 20,692 unigenes were annotated from the seven public databases and around 18,536 protein-coding genes have been predicted from the analysis of D. provosti transcriptome. About 5036 orthologous cutlers were identified among four species and 494 unique clusters were identified from D. provosti larvae including the visual perception. Furthermore, to reveal the molecular difference between D. provosti and the Colorado potato beetle Leptinotarsa decemlineata, a comparison between CDS of the two beetles was conducted and 6627 orthologous gene pairs were identified. Based on the ratio of nonsynonymous and synonymous substitutions, 93 orthologous gene pairs were found evolving under positive selection. Interestingly, our results also show that there are 4 orthologous gene pairs of the 93 gene pairs were associated with the “mTOR signaling pathway”, which are predicted to be involved in the molecular mechanism of D. provosti adaptation to the underwater environment. This study will provide us with an important scientific basis for building effective prevention and control system of the aquatic leaf beetle Donacia provosti.

Pantothenate mediates the coordination of whitefly and symbiont fitness

Intracellular symbionts in insects often have reduced genomes. Host acquisition of genes from bacteria is an important adaptation that supports symbionts. However, the function of horizontally transferred genes in insect symbiosis remains largely unclear. The primary symbiont Portiera housed in bacteriocytes lacks pantothenate synthesis genes: panB and panC, which is presumably complemented by a fused gene panB-panC (hereafter panBC) horizontally transferred from bacteria in Bemisia tabaci MEAM1. We found panBC in many laboratory cultures, and species of B. tabaci shares a common evolutionary origin. We demonstrated that complementation with whitefly panBC rescued E. coli pantothenate gene knockout mutants. Portiera elimination decreased the pantothenate level and PanBC abundance in bacteriocytes, and reduced whitefly survival and fecundity. Silencing PanBC decreased the Portiera titer, reduced the pantothenate level, and decreased whitefly survival and fecundity. Supplementation with pantothenate restored the symbiont titer, PanBC level, and fitness of RNAi whiteflies. These data suggest that pantothenate synthesis requires cooperation and coordination of whitefly PanBC expression and Portiera. This host-symbiont co-regulation was mediated by the pantothenate level. Our findings demonstrated that pantothenate production, by the cooperation of a horizontally acquired, fused bacteria gene and Portiera, facilitates the coordination of whitefly and symbiont fitness. Thus, this study extends our understanding on the basis of complex host-symbiont interactions.

Glucosylation prevents plant defense activation in phloem-feeding insects

The metabolic adaptations by which phloem-feeding insects counteract plant defense compounds are poorly known. Two-component plant defenses, such as glucosinolates, consist of a glucosylated protoxin that is activated by a glycoside hydrolase upon plant damage. Phloem-feeding herbivores are not generally believed to be negatively impacted by two-component defenses due to their slender piercing-sucking mouthparts, which minimize plant damage. However, here we document that glucosinolates are indeed activated during feeding by the whitefly Bemisia tabaci. This phloem feeder was also found to detoxify the majority of the glucosinolates it ingests by the stereoselective addition of glucose moieties, which prevents hydrolytic activation of these defense compounds. Glucosylation of glucosinolates in B. tabaci was accomplished via a transglucosidation mechanism, and two glycoside hydrolase family 13 (GH13) enzymes were shown to catalyze these reactions. This detoxification reaction was also found in a range of other phloem-feeding herbivores.

Comparative transcriptome analysis of three invasive leafminer flies provides insights into interspecific competition

Liriomyza spp. (Diptera: Agromyzidae) represent a group of economically-significant highly polyphagous pests of plants grown in field and greenhouse conditions. Liriomyza spp. share similar biological and morphological characteristics, and complex interspecific interactions have been documented among these species in various geographical regions. Where the displacement of one of these species by the other has been studied, no unique mechanisms have been identified as causing it. The impact of competitive factors (such as, insecticide tolerance, thermotolerance, and adaptability to cropping systems) may be unique to specific geographic regions of Liriomyza spp., but more research is needed to confirm these hypotheses. In this study, RNA-seq was used to determine the transcriptomes of three closely-related leafminers, e.g. L. sativae, L. trifolii, and L. huidobrensis. Over 20 Gb of clean reads were generated and assembled into unique transcriptomes, and 38,747 unigenes were annotated in different databases. In pairwise comparisons, L. trifolii and L. sativae had more up-regulated genes than L. huidobrensis. With respect to common differentially-expressed genes (Co-DEGs), the three leafminers exhibited distinct groups of highly-expressed gene clusters. When genes related to competitive factors were compared, expression patterns in L. trifolii and L. sativae were more closely related to each other than to L. huidobrensis. The data suggest that DEGs involved in competitive factors may play a key role in competition and displacement of leafminers. The divergent genes identified in this study will be valuable in revealing possible mechanisms of invasion, displacement and interspecific competition in Liriomyza spp.

Biotin provisioning by horizontally transferred genes from bacteria confers animal fitness benefits

Insect symbionts are widespread in nature and lateral gene transfer is prevalent in insect symbiosis. However, the function of horizontally transferred genes (HTGs) in insect symbiosis remains speculative, including the mechanism that enables insects to feed on plant phloem deficient in B vitamins. Previously, we found there is redundancy in biotin synthesis pathways from both whitefly Bemisia tabaci and symbiotic Hamiltonella due to the presence of whitefly HTGs. Here, we demonstrate that elimination of Hamiltonella decreased biotin levels but elevated the expression of horizontally transferred biotin genes in whiteflies. HTGs proteins exhibit specific expression patterns in specialized insect cells called bacteriocytes housing symbionts. Complementation with whitefly HTGs rescued E. coli biotin gene knockout mutants. Furthermore, silencing whitefly HTGs in Hamiltonella-infected whiteflies reduced biotin levels and hindered adult survival and fecundity, which was partially rescued by biotin supplementation. Each of horizontally transferred biotin genes are conserved in various laboratory cultures and species of whiteflies with geographically diverse distributions, which shares an evolutionary origin. We provide the first experimental evidence that biotin synthesized through acquired HTGs is important in whiteflies and may be as well in other animals. Our findings suggest that B vitamin provisioning in animal-microbe symbiosis frequently evolved from bacterial symbionts to animal hosts through horizontal gene transfer events. This study will also shed light on how the animal genomes evolve through functional transfer of genes with bacterial origin in the wider contexts of microbial ecology.

Transcriptome analysis and comparison reveal divergence between the Mediterranean and the greenhouse whiteflies

Both the Mediterranean (MED) species of the Bemisia tabaci whitefly complex and the greenhouse whitefly (Trialeurodes vaporariorum, TV) are important agricultural pests. The two species of whiteflies differ in many aspects such as morphology, geographical distribution, host plant range, plant virus transmission, and resistance to insecticides. However, the molecular basis underlying their differences remains largely unknown. In this study, we analyzed the genetic divergences between the transcriptomes of MED and TV. In total, 2,944 pairs of orthologous genes were identified. The average identity of amino acid sequences between the two species is 93.6%. The average nonsynonymous (Ka) and synonymous (Ks) substitution rates and the ratio of Ka/Ks of the orthologous genes are 0.0389, 2.23 and 0.0204, respectively. The low average Ka/Ks ratio indicates that orthologous genes tend to be under strong purified selection. The most divergent gene classes are related to the metabolisms of xenobiotics, cofactors, vitamins and amino acids, and this divergence may underlie the different biological characteristics between the two species of whiteflies. Genes of differential expression between the two species are enriched in carbohydrate metabolism and regulation of autophagy. These findings provide molecular clues to uncover the biological and molecular differences between the two species of whiteflies.

Evolutionary Patterns of Sex-Biased Genes in Three Species of Haplodiploid Insects

Females and males often differ obviously in morphology and behavior, and the differences between sexes are the result of natural selection and/or sexual selection. To a great extent, the differences between the two sexes are the result of differential gene expression. In haplodiploid insects, this phenomenon is obvious, since males develop from unfertilized zygotes and females develop from fertilized zygotes. Whiteflies of the Bemisia tabaci species complex are typical haplodiploid insects, and some species of this complex are important pests of many crops worldwide. Here, we report the transcriptome profiles of males and females in three species of this whitefly complex. Between-species comparisons revealed that non-sex-biased genes display higher variation than male-biased or female-biased genes. Sex-biased genes evolve at a slow rate in protein coding sequences and gene expression and have a pattern of evolution that differs from those of social haplodiploid insects and diploid animals. Genes with high evolutionary rates are more related to non-sex-biased traits-such as nutrition, immune system, and detoxification-than to sex-biased traits, indicating that the evolution of protein coding sequences and gene expression has been mainly driven by non-sex-biased traits.

Examination of Genomic and Transcriptomic Alterations in a Morphologically Stable Line, MU1, Generated by Intergeneric Pollination

Interspecific hybridization creates genetic variation useful for crop improvement. However, whether pollen from a different genus affects the genomic stability and/or transcriptome of the recipient species during intergeneric pollination has not been investigated. Here, we crossed japonica rice cv. Z12 with the maize accession B73 (pollen donor) and obtained a morphologically stable line, MU1, exhibiting moderate dwarfism, higher tiller number, and increased grain weight compared with Z12. To reveal the genetic basis of these morphological changes in MU1, we performed whole-genome resequencing of MU1 and Z12. Compared with Z12, MU1 showed 107,250 single nucleotide polymorphisms (SNPs) and 23,278 insertion/deletions (InDels). Additionally, 5'-upstream regulatory regions (5'UTRs) of 429 and 309 differentially expressed genes (DEGs) in MU1 contained SNPs and InDels, respectively, suggesting that a subset of these DEGs account for the variation in 5'UTRs. Transcriptome analysis revealed 2190 DEGs in MU1 compared with Z12. Genes up-regulated in MU1 were mainly involved in photosynthesis, generation of precursor metabolites, and energy and cellular biosynthetic processes; whereas those down-regulated in MU1 were involved in plant hormone signal transduction pathway and response to stimuli and stress processes. Quantitative PCR (qPCR) further identified the expression levels of the up- or down-regulated gene in plant hormone signal transduction pathway. The expression level changes of plant hormone signal transduction pathway may be significant for plant growth and development. These findings suggest that mutations caused by intergeneric pollination could be the important reason for changes of MU1 in agronomic traits.

Impacts of Salt Stress on Locomotor and Transcriptomic Responses in the Intertidal Gastropod Batillaria attramentaria

Salinity is one of the most crucial environmental factors that structures biogeographic boundaries of aquatic organisms, affecting distribution, abundance, and behavior. However, the association between behavior and gene regulation underlying acclimation to changes in salinity remains poorly understood. In this study, we investigated the effects of salinity stress on behavior (movement distance) and patterns of gene expression (using RNA sequencing) of the intertidal gastropod Batillaria attramentaria. We examined responses to short-term (1-hour) and long-term (30-day) acclimation to a range of salinities (43, 33 [control], 23, 13, and 3 psu). We found that the intertidal B. attramentaria is able to tolerate a broad range of salinity from 13 to 43 psu but not the acute low salinity of 3 psu. Behavioral experiments showed that salt stress significantly influenced snails' movement, with lower salinity resulting in shorter movement distance. Transcriptomic analyses revealed critical metabolic pathways and genes potentially involved in acclimation to salinity stress, including ionic and osmotic regulation, signal and hormonal transduction pathways, water exchange, cell protection, and gene regulation or epigenetic modification. In general, our study presents a robust, integrative laboratory-based approach to investigate the effects of salt stress on a nonmodel gastropod facing detrimental consequences of environmental change. The current genetic results provide a wealth of reference data for further research on mechanisms of ionic and osmotic regulation and adaptive evolution of this coastal gastropod.

Two types of lysozymes from the whitefly Bemisia tabaci: Molecular characterization and functional diversification

Lysozyme is well-known as an immune effector in the immune system. Here we identified three genes including one c-type lysozyme, Btlysc, and two i-type lysozymes, Btlysi1 and Btlysi2, from the whitefly Bemisia tabaci. All three lysozymes were constitutively expressed in different tissues and developmental stages, but the two types of lysozymes showed different expression patterns. The expression levels of Btlysi1 and Btlysi2 were dramatically induced after the whitefly fed with different host plants while the expression level of Btlysc kept unchanged. After fungal infection and begomovirus acquisition, Btlysc expression was significantly upregulated while Btlysi1 and Btlysi2 expression were basically not induced. Furthermore, we found that Btlysc showed muramidase and antibacterial activities. Altogether, our results suggest that the two types of lysozymes act in two different ways in B. tabaci, that is, Btlysc is involved in the whitefly immune system while Btlysi1 and Btlysi2 may play a role in digestion or nutrition absorption.

The first transcriptomes from field-collected individual whiteflies ( Bemisia tabaci, Hemiptera: Aleyrodidae): a case study of the endosymbiont composition

Background: Bemisia tabaci species ( B. tabaci), or whiteflies, are the world's most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations. Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes. Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop. Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

The first transcriptomes from field-collected individual whiteflies ( Bemisia tabaci, Hemiptera: Aleyrodidae): a case study of the endosymbiont composition

Background: Bemisia tabaci species ( B. tabaci), or whiteflies, are the world’s most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations.

Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes.

Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop.

Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

The first transcriptomes from field-collected individual whiteflies (Bemisia tabaci, Hemiptera: Aleyrodidae)

Background: Bemisia tabaci species (B. tabaci), or whiteflies, are the world’s most devastating insect pests. They cause billions of dollars (US) of damage each year, and are leaving farmers in the developing world food insecure. Currently, all publically available transcriptome data for B. tabaci are generated from pooled samples, which can lead to high heterozygosity and skewed representation of the genetic diversity. The ability to extract enough RNA from a single whitefly has remained elusive due to their small size and technological limitations. Methods: In this study, we optimised a single whitefly RNA extraction procedure, and sequenced the transcriptome of four individual adult Sub-Saharan Africa 1 (SSA1) B. tabaci. Transcriptome sequencing resulted in 39-42 million raw reads. De novo assembly of trimmed reads yielded between 65,000-162,000 Contigs across B. tabaci transcriptomes. Results: Bayesian phylogenetic analysis of mitochondrion cytochrome I oxidase (mtCOI) grouped the four whiteflies within the SSA1 clade. BLASTn searches on the four transcriptomes identified five endosymbionts; the primary endosymbiont Portiera aleyrodidarum and four secondary endosymbionts: Arsenophonus, Wolbachia, Rickettsia, and Cardinium spp. that were predominant across all four SSA1 B. tabaci samples with prevalence levels of between 54.1 to 75%. Amino acid alignments of the NusG gene of P. aleyrodidarum for the SSA1 B. tabaci transcriptomes of samples WF2 and WF2b revealed an eleven amino acid residue deletion that was absent in samples WF1 and WF2a. Comparison of the protein structure of the NusG protein from P. aleyrodidarum in SSA1 with known NusG structures showed the deletion resulted in a shorter D loop. Conclusions: The use of field-collected specimens means time and money will be saved in future studies using single whitefly transcriptomes in monitoring vector and viral interactions. Our method is applicable to any small organism where RNA quantity has limited transcriptome studies.

African ancestry of New World, Bemisia tabaci-whitefly species

Bemisia tabaci whitefly species are some of the world’s most devastating agricultural pests and plant-virus disease vectors. Elucidation of the phylogenetic relationships in the group is the basis for understanding their evolution, biogeography, gene-functions and development of novel control technologies. We report here the discovery of five new Sub-Saharan Africa (SSA) B. tabaci putative species, using the partial mitochondrial cytochrome oxidase 1 gene: SSA9, SSA10, SSA11, SSA12 and SSA13. Two of them, SSA10 and SSA11 clustered with the New World species and shared 84.8‒86.5% sequence identities. SSA10 and SSA11 provide new evidence for a close evolutionary link between the Old and New World species. Re-analysis of the evolutionary history of B. tabaci species group indicates that the new African species (SSA10 and SSA11) diverged from the New World clade c. 25 million years ago. The new putative species enable us to: (i) re-evaluate current models of B. tabaci evolution, (ii) recognise increased diversity within this cryptic species group and (iii) re-estimate divergence dates in evolutionary time.

Towards an understanding of the molecular basis of effective RNAi against a global insect pest, the whitefly Bemisia tabaci

In planta RNAi against essential insect genes offers a promising route to control insect crop pests, but is constrained for many insect groups, notably phloem sap-feeding hemipterans, by poor RNAi efficacy. This study conducted on the phloem-feeding whitefly Bemisia tabaci reared on tomato plants investigated the causes of low RNAi efficacy and routes to ameliorate the problem. Experiments using tomato transgenic lines containing ds-GFP (green fluorescent protein) revealed that full-length dsRNA is phloem-mobile, ingested by the insects, and degraded in the insect. We identified B. tabaci homologs of nuclease genes (dsRNases) in other insects that degrade dsRNA, and demonstrated that degradation of ds-GFP in B. tabaci is suppressed by administration of dsRNA against these genes. dsRNA against the nuclease genes was co-administered with dsRNA against two insect genes, an aquaporin AQP1 and sucrase SUC1, that are predicted to protect B. tabaci against osmotic collapse. When dsRNA constructs for AQP1, SUC1, dsRNase1 and dsRNase2 were stacked, insect mortality was significantly elevated to 50% over 6 days on artificial diets. This effect was accompanied by significant reduction in gene expression of the target genes in surviving diet-fed insects. This study offers proof-of-principle that the efficacy of RNAi against insect pests can be enhanced by using dsRNA to suppress the activity of RNAi-suppressing nuclease genes, especially where multiple genes with related physiological function but different molecular function are targeted.

Genetics and Genomics of Cotton Leaf Curl Disease, Its Viral Causal Agents and Whitefly Vector: A Way Forward to Sustain Cotton Fiber Security

Cotton leaf curl disease (CLCuD) after its first epidemic in 1912 in Nigeria, has spread to different cotton growing countries including United States, Pakistan, India, and China. The disease is of viral origin-transmitted by the whitefly Bemisia tabaci, which is difficult to control because of the prevalence of multiple virulent viral strains or related species. The problem is further complicated as the CLCuD causing virus complex has a higher recombination rate. The availability of alternate host crops like tomato, okra, etc., and practicing mixed type farming system have further exaggerated the situation by adding synergy to the evolution of new viral strains and vectors. Efforts to control this disease using host plant resistance remained successful using two gene based-resistance that was broken by the evolution of new resistance breaking strain called Burewala virus. Development of transgenic cotton using both pathogen and non-pathogenic derived approaches are in progress. In future, screening for new forms of host resistance, use of DNA markers for the rapid incorporation of resistance into adapted cultivars overlaid with transgenics and using genome editing by CRISPR/Cas system would be instrumental in adding multiple layers of defense to control the disease-thus cotton fiber production will be sustained.

The Involvement of Heat Shock Proteins in the Establishment of Tomato Yellow Leaf Curl Virus Infection

Tomato yellow leaf curl virus (TYLCV), a begomovirus, induces protein aggregation in infected tomatoes and in its whitefly vector Bemisia tabaci. The interactions between TYLCV and HSP70 and HSP90 in plants and vectors are necessity for virus infection to proceed. In infected host cells, HSP70 and HSP90 are redistributed from a soluble to an aggregated state. These aggregates contain, together with viral DNA/proteins and virions, HSPs and components of the protein quality control system such as ubiquitin, 26S proteasome subunits, and the autophagy protein ATG8. TYLCV CP can form complexes with HSPs in tomato and whitefly. Nonetheless, HSP70 and HSP90 play different roles in the viral cell cycle in the plant host. In the infected host cell, HSP70, but not HSP90, participates in the translocation of CP from the cytoplasm into the nucleus. Viral amounts decrease when HSP70 is inhibited, but increase when HSP90 is downregulated. In the whitefly vector, HSP70 impairs the circulative transmission of TYLCV; its inhibition increases transmission. Hence, the efficiency of virus acquisition by whiteflies depends on the functionality of both plant chaperones and their cross-talk with other protein mechanisms controlling virus-induced aggregation.

Whitefly interactions with plants

Whiteflies are important pests of many crops worldwide. They are polyphagous and effectively feed on phloem sap using mouthparts modified into long, flexible stylets. Plants respond to whitefly attack by activating defense genes leading to production of toxic compounds. To reach plant phloem and survive on host plants, whiteflies secret effectors in the saliva to regulate plant responses and activate detoxification system to cope with plant defenses. Additionally, whitefly-transmitted viruses may exert substantial effects on host plants and in turn the performance of whiteflies. Understanding the interactions between whiteflies and host plants will promote the development of novel strategies for controlling whiteflies. Here, we summarize the genetics, molecular genetics and genomics of the whitefly's interactions with plants.

The draft genome of whitefly Bemisia tabaci MEAM1, a global crop pest, provides novel insights into virus transmission, host adaptation, and insecticide resistance

Background

The whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) is among the 100 worst invasive species in the world. As one of the most important crop pests and virus vectors, B. tabaci causes substantial crop losses and poses a serious threat to global food security.

Results

We report the 615-Mb high-quality genome sequence of B. tabaci Middle East-Asia Minor 1 (MEAM1), the first genome sequence in the Aleyrodidae family, which contains 15,664 protein-coding genes. The B. tabaci genome is highly divergent from other sequenced hemipteran genomes, sharing no detectable synteny. A number of known detoxification gene families, including cytochrome P450s and UDP-glucuronosyltransferases, are significantly expanded in B. tabaci. Other expanded gene families, including cathepsins, large clusters of tandemly duplicated B. tabaci-specific genes, and phosphatidylethanolamine-binding proteins (PEBPs), were found to be associated with virus acquisition and transmission and/or insecticide resistance, likely contributing to the global invasiveness and efficient virus transmission capacity of B. tabaci. The presence of 142 horizontally transferred genes from bacteria or fungi in the B. tabaci genome, including genes encoding hopanoid/sterol synthesis and xenobiotic detoxification enzymes that are not present in other insects, offers novel insights into the unique biological adaptations of this insect such as polyphagy and insecticide resistance. Interestingly, two adjacent bacterial pantothenate biosynthesis genes, panB and panC, have been co-transferred into B. tabaci and fused into a single gene that has acquired introns during its evolution.

Conclusions

The B. tabaci genome contains numerous genetic novelties, including expansions in gene families associated with insecticide resistance, detoxification and virus transmission, as well as numerous horizontally transferred genes from bacteria and fungi. We believe these novelties likely have shaped B. tabaci as a highly invasive polyphagous crop pest and efficient vector of plant viruses. The genome serves as a reference for resolving the B. tabaci cryptic species complex, understanding fundamental biological novelties, and providing valuable genetic information to assist the development of novel strategies for controlling whiteflies and the viruses they transmit.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0321-y) contains supplementary material, which is available to authorized users.

Transcriptomic Analysis of the Endangered Neritid Species Clithon retropictus: De Novo Assembly, Functional Annotation, and Marker Discovery

An aquatic gastropod belonging to the family Neritidae, Clithon retropictus is listed as an endangered class II species in South Korea. The lack of information on its genomic background limits the ability to obtain functional data resources and inhibits informed conservation planning for this species. In the present study, the transcriptomic sequencing and de novo assembly of C. retropictus generated a total of 241,696,750 high-quality reads. These assembled to 282,838 unigenes with mean and N50 lengths of 736.9 and 1201 base pairs, respectively. Of these, 125,616 unigenes were subjected to annotation analysis with known proteins in Protostome DB, COG, GO, and KEGG protein databases (BLASTX; E ≤ 0.00001) and with known nucleotides in the Unigene database (BLASTN; E ≤ 0.00001). The GO analysis indicated that cellular process, cell, and catalytic activity are the predominant GO terms in the biological process, cellular component, and molecular function categories, respectively. In addition, 2093 unigenes were distributed in 107 different KEGG pathways. Furthermore, 49,280 simple sequence repeats were identified in the unigenes (>1 kilobase sequences). This is the first report on the identification of transcriptomic and microsatellite resources for C. retropictus, which opens up the possibility of exploring traits related to the adaptation and acclimatization of this species.

Transcription analysis of neonicotinoid resistance in Mediterranean (MED) populations of B. tabaci reveal novel cytochrome P450s, but no nAChR mutations associated with the phenotype

Background

Bemisia tabaci is one of the most damaging agricultural pests world-wide. Although its control is based on insecticides, B. tabaci has developed resistance against almost all classes of insecticides, including neonicotinoids.

Results

We employed an RNA-seq approach to generate genome wide expression data and identify genes associated with neonicotinoid resistance in Mediterranean (MED) B. tabaci (Q1 biotype). Twelve libraries from insecticide resistant and susceptible whitefly populations were sequenced on an Illumina Next-generation sequencing platform, and genomic sequence information of approximately 73 Gbp was generated.

A reference transcriptome was built by de novo assembly and functionally annotated. A total of 146 P450s, 18 GSTs and 23 CCEs enzymes (unigenes) potentially involved in the detoxification of xenobiotics were identified, along with 78 contigs encoding putative target proteins of six different insecticide classes. Ten unigenes encoding nicotinic Acetylcholine Receptors (nAChR), the target of neoinicotinoids, were identified and phylogenetically classified. No nAChR polymorphism potentially related with the resistant phenotypes, was observed among the studied strains.

DE analysis revealed that among the 550 differentially (logFC > 1) over-transcribed unigenes, 52 detoxification enzymes were over expressed including unigenes with orthologues in P450s, GSTs, CCE and UDP-glucuronosyltransferases.

Eight P450 unigenes belonging to clades CYP2, CYP3 and CYP4 were highly up-regulated (logFC > 2) including CYP6CM1, a gene already known to confer imidacloprid resistance in B. tabaci. Using quantitative qPCRs, a larger screening of field MED B. tabaci from Crete with known neonicotinoid phenotype was performed to associate expression levels of P450s with resistance levels. Expression levels of five P450s, including CYP6CM1, were found associated with neonicotinoid resistance. However, a significant correlation was found only in CYP303 and CYP6CX3, with imidacloprid and acetamiprid respectively.

Conclusion

Our work has generated new toxicological data and genomic resources which will significantly enrich the available dataset and substantially facilitate the molecular studies in MED B. tabaci. No evidence of target site neonicotinoid resistance has been found. Eight P450 unigenes, including CYP6CM1, were found significantly over-expressed in resistant B. tabaci. This study suggests at least two novel P450s (CYP303 and CYP6CX3) as candidates for their functional characterization as detoxification mechanisms of neonicotinoid resistance in B. tabaci.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2161-5) contains supplementary material, which is available to authorized users.

De Novo Assembly and Characterization of the Invasive Northern Pacific Seastar Transcriptome

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI’s non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

Metabolic Coevolution in the Bacterial Symbiosis of Whiteflies and Related Plant Sap-Feeding Insects

Genomic decay is a common feature of intracellular bacteria that have entered into symbiosis with plant sap-feeding insects. This study of the whitefly Bemisia tabaci and two bacteria (Portiera aleyrodidarum and Hamiltonella defensa) cohoused in each host cell investigated whether the decay of Portiera metabolism genes is complemented by host and Hamiltonella genes, and compared the metabolic traits of the whitefly symbiosis with other sap-feeding insects (aphids, psyllids, and mealybugs). Parallel genomic and transcriptomic analysis revealed that the host genome contributes multiple metabolic reactions that complement or duplicate Portiera function, and that Hamiltonella may contribute multiple cofactors and one essential amino acid, lysine. Homologs of the Bemisia metabolism genes of insect origin have also been implicated in essential amino acid synthesis in other sap-feeding insect hosts, indicative of parallel coevolution of shared metabolic pathways across multiple symbioses. Further metabolism genes coded in the Bemisia genome are of bacterial origin, but phylogenetically distinct from Portiera, Hamiltonella and horizontally transferred genes identified in other sap-feeding insects. Overall, 75% of the metabolism genes of bacterial origin are functionally unique to one symbiosis, indicating that the evolutionary history of metabolic integration in these symbioses is strongly contingent on the pattern of horizontally acquired genes. Our analysis, further, shows that bacteria with genomic decay enable host acquisition of complex metabolic pathways by multiple independent horizontal gene transfers from exogenous bacteria. Specifically, each horizontally acquired gene can function with other genes in the pathway coded by the symbiont, while facilitating the decay of the symbiont gene coding the same reaction.

Transcriptome Sequencing of Lima Bean (Phaseolus lunatus) to Identify Putative Positive Selection in Phaseolus and Legumes

The identification of genes under positive selection is a central goal of evolutionary biology. Many legume species, including Phaseolus vulgaris (common bean) and Phaseolus lunatus (lima bean), have important ecological and economic value. In this study, we sequenced and assembled the transcriptome of one Phaseolus species, lima bean. A comparison with the genomes of six other legume species, including the common bean, Medicago, lotus, soybean, chickpea, and pigeonpea, revealed 15 and 4 orthologous groups with signatures of positive selection among the two Phaseolus species and among the seven legume species, respectively. Characterization of these positively selected genes using Non redundant (nr) annotation, gene ontology (GO) classification, GO term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes are mostly involved in thylakoids, photosynthesis and metabolism. This study identified genes that may be related to the divergence of the Phaseolus and legume species. These detected genes are particularly good candidates for subsequent functional studies.

Transcriptional responses of invasive and indigenous whiteflies to different host plants reveal their disparate capacity of adaptation

The whitefly Bemisia tabaci contains more than 35 cryptic species. The higher adaptability of Middle East-Asia Minor 1 (MEAM1) cryptic species has been recognized as one important factor for its invasion and displacement of other indigenous species worldwide. Here we compared the performance of the invasive MEAM1 and the indigenous Asia II 3 whitefly species following host plant transfer from a suitable host (cotton) to an unsuitable host (tobacco) and analyzed their transcriptional responses. After transfer to tobacco for 24 h, MEAM1 performed much better than Asia II 3. Transcriptional analysis showed that the patterns of gene regulation were very different with most of the genes up-regulated in MEAM1 but down-regulated in Asia II 3. Whereas carbohydrate and energy metabolisms were repressed in Asia II 3, the gene expression and protein metabolisms were activated in MEAM1. Compared to the constitutive high expression of detoxification genes in MEAM1, most of the detoxification genes were down-regulated in Asia II 3. Enzymatic activities of P450, GST and esterase further verified that the detoxification of MEAM1 was much higher than that of Asia II 3. These results reveal obvious differences in responses of MEAM1 and Asia II 3 to host transfer.

Whitefly genome expression reveals host-symbiont interaction in amino acid biosynthesis

Background

Whitefly (Bemisia tabaci) complex is a serious insect pest of several crop plants worldwide. It comprises several morphologically indistinguishable species, however very little is known about their genetic divergence and biosynthetic pathways. In the present study, we performed transcriptome sequencing of Asia 1 species of B. tabaci complex and analyzed the interaction of host-symbiont genes in amino acid biosynthetic pathways.

Methodology/Principal Findings

We obtained about 83 million reads using Illumina sequencing that assembled into 72716 unitigs. A total of 21129 unitigs were annotated at stringent parameters. Annotated unitigs were mapped to 52847 gene ontology (GO) terms and 131 Kyoto encyclopedia of genes and genomes (KEGG) pathways. Expression analysis of the genes involved in amino acid biosynthesis pathways revealed the complementation between whitefly and its symbiont partner Candidatus Portiera aleyrodidarum. Most of the non-essential amino acids and intermediates of essential amino acid pathways were supplied by the host insect to its symbiont. The symbiont expressed the pathways for the essential amino acids arginine, threonine and tryptophan and the immediate precursors of valine, leucine, isoleucine and phenyl-alanine. High level expression of the amino acid transporters in the whitefly suggested the molecular mechanisms for the exchange of amino acids between the host and the symbiont.

Conclusions/Significance

Our study provides a comprehensive transcriptome data for Asia 1 species of B. tabaci complex that focusses light on integration of host and symbiont genes in amino acid biosynthesis pathways.

Transcriptome analysis between invasive Pomacea canaliculata and indigenous Cipangopaludina cahayensis reveals genomic divergence and diagnostic microsatellite/SSR markers

BACKGROUND

Pomacea canaliculata is an important invasive species worldwide. However, little is known about the molecular mechanisms behind species displacement, adaptational abilities, and pesticide resistance, partly because of the lack of genomic information that is available for this species. Here, the transcriptome sequences for the invasive golden apple snail P. canaliculata and the native mudsnail Cipangopaludina cahayensis were obtained by next-generation-sequencing and used to compare genomic divergence and identify molecular markers.

RESULTS

More than 46 million high quality sequencing reads were generated from P. canaliculata and C. cahayensis using Illumina paired-end sequencing technology. Our analysis indicated that 11,312 unigenes from P. canaliculata and C. cahayensis showed significant similarities to known proteins families, among which a total of 4,320 specific protein families were identified. KEGG pathway enrichment was analyzed for the unique unigenes with 17 pathways (p-value < 10(-5)) in P. canaliculata relating predominantly to lysosomes and vitamin digestion and absorption, and with 12 identified in C. cahayensis, including cancer and toxoplasmosis pathways, respectively. Our analysis also indicated that the comparatively high number of P450 genes in the P. canaliculata transcriptome may be associated with the pesticide resistance in this species. Additionally, 16,717 simple sequence repeats derived from expressed sequence tags (EST-SSRs) were identified from the 14,722 unigenes in P. canaliculata and 100 of them were examined by PCR, revealing a species-specific molecular marker that could distinguish between the morphologically similar P. canaliculata and C. cahayensis snails.

CONCLUSIONS

Here, we present the genomic resources of P. canaliculata and C. cahayensis. Differentially expressed genes in the transcriptome of P. canaliculata compared with C. cahayensis corresponded to critical metabolic pathways, and genes specifically related to environmental stress response were detected. The CYP4 family of P450 cytochromes that may be important factors in pesticide metabolism in P. canaliculata was identified. Overall, these findings will provide valuable genetic data for the further characterization of the molecular mechanisms that support the invasive and adaptive abilities of P. canaliculata.

Cloning of a putative extracellular Cu/Zn superoxide dismutase and functional differences of superoxide dismutases in invasive and indigenous whiteflies

Superoxide dismutases (SODs) are a group of important antioxidant defense enzymes. In this study, a putative extracellular Cu/Zn superoxide dismutase (ecCuZnSOD) complementary DNA was cloned and characterized from the whitefly, Bemisia tabaci. Quantitative polymerase chain reaction analysis showed that the expression level of BtecCuZnSOD was more than 10-fold higher in the invasive Middle East Asia Minor 1 (MEAM1) than in the native Asia II 3 species of the B. tabaci species complex. After exposure to low temperature (4 °C), the expression of Bt-ecCuZnSOD gene was significantly up-regulated in MEAM1 but not in Asia II 3. Furthermore, the expression level of B. tabaci intracellular CuZnSOD (Bt-icCuZnSOD), Bt-ecCuZnSOD and mitochondrial MnSOD (Bt-mMnSOD) was compared after transferring MEAM1 and Asia II 3 whiteflies from favorable (cotton) to unfavorable host plants (tobacco). On cotton, both CuZnSOD genes were expressed at a higher level in MEAM1 compared with Asia II 3. Interestingly, after transferring onto tobacco, the expression of Bt-ecCuZnSOD was significantly induced in Asia II 3 but not in MEAM1. On the other hand, while Bt-mMnSOD was expressed equally in both species on cotton, Bt-mMnSOD messenger RNA was up-regulated in MEAM1 on tobacco. Consistently, enzymatic activity assays of CuZnSOD and MnSOD demonstrated that CuZnSOD might play an important protective role against oxidative stress in Asia II 3, whereas MnSOD activation was critical for MEAM1 whiteflies during host adaptation. Taken together, our results suggest that the successful invasion of MEAM1 is correlated with its constitutive high activity of CuZnSOD and inducible expression of MnSOD under stress conditions.

Transcriptome analysis of host-associated differentiation in Bemisia tabaci (Hemiptera: Aleyrodidae)

Host-associated differentiation is one of the driving forces behind the diversification of phytophagous insects. In this study, host induced transcriptomic differences were investigated in the sweetpotato whitefly Bemisia tabaci, an invasive agricultural pest worldwide. Comparative transcriptomic analyses using coding sequence (CDS), 5′ and 3′ untranslated regions (UTR) showed that sequence divergences between the original host plant, cabbage, and the derived hosts, including cotton, cucumber and tomato, were 0.11–0.14%, 0.19–0.26%, and 0.15–0.21%, respectively. In comparison to the derived hosts, 418 female and 303 male transcripts, respectively, were up-regulated in the original cabbage strain. Among them, 17 transcripts were consistently up-regulated in both female and male whiteflies originated from the cabbage host. Specifically, two ESTs annotated as Cathepsin B or Cathepsin B-like genes were significantly up-regulated in the original cabbage strain, representing a transcriptomic response to the dietary challenges imposed by the host shifting. Results from our transcriptome analysis, in conjunction with previous reports documenting the minor changes in their reproductive capacity, insecticide susceptibility, symbiotic composition and feeding behavior, suggest that the impact of host-associated differentiation in whiteflies is limited. Furthermore, it is unlikely the major factor contributing to their rapid range expansion/invasiveness.

Transient receptor potential is essential for high temperature tolerance in invasive Bemisia tabaci Middle East Asia minor 1 cryptic species

Temperature is an important factor in affecting population dynamics and diffusion distribution of organisms. Alien species can successfully invade and colonize to various temperature environments, and one of important reasons is that alien species have a strong resistance to stress temperature. Recently, researchers have focused on the mechanisms of temperature sensing to determine the sensing and regulation mechanisms of temperature adaptation. The transient receptor potential (TRP) is one of the key components of an organism’s temperature perception system. TRP plays important roles in perceiving temperature, such as avoiding high temperature, low temperature and choosing the optimum temperature. To assess high temperature sensation and the heat resistance role of the TRP gene, we used 3′ and 5′ rapid-amplification of cDNA ends to isolate the full-length cDNA sequence of the TRP gene from Bemisia tabaci (Gennadius) MEAM1 (Middle East Asia Minor 1), examined the mRNA expression profile under various temperature conditions, and identified the heat tolerance function. This is the first study to characterize the TRP gene of invasive B. tabaci MEAM1 (MEAM1 BtTRP). The full-length cDNA of MEAM1 BtTRP was 3871 bp, and the open reading frames of BtTRP was 3501 bp, encoding 1166 amino acids. Additionally, the BtTRP mRNA expression level was significantly increased at 35°C. Furthermore, compared with control treatments, the survival rate of B. tabaci MEAM1 adults was significantly decreased under high temperature stress conditions after feeding with dsRNA BtTRP. Collectively, these results showed that MEAM1 BtTRP is a key element in sensing high temperature and plays an essential role in B. tabaci MEAM1 heat tolerance ability. Our data improved our understanding of the mechanism of temperature sensation in B. tabaci MEAM1 at the molecular level and could contribute to the understanding of the thermal biology of B. tabaci MEAM1 within the context of global climate change.

Developing conversed microsatellite markers and their implications in evolutionary analysis of the Bemisia tabaci complex

The study of population genetics among the Bemisia tabaci complex is limited due to the lack of conserved molecular markers. In this study, 358, 433 and 322 new polynucleotide microsatellites are separately identified from the transcriptome sequences of three cryptic species of the B. tabaci complex. The cross species transferability of 57 microsatellites was then experimentally validated. The results indicate that these markers are conserved and have high inter-taxon transferability. Thirteen markers were employed to assess the genetic relationships among six cryptic species of the B. tabaci complex. To our surprise, the inferred phylogeny was consistent with that of mitochondrial COI sequences, indicating that microsatellites have the potential to distinguish species of the B. tabaci complex. Our results demonstrate that development of microsatellites from transcriptome data is a fast and cost-effective approach. These markers can be used to analyze the population genetics and evolutionary patterns of the B. tabaci complex.

Transcriptome analysis of the oil-rich tea plant, Camellia oleifera, reveals candidate genes related to lipid metabolism

Background

Rapidly driven by the need for developing sustainable sources of nutritionally important fatty acids and the rising concerns about environmental impacts after using fossil oil, oil-plants have received increasing awareness nowadays. As an important oil-rich plant in China, Camellia oleifera has played a vital role in providing nutritional applications, biofuel productions and chemical feedstocks. However, the lack of C. oleifera genome sequences and little genetic information have largely hampered the urgent needs for efficient utilization of the abundant germplasms towards modern breeding efforts of this woody oil-plant.

Results

Here, using the 454 GS-FLX sequencing platform, we generated approximately 600,000 RNA-Seq reads from four tissues of C. oleifera. These reads were trimmed and assembled into 104,842 non-redundant putative transcripts with a total length of ∼38.9 Mb, representing more than 218-fold of all the C. oleifera sequences currently deposited in the GenBank (as of March 2014). Based on the BLAST similarity searches, nearly 42.6% transcripts could be annotated with known genes, conserved domains, or Gene Ontology (GO) terms. Comparisons with the cultivated tea tree, C. sinensis, identified 3,022 pairs of orthologs, of which 211 exhibited the evidence under positive selection. Pathway analysis detected the majority of genes potentially related to lipid metabolism. Evolutionary analysis of omega-6 fatty acid desaturase (FAD2) genes among 20 oil-plants unexpectedly suggests that a parallel evolution may occur between C. oleifera and Olea oleifera. Additionally, more than 2,300 simple sequence repeats (SSRs) and 20,200 single-nucleotide polymorphisms (SNPs) were detected in the C. oleifera transcriptome.

Conclusions

The generated transcriptome represents a considerable increase in the number of sequences deposited in the public databases, providing an unprecedented opportunity to discover all related-genes associated with lipid metabolic pathway in C. oleifera. It will greatly enhance the generation of new varieties of C. oleifera with increased yields and high quality.

Transcriptomic analyses reveal the adaptive features and biological differences of guts from two invasive whitefly species

Background

The gut of phloem feeding insects is critical for nutrition uptake and xenobiotics degradation. However, partly due to its tiny size, genomic information for the gut of phloem feeding insects is limited.

Results

In this study, the gut transcriptomes of two species of invasive whiteflies in the Bemisia tabaci complex, Middle East Asia Minor 1 (MEAM1) and Mediterranean (MED), were analyzed using the Illumina sequencing. A total of 12,879 MEAM1 transcripts and 11,246 MED transcripts were annotated with a significant Blastx hit. In addition, 7,000 and 5,771 gut specific genes were respectively identified for MEAM1 and MED. Functional analyses on these gut specific genes demonstrated the important roles of gut in metabolism of insecticides and secondary plant chemicals. To reveal the molecular difference between guts of MEAM1 and MED, a comparison between gut transcriptomes of the two species was conducted and 3,910 pairs of orthologous genes were identified. Based on the ratio of nonsynonymous and synonymous substitutions, 15 genes were found evolving under positive selection. Many of those genes are predicted to be involved in metabolism and insecticide resistance. Furthermore, many genes related to detoxification were expressed at an elevated level in the gut of MED compared to MEAM1, which might be responsible for the MED’s higher resistance to insecticides and environmental stresses.

Conclusion

The sequencing of MED and MEAM1 gut transcriptomes and extensive comparisons of MEAM1 and MED gut transcripts provide substantial sequence information for revealing the role of gut in whiteflies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-370) contains supplementary material, which is available to authorized users.

Disruption of the homogentisate solanesyltransferase gene results in albino and dwarf phenotypes and root, trichome and stomata defects in Arabidopsis thaliana

Homogentisate solanesyltransferase (HST) plays an important role in plastoquinone (PQ) biosynthesis and acts as the electron acceptor in the carotenoids and abscisic acid (ABA) biosynthesis pathways. We isolated and identified a T-DNA insertion mutant of the HST gene that displayed the albino and dwarf phenotypes. PCR analyses and functional complementation also confirmed that the mutant phenotypes were caused by disruption of the HST gene. The mutants also had some developmental defects, including trichome development and stomata closure defects. Chloroplast development was also arrested and chlorophyll (Chl) was almost absent. Developmental defects in the chloroplasts were consistent with the SDS-PAGE result and the RNAi transgenic phenotype. Exogenous gibberellin (GA) could partially rescue the dwarf phenotype and the root development defects and exogenous ABA could rescue the stomata closure defects. Further analysis showed that ABA and GA levels were both very low in the pds2-1 mutants, which suggested that biosynthesis inhibition by GAs and ABA contributed to the pds2-1 mutants' phenotypes. An early flowering phenotype was found in pds2-1 mutants, which showed that disruption of the HST gene promoted flowering by partially regulating plant hormones. RNA-sequencing showed that disruption of the HST gene resulted in expression changes to many of the genes involved in flowering time regulation and in the biosynthesis of PQ, Chl, GAs, ABA and carotenoids. These results suggest that HST is essential for chloroplast development, hormone biosynthesis, pigment accumulation and plant development.

Global haplotype analysis of the whitefly Bemisia tabaci cryptic species Asia I in Asia

The whitefly, Bemisia tabaci (Hemiptera: Aleyrodidiae), is a cryptic species complex comprising a minimum of 24 cryptic species. Some members of this complex are important agricultural pests, causing considerable damage to vegetable as well as ornamental and horticultural crops. Asia I, one of the cryptic species of B. tabaci, is widely distributed in Asia. One hundred and sixty mitochondrial cytochrome oxidase I (COI) sequences from eight countries have been analyzed to investigate the geographic origin and current genetic structure of this cryptic species. Sixty different haplotypes were identified, with levels of genetic distances ranging from 0.001 to 0.021. A sign of possible genetic differentiation emerges from the differential distribution of dominant haplotypes in Indonesia and India compared to China. A possible ancient separation between Asia I in India and Indonesia and secondary contact in China has been hypothesized.

Parallel histories of horizontal gene transfer facilitated extreme reduction of endosymbiont genomes in sap-feeding insects

Bacteria confined to intracellular environments experience extensive genome reduction. In extreme cases, insect endosymbionts have evolved genomes that are so gene-poor that they blur the distinction between bacteria and endosymbiotically derived organelles such as mitochondria and plastids. To understand the host's role in this extreme gene loss, we analyzed gene content and expression in the nuclear genome of the psyllid Pachypsylla venusta, a sap-feeding insect that harbors an ancient endosymbiont (Carsonella) with one of the most reduced bacterial genomes ever identified. Carsonella retains many genes required for synthesis of essential amino acids that are scarce in plant sap, but most of these biosynthetic pathways have been disrupted by gene loss. Host genes that are upregulated in psyllid cells housing Carsonella appear to compensate for endosymbiont gene losses, resulting in highly integrated metabolic pathways that mirror those observed in other sap-feeding insects. The host contribution to these pathways is mediated by a combination of native eukaryotic genes and bacterial genes that were horizontally transferred from multiple donor lineages early in the evolution of psyllids, including one gene that appears to have been directly acquired from Carsonella. By comparing the psyllid genome to a recent analysis of mealybugs, we found that a remarkably similar set of functional pathways have been shaped by independent transfers of bacterial genes to the two hosts. These results show that horizontal gene transfer is an important and recurring mechanism driving coevolution between insects and their bacterial endosymbionts and highlight interesting similarities and contrasts with the evolutionary history of mitochondria and plastids.

siRNA machinery in whitefly (Bemisia tabaci)

BACKGROUND

RNA interference has been emerged as an utmost tool for the control of sap sucking insect pests. Systemic response is necessary to control them in field condition. Whitefly is observed to be more prone to siRNA in recent studies, however the siRNA machinery and mechanism is not well established.

METHODOLOGY/PRINCIPAL FINDINGS

To identify the core siRNA machinery, we curated transcriptome data of whitefly from NCBI database. Partial mRNA sequences encoding Dicer2, R2D2, Argonaute2 and Sid1 were identified by tblastn search of homologous sequences from Aphis glycines and Tribolium castaneum. Complete encoding sequences were obtained by RACE, protein sequences derived by Expasy translate tool and confirmed by blastp analysis. Conserved domain search and Prosite-Scan showed similar domain architecture as reported in homologs from related insects. We found helicase, PAZ, RNaseIIIa, RNaseIIIb and double-stranded RNA-binding fold (DSRBF) in Dicer2; DsRBD in R2D2; and PAZ and PIWI domains in Argonaute2. Eleven transmembrane domains were detected in Sid1. Sequence homology and phylogenetic analysis revealed that RNAi machinery of whitefly is close to Aphids. Real-time PCR analysis showed similar expression of these genes in different developmental stages as reported in A. glycines and T. castaneum. Further, the expression level of above genes was quite similar to the housekeeping gene actin.

CONCLUSIONS/SIGNIFICANCE

Availability of core siRNA machinery including the Sid1 and their universal expression in reasonable quantity indicated significant response of whitefly towards siRNA. Present report opens the way for controlling whitefly, one of the most destructive crop insect pest.

High Variation in Single Nucleotide Polymorphisms (SNPs) and Insertions/Deletions (Indels) in the Highly Invasive Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Middle East-Asia Minor 1 (MEAM1)

Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) Middle East-Asia Minor 1 (MEAM1) is invasive and adaptive to varied environments throughout the world. The adaptability is closely related to genomic variation such as single nucleotide polymorphisms (SNPs) and insertions/deletions (indels). In order to elucidate the feature of SNPs and indels in MEAM1, and reveal the association between SNPs/indels and adaptive capacity to various environments, a computational approach with QualitySNP was used to identify reliable SNPs and indels on the basis of 9110-expressed sequence tags of MEAM1 present in the NCBI database. There were 575 SNPs detected with a density of 10.1 SNPs/kb and 6.4 SNPs/contig. Also, 237 transitions (39.3%) and 366 transversions (60.7%) were obtained, where the ratio of transitions to transversions was 0.65:1. In addition, 581 indels with a density of 14.1 indels/kb and 9.2 indels/contig were detected. Collectively, it showed that invasive MEAM1 has high SNPs density, and higher SNPs percentage than non-invasive B. tabaci species. A high SNPs density/percentage in MEAM1 yielded a high genomic variation that might have allowed it to adapt to varied environments, which provides some support to understand the invasive nature of MEAM1 at the genomic level. High levels of genomic variation are implicated in the level of adaptive capacity and invasive species are thought to exhibit higher levels of adaptive capacity than non-invasive species.

Genomewide variation in an introgression line of rice-Zizania revealed by whole-genome re-sequencing

BACKGROUND

Hybridization between genetically diverged organisms is known as an important avenue that drives plant genome evolution. The possible outcomes of hybridization would be the occurrences of genetic instabilities in the resultant hybrids. It remained under-investigated however whether pollination by alien pollens of a closely related but sexually "incompatible" species could evoke genomic changes and to what extent it may result in phenotypic novelties in the derived progenies.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we have re-sequenced the genomes of Oryza sativa ssp. japonica cv. Matsumae and one of its derived introgressant RZ35 that was obtained from an introgressive hybridization between Matsumae and Zizanialatifolia Griseb. in general, 131 millions 90 base pair (bp) paired-end reads were generated which covered 13.2 and 21.9 folds of the Matsumae and RZ35 genomes, respectively. Relative to Matsumae, a total of 41,724 homozygous single nucleotide polymorphisms (SNPs) and 17,839 homozygous insertions/deletions (indels) were identified in RZ35, of which 3,797 SNPs were nonsynonymous mutations. Furthermore, rampant mobilization of transposable elements (TEs) was found in the RZ35 genome. The results of pathogen inoculation revealed that RZ35 exhibited enhanced resistance to blast relative to Matsumae. Notably, one nonsynonymous mutation was found in the known blast resistance gene Pid3/Pi25 and real-time quantitative (q) RT-PCR analysis revealed constitutive up-regulation of its expression, suggesting both altered function and expression of Pid3/Pi25 may be responsible for the enhanced resistance to rice blast by RZ35.

CONCLUSIONS/SIGNIFICANCE

Our results demonstrate that introgressive hybridization by Zizania has provoked genomewide, extensive genomic changes in the rice genome, and some of which have resulted in important phenotypic novelties. These findings suggest that introgressive hybridization by alien pollens of even a sexually incompatible species may represent a potent means to generate novel genetic diversities, and which may have played relevant roles in plant evolution and can be manipulated for crop improvements.

Signatures of rapid evolution in urban and rural transcriptomes of white-footed mice (Peromyscus leucopus) in the New York metropolitan area

Urbanization is a major cause of ecological degradation around the world, and human settlement in large cities is accelerating. New York City (NYC) is one of the oldest and most urbanized cities in North America, but still maintains 20% vegetation cover and substantial populations of some native wildlife. The white-footed mouse, Peromyscusleucopus, is a common resident of NYC's forest fragments and an emerging model system for examining the evolutionary consequences of urbanization. In this study, we developed transcriptomic resources for urban P. leucopus to examine evolutionary changes in protein-coding regions for an exemplar "urban adapter." We used Roche 454 GS FLX+ high throughput sequencing to derive transcriptomes from multiple tissues from individuals across both urban and rural populations. From these data, we identified 31,015 SNPs and several candidate genes potentially experiencing positive selection in urban populations of P. leucopus. These candidate genes are involved in xenobiotic metabolism, innate immune response, demethylation activity, and other important biological phenomena in novel urban environments. This study is one of the first to report candidate genes exhibiting signatures of directional selection in divergent urban ecosystems.

Silencing the ecdysone synthesis and signaling pathway genes disrupts nymphal development in the whitefly

Sap-sucking insects are important pests in agriculture and good models to study insect biology. The role of ecdysone pathway genes in the life history of this group of insects is largely unknown likely due to a lack of efficient gene silencing methods allowing functional genetic analyses. Here, we developed a new and high throughput method to silence whitefly genes using a leaf-mediated dsRNA feeding method. We have applied this method to explore the roles of genes within the molting hormone-ecdysone synthesis and signaling pathway for the survival, reproduction and development of whiteflies. Silencing of genes in the ecdysone pathway had a limited effect on the survival and fecundity of adult whiteflies. However, gene silencing reduced survival and delayed development of the whitefly during nymphal stages. These data suggest that the silencing method developed here provides a useful tool for functional gene discovery studies of sap-sucking insects, and further indicate the potential of regulating the ecdysone pathway in whitefly control.

Analysis of the transcriptional differences between indigenous and invasive whiteflies reveals possible mechanisms of whitefly invasion

BACKGROUND

The whitefly Bemisa tabaci is a species complex of more than 31 cryptic species which include some of the most destructive invasive pests of crops worldwide. Among them, Middle East-Asia Minor 1 (MEAM1) and Mediterranean have invaded many countries and displaced the native whitefly species. The successful invasion of the two species is largely due to their wide range of host plants, high resistance to insecticides and remarkable tolerance to environmental stresses. However, the molecular differences between invasive and indigenous whiteflies remain largely unknown.

METHODOLOGY/PRINCIPAL FINDINGS

Here the global transcriptional difference between the two invasive whitefly species (MEAM1, MED) and one indigenous whitefly species (Asia II 3) were analyzed using the Illumina sequencing. Our analysis indicated that 2,422 genes between MEAM1 and MED; 3,073 genes between MEAM1 and Asia II 3; and 3,644 genes between MED and Asia II 3 were differentially expressed. Gene Ontology enrichment analysis revealed that the differently expressed genes between the invasive and indigenous whiteflies were significantly enriched in the term of 'oxidoreductase activity'. Pathway enrichment analysis showed that carbohydrate, amino acid and glycerolipid metabolisms were more active in MEAM1 and MED than in Asia II 3, which may contribute to their differences in biological characteristics. Our analysis also illustrated that the majority of genes involved in 'drug metabolic pathway' were expressed at a higher level in MEAM1 and MED than in Asia II 3. Taken together, these results revealed that the genes related to basic metabolism and detoxification were expressed at an elevated level in the invasive whiteflies, which might be responsible for their higher resistance to insecticides and environmental stresses.

CONCLUSIONS/SIGNIFICANCE

The extensive comparison of MEAM1, MED and Asia II 3 gene expression may serve as an invaluable resource for revealing the molecular mechanisms underlying their biological differences and the whitefly invasion.