The First Mitochondrial Genome of the Sepsid Fly Nemopoda mamaevi Ozerov, 1997 (Diptera: Sciomyzoidea: Sepsidae), with Mitochondrial Genome Phylogeny of Cyclorrhapha

Sepsid flies (Diptera: Sepsidae) are important model insects for sexual selection research. In order to develop mitochondrial (mt) genome data for this significant group, we sequenced the first complete mt genome of the sepsid fly Nemopoda mamaevi Ozerov, 1997. The circular 15,878 bp mt genome is typical of Diptera, containing all 37 genes usually present in bilaterian animals. We discovered inaccurate annotations of fly mt genomes previously deposited on GenBank and thus re-annotated all published mt genomes of Cyclorrhapha. These re-annotations were based on comparative analysis of homologous genes, and provide a statistical analysis of start and stop codon positions. We further detected two 18 bp of conserved intergenic sequences from tRNA^Glu-tRNA^Phe and ND1-tRNA^Ser(UCN) across Cyclorrhapha, which are the mtTERM binding site motifs. Additionally, we compared automated annotation software MITOS with hand annotation method. Phylogenetic trees based on the mt genome data from Cyclorrhapha were inferred by Maximum-likelihood and Bayesian methods, strongly supported a close relationship between Sepsidae and the Tephritoidea.     

A benchmark and an algorithm for detecting germline transposon insertions and measuring de novo transposon insertion frequencies

Transposons are genomic parasites, and their new insertions can cause instability and spur the evolution of their host genomes. Rapid accumulation of short-read whole-genome sequencing data provides a great opportunity for studying new transposon insertions and their impacts on the host genome. Although many algorithms are available for detecting transposon insertions, the task remains challenging and existing tools are not designed for identifying de novo insertions. Here, we present a new benchmark fly dataset based on PacBio long-read sequencing and a new method TEMP2 for detecting germline insertions and measuring de novo ‘singleton’ insertion frequencies in eukaryotic genomes. TEMP2 achieves high sensitivity and precision for detecting germline insertions when compared with existing tools using both simulated data in fly and experimental data in fly and human. Furthermore, TEMP2 can accurately assess the frequencies of de novo transposon insertions even with high levels of chimeric reads in simulated datasets; such chimeric reads often occur during the construction of short-read sequencing libraries. By applying TEMP2 to published data on hybrid dysgenic flies inflicted by de-repressed P-elements, we confirmed the continuous new insertions of P-elements in dysgenic offspring before they regain piRNAs for P-element repression. TEMP2 is freely available at Github: https://github.com/weng-lab/TEMP2.     

Interspecific Transfer of Bacterial Endosymbionts between Tsetse Fly Species: Infection Establishment and Effect on Host Fitness

Tsetse flies (Glossina spp.) can harbor up to three distinct species of endosymbiotic bacteria that exhibit unique modes of transmission and evolutionary histories with their host. Two mutualist enterics, Wigglesworthia and Sodalis, are transmitted maternally to tsetse flies' intrauterine larvae. The third symbiont, from the genus Wolbachia, parasitizes developing oocytes. In this study, we determined that Sodalis isolates from several tsetse fly species are virtually identical based on a phylogenetic analysis of their ftsZ gene sequences. Furthermore, restriction fragment-length polymorphism analysis revealed little variation in the genomes of Sodalis isolates from tsetse fly species within different subgenera (Glossina fuscipes fuscipes and Glossina morsitans morsitans). We also examined the impact on host fitness of transinfecting G. fuscipes fuscipes and G. morsitans morsitans flies with reciprocal Sodalis strains. Tsetse flies cleared of their native Sodalis symbionts were successfully repopulated with the Sodalis species isolated from a different tsetse fly species. These transinfected flies effectively transmitted the novel symbionts to their offspring and experienced no detrimental fitness effects compared to their wild-type counterparts, as measured by longevity and fecundity. Quantitative PCR analysis revealed that transinfected flies maintained their Sodalis populations at densities comparable to those in flies harboring native symbionts. Our ability to transinfect tsetse flies is indicative of Sodalis ' recent evolutionary history with its tsetse fly host and demonstrates that this procedure may be used as a means of streamlining future paratransgenesis experiments.     

