Genomes and transcriptomes help unravel the complex life cycle of the blastoclad fungus, Coelomomyces lativittatus, an obligate parasite of mosquitoes and microcrustaceans.

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

Extended in vivo transcriptomes of two ascoviruses with different tissue tropisms reveal alternative mechanisms for enhancing virus reproduction in hemolymph.

Ascoviruses are large dsDNA viruses characterized by the extraordinary changes they induce in cellular pathogenesis and architecture whereby after nuclear lysis and extensive hypertrophy, each cell is cleaved into numerous vesicles for virion reproduction. However, the level of viral replication and transcription in vesicles compared to other host tissues remains uncertain. Therefore, we applied RNA-Sequencing to compare the temporal transcriptome of Spodoptera frugiperda ascovirus (SfAV) and Trichoplusia ni ascovirus (TnAV) at 7, 14, and 21 days post-infection (dpi). We found most transcription occurred in viral vesicles, not in initial tissues infected, a remarkably novel reproduction mechanism compared to all other viruses and most other intracellular pathogens. Specifically, the highest level of viral gene expression occurred in hemolymph, for TnAV at 7 dpi, and SfAV at 14 dpi. Moreover, we found that host immune genes were partially down-regulated in hemolymph, where most viral replication occurred in highly dense accumulations of vesicles.

Ecdysis triggering hormone receptors regulate male courtship behavior via antennal lobe interneurons in Drosophila.

Ecdysis triggering hormone receptors (ETHR) regulate the behavioral sequence necessary for cuticle shedding. Recent reports have documented functions for ETHR signaling in adult Drosophila melanogaster. In this study, we report that ETHR silencing in local interneurons of the antennal lobes and fruitless neurons leads to sharply increased rates of male-male courtship. RNAseq analysis of ETHR knockdown flies reveals differential expression of genes involved in axon guidance, courtship behavior and chemosensory functions. Our findings indicate an important role for ETHR in regulation of Drosophila courtship behavior through chemosensory processing in the antennal lobe.

Highly Effective Broad Spectrum Chimeric Larvicide That Targets Vector Mosquitoes Using a Lipophilic Protein.

Two mosquitocidal bacteria, Bacillus thuringiensis subsp. israelensis (Bti) and Lysinibacillus sphaericus (Ls) are the active ingredients of commercial larvicides used widely to control vector mosquitoes. Bti's efficacy is due to synergistic interactions among four proteins, Cry4Aa, Cry4Ba, Cry11Aa, and Cyt1Aa, whereas Ls's activity is caused by Bin, a heterodimer consisting of BinA, the toxin, and BinB, a midgut-binding protein. Cyt1Aa is lipophilic and synergizes Bti Cry proteins by increasing midgut binding. We fused Bti's Cyt1Aa to Ls's BinA yielding a broad-spectrum chimeric protein highly mosquitocidal to important vector species including Anopheles gambiae, Culex quinquefasciatus, and Aedes aegypti, the latter an important Zika and Dengue virus vector insensitive to Ls Bin. Aside from its vector control potential, our bioassay data, in contrast to numerous other reports, provide strong evidence that BinA does not require conformational interactions with BinB or microvillar membrane lipids to bind to its intracellular target and kill mosquitoes.

De novo phasing with X-ray laser reveals mosquito larvicide BinAB structure.

BinAB is a naturally occurring paracrystalline larvicide distributed worldwide to combat the devastating diseases borne by mosquitoes. These crystals are composed of homologous molecules, BinA and BinB, which play distinct roles in the multi-step intoxication process, transforming from harmless, robust crystals, to soluble protoxin heterodimers, to internalized mature toxin, and finally to toxic oligomeric pores. The small size of the crystals-50 unit cells per edge, on average-has impeded structural characterization by conventional means. Here we report the structure of Lysinibacillus sphaericus BinAB solved de novo by serial-femtosecond crystallography at an X-ray free-electron laser. The structure reveals tyrosine- and carboxylate-mediated contacts acting as pH switches to release soluble protoxin in the alkaline larval midgut. An enormous heterodimeric interface appears to be responsible for anchoring BinA to receptor-bound BinB for co-internalization. Remarkably, this interface is largely composed of propeptides, suggesting that proteolytic maturation would trigger dissociation of the heterodimer and progression to pore formation.

Effect of Promoters and Plasmid Copy Number on Cyt1A Synthesis and Crystal Assembly in Bacillus thuringiensis.

