A novel highly divergent protein family identified from a viviparous insect by RNA-seq analysis: a potential target for tsetse fly-specific abortifacients

Live-bearing cockroach genome reveals convergent evolutionary mechanisms linked to viviparity in insects and beyond

Live birth (viviparity) has arisen repeatedly and independently among animals. We sequenced the genome and transcriptome of the viviparous Pacific beetle-mimic cockroach and performed comparative analyses with two other viviparous insect lineages, tsetse flies and aphids, to unravel the basis underlying the transition to viviparity in insects. We identified pathways undergoing adaptive evolution for insects, involved in urogenital remodeling, tracheal system, heart development, and nutrient metabolism. Transcriptomic analysis of cockroach and tsetse flies revealed that uterine remodeling and nutrient production are increased and the immune response is altered during pregnancy, facilitating structural and physiological changes to accommodate and nourish the progeny. These patterns of convergent evolution of viviparity among insects, together with similar adaptive mechanisms identified among vertebrates, highlight that the transition to viviparity requires changes in urogenital remodeling, enhanced tracheal and heart development (corresponding to angiogenesis in vertebrates), altered nutrient metabolism, and shifted immunity in animal systems.

Reduced male fertility of an Antarctic mite following extreme heat stress could prompt localized population declines

Climate change is leading to substantial global thermal changes, which are particularly pronounced in polar regions. Therefore, it is important to examine the impact of heat stress on the reproduction of polar terrestrial arthropods, specifically, how brief extreme events may alter survival. We observed that sublethal heat stress reduces male fecundity in an Antarctic mite, yielding females that produced fewer viable eggs. Females and males collected from microhabitats with high temperatures showed a similar reduction in fertility. This impact is temporary, as indicated by recovery of male fecundity following return to cooler, stable conditions. The diminished fecundity is likely due to a drastic reduction in the expression of male-associated factors that occur in tandem with a substantial increase in the expression of heat shock proteins. Cross-mating between mites from different sites confirmed that heat-exposed populations have impaired male fertility. However, the negative impacts are transient as the effect on fertility declines with recovery time under less stressful conditions. Modeling indicated that heat stress is likely to reduce population growth and that short bouts of non-lethal heat stress could have substantial reproductive effects on local populations of Antarctic arthropods.

Lipid metabolism dysfunction following symbiont elimination is linked to altered Kennedy pathway homeostasis

Lipid metabolism is critical for insect reproduction, especially for species that invest heavily in the early developmental stages of their offspring. The role of symbiotic bacteria during this process is understudied but likely essential. We examined the role of lipid metabolism during the interaction between the viviparous tsetse fly (Glossina morsitans morsitans) and its obligate endosymbiotic bacteria (Wigglesworthia glossinidia) during tsetse pregnancy. We observed increased CTP:phosphocholine cytidylyltransferase (cct1) expression during pregnancy, which is critical for phosphatidylcholine biosynthesis in the Kennedy pathway. Experimental removal of Wigglesworthia impaired lipid metabolism via disruption of the Kennedy pathway, yielding obese mothers whose developing progeny starve. Functional validation via experimental cct1 suppression revealed a phenotype similar to females lacking obligate Wigglesworthia symbionts. These results indicate that, in Glossina, symbiont-derived factors, likely B vitamins, are critical for the proper function of both lipid biosynthesis and lipolysis to maintain tsetse fly fecundity.

The genome of sheep ked (Melophagus ovinus) reveals potential mechanisms underlying reproduction and narrower ecological niches

BACKGROUND

Melophagus ovinus is considered to be of great veterinary health significance. However, little is known about the information on genetic mechanisms of the specific biological characteristics and novel methods for controlling M. ovinus.

RESULTS

In total, the de novo genome assembly of M. ovinus was 188.421 Mb in size (330 scaffolds, N50 Length: 10.666 Mb), with a mean GC content of 27.74%. A total of 13,372 protein-coding genes were functionally annotated. Phylogenetic analysis indicated that the diversification of M. ovinus and Glossina fuscipes took place 72.76 Mya within the Late Cretaceous. Gene family expansion and contraction analysis revealed that M. ovinus has 65 rapidly-evolving families (26 expansion and 39 contractions) mainly involved DNA metabolic activity, transposases activity, odorant receptor 59a/67d-like, IMD domain-containing protein, and cuticle protein, etc. The universal and tightly conserved list of milk protein orthologues has been assembled from the genome of M. ovinus. Contractions and losses of sensory receptors and vision-associated Rhodopsin genes were significant in M. ovinus, which indicate that the M. ovinus has narrower ecological niches.

CONCLUSIONS

We sequenced, assembled, and annotated the whole genome sequence of M. ovinus, and launches into the preliminary genetic mechanisms analysis of the adaptive evolution characteristics of M. ovinus. These resources will provide insights to understand the biological underpinnings of this parasite and the disease control strategies.

Warm Blood Meal Increases Digestion Rate and Milk Protein Production to Maximize Reproductive Output for the Tsetse Fly, Glossina morsitans

The ingestion of blood represents a significant burden that immediately increases water, oxidative, and thermal stress, but provides a significant nutrient source to generate resources necessary for the development of progeny. Thermal stress has been assumed to solely be a negative byproduct that has to be alleviated to prevent stress. Here, we examined if the short thermal bouts incurred during a warm blood meal are beneficial to reproduction. To do so, we examined the duration of pregnancy and milk gland protein expression in the tsetse fly, Glossina morsitans, that consumed a warm or cool blood meal. We noted that an optimal temperature for blood ingestion yielded a reduction in the duration of pregnancy. This decline in the duration of pregnancy is due to increased rate of blood digestion when consuming warm blood. This increased digestion likely provided more energy that leads to increased expression of transcript for milk-associated proteins. The shorter duration of pregnancy is predicted to yield an increase in population growth compared to those that consume cool or above host temperatures. These studies provide evidence that consumption of a warm blood meal is likely beneficial for specific aspects of vector biology.

Functional disruption of transferrin expression alters reproductive physiology in Anopheles culicifacies

BACKGROUND