The COG database: an updated version includes eukaryotes

Background

The availability of multiple, essentially complete genome sequences of prokaryotes and eukaryotes spurred both the demand and the opportunity for the construction of an evolutionary classification of genes from these genomes. Such a classification system based on orthologous relationships between genes appears to be a natural framework for comparative genomics and should facilitate both functional annotation of genomes and large-scale evolutionary studies.

Results

We describe here a major update of the previously developed system for delineation of Clusters of Orthologous Groups of proteins (COGs) from the sequenced genomes of prokaryotes and unicellular eukaryotes and the construction of clusters of predicted orthologs for 7 eukaryotic genomes, which we named KOGs after eukaryotic orthologous groups. The COG collection currently consists of 138,458 proteins, which form 4873 COGs and comprise 75% of the 185,505 (predicted) proteins encoded in 66 genomes of unicellular organisms. The eukaryotic orthologous groups (KOGs) include proteins from 7 eukaryotic genomes: three animals (the nematode Caenorhabditis elegans, the fruit fly Drosophila melanogaster and Homo sapiens), one plant, Arabidopsis thaliana, two fungi (Saccharomyces cerevisiae and Schizosaccharomyces pombe), and the intracellular microsporidian parasite Encephalitozoon cuniculi. The current KOG set consists of 4852 clusters of orthologs, which include 59,838 proteins, or ~54% of the analyzed eukaryotic 110,655 gene products. Compared to the coverage of the prokaryotic genomes with COGs, a considerably smaller fraction of eukaryotic genes could be included into the KOGs; addition of new eukaryotic genomes is expected to result in substantial increase in the coverage of eukaryotic genomes with KOGs. Examination of the phyletic patterns of KOGs reveals a conserved core represented in all analyzed species and consisting of ~20% of the KOG set. This conserved portion of the KOG set is much greater than the ubiquitous portion of the COG set (~1% of the COGs). In part, this difference is probably due to the small number of included eukaryotic genomes, but it could also reflect the relative compactness of eukaryotes as a clade and the greater evolutionary stability of eukaryotic genomes.

Conclusion

The updated collection of orthologous protein sets for prokaryotes and eukaryotes is expected to be a useful platform for functional annotation of newly sequenced genomes, including those of complex eukaryotes, and genome-wide evolutionary studies.     

A simplified and efficient germline-specific CRISPR/Cas9 system for Drosophila genomic engineering

The type II CRISPR/Cas9 system (clustered regularly interspaced short palindromic repeats/CRISPR-associated) has recently emerged as an efficient and simple tool for site-specific engineering of eukaryotic genomes. To improve its applications in Drosophila genome engineering, we simplified the standard two-component CRISPR/Cas9 system by generating a stable transgenic fly line expressing the Cas9 endonuclease in the germline (Vasa-Cas9 line). By injecting vectors expressing engineered target-specific guide RNAs into Vasa-Cas9 fly embryos, mutations were generated from site-specific DNA cleavages and efficiently transmitted into progenies. Because Cas9 endonuclease is the universal component of the type II CRISPR/Cas9 system, site-specific genomic engineering based on this improved platform can be achieved with lower complexity and toxicity, greater consistency, and excellent versatility.     

On the persistence of P element in cultured lineages of Drosophila melanogaster

P transposon is known to have invaded the Drosophila melanogaster genome in the 1950s as a result of horizontal transmission from D. willistoni. Part of the evidence supporting the timing of its invasion comes from analyses of cultured drosophila lineages originating from wild flies cultivated long time in laboratory before analysis. Such analyses have shown that P element was absent from the genomes of cultured lineages established from wild flies caught from the wild before the 1950s. Although the hypothesis of P element transmission has obtained multiple lines of evidence and is beyond doubt today, we decided to test whether analysis of cultured lineages can provide some temporal information on the P element population dynamics. In the present work we demonstrate that P element present the in wild-caught flies may be lost in the cultured fly lineages after some generations. This result is in accordance with the results of at least one published work and suggests that analysis of the cultured fly lineages may sometimes be unreliable in establishing historical trends in P element population dynamics, as the transposon may be occasionally lost, perhaps in the highly inbred lineages in which not all founding females carry it.     