Cyt1Aa is a major mosquitocidal protein synthesized during sporulation of Bacillus thuringiensis subsp. israelensis, composing more than 50% of its parasporal body. This high level of synthesis is due to several factors including three strong sporulation-dependent promoters, a strong transcription termination sequence, and an associated 20-kDa helper protein. Cyt1Aa's toxicity is low compared to the Cry proteins of this species, namely, Cry4Aa, Cry4Ba, and Cry11Aa, but it nevertheless plays an important role in the biology of B. thuringiensis subsp. israelensis in that it synergizes their mosquitocidal toxicity and suppresses the evolution of resistance. In the present study, the effects of using different cyt1Aa promoter combinations and plasmid copy number on synthesis of Cyt1Aa were evaluated. Using the 4Q7 (plasmid-cured) strain of B. thuringiensis subsp. israelensis as an experimental host, a plasmid copy number of two or three yielded no Cyt1Aa, whereas a copy number of four yielded only small crystals, even when expression was driven by one of the wild-type promoters. However, using all three wild-type promoters and a plasmid copy number of 20 yielded Cyt1A crystals tenfold larger than those produced by one promoter and a plasmid copy number of four. High levels of Cyt1Aa synthesis resulted in significantly fewer spores per unit medium and imperfectly formed crystals. Similar results were obtained when Cyt1Aa synthesis was evaluated using the same expression constructs in a mutant strain of B. thuringiensis subsp. israelensis that lacks the cyt1Aa gene.

A resurrected mammalian hAT transposable element and a closely related insect element are highly active in human cell culture.

Chromosome structure and function are influenced by transposable elements, which are mobile DNA segments that can move from place to place. hAT elements are a superfamily of DNA cut and paste elements that move by excision and integration. We have characterized two hAT elements, TcBuster and Space Invaders (SPIN), that are members of a recently described subfamily of hAT elements called Buster elements. We show that TcBuster, from the red flour beetle Tribolium castaneum, is highly active in human cells. SPIN elements are currently inactive elements that were recently highly active in multiple vertebrate genomes, and the high level of sequence similarity across widely diverged species and patchy phylogenetic distribution suggest that they may have moved between genomes by horizontal transfer. We have generated an intact version of this element, SPIN(ON), which is highly active in human cells. In vitro analysis of TcBuster and SPIN(ON) shows that no proteins other than transposase are essential for recombination, a property that may contribute to the ability of SPIN to successfully invade multiple organisms. We also analyze the target site preferences of de novo insertions in the human genome of TcBuster and SPIN(ON) and compare them with the preferences of Sleeping Beauty and piggyBac, showing that each superfamily has a distinctive pattern of insertion. The high-frequency transposition of both TcBuster and SPIN(ON) suggests that these transposon systems offer powerful tools for genome engineering. Finally, we describe a Saccharomyces cerevisiae assay for TcBuster that will provide a means for isolation of hyperactive and other interesting classes of transposase mutants.

Comparative analysis of the recently discovered hAT transposon TcBuster in human cells.

BACKGROUND: Transposons are useful tools for creating transgenic organisms, insertional mutagenesis, and genome engineering. TcBuster, a novel hAT-family transposon system derived from the red flour beetle Tribolium castaneum, was shown to be highly active in previous studies in insect embryoes. METHODOLOGY/PRINCIPAL FINDINGS: We tested TcBuster for its activity in human embryonic kidney 293 (HEK-293) cells. Excision footprints obtained from HEK-293 cells contained small insertions and deletions consistent with a hAT-type repair mechanism of hairpin formation and non-homologous end-joining. Genome-wide analysis of 23,417 piggyBac, 30,303 Sleeping Beauty, and 27,985 TcBuster integrations in HEK-293 cells revealed a uniquely different integration pattern when compared to other transposon systems with regards to genomic elements. TcBuster experimental conditions were optimized to assay TcBuster activity in HEK-293 cells by colony assay selection for a neomycin-containing transposon. Increasing transposon plasmid increased the number of colonies, whereas gene transfer activity dependent on codon-optimized transposase plasmid peaked at 100 ng with decreased colonies at the highest doses of transposase DNA. Expression of the related human proteins Buster1, Buster3, and SCAND3 in HEK-293 cells did not result in genomic integration of the TcBuster transposon. TcBuster, Tol2, and piggyBac were compared directly at different ratios of transposon to transposase and found to be approximately comparable while having their own ratio preferences. CONCLUSIONS/SIGNIFICANCE: TcBuster was found to be highly active in mammalian HEK-293 cells and represents a promising tool for mammalian genome engineering.