Iron metabolism is crucial to maintain optimal physiological homeostasis of every organism and any alteration of the iron concentration (i.e. deficit or excess) can have adverse consequences. Transferrins are glycoproteins that play important role in iron transportation and have been widely characterized in vertebrates and insects, but poorly studied in blood-feeding mosquitoes.

RESULTS

We characterized a 2102 bp long transcript AcTrf1a with complete CDS of 1872bp, and 226bp UTR region, encoding putative transferrin homolog protein from mosquito An. culicifacies. A detailed in silico analysis predicts AcTrf1a encodes 624 amino acid (aa) long polypeptide that carries transferrin domain. AcTrf1a also showed a putative N-linked glycosylation site, a characteristic feature of most of the mammalian transferrins and certain non-blood feeding insects. Structure modelling prediction confirms the presence of an iron-binding site at the N-terminal lobe of the transferrin. Our spatial and temporal expression analysis under altered pathophysiological conditions showed that AcTrf1a is abundantly expressed in the fat-body, ovary, and its response is significantly altered (enhanced) after blood meal uptake, and exogenous bacterial challenge. Additionally, non-heme iron supplementation of FeCl3 at 1 mM concentration not only augmented the AcTrf1a transcript expression in fat-body but also enhanced the reproductive fecundity of gravid adult female mosquitoes. RNAi-mediated knockdown of AcTrf1a causes a significant reduction in fecundity, confirming the important role of transferrin in oocyte maturation.

CONCLUSION

All together our results advocate that detailed characterization of newly identified AcTrf1a transcript may help to select it as a unique target to impair the mosquito reproductive outcome.

Expression Profiles of Digestive Genes in the Gut and Salivary Glands of Tarnished Plant Bug (Hemiptera: Miridae)

Host plant preference of agricultural pests may shift throughout the growing season, allowing the pests to persist on wild hosts when crops are not available. Lygus Hahn (Hemiptera: Miridae) bugs are severe pests of cotton during flowering and fruiting stages, but can persist on alternative crops, or on weed species. Diversity of digestive enzymes produced by salivary glands and gut tissues play a pivotal role in an organism's ability to utilize various food sources. Polyphagous insects produce an array of enzymes that can process carbohydrates, lipids, and proteins. In this study, the digestive enzyme repertoire of the tarnished plant bug, Lygus lineolaris (Palisot de Beauvois), was identified by high-throughput sequencing followed by cDNA cloning and sequencing. This study identified 87 digestive genes, including 30 polygalacturonases (PG), one β-galactosidase, three α-glucosidases, six β-glucosidases, 28 trypsin-like proteases, three serine proteases, one apyrase-like protease, one cysteine protease, 12 lipases, and two transcripts with low similarity to a xylanase A-like genes. RNA-Seq expression profiles of these digestive genes in adult tarnished plant bugs revealed that 57 and 12 genes were differentially expressed in the salivary gland and gut (≥5-fold, P ≤ 0.01), respectively. All polygalacturonase genes, most proteases, and two xylanase-like genes were differentially expressed in salivary glands, while most of the carbohydrate and lipid processing enzymes were differentially expressed in the gut. Seven of the proteases (KF208689, KF208697, KF208698, KF208699, KF208700, KF208701, and KF208702) were not detected in either the gut or salivary glands.

The genome of the stable fly, Stomoxys calcitrans, reveals potential mechanisms underlying reproduction, host interactions, and novel targets for pest control

Background

The stable fly, Stomoxys calcitrans, is a major blood-feeding pest of livestock that has near worldwide distribution, causing an annual cost of over $2 billion for control and product loss in the USA alone. Control of these flies has been limited to increased sanitary management practices and insecticide application for suppressing larval stages. Few genetic and molecular resources are available to help in developing novel methods for controlling stable flies.

Results

This study examines stable fly biology by utilizing a combination of high-quality genome sequencing and RNA-Seq analyses targeting multiple developmental stages and tissues. In conjunction, 1600 genes were manually curated to characterize genetic features related to stable fly reproduction, vector host interactions, host-microbe dynamics, and putative targets for control. Most notable was characterization of genes associated with reproduction and identification of expanded gene families with functional associations to vision, chemosensation, immunity, and metabolic detoxification pathways.

Conclusions

The combined sequencing, assembly, and curation of the male stable fly genome followed by RNA-Seq and downstream analyses provide insights necessary to understand the biology of this important pest. These resources and new data will provide the groundwork for expanding the tools available to control stable fly infestations. The close relationship of Stomoxys to other blood-feeding (horn flies and Glossina) and non-blood-feeding flies (house flies, medflies, Drosophila) will facilitate understanding of the evolutionary processes associated with development of blood feeding among the Cyclorrhapha.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-00975-9.

Chromosome-level de novo genome assembly of Sarcophaga peregrina provides insights into the evolutionary adaptation of flesh flies

Sarcophaga peregrina is considered to be of great ecological, medical and forensic significance, and has unusual biological characteristics such as an ovoviviparous reproductive pattern and adaptation to feed on carrion. The availability of a high-quality genome will help to further reveal the mechanisms underlying these charcateristics. Here we present a de novo-assembled genome at chromosome scale for S. peregrina. The final assembled genome was 560.31 Mb with contig N50 of 3.84 Mb. Hi-C scaffolding reliably anchored six pseudochromosomes, accounting for 97.76% of the assembled genome. Moreover, 45.70% of repeat elements were identified in the genome. A total of 14,476 protein-coding genes were functionally annotated, accounting for 92.14% of all predicted genes. Phylogenetic analysis indicated that S. peregrina and S. bullata diverged ~ 7.14 million years ago. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid in clarifying its ovoviviparous reproduction and carrion-feeding adaptations, such as lipid metabolism, olfactory receptor activity, antioxidant enzymes, proteolysis and serine-type endopeptidase activity. Protein-coding genes associated with ovoviparity, such as yolk proteins, transferrin and acid sphingomyelinase, were identified. This study provides a valuable genomic resource for S. peregrina, and sheds insight into further revealing the underlying molecular mechanisms of adaptive evolution.

Survey for positively selected coding regions in the genome of the hematophagous tsetse fly Glossina morsitans identifies candidate genes associated with feeding habits and embryonic development