“Wigglesworthia morsitans” Folate (Vitamin B9) Biosynthesis Contributes to Tsetse Host Fitness

Closely related ancient endosymbionts may retain minor genomic distinctions through evolutionary time, yet the biological relevance of these small pockets of unique loci remains unknown. The tsetse fly (Diptera: Glossinidae), the sole vector of lethal African trypanosomes (Trypanosoma spp.), maintains an ancient and obligate mutualism with species belonging to the gammaproteobacterium Wigglesworthia. Extensive concordant evolution with associated Wigglesworthia species has occurred through tsetse species radiation. Accordingly, the retention of unique symbiont loci between Wigglesworthia genomes may prove instrumental toward host species-specific biological traits. Genome distinctions between “Wigglesworthia morsitans” (harbored within Glossina morsitans bacteriomes) and the basal species Wigglesworthia glossinidia (harbored within Glossina brevipalpis bacteriomes) include the retention of chorismate and downstream folate (vitamin B₉) biosynthesis capabilities, contributing to distinct symbiont metabolomes. Here, we demonstrate that these W. morsitans pathways remain functionally intact, with folate likely being systemically disseminated through a synchronously expressed tsetse folate transporter within bacteriomes. The folate produced by W. morsitans is demonstrated to be pivotal for G. morsitans sexual maturation and reproduction. Modest differences between ancient symbiont genomes may still play key roles in the evolution of their host species, particularly if loci are involved in shaping host physiology and ecology. Enhanced knowledge of the Wigglesworthia-tsetse mutualism may also provide novel and specific avenues for vector control.     

Comparative Analysis of Codon Usage Bias and Codon Context Patterns between Dipteran and Hymenopteran Sequenced Genomes

Background

Codon bias is a phenomenon of non-uniform usage of codons whereas codon context generally refers to sequential pair of codons in a gene. Although genome sequencing of multiple species of dipteran and hymenopteran insects have been completed only a few of these species have been analyzed for codon usage bias.

Methods and Principal Findings

Here, we use bioinformatics approaches to analyze codon usage bias and codon context patterns in a genome-wide manner among 15 dipteran and 7 hymenopteran insect species. Results show that GAA is the most frequent codon in the dipteran species whereas GAG is the most frequent codon in the hymenopteran species. Data reveals that codons ending with C or G are frequently used in the dipteran genomes whereas codons ending with A or T are frequently used in the hymenopteran genomes. Synonymous codon usage orders (SCUO) vary within genomes in a pattern that seems to be distinct for each species. Based on comparison of 30 one-to-one orthologous genes among 17 species, the fruit fly Drosophila willistoni shows the least codon usage bias whereas the honey bee (Apis mellifera) shows the highest bias. Analysis of codon context patterns of these insects shows that specific codons are frequently used as the 3′- and 5′-context of start and stop codons, respectively.

Conclusions

Codon bias pattern is distinct between dipteran and hymenopteran insects. While codon bias is favored by high GC content of dipteran genomes, high AT content of genes favors biased usage of synonymous codons in the hymenopteran insects. Also, codon context patterns vary among these species largely according to their phylogeny.     

Stabilization of mineralized and sclerotized puparial cuticle of muscid flies

Calcium, magnesium and phosphorus are the major mineral elements in puparial exuviae of the face fly, Musca autumnalis, house fly, M. domestica and stable fly, Stomoxys calcitrans, but they are 20–50 times more prevalent in face fly than in the other two species that sclerotize the puparium. Carbon and nitrogen are approx. 5 times more abundant in house fly puparia than in face fly puparia. Face fly puparia contain two and three-fold less total amino acids than the house fly and stable fly, respectively. β-Alanine is a major amino acid in puparial cuticle of the house fly and stable fly, but it is absent in the face fly. There is no significant difference in glucosamine (chitin) content between the three species. Dopamine is the major catechol detected in face fly puparial cuticle while N-β-alanyldopamine (NBAD) is 10 to 15 times more prevalent than other catechols such as dopamine, N-acetyldopamine (NADA), 3,4-dihydroxyphenylalanine (DOPA) and 3,4-dihydroxyphenylacetic acid (DOPAC) in house fly and stable fly puparial cuticles. The latter two species have 75 to nearly 200 times higher levels of extractable catechols than the face fly. At the onset of pupariation, dopamine and NBAD attain nearly equivalent titres in puparial cuticles of face fly and house fly, respectively. Dopamine subsequently decreases more than 40-fold in the face fly as the cuticle becomes stabilized, while NBAD continues to accumulate in the house fly. The house fly covalently incorporates about 150 times more catechols in the puparium than does the face fly. The force required to fracture house fly and stable fly puparia is about three-fold greater than that required to fracture face fly puparia of comparable thickness. However, the face fly puparium attains a strength comparable to those of house fly and stable fly puparia by significantly increasing its thickness. These results demonstrate that dipterans use both catecholamines and minerals for stabilization of puparial cuticle with the house fly and stable fly relying primarily on sclerotization and the face fly on mineralization.     