The mosquito Aedes aegypti has a large genome size and high transposable element load but contains a low proportion of transposon-specific piRNAs.

BACKGROUND: The piRNA pathway has been shown in model organisms to be involved in silencing of transposons thereby providing genome stability. In D. melanogaster the majority of piRNAs map to these sequences. The medically important mosquito species Aedes aegypti has a large genome size, a high transposon load which includes Miniature Inverted repeat Transposable Elements (MITES) and an expansion of the piRNA biogenesis genes. Studies of transgenic lines of Ae. aegypti have indicated that introduced transposons are poorly remobilized and we sought to explore the basis of this. We wished to analyze the piRNA profile of Ae. aegypti and thereby determine if it is responsible for transposon silencing in this mosquito. RESULTS: Estimated piRNA sequence diversity was comparable between Ae. aegypti and D. melanogaster, but surprisingly only 19% of mosquito piRNAs mapped to transposons compared to 51% for D. melanogaster. Ae. aegypti piRNA clusters made up a larger percentage of the total genome than those of D. melanogaster but did not contain significantly higher percentages of transposon derived sequences than other regions of the genome. Ae. aegypti contains a number of protein coding genes that may be sources of piRNA biogenesis with two, traffic jam and maelstrom, implicated in this process in model organisms. Several genes of viral origin were also targeted by piRNAs. Examination of six mosquito libraries that had previously been transformed with transposon derived sequence revealed that new piRNA sequences had been generated to the transformed sequences, suggesting that they may have stimulated a transposon inactivation mechanism. CONCLUSIONS: Ae. aegypti has a large piRNA complement that maps to transposons but primarily gene sequences, including many viral-derived sequences. This, together the more uniform distribution of piRNA clusters throughout its genome, suggest that some aspects of the piRNA system differ between Ae. aegypti and D. melanogaster.

DNA sequence requirements for hobo transposable element transposition in Drosophila melanogaster.

We have conducted a structure and functional analysis of the hobo transposable element of Drosophila melanogaster. A minimum of 141 bp of the left (L) end and 65 bp of the right (R) end of the hobo were shown to contain sequences sufficient for transposition. Both ends of hobo contain multiple copies of the motifs GGGTG and GTGGC and we show that the frequency of hobo transposition increases as a function of the copy number of these motifs. The R end of hobo contains a unique 12 bp internal inverted repeat that is identical to the hobo terminal inverted repeats. We show that this internal inverted repeat suppresses transposition activity in a hobo element containing an intact L end and only 475 bp of the R end. In addition to establishing cis-sequences requirements for transposition, we analyzed trans-sequence effects of the hobo transposase. We show a hobo transposase lacking the first 49 amino acids catalyzed hobo transposition at a higher frequency than the full-length transposase suggesting that, similar to the related Ac transposase, residues at the amino end of the transposase reduce transposition. Finally, we compared target site sequences of hobo with those of the related Hermes element and found both transposons have strong preferences for the same insertion sites.

DNA binding activities of the Herves transposase from the mosquito Anopheles gambiae.

BACKGROUND: Determining the mechanisms by which transposable elements move within a genome increases our understanding of how they can shape genome evolution. Class 2 transposable elements transpose via a 'cut-and-paste' mechanism mediated by a transposase that binds to sites at or near the ends of the transposon. Herves is a member of the hAT superfamily of class 2 transposons and was isolated from Anopheles gambiae, a medically important mosquito species that is the major vector of malaria in sub-Saharan Africa. Herves is transpositionally active and intact copies of it are found in field populations of A gambiae. In this study we report the binding activities of the Herves transposase to the sequences at the ends of the Herves transposon and compare these to other sequences recognized by hAT transposases isolated from other organisms. RESULTS: We identified the specific DNA-binding sites of the Herves transposase. Active Herves transposase was purified using an Escherichia coli expression system and bound in a site-specific manner to the subterminal and terminal sequences of the left and right ends of the element, respectively, and also interacted with the right but not the left terminal inverted repeat. We identified a common subterminal DNA-binding motif (CG/AATTCAT) that is critical and sufficient for Herves transposase binding. CONCLUSIONS: The Herves transposase binds specifically to a short motif located at both ends of the transposon but shows differential binding with respect to the left and right terminal inverted repeats. Despite similarities in the overall structures of hAT transposases, the regions to which they bind in their respective transposons differ in sequence ensuring the specificity of these enzymes to their respective transposon. The asymmetry with which the Herves terminal inverted repeats are bound by the transposase may indicate that these differ in their interactions with the enzyme.