Tsetse flies are responsible for the transmission of Trypanossoma sp. to vertebrate animals in Africa causing huge health issues and economic loss. The availability of the genome sequence of Glossina morsitans enabled the discovery of several genes related to medically important phenotypes and novel physiological features. However, a genome-wide scan for coding regions that underwent positive selection is still missing, which is surprising given the evolution of traits associated with the hematophagy in this lineage. In this study, we employed an experimental design that controlled for the rate of false positives and we performed a scan of 3,318 G. morsitans genes. We found 145 genes with significant historical signal of positive selection. These genes were categorized into 18 functional classes after careful manual annotation. Based on their attributed functions, we identified candidate genes related with feeding habits and embryonic development. When our results were contrasted with gene expression data, we confirmed that most genes that underwent adaptive molecular evolution were frequently expressed in organs associated with key physiological evolutionary innovations in the G. morsitans lineage, namely, the salivary gland, the midgut, fat body tissue, and in the spermatophore.

Molecular mechanisms underlying milk production and viviparity in the cockroach, Diploptera punctata

Viviparous reproduction is characterized by maternal retention of developing offspring within the reproductive tract during gestation, culminating in live birth. In some cases, a mother will provide nutrition beyond that present in the yolk; this is known as matrotrophic viviparity. While this phenomenon is best associated with mammals, it is observed in insects such as the viviparous cockroach, Diploptera punctata. Female D. punctata carry developing embryos in the brood sac, a reproductive organ that acts as both a uterus and a placenta by protecting and providing a nutritive secretion to the intrauterine developing progeny. While the basic physiology of D. punctata pregnancy has been characterized, little is known about the molecular mechanisms underlying this phenomenon. This study combined RNA-seq analysis, RNA interference, and other assays to characterize molecular and physiological changes associated with D. punctata reproduction. A comparison of four stages of the female reproductive cycle and males revealed unique gene expression profiles corresponding to each stage and between sexes. Differentially regulated transcripts of interest include the previously identified family of milk proteins and transcripts associated with juvenile hormone metabolism. RNA interference and methoprene application experiments established the potential impacts of bothbreakdown and synthesis reduction of juvenile hormone in maintaining pregnancy in D. punctata. These studies provide the comprehensive molecular mechanisms associated with cockroach viviparity, which will be a critical resource for comparative purposes among viviparity in insect systems.

Species-specific transcriptional profiles of the gut and gut microbiome of Ceratitis quilicii and Ceratitis rosa sensu stricto

The fruit fly species, Ceratitis rosa sensu stricto and Ceratitis quilicii, are sibling species restricted to the lowland and highland regions, respectively. Until recently, these sibling species were considered as allopatric populations of C. rosa with distinct bionomics. We used deep Next Generation Sequencing (NGS) technology on intact guts of individuals from the two sibling species to compare their transcriptional profiles and simultaneously understand gut microbiome and host molecular processes and identify distinguishing genetic differences between the two species. Since the genomes of both species had not been published previously, the transcriptomes were assembled de novo into transcripts. Microbe-specific transcript orthologs were separated from the assembly by filtering searches of the transcripts against microbe databases using OrthoMCL. We then used differential expression analysis of host-specific transcripts (i.e. those remaining after the microbe-specific transcripts had been removed) and microbe-specific transcripts from the two-sibling species to identify defining species-specific transcripts that were present in only one fruit fly species or the other, but not in both. In C. quilicii females, bacterial transcripts of Pectobacterium spp., Enterobacterium buttiauxella, Enterobacter cloacae and Klebsiella variicola were upregulated compared to the C. rosa s.s. females. Comparison of expression levels of the host transcripts revealed a heavier investment by C. quilicii (compared with C. rosa s.s.) in: immunity; energy production; cell proliferation; insecticide resistance; reproduction and proliferation; and redox reactions that are usually associated with responses to stress and degradation of fruit metabolites.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes

BACKGROUND

Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity.

RESULTS

Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges.

CONCLUSIONS

Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.

Comparative genomic analysis of six Glossina genomes, vectors of African trypanosomes

BACKGROUND

Tsetse flies (Glossina sp.) are the vectors of human and animal trypanosomiasis throughout sub-Saharan Africa. Tsetse flies are distinguished from other Diptera by unique adaptations, including lactation and the birthing of live young (obligate viviparity), a vertebrate blood-specific diet by both sexes, and obligate bacterial symbiosis. This work describes the comparative analysis of six Glossina genomes representing three sub-genera: Morsitans (G. morsitans morsitans, G. pallidipes, G. austeni), Palpalis (G. palpalis, G. fuscipes), and Fusca (G. brevipalpis) which represent different habitats, host preferences, and vectorial capacity.

RESULTS

Genomic analyses validate established evolutionary relationships and sub-genera. Syntenic analysis of Glossina relative to Drosophila melanogaster shows reduced structural conservation across the sex-linked X chromosome. Sex-linked scaffolds show increased rates of female-specific gene expression and lower evolutionary rates relative to autosome associated genes. Tsetse-specific genes are enriched in protease, odorant-binding, and helicase activities. Lactation-associated genes are conserved across all Glossina species while male seminal proteins are rapidly evolving. Olfactory and gustatory genes are reduced across the genus relative to other insects. Vision-associated Rhodopsin genes show conservation of motion detection/tracking functions and variance in the Rhodopsin detecting colors in the blue wavelength ranges.

CONCLUSIONS

Expanded genomic discoveries reveal the genetics underlying Glossina biology and provide a rich body of knowledge for basic science and disease control. They also provide insight into the evolutionary biology underlying novel adaptations and are relevant to applied aspects of vector control such as trap design and discovery of novel pest and disease control strategies.

Gene expression in reproductive organs of tsetse females - initial data in an approach to reduce fecundity

BACKGROUND

Tsetse flies are vectors of African trypanosomes, and their vectorial capacity results in a major public health emergency and vast economic losses in sub-Saharan Africa. Given the limited ability of trypanosome prevention and eradication, tsetse vectors remain major targets of control efforts. Larvae of all three instars are developed in mothers' uteri, nourished through milk, and 'larviposited' shortly before pupation. The past few years have witnessed the emergence of approaches based on knockdown of genes involved in milk production, resulting in a significant reduction of fecundity.