Positive Selection of Iris, a Retroviral Envelope–Derived Host Gene in Drosophila melanogaster

Eukaryotic genomes can usurp enzymatic functions encoded by mobile elements for their own use. A particularly interesting kind of acquisition involves the domestication of retroviral envelope genes, which confer infectious membrane-fusion ability to retroviruses. So far, these examples have been limited to vertebrate genomes, including primates where the domesticated envelope is under purifying selection to assist placental function. Here, we show that in Drosophila genomes, a previously unannotated gene (CG4715, renamed Iris) was domesticated from a novel, active Kanga lineage of insect retroviruses at least 25 million years ago, and has since been maintained as a host gene that is expressed in all adult tissues. Iris and the envelope genes from Kanga retroviruses are homologous to those found in insect baculoviruses and gypsy and roo insect retroviruses. Two separate envelope domestications from the Kanga and roo retroviruses have taken place, in fruit fly and mosquito genomes, respectively. Whereas retroviral envelopes are proteolytically cleaved into the ligand-interaction and membrane-fusion domains, Iris appears to lack this cleavage site. In the takahashii/suzukii species groups of Drosophila, we find that Iris has tandemly duplicated to give rise to two genes (Iris-A and Iris-B). Iris-B has significantly diverged from the Iris-A lineage, primarily because of the “invention” of an intron de novo in what was previously exonic sequence. Unlike domesticated retroviral envelope genes in mammals, we find that Iris has been subject to strong positive selection between Drosophila species. The rapid, adaptive evolution of Iris is sufficient to unambiguously distinguish the phylogenies of three closely related sibling species of Drosophila (D. simulans, D. sechellia, and D. mauritiana), a discriminative power previously described only for a putative “speciation gene.” Iris represents the first instance of a retroviral envelope–derived host gene outside vertebrates. It is also the first example of a retroviral envelope gene that has been found to be subject to positive selection following its domestication. The unusual selective pressures acting on Iris suggest that it is an active participant in an ongoing genetic conflict. We propose a model in which Iris has “switched sides,” having been recruited by host genomes to combat baculoviruses and retroviruses, which employ homologous envelope genes to mediate infection.     

RAG1 core and V(D)J recombination signal sequences were derived from Transib transposons

The V(D)J recombination reaction in jawed vertebrates is catalyzed by the RAG1 and RAG2 proteins, which are believed to have emerged approximately 500 million years ago from transposon-encoded proteins. Yet no transposase sequence similar to RAG1 or RAG2 has been found. Here we show that the approximately 600-amino acid “core” region of RAG1 required for its catalytic activity is significantly similar to the transposase encoded by DNA transposons that belong to the Transib superfamily. This superfamily was discovered recently based on computational analysis of the fruit fly and African malaria mosquito genomes. Transib transposons also are present in the genomes of sea urchin, yellow fever mosquito, silkworm, dog hookworm, hydra, and soybean rust. We demonstrate that recombination signal sequences (RSSs) were derived from terminal inverted repeats of an ancient Transib transposon. Furthermore, the critical DDE catalytic triad of RAG1 is shared with the Transib transposase as part of conserved motifs. We also studied several divergent proteins encoded by the sea urchin and lancelet genomes that are 25%−30% identical to the RAG1 N-terminal domain and the RAG1 core. Our results provide the first direct evidence linking RAG1 and RSSs to a specific superfamily of DNA transposons and indicate that the V(D)J machinery evolved from transposons. We propose that only the RAG1 core was derived from the Transib transposase, whereas the N-terminal domain was assembled from separate proteins of unknown function that may still be active in sea urchin, lancelet, hydra, and starlet sea anemone. We also suggest that the RAG2 protein was not encoded by ancient Transib transposons but emerged in jawed vertebrates as a counterpart of RAG1 necessary for the V(D)J recombination reaction.     

