Cassava (Manihot esculenta Crantz) production can be severely affected by cassava mosaic disease (CMD) caused by viruses in the genus Begomovirus of the family Geminiviridae. Eight begomoviruses associated with CMD have been recorded on the African continent and two on the Indian subcontinent (Brown et al. 2015). Duraisamy et al. (2013) reported that two cassava geminiviruses on the Indian subcontinent, i.e. Indian cassava mosaic virus (ICMV) and Sri Lankan cassava mosaic virus were transmitted by the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae), but the genetic group(s) or cryptic species of the whitefly responsible for the transmission were not determined, as B. tabaci is now known as a complex consisting of >35 morphologically indistinguishable species (Liu et al. 2012). Cassava cultivation in Southeast Asia started in early 1900s, and today it has become one of the major crops in many countries in this region including Cambodia (FAOSTAT). However, cassava in Southeast Asia has not prev...

The whiteflies under the name Bemisia tabaci (Gennadius) (Aleyrodidae: Hemiptera) are species complex of at least 31 cryptic species some of which are globally invasive agricultural pests. Previously, the mitochondrial genome (mitogenome) of the indigenous New World B. tabaci species was sequenced and major differences of gene order from the postulated whitefly ancestral gene order were found. However, the sequence and gene order of mitogenomes in other B. tabaci species are unknown. In addition, the sequence divergences and gene expression profiles of mitogenomes in the B. tabaci species complex remain completely unexplored. In this study, we obtained the complete mitogenome (15,632 bp) of the invasive Mediterranean (MED), which has been identified as the type species of the B. tabaci complex. It encodes 37 genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNAs and 22 transfer RNAs (tRNA). Comparative analyses of the mitogenomes from MED and New World (previously published) species reveal that there are no gene arrangements. Based on the Illumina sequencing data, the gene expression profile of the MED mitogenome was analyzed. We found that a number of genes were polyadenylated and the partial stop codons in cox1, cox2 and nd5 are completed via polyadenylation that changed T to the TAA stop codon. In addition, combining the transcriptome with the sequence alignment data, the possible termination site of some PCGs were defined. Our analyses also revealed that atp6 and atp8, nd4 and nd4l, nd6 and cytb were found on the same cistronic transcripts, whereas the other mature mitochondrial transcripts were monocistronic. Furthermore, RT-PCR analyses of the mitochondrial PCGs expression in different developmental stages revealed that the expression level of individual mitochondrial genes varied in each developmental stage of nymph, pupa and adult. Interestingly, mRNA levels showed significant differences among genes located in the same transcription unit suggesting that mitochondrial mRNA abundance is heavily modulated by post-transcriptional regulation. This work provides novel insights into the mitogenome evolution of B. tabaci species and demonstrates that utilizing RNA-seq data to obtain the mitogenome and analyze mitochondrial gene expression characteristics is practical.

The 37 currently recognized Bemisia tabaci cryptic species are economically important species and contain both primary and secondary endosymbionts, but their diversity has never been mapped systematically across the group. To achieve this, PacBio sequencing of full-length bacterial 16S rRNA gene amplicons was carried out on 21 globally collected species in the B. tabaci complex, and two samples from B. afer were used here as outgroups. The microbial diversity was first explored across the major lineages of the whole group and 15 new putative bacterial sequences were observed. Extensive comparison of our results with previous endosymbiont diversity surveys which used PCR or multiplex 454 pyrosequencing platforms showed that the bacterial diversity was underestimated. To validate these new putative bacteria, one of them (Halomonas) was first confirmed to be present in MED B. tabaci using Hiseq2500 and FISH technologies. These results confirmed PacBio is a reliable and informative venue to reveal the bacterial diversity of insects. In addition, many new secondary endosymbiotic strains of Rickettsia and Arsenophonus were found, increasing the known diversity in these groups. For the previously described primary endosymbionts, one Portiera Operational Taxonomic Units (OTU) was shared by all B. tabaci species. The congruence of the B. tabaci-host and Portiera phylogenetic trees provides strong support for the hypothesis that primary endosymbionts co-speciated with their hosts. Likewise, a comparison of bacterial alpha diversities, Principal Coordinate Analysis, indistinct endosymbiotic communities harbored by different species and the co-divergence analyses suggest a lack of association between overall microbial diversity with cryptic species, further indicate that the secondary endosymbiont-mediated speciation is unlikely to have occurred in the B. tabaci species group.

Summary: The bacterium Rickettsia is found widely in phytophagous insects and often exerts profound effects on the phenotype and fitness of its hosts. Here, we decrypt a new, independent, phylogenetically ancient Torix Rickettsia endosymbiont found constantly in a laboratory line of an economically important insect Asia II 7, a putative species of the Bemisia tabaci whitefly complex (Hemiptera: Aleyrodidae), and occasionally in field whitefly populations. This new Rickettsia distributes throughout the body of its whitefly host. Genetically, compared to Rickettsia_ bellii_MEAM1 found earlier in whiteflies, the new Rickettsia species has more gene families and pathways, which may be important factors in shaping specific symbiotic relationships. We propose the name ‘ Candidatus Rickettsia_Torix_ Bemisia_tabaci (RiTBt)’ for this new endosymbiont associated with whiteflies. Comparative genomic analyses indicate that RiTBi may be a relatively recent intruder in whiteflies given its low abundance in the field and relatively larger genome compared to Rickettsia _bellii_MEAM1.