RESULTS

In order to identify further genes applicable in the control of tsetse flies, we determined the expression of protein-coding genes in ovaries and uteri from both virgin and heavily pregnant Glossina morsitans morsitans females. Comparison of expression profiles allowed us to identify candidate genes with increased expression in pregnant individuals. Lists with the highest increases include genes involved in oocyte and embryonic development, or nourishment. Maximum ovarian fold change does not exceed 700, while the highest uterine fold change reaches to more than 4000. Relatively high fold changes of two neuropeptide receptors (for corazonin and myosuppressin) propose the corresponding genes alternative targets.

CONCLUSIONS

Given the higher fold changes in the uterus, targeting gene expression in this tissue may result in a more evident reduction of fecundity. However, ovaries should not be neglected, as manifested by several genes with top fold changes involved in early developmental stages. Apart from focusing on the highest fold changes, neuropeptide receptors with moderate increases in expression should be also verified as targets, given their roles in mediating the tissue control. However, this data needs to be considered initial, and the potential of these genes in affecting female fecundity needs to be verified experimentally.

Enhancing vector refractoriness to trypanosome infection: achievements, challenges and perspectives

With the absence of effective prophylactic vaccines and drugs against African trypanosomosis, control of this group of zoonotic neglected tropical diseases depends the control of the tsetse fly vector. When applied in an area-wide insect pest management approach, the sterile insect technique (SIT) is effective in eliminating single tsetse species from isolated populations. The need to enhance the effectiveness of SIT led to the concept of investigating tsetse-trypanosome interactions by a consortium of researchers in a five-year (2013-2018) Coordinated Research Project (CRP) organized by the Joint Division of FAO/IAEA. The goal of this CRP was to elucidate tsetse-symbiome-pathogen molecular interactions to improve SIT and SIT-compatible interventions for trypanosomoses control by enhancing vector refractoriness. This would allow extension of SIT into areas with potential disease transmission. This paper highlights the CRP's major achievements and discusses the science-based perspectives for successful mitigation or eradication of African trypanosomosis.

Dehydration prompts increased activity and blood feeding by mosquitoes

Current insights into the mosquito dehydration response rely on studies that examine specific responses but ultimately fail to provide an encompassing view of mosquito biology. Here, we examined underlying changes in the biology of mosquitoes associated with dehydration. Specifically, we show that dehydration increases blood feeding in the northern house mosquito, Culex pipiens, which was the result of both higher activity and a greater tendency to land on a host. Similar observations were noted for Aedes aegypti and Anopheles quadrimaculatus. RNA-seq and metabolome analyses in C. pipiens following dehydration revealed that factors associated with carbohydrate metabolism are altered, specifically the breakdown of trehalose. Suppression of trehalose breakdown in C. pipiens by RNA interference reduced phenotypes associated with lower hydration levels. Lastly, mesocosm studies for C. pipiens confirmed that dehydrated mosquitoes were more likely to host feed under ecologically relevant conditions. Disease modeling indicates dehydration bouts will likely enhance viral transmission. This dehydration-induced increase in blood feeding is therefore likely to occur regularly and intensify during periods when availability of water is low.

Rapid autophagic regression of the milk gland during involution is critical for maximizing tsetse viviparous reproductive output

Tsetse flies are important vectors of human and animal trypanosomiasis. Ability to reduce tsetse populations is an effective means of disease control. Lactation is an essential component of tsetse’s viviparous reproductive physiology and requires a dramatic increase in the expression and synthesis of milk proteins by the milk gland organ in order to nurture larval growth. In between each gonotrophic cycle, tsetse ceases milk production and milk gland tubules undergo a nearly two-fold reduction in width (involution). In this study, we examined the role autophagy plays during tsetse fly milk gland involution and reproductive output. Autophagy genes show elevated expression in tissues associated with lactation, immediately before or within two hours post-parturition, and decline at 24-48h post-parturition. This expression pattern is inversely correlated with that of the milk gland proteins (lactation-specific protein coding genes) and the autophagy inhibitor fk506-bp1. Increased expression of Drosophila inhibitor of apoptosis 1, diap1, was also observed in the milk gland during involution, when it likely prevents apoptosis of milk gland cells. RNAi-mediated knockdown of autophagy related gene 8a (atg8a) prevented rapid milk gland autophagy during involution, prolonging gestation, and reducing fecundity in the subsequent gonotrophic cycle. The resultant inhibition of autophagy reduced the recovery of stored lipids during the dry (non-lactating) periods by 15–20%. Ecdysone application, similar to levels that occur immediately before birth, induced autophagy, and increased milk gland involution even before abortion. This suggests that the ecdysteroid peak immediately preceding parturition likely triggers milk gland autophagy. Population modeling reveals that a delay in involution would yield a negative population growth rate. This study indicates that milk gland autophagy during involution is critical to restore nutrient reserves and allow efficient transition between pregnancy cycles. Targeting post-birth phases of reproduction could be utilized as a novel mechanism to suppress tsetse populations and reduce trypanosomiasis.

Tsetse flies are vectors for trypanosomes that cause both African sleeping sickness in humans and Nagana in animals. The reduction of tsetse populations is the most efficient way to reduce the prevalence of this economically important disease with current control methods including pesticide application, traps, and sterile insect techniques. Tsetse pregnancy and milk production represent a species-specific target for population control and milk gland transition during each larval growth cycle could represent a novel target for tsetse control. Within one day after birth, the milk gland organ, essential for provisioning nutrients to the intrauterine larva, undergoes involution marked by an ecdysone driven increase in autophagy that allows breakdown of this gland. Inhibiting the process of autophagy prevents the timely transition from the lactation phase to the dry phase, triggering a delay in subsequent pregnancy cycle. This misregulation of milk gland involution leads to an overall decrease in the number of offspring that each female can produce per lifetime. This study has determined the molecular components of this process, and reveals new targets of interference for vector control.

Symbiont-induced odorant binding proteins mediate insect host hematopoiesis