In vivo zinc toxicity phenotypes provide a sensitized background that suggests zinc transport activities for most of the Drosophila Zip and ZnT genes

Members of the ZIP (SLC39A) and ZnT (SLC30A) families of transmembrane domain proteins are predicted to transport the essential transition metal zinc across membranes, regulating cellular zinc content and distribution via uptake and efflux at the outer plasma and organellar membranes. Twenty-four ZIP and ZnT proteins are encoded in mammalian genomes, raising questions of whether all actually transport zinc, whether several function together in the same tissues/cell types, and how the activity of these transporters is coordinated. To address these questions, we have taken advantage of the ability to manipulate several genes simultaneously in targeted cell types in Drosophila. Previously we reported zinc toxicity phenotypes caused by combining overexpression of a zinc uptake gene, dZip42C.1, with suppression of a zinc efflux gene, dZnT63C. Here we show that these phenotypes can be used as a sensitized in vivo system to detect subtle alterations in zinc transport activity that would be buffered in healthy cells. Using two adult tissues, the fly eye and midline (thorax/abdomen), we find that when overexpressed, most of the 17 Drosophila Zip and ZnT genes modify the zinc toxicity phenotypes in a manner consistent with their predicted zinc transport activity. In most cases, we can reconcile that activity with the cellular localization of an enhanced green fluorescent protein tagged version of the protein. Additionally, targeted suppression of each gene by RNA interference reveals several of the fly Zip and ZnT genes are required in the eye, indicating that numerous independent zinc transport genes are acting together in a single tissue.     

Presence of the endosymbiont Wolbachia among some fruit flies (Diptera: Tephritidae) from Iran: A multilocus sequence typing approach

Wolbachia is a widespread endosymbiont of insects with a diverse range of biological effects on its hosts. We studied the prevalence of Wolbachia in some important species of tephritids in Iran. Among different populations of five fruit fly species, Dacus ciliatus (cucurbit fly), Rhagoletis cerasi (cherry fruit fly), Ceratitis capitata (Mediterranean fruit fly), Myiopardalis pardalina (melon fly) and Carypomya vesuviana (jujube fly), two species, R. cerasi and C. vesuviana, showed infection with separate Wolbachia strains, namely wCer6 and wVes1, respectively. C. vesuviana is introduced here as a novel host for Wolbachia. Genotyping of Wolbachia strains in 12 populations of five fruit fly species, using multilocus sequence typing (MLST) and the wsp gene sequence showed the occurrence of two new strains as well as a new strain type (ST) belonging to the A supergroup. On the basis of the results of this study, 12 barcodes under five species of Iranian tephritids have been added to the database of DNA barcodes. Inter- and intra-specific differences among COI sequences showed a clear gap in barcoding among most fruit flies.     

Targeting the tsetse-trypanosome interplay using genetically engineered Sodalis glossinidius

Sodalis glossinidius, a secondary bacterial symbiont of the tsetse fly, is currently considered as a potential delivery system for anti-trypanosomal components interfering with African trypanosome transmission (i.e. paratransgenesis). Nanobodies (Nbs) have been proposed as potential candidates to target the parasite during development in the tsetse fly. In this study, we have generated an immune Nb-library and developed a panning strategy to select Nbs against the Trypanosoma brucei brucei procyclic developmental stage present in the tsetse fly midgut. Selected Nbs were expressed, purified, assessed for binding and tested for their impact on the survival and growth of in vitro cultured procyclic T. b. brucei parasites. Next, we engineered S. glossinidius to express the selected Nbs and validated their ability to block T. brucei development in the tsetse fly midgut. Genetically engineered S. glossinidius expressing Nb_88 significantly compromised parasite development in the tsetse fly midgut both at the level of infection rate and parasite load. Interestingly, expression of Nb_19 by S. glossinidius resulted in a significantly enhanced midgut establishment. These data are the first to show in situ delivery by S. glossinidius of effector molecules that can target the trypanosome-tsetse fly crosstalk, interfering with parasite development in the fly. These proof-of-principle data represent a major step forward in the development of a control strategy based on paratransgenic tsetse flies. Finally, S. glossinidius-based Nb delivery can also be applied as a powerful laboratory tool to unravel the molecular determinants of the parasite-vector association.     