Phylogenetic and functional characterization of the hAT transposon superfamily.

Transposons are found in virtually all organisms and play fundamental roles in genome evolution. They can also acquire new functions in the host organism and some have been developed as incisive genetic tools for transformation and mutagenesis. The hAT transposon superfamily contains members from the plant and animal kingdoms, some of which are active when introduced into new host organisms. We have identified two new active hAT transposons, AeBuster1, from the mosquito Aedes aegypti and TcBuster from the red flour beetle Tribolium castaneum. Activity of both transposons is illustrated by excision and transposition assays performed in Drosophila melanogaster and Ae. aegypti and by in vitro strand transfer assays. These two active insect transposons are more closely related to the Buster sequences identified in humans than they are to the previously identified active hAT transposons, Ac, Tam3, Tol2, hobo, and Hermes. We therefore reexamined the structural and functional relationships of hAT and hAT-like transposase sequences extracted from genome databases and found that the hAT superfamily is divided into at least two families. This division is supported by a difference in target-site selections generated by active transposons of each family. We name these families the Ac and Buster families after the first identified transposon or transposon-like sequence in each. We find that the recently discovered SPIN transposons of mammals are located within the family of Buster elements.

Transpositionally active episomal hAT elements.

BACKGROUND: hAT elements and V(D)J recombination may have evolved from a common ancestral transposable element system. Extrachromosomal, circular forms of transposable elements (referred to here as episomal forms) have been reported yet their biological significance remains unknown. V(D)J signal joints, which resemble episomal transposable elements, have been considered non-recombinogenic products of V(D)J recombination and a safe way to dispose of excised chromosomal sequences. V(D)J signal joints can, however, participate in recombination reactions and the purpose of this study was to determine if hobo and Hermes episomal elements are also recombinogenic. RESULTS: Up to 50% of hobo/Hermes episomes contained two intact, inverted-terminal repeats and 86% of these contained from 1-1000 bp of intercalary DNA. Episomal hobo/Hermes elements were recovered from Musca domestica (a natural host of Hermes), Drosophila melanogaster (a natural host of hobo) and transgenic Drosophila melanogaster and Aedes aegypti (with autonomous Hermes elements). Episomal Hermes elements were recovered from unfertilized eggs of M. domestica and D. melanogaster demonstrating their potential for extrachromosomal, maternal transmission. Reintegration of episomal Hermes elements was observed in vitro and in vivo and the presence of Hermes episomes resulted in lower rates of canonical Hermes transposition in vivo. CONCLUSION: Episomal hobo/Hermes elements are common products of element excision and can be maternally transmitted. Episomal forms of Hermes are capable of integration and also of influencing the transposition of canonical elements suggesting biological roles for these extrachromosomal elements in element transmission and regulation.

Gene vector and transposable element behavior in mosquitoes.

The development of efficient germ-line transformation technologies for mosquitoes has increased the ability of entomologists to find, isolate and analyze genes. The utility of the currently available systems will be determined by a number of factors including the behavior of the gene vectors during the initial integration event and their behavior after chromosomal integration. Post-integration behavior will determine whether the transposable elements being employed currently as primary gene vectors will be useful as gene-tagging and enhancer-trapping agents. The post-integration behavior of existing insect vectors has not been extensively examined. Mos1 is useful as a primary germ-line transformation vector in insects but is inefficiently remobilized in Drosophila melanogaster and Aedes aegypti. Hermes transforms D. melanogaster efficiently and can be remobilized in this species. This element is also useful for creating transgenic A. aegypti, but its mode of integration in mosquitoes results in the insertion of flanking plasmid DNA. Hermes can be remobilized in the soma of A. aegypti and transposes using a common cut-and-paste mechanism; however, the element does not remobilize in the germ line. piggyBac can be used to create transgenic mosquitoes and occasionally integrates using a mechanism other than a simple cut-and-paste mechanism. Preliminary data suggest that remobilization is infrequent. Minos also functions in mosquitoes and, like the other gene vectors, appears to remobilize inefficiently following integration. These results have implications for future gene vector development efforts and applications.