Symbiotic bacteria assist in maintaining homeostasis of the animal immune system. However, the molecular mechanisms that underlie symbiont-mediated host immunity are largely unknown. Tsetse flies (Glossina spp.) house maternally transmitted symbionts that regulate the development and function of their host's immune system. Herein we demonstrate that the obligate mutualist, Wigglesworthia, up-regulates expression of odorant binding protein six in the gut of intrauterine tsetse larvae. This process is necessary and sufficient to induce systemic expression of the hematopoietic RUNX transcription factor lozenge and the subsequent production of crystal cells, which actuate the melanotic immune response in adult tsetse. Larval Drosophila's indigenous microbiota, which is acquired from the environment, regulates an orthologous hematopoietic pathway in their host. These findings provide insight into the molecular mechanisms that underlie enteric symbiont-stimulated systemic immune system development, and indicate that these processes are evolutionarily conserved despite the divergent nature of host-symbiont interactions in these model systems.

Disease vectors in the era of next generation sequencing

Almost 20 % of all infectious human diseases are vector borne and, together, are responsible for over one million deaths per annum. Over the past decade, the decreasing costs of massively parallel sequencing technologies have facilitated the agnostic interrogation of insect vector genomes, giving medical entomologists access to an ever-expanding volume of high-quality genomic and transcriptomic data. In this review, we highlight how genomics resources have provided new insights into the physiology, behavior, and evolution of human disease vectors within the context of the global health landscape.

Duplication and Diversification of Dipteran Argonaute Genes, and the Evolutionary Divergence of Piwi and Aubergine

Genetic studies of Drosophila melanogaster have provided a paradigm for RNA interference (RNAi) in arthropods, in which the microRNA and antiviral pathways are each mediated by a single Argonaute (Ago1 and Ago2) and germline suppression of transposable elements is mediated by a trio of Piwi-subfamily Argonaute proteins (Ago3, Aub, and Piwi). Without a suitable evolutionary context, deviations from this can be interpreted as derived or idiosyncratic. Here we analyze the evolution of Argonaute genes across the genomes and transcriptomes of 86 Dipteran species, showing that variation in copy number can occur rapidly, and that there is constant flux in some RNAi mechanisms. The lability of the RNAi pathways is illustrated by the divergence of Aub and Piwi (182-156 Ma), independent origins of multiple Piwi-family genes in Aedes mosquitoes (less than 25Ma), and the recent duplications of Ago2 and Ago3 in the tsetse fly Glossina morsitans. In each case the tissue specificity of these genes has altered, suggesting functional divergence or innovation, and consistent with the action of dynamic selection pressures across the Argonaute gene family. We find there are large differences in evolutionary rates and gene turnover between pathways, and that paralogs of Ago2, Ago3, and Piwi/Aub show contrasting rates of evolution after duplication. This suggests that Argonautes undergo frequent evolutionary expansions that facilitate functional divergence.

The Spermatophore in Glossina morsitans morsitans: Insights into Male Contributions to Reproduction

Male Seminal Fluid Proteins (SFPs) transferred during copulation modulate female reproductive physiology and behavior, impacting sperm storage/use, ovulation, oviposition, and remating receptivity. These capabilities make them ideal targets for developing novel methods of insect disease vector control. Little is known about the nature of SFPs in the viviparous tsetse flies (Diptera: Glossinidae), vectors of Human and Animal African trypanosomiasis. In tsetse, male ejaculate is assembled into a capsule-like spermatophore structure visible post-copulation in the female uterus. We applied high-throughput approaches to uncover the composition of the spermatophore in Glossina morsitans morsitans. We found that both male accessory glands and testes contribute to its formation. The male accessory glands produce a small number of abundant novel proteins with yet unknown functions, in addition to enzyme inhibitors and peptidase regulators. The testes contribute sperm in addition to a diverse array of less abundant proteins associated with binding, oxidoreductase/transferase activities, cytoskeletal and lipid/carbohydrate transporter functions. Proteins encoded by female-biased genes are also found in the spermatophore. About half of the proteins display sequence conservation relative to other Diptera, and low similarity to SFPs from other studied species, possibly reflecting both their fast evolutionary pace and the divergent nature of tsetse's viviparous biology.

A comparative analysis of reproductive biology of insect vectors of human disease

Studying the reproductive strategies of insect species that transmit diseases to humans can identify new exploitable targets for the development of vector control methods. Here we describe shared characteristics and individual features of the reproductive biology of three major disease vectors: Anopheles gambiae, Aedes aegypti and Glossina morsitans. Current studies are identifying i) species-specific molecular cascades that determine female monandrous behavior, ii) core aspects of egg development that could be disrupted for controlling natural populations, and iii) the increasingly apparent role of resident microbiota in shaping reproductive success and disease transmission potential. The recent completion of multiple genome sequencing projects is allowing comparative genomics studies that not only increase our knowledge of reproductive processes but also facilitate the identification of novel targets for vector control.

Adenotrophic viviparity in tsetse flies: potential for population control and as an insect model for lactation

Tsetse flies (Glossina spp.), vectors of African trypanosomes, are distinguished by their specialized reproductive biology, defined by adenotrophic viviparity (maternal nourishment of progeny by glandular secretions followed by live birth). This trait has evolved infrequently among insects and requires unique reproductive mechanisms. A key event in Glossina reproduction involves the transition between periods of lactation and nonlactation (dry periods). Increased lipolysis, nutrient transfer to the milk gland, and milk-specific protein production characterize lactation, which terminates at the birth of the progeny and is followed by a period of involution. The dry stage coincides with embryogenesis of the progeny, during which lipid reserves accumulate in preparation for the next round of lactation. The obligate bacterial symbiont Wigglesworthia glossinidia is critical to tsetse reproduction and likely provides B vitamins required for metabolic processes underlying lactation and/or progeny development. Here we describe findings that utilized transcriptomics, physiological assays, and RNA interference-based functional analysis to understand different components of adenotrophic viviparity in tsetse flies.

Trypanosome Transmission Dynamics in Tsetse