Comparison of the thermal performance between a population of the olive fruit fly and its co-adapted parasitoids

Long-term separation of a host from its native parasitoids may result in divergent thermal adaptation between host and parasitoid. The olive fruit fly, Bactrocera oleae (Rossi), most likely originated from Sub-Saharan Africa, but has since had a long invasion history in cultivated olives that spans geographical barriers and continents. This study compared three major thermal performance profiles (development, survival, and reproduction) across a wide range of temperatures (10–34 °C) among a Californian population of the olive fruit fly and two African parasitoids, Psyttalia lounsburyi (Silvestri) and Psyttalia humilis (Silvestri), believed to have co-adapted with the fruit fly in its native range. Temperature ranges for the development and survival were 10–30 °C for the fly, 10–28 °C for P. lounsburyi, and 14–32 °C for P. humilis. There was no difference in any thermal performance measured between two P. humilis populations (Kenya and Namibia) tested. The most suitable temperature ranges for reproduction were 22–30 °C for the fly, 18–32 °C for P. humilis, and 18–26 °C for P. lounsburyi. The results showed slight differences in the thermal profiles among olive fruit fly and both parasitoids species, with P. humilis being more heat tolerant whereas P. lounsburyi was less heat tolerant than the fruit fly. The results are discussed with respect to thermal co-adaptation and classical biological control of the olive fruit fly.     

Tsetse EP Protein Protects the Fly Midgut from Trypanosome Establishment

African trypanosomes undergo a complex developmental process in their tsetse fly vector before transmission back to a vertebrate host. Typically, 90% of fly infections fail, most during initial establishment of the parasite in the fly midgut. The specific mechanism(s) underpinning this failure are unknown. We have previously shown that a Glossina-specific, immunoresponsive molecule, tsetse EP protein, is up regulated by the fly in response to gram-negative microbial challenge. Here we show by knockdown using RNA interference that this tsetse EP protein acts as a powerful antagonist of establishment in the fly midgut for both Trypanosoma brucei brucei and T. congolense. We demonstrate that this phenomenon exists in two species of tsetse, Glossina morsitans morsitans and G. palpalis palpalis, suggesting tsetse EP protein may be a major determinant of vector competence in all Glossina species. Tsetse EP protein levels also decline in response to starvation of the fly, providing a possible explanation for increased susceptibility of starved flies to trypanosome infection. As starvation is a common field event, this fact may be of considerable importance in the epidemiology of African trypanosomiasis.     

Cluster J Mycobacteriophages: Intron Splicing in Capsid and Tail Genes

Bacteriophages isolated on Mycobacterium smegmatis mc²155 represent many distinct genomes sharing little or no DNA sequence similarity. The genomes are architecturally mosaic and are replete with genes of unknown function. A new group of genomes sharing substantial nucleotide sequences constitute Cluster J. The six mycobacteriophages forming Cluster J are morphologically members of the Siphoviridae, but have unusually long genomes ranging from 106.3 to 117 kbp. Reconstruction of the capsid by cryo-electron microscopy of mycobacteriophage BAKA reveals an icosahedral structure with a triangulation number of 13. All six phages are temperate and homoimmune, and prophage establishment involves integration into a tRNA-Leu gene not previously identified as a mycobacterial attB site for phage integration. The Cluster J genomes provide two examples of intron splicing within the virion structural genes, one in a major capsid subunit gene, and one in a tail gene. These genomes also contain numerous free-standing HNH homing endonuclease, and comparative analysis reveals how these could contribute to genome mosaicism. The unusual Cluster J genomes provide new insights into phage genome architecture, gene function, capsid structure, gene mobility, intron splicing, and evolution.