Tsetse flies (Diptera:Glossinidae) are vectors of African trypanosomes. Tsetse undergo viviparous reproductive biology, and depend on their obligate endosymbiont (genus Wigglesworthia) for the maintenance of fecundity and immune system development. Trypanosomes establish infections in the midgut and salivary glands of the fly. Tsetse's resistance to trypanosome infection increases as a function of age. Among the factors that mediate resistance to parasites are antimicrobial peptides (AMPs) produced by the Immune deficiency (Imd) signaling pathway, peptidoglycan recognition protein (PGRP) LB, tsetse-EP protein and the integrity of the midgut peritrophic matrix (PM) barrier. The presence of obligate Wigglesworthia during larval development is essential for adult immune system maturation and PM development. Thus, Wigglesworthia prominently influences the vector competency of it's tsetse host.

Vitamin B6 generated by obligate symbionts is critical for maintaining proline homeostasis and fecundity in tsetse flies

The viviparous tsetse fly utilizes proline as a hemolymph-borne energy source. In tsetse, biosynthesis of proline from alanine involves the enzyme alanine-glyoxylate aminotransferase (AGAT), which requires pyridoxal phosphate (vitamin B6) as a cofactor. This vitamin can be synthesized by tsetse's obligate symbiont, Wigglesworthia glossinidia. In this study, we examined the role of Wigglesworthia-produced vitamin B6 for maintenance of proline homeostasis, specifically during the energetically expensive lactation period of the tsetse's reproductive cycle. We found that expression of agat, as well as genes involved in vitamin B6 metabolism in both host and symbiont, increases in lactating flies. Removal of symbionts via antibiotic treatment of flies (aposymbiotic) led to hypoprolinemia, reduced levels of vitamin B6 in lactating females, and decreased fecundity. Proline homeostasis and fecundity recovered partially when aposymbiotic tsetse were fed a diet supplemented with either yeast or Wigglesworthia extracts. RNA interference-mediated knockdown of agat in wild-type flies reduced hemolymph proline levels to that of aposymbiotic females. Aposymbiotic flies treated with agat short interfering RNA (siRNA) remained hypoprolinemic even upon dietary supplementation with microbial extracts or B vitamins. Flies infected with parasitic African trypanosomes display lower hemolymph proline levels, suggesting that the reduced fecundity observed in parasitized flies could result from parasite interference with proline homeostasis. This interference could be manifested by competition between tsetse and trypanosomes for vitamins, proline, or other factors involved in their synthesis. Collectively, these results indicate that the presence of Wigglesworthia in tsetse is critical for the maintenance of proline homeostasis through vitamin B6 production.

Genome sequence of the tsetse fly (Glossina morsitans): vector of African trypanosomiasis

Tsetse flies are the sole vectors of human African trypanosomiasis throughout sub-Saharan Africa. Both sexes of adult tsetse feed exclusively on blood and contribute to disease transmission. Notable differences between tsetse and other disease vectors include obligate microbial symbioses, viviparous reproduction, and lactation. Here, we describe the sequence and annotation of the 366-megabase Glossina morsitans morsitans genome. Analysis of the genome and the 12,308 predicted protein-encoding genes led to multiple discoveries, including chromosomal integrations of bacterial (Wolbachia) genome sequences, a family of lactation-specific proteins, reduced complement of host pathogen recognition proteins, and reduced olfaction/chemosensory associated genes. These genome data provide a foundation for research into trypanosomiasis prevention and yield important insights with broad implications for multiple aspects of tsetse biology.

Presence of extensive Wolbachia symbiont insertions discovered in the genome of its host Glossina morsitans morsitans

Tsetse flies (Glossina spp.) are the cyclical vectors of Trypanosoma spp., which are unicellular parasites responsible for multiple diseases, including nagana in livestock and sleeping sickness in humans in Africa. Glossina species, including Glossina morsitans morsitans (Gmm), for which the Whole Genome Sequence (WGS) is now available, have established symbiotic associations with three endosymbionts: Wigglesworthia glossinidia, Sodalis glossinidius and Wolbachia pipientis (Wolbachia). The presence of Wolbachia in both natural and laboratory populations of Glossina species, including the presence of horizontal gene transfer (HGT) events in a laboratory colony of Gmm, has already been shown. We herein report on the draft genome sequence of the cytoplasmic Wolbachia endosymbiont (cytWol) associated with Gmm. By in silico and molecular and cytogenetic analysis, we discovered and validated the presence of multiple insertions of Wolbachia (chrWol) in the host Gmm genome. We identified at least two large insertions of chrWol, 527,507 and 484,123 bp in size, from Gmm WGS data. Southern hybridizations confirmed the presence of Wolbachia insertions in Gmm genome, and FISH revealed multiple insertions located on the two sex chromosomes (X and Y), as well as on the supernumerary B-chromosomes. We compare the chrWol insertions to the cytWol draft genome in an attempt to clarify the evolutionary history of the HGT events. We discuss our findings in light of the evolution of Wolbachia infections in the tsetse fly and their potential impacts on the control of tsetse populations and trypanosomiasis.

African trypanosomes are transmitted to man and animals by tsetse fly, a blood sucking insect. Tsetse flies include all Glossina species with the genome of Glossina morsitans morsitans (Gmm) being sequenced under the International Glossina Genome Initiative. The endosymbionts Wigglesworthia glossinidia, Sodalis glossinidius and Wolbachia pipientis (Wolbachia) have been found to establish symbiotic associations with Gmm. Wolbachia is known to be present in natural and laboratory populations of Glossina species. In this study we report the genome sequence of the Wolbachia strain that is associated with Gmm. With the aid of in silico and molecular and cytogenetic analyses, multiple insertions of the Wolbachia genome were revealed and confirmed in Gmm chromosome. Comparison of the cytoplasmic Wolbachia draft genome and the chromosomal insertions enabled us to infer the evolutionary history of the Wolbachia horizontal transfer events. These findings are discussed in relation to their impact on the development of Wolbachia-based strategies for the control of tsetse flies and trypanosomiasis.

Insights into the trypanosome-host interactions revealed through transcriptomic analysis of parasitized tsetse fly salivary glands

The agents of sleeping sickness disease, Trypanosoma brucei complex parasites, are transmitted to mammalian hosts through the bite of an infected tsetse. Information on tsetse-trypanosome interactions in the salivary gland (SG) tissue, and on mammalian infective metacyclic (MC) parasites present in the SG, is sparse. We performed RNA-seq analyses from uninfected and T. b. brucei infected SGs of Glossina morsitans morsitans. Comparison of the SG transcriptomes to a whole body fly transcriptome revealed that only 2.7% of the contigs are differentially expressed during SG infection, and that only 263 contigs (0.6%) are preferentially expressed in the SGs (SG-enriched). The expression of only 37 contigs (0.08%) and 27 SG-enriched contigs (10%) were suppressed in infected SG. These suppressed contigs accounted for over 55% of the SG transcriptome, and included the most abundant putative secreted proteins with anti-hemostatic functions present in saliva. In contrast, expression of putative host proteins associated with immunity, stress, cell division and tissue remodeling were enriched in infected SG suggesting that parasite infections induce host immune and stress response(s) that likely results in tissue renewal. We also performed RNA-seq analysis from mouse blood infected with the same parasite strain, and compared the transcriptome of bloodstream form (BSF) cells with that of parasites obtained from the infected SG. Over 30% of parasite transcripts are differentially regulated between the two stages, and reflect parasite adaptations to varying host nutritional and immune ecology. These differences are associated with the switch from an amino acid based metabolism in the SG to one based on glucose utilization in the blood, and with surface coat modifications that enable parasite survival in the different hosts. This study provides a foundation on the molecular aspects of the trypanosome dialogue with its tsetse and mammalian hosts, necessary for future functional investigations.

Tsetse flies transmit the causative agents of African sleeping sickness and nagana in sub-Saharan Africa. The parasites are acquired when tsetse flies feed on an infected host, undergo multiplication in the fly gut and migrate to the salivary glands (SG). The cycle resumes once this infected fly transmits the parasites in conjunction with saliva to another host when feeding. We compared gene expression changes between parasitized and uninfected tsetse SG. We also assessed changes in parasite gene expression in the tsetse SG in relation to those present within vertebrate blood. We found that parasite infections increase expression of host proteins associated with stress and cell division, indicative of extensive cellular damage in SG. We also found that parasite infections reduce expression of the most highly expressed SG-specific secreted proteins, suggesting modification of saliva composition. The parasite transcriptome reveals changes in specific cell surface proteins and in metabolism related to glucose-amino acid utilization in the different host environments. This study provides information for critical understanding of tsetse-trypanosome interactions, and transcriptional changes that likely enable the parasite to persist in the varying environment of its insect and vertebrate hosts.

Amelioration of reproduction-associated oxidative stress in a viviparous insect is critical to prevent reproductive senescence

Impact of reproductive processes upon female health has yielded conflicting results; particularly in relation to the role of reproduction-associated stress. We used the viviparous tsetse fly to determine if lactation, birth and involution lead to damage from oxidative stress (OS) that impairs subsequent reproductive cycles. Tsetse females carry an intrauterine larva to full term at each pregnancy cycle, and lactate to nourish them with milk secretions produced by the accessory gland ( = milk gland) organ. Unlike most K-strategists, tsetse females lack an apparent period of reproductive senescence allowing the production of 8-10 progeny over their entire life span. In a lactating female, over 47% of the maternal transcriptome is associated with the generation of milk proteins. The resulting single larval offspring weighs as much as the mother at birth. In studying this process we noted an increase in specific antioxidant enzyme (AOE) transcripts and enzymatic activity at critical times during lactation, birth and involution in the milk gland/fat body organ and the uterus. Suppression of superoxide dismutase (sod) decreased fecundity in subsequent reproductive cycles in young mothers and nearly abolished fecundity in geriatric females. Loss of fecundity was in part due to the inability of the mother to produce adequate milk to support larval growth. Longevity was also impaired after sod knockdown. Generation of OS in virgin females through exogenous treatment with hydrogen peroxide at times corresponding to pregnancy intervals reduced survival, which was exacerbated by sod knockdown. AOE expression may prevent oxidative damage associated with the generation of nutrients by the milk gland, parturition and milk gland breakdown. Our results indicate that prevention of OS is essential for females to meet the growing nutritional demands of juveniles during pregnancy and to repair the damage that occurs at birth. This process is particularly important for females to remain fecund during the latter portion of their lifetime.

Aquaporins are critical for provision of water during lactation and intrauterine progeny hydration to maintain tsetse fly reproductive success

Tsetse flies undergo drastic fluctuations in their water content throughout their adult life history due to events such as blood feeding, dehydration and lactation, an essential feature of the viviparous reproductive biology of tsetse. Aquaporins (AQPs) are transmembrane proteins that allow water and other solutes to permeate through cellular membranes. Here we identify tsetse aquaporin (AQP) genes, examine their expression patterns under different physiological conditions (blood feeding, lactation and stress response) and perform functional analysis of three specific genes utilizing RNA interference (RNAi) gene silencing. Ten putative aquaporins were identified in the Glossina morsitans morsitans (Gmm) genome, two more than has been previously documented in any other insect. All organs, tissues, and body parts examined had distinct AQP expression patterns. Two AQP genes, gmmdripa and gmmdripb ( = gmmaqp1a and gmmaqp1b) are highly expressed in the milk gland/fat body tissues. The whole-body transcript levels of these two genes vary over the course of pregnancy. A set of three AQPs (gmmaqp5, gmmaqp2a, and gmmaqp4b) are expressed highly in the Malpighian tubules. Knockdown of gmmdripa and gmmdripb reduced the efficiency of water loss following a blood meal, increased dehydration tolerance and reduced heat tolerance of adult females. Knockdown of gmmdripa extended pregnancy length, and gmmdripb knockdown resulted in extended pregnancy duration and reduced progeny production. We found that knockdown of AQPs increased tsetse milk osmolality and reduced the water content in developing larva. Combined knockdown of gmmdripa, gmmdripb and gmmaqp5 extended pregnancy by 4–6 d, reduced pupal production by nearly 50%, increased milk osmolality by 20–25% and led to dehydration of feeding larvae. Based on these results, we conclude that gmmDripA and gmmDripB are critical for diuresis, stress tolerance and intrauterine lactation through the regulation of water and/or other uncharged solutes.

Glossina sp. are responsible for transmission of African trypanosomes, the causative agents of sleeping sickness in humans and Nagana in cattle. Blood feeding and nutrient provisioning through lactation during intrauterine progeny development are periods when considerable water movement occurs within tsetse flies. With the completion of the tsetse fly genome, we sought to characterize the role of aquaporins in relation water homeostasis during blood feeding, stress tolerance and the lactation cycle. We provide evidence that specific AQPs are 1. critical during diuresis following a bloodmeal, 2. important in the regulation of dehydration resistance and heat tolerance and 3. crucial in the allocation of water within tsetse milk that is necessary for progeny hydration. Specifically, we discovered a novel tsetse AQP that is imperative to lactation and may represent a potential target for population control of this disease vector.

The homeodomain protein ladybird late regulates synthesis of milk proteins during pregnancy in the tsetse fly (Glossina morsitans)

Regulation of tissue and development specific gene expression patterns underlies the functional specialization of organs in multi-cellular organisms. In the viviparous tsetse fly (Glossina), the female accessory gland is specialized to generate nutrients in the form of a milk-like secretion to support growth of intrauterine larva. Multiple milk protein genes are expressed specifically in the female accessory gland and are tightly linked with larval development. Disruption of milk protein synthesis deprives developing larvae of nutrients and results in extended larval development and/or in abortion. The ability to cause such a disruption could be utilized as a tsetse control strategy. Here we identify and delineate the regulatory sequence of a major milk protein gene (milk gland protein 1:mgp1) by utilizing a combination of molecular techniques in tsetse, Drosophila transgenics, transcriptomics and in silico sequence analyses. The function of this promoter is conserved between tsetse and Drosophila. In transgenic Drosophila the mgp1 promoter directs reporter gene expression in a tissue and stage specific manner orthologous to that of Glossina. Analysis of the minimal required regulatory region of mgp1, and the regulatory regions of other Glossina milk proteins identified putative homeodomain protein binding sites as the sole common feature. Annotation and expression analysis of Glossina homeodomain proteins identified ladybird late (lbl) as being accessory gland/fat body specific and differentially expressed between lactating/non-lactating flies. Knockdown of lbl in tsetse resulted in a significant reduction in transcript abundance of multiple milk protein genes and in a significant loss of fecundity. The role of Lbl in adult reproductive physiology is previously unknown. These results suggest that Lbl is part of a conserved reproductive regulatory system that could have implications beyond tsetse to other vector insects such as mosquitoes. This system is critical for tsetse fecundity and provides a potential target for development of a reproductive inhibitor.

Female tsetse flies (Diptera: Glossina) harbor and give birth to live young. To do this, they nourish their intrauterine larvae with milk secretions. This work focuses upon understanding the regulation of tsetse milk proteins, which are essential for fecundity and are expressed in a temporally and spatially specific manner by pregnant females. We identified the minimal upstream regulatory DNA sequence of the major milk protein gene mgp1, which confers tissue specific expression in the female accessory glands of reproductively active flies. This regulatory sequence functions similarly in transgenic fruit flies (Drosophila melanogaster) and drives expression of reporter gene products in the adult female accessory gland. Comparison of this regulatory sequence with sequences from other characterized milk proteins indicates that conserved homeodomain transcription factors may be responsible for regulating these genes. Analysis of Glossina homeodomain proteins identified an accessory gland/fat body specific factor, Ladybird late (lbl), which appears to regulate the expression of multiple milk proteins. Reduction of lbl levels interferes with milk protein gene expression, which in turn reduces Glossina fecundity. These results suggest that milk proteins in Glossina are regulated by a conserved regulatory system mediated in part by the homeodomain transcription factor lbl. Components of this system could provide a target for development of a tsetse reproductive inhibitor.

Odorant and gustatory receptors in the tsetse fly Glossina morsitans morsitans

Tsetse flies use olfactory and gustatory responses, through odorant and gustatory receptors (ORs and GRs), to interact with their environment. Glossina morsitans morsitans genome ORs and GRs were annotated using homologs of these genes in Drosophila melanogaster and an ab initio approach based on OR and GR specific motifs in G. m. morsitans gene models coupled to gene ontology (GO). Phylogenetic relationships among the ORs or GRs and the homologs were determined using Maximum Likelihood estimates. Relative expression levels among the G. m. morsitans ORs or GRs were established using RNA-seq data derived from adult female fly. Overall, 46 and 14 putative G. m. morsitans ORs and GRs respectively were recovered. These were reduced by 12 and 59 ORs and GRs respectively compared to D. melanogaster. Six of the ORs were homologous to a single D. melanogaster OR (DmOr67d) associated with mating deterrence in females. Sweet taste GRs, present in all the other Diptera, were not recovered in G. m. morsitans. The GRs associated with detection of CO₂ were conserved in G. m. morsitans relative to D. melanogaster. RNA-sequence data analysis revealed expression of GmmOR15 locus represented over 90% of expression profiles for the ORs. The G. m. morsitans ORs or GRs were phylogenetically closer to those in D. melanogaster than to other insects assessed. We found the chemoreceptor repertoire in G. m. morsitans smaller than other Diptera, and we postulate that this may be related to the restricted diet of blood-meal for both sexes of tsetse flies. However, the clade of some specific receptors has been expanded, indicative of their potential importance in chemoreception in the tsetse.

Tsetse flies navigate their environments using chemosensory receptors, which permit them to locate hosts, mating partners, and resting and larviposition sites. The genome of G. m. morsitans was interrogated for coding genes of odorant receptors (ORs) and gustatory receptors (GRs) that express in antennae and maxillary palp, and detect the volatile and soluble chemical signals. The signals are then transmitted to the central nervous system and translated to phenotypes. Majority of these genes in G. m. morsitans were spread across different scaffolds, but a few were found to occur in clusters, which suggested possible co-regulation of their expression. The number of ORs and GRs were much reduced in the G. m. morsitans genome, including the apparent loss of receptors for sugar when compared to selected Diptera. There was also an apparent numerical expansion of some receptors, presumably to maximize on their restricted blood-meal diet. The annotation of the chemoreceptor package of G. m. morsitans provides a resource for investigating key activities of tsetse flies that could be exploited for their control.