Morphology and husbandry of the large milkweed bug, Oncopeltus fasciatus

Theme and variation in the evolution of insect sex determination

The development of dimorphic adult sexes is a critical process for most animals, one that is subject to intense selection. Work in vertebrate and insect model species has revealed that sex determination mechanisms vary widely among animal groups. However, this variation is not uniform, with a limited number of conserved factors. Therefore, sex determination offers an excellent context to consider themes and variations in gene network evolution. Here we review the literature describing sex determination in diverse insects. We have screened public genomic sequence databases for orthologs and duplicates of 25 genes involved in insect sex determination, identifying patterns of presence and absence. These genes and a 3.5 reference set of 43 others were used to infer phylogenies and compared to accepted organismal relationships to examine patterns of congruence and divergence. The function of candidate genes for roles in sex determination (virilizer, female-lethal-2-d, transformer-2) and sex chromosome dosage compensation (male specific lethal-1, msl-2, msl-3) were tested using RNA interference in the milkweed bug, Oncopeltus fasciatus. None of these candidate genes exhibited conserved roles in these processes. Amidst this variation we wish to highlight the following themes for the evolution of sex determination: (1) Unique features within taxa influence network evolution. (2) Their position in the network influences a component's evolution. Our analyses also suggest an inverse association of protein sequence conservation with functional conservation.

Distinct developmental mechanisms influence sexual dimorphisms in the milkweed bug Oncopeltus fasciatus

Sexual dimorphism is common in animals. The most complete model of sex determination comes from Drosophila melanogaster, where the relative dosage of autosomes and X chromosomes leads indirectly to sex-specific transcripts of doublesex (dsx). Female Dsx interacts with a mediator complex protein encoded by intersex (ix) to activate female development. In males, the transcription factor encoded by fruitless (fru) promotes male-specific behaviour. The genetics of sex determination have been examined in a small number of other insects, yet several questions remain about the plesiomorphic state. Is dsx required for female and male development? Is fru conserved in male behaviour or morphology? Are other components such as ix functionally conserved? To address these questions, we report expression and functional tests of dsx, ix and fru in the hemipteran Oncopeltus fasciatus, characterizing three sexual dimorphisms. dsx prevents ix phenotypes in all sexes and dimorphic traits in the milkweed bug. ix and fru are expressed across the body, in females and males. fru and ix also affect the genitalia of both sexes, but have effects limited to different dimorphic structures in different sexes. These results reveal roles for ix and fru distinct from other insects, and demonstrate distinct development mechanisms in different sexually dimorphic structures.

High-Efficiency CRISPR/Cas9 Mutagenesis of the white Gene in the Milkweed Bug Oncopeltus fasciatus

In this manuscript, we report that clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 is highly efficient in the hemipteran Oncopeltus fasciatus. The white gene is well characterized in Drosophila where mutation causes loss of eye pigmentation; white is a reliable marker for transgenesis and other genetic manipulations. Accordingly, white has been targeted in a number of nonmodel insects to establish tools for genetic studies. Here, we generated mutations in the Of-white (Of-w) locus using CRISPR/Cas9. We found that Of-w is required for pigmentation throughout the body of Oncopeltus, not just the ommatidia. High rates of somatic mosaicism were observed in the injected generation, reflecting biallelic mutations, and a high rate of germline mutation was evidenced by the large proportion of heterozygous G1s. However, Of-w mutations are homozygous lethal; G2 homozygotes lacked pigment dispersion throughout the body and did not hatch, precluding the establishment of a stable mutant line. Embryonic and parental RNA interference (RNAi) were subsequently performed to rule out off-target mutations producing the observed phenotype and to evaluate the efficacy of RNAi in ablating gene function compared to a loss-of-function mutation. RNAi knockdowns phenocopied Of-w homozygotes, with an unusual accumulation of orange granules observed in unhatched embryos. This is, to our knowledge, the first CRISPR/Cas9-targeted mutation generated in Oncopeltus. While we were unable to establish white as a useful visible marker for Oncopeltus, these findings are instructive for the selection of visible markers in nonmodel species and reveal an unusual role for an ortholog of a classic Drosophila gene.

Codon and Amino Acid Usage Are Shaped by Selection Across Divergent Model Organisms of the Pancrustacea

In protein-coding genes, synonymous codon usage and amino acid composition correlate to expression in some eukaryotes, and may result from translational selection. Here, we studied large-scale RNA-seq data from three divergent arthropod models, including cricket (Gryllus bimaculatus), milkweed bug (Oncopeltus fasciatus), and the amphipod crustacean Parhyale hawaiensis, and tested for optimization of codon and amino acid usage relative to expression level. We report strong signals of AT3 optimal codons (those favored in highly expressed genes) in G. bimaculatus and O. fasciatus, whereas weaker signs of GC3 optimal codons were found in P. hawaiensis, suggesting selection on codon usage in all three organisms. Further, in G. bimaculatus and O. fasciatus, high expression was associated with lowered frequency of amino acids with large size/complexity (S/C) scores in favor of those with intermediate S/C values; thus, selection may favor smaller amino acids while retaining those of moderate size for protein stability or conformation. In P. hawaiensis, highly transcribed genes had elevated frequency of amino acids with large and small S/C scores, suggesting a complex dynamic in this crustacean. In all species, the highly transcribed genes appeared to favor short proteins, high optimal codon usage, specific amino acids, and were preferentially involved in cell-cycling and protein synthesis. Together, based on examination of 1,680,067, 1,667,783, and 1,326,896 codon sites in G. bimaculatus, O. fasciatus, and P. hawaiensis, respectively, we conclude that translational selection shapes codon and amino acid usage in these three Pancrustacean arthropods.

Dynamic BMP signaling polarized by Toll patterns the dorsoventral axis in a hemimetabolous insect

Toll-dependent patterning of the dorsoventral axis in Drosophila represents one of the best understood gene regulatory networks. However, its evolutionary origin has remained elusive. Outside the insects Toll is not known for a patterning function, but rather for a role in pathogen defense. Here, we show that in the milkweed bug Oncopeltus fasciatus, whose lineage split from Drosophila's more than 350 million years ago, Toll is only required to polarize a dynamic BMP signaling network. A theoretical model reveals that this network has self-regulatory properties and that shallow Toll signaling gradients are sufficient to initiate axis formation. Such gradients can account for the experimentally observed twinning of insect embryos upon egg fragmentation and might have evolved from a state of uniform Toll activity associated with protecting insect eggs against pathogens.

Dual mode of embryonic development is highlighted by expression and function of Nasonia pair-rule genes

Embryonic anterior-posterior patterning is well understood in Drosophila, which uses 'long germ' embryogenesis, in which all segments are patterned before cellularization. In contrast, most insects use 'short germ' embryogenesis, wherein only head and thorax are patterned in a syncytial environment while the remainder of the embryo is generated after cellularization. We use the wasp Nasonia (Nv) to address how the transition from short to long germ embryogenesis occurred. Maternal and gap gene expression in Nasonia suggest long germ embryogenesis. However, the Nasonia pair-rule genes even-skipped, odd-skipped, runt and hairy are all expressed as early blastoderm pair-rule stripes and late-forming posterior stripes. Knockdown of Nv eve, odd or h causes loss of alternate segments at the anterior and complete loss of abdominal segments. We propose that Nasonia uses a mixed mode of segmentation wherein pair-rule genes pattern the embryo in a manner resembling Drosophila at the anterior and ancestral Tribolium at the posterior. DOI: http://dx.doi.org/10.7554/eLife.01440.001.

ASGARD: an open-access database of annotated transcriptomes for emerging model arthropod species

The increased throughput and decreased cost of next-generation sequencing (NGS) have shifted the bottleneck genomic research from sequencing to annotation, analysis and accessibility. This is particularly challenging for research communities working on organisms that lack the basic infrastructure of a sequenced genome, or an efficient way to utilize whatever sequence data may be available. Here we present a new database, the Assembled Searchable Giant Arthropod Read Database (ASGARD). This database is a repository and search engine for transcriptomic data from arthropods that are of high interest to multiple research communities but currently lack sequenced genomes. We demonstrate the functionality and utility of ASGARD using de novo assembled transcriptomes from the milkweed bug Oncopeltus fasciatus, the cricket Gryllus bimaculatus and the amphipod crustacean Parhyale hawaiensis. We have annotated these transcriptomes to assign putative orthology, coding region determination, protein domain identification and Gene Ontology (GO) term annotation to all possible assembly products. ASGARD allows users to search all assemblies by orthology annotation, GO term annotation or Basic Local Alignment Search Tool. User-friendly features of ASGARD include search term auto-completion suggestions based on database content, the ability to download assembly product sequences in FASTA format, direct links to NCBI data for predicted orthologs and graphical representation of the location of protein domains and matches to similar sequences from the NCBI non-redundant database. ASGARD will be a useful repository for transcriptome data from future NGS studies on these and other emerging model arthropods, regardless of sequencing platform, assembly or annotation status. This database thus provides easy, one-stop access to multi-species annotated transcriptome information. We anticipate that this database will be useful for members of multiple research communities, including developmental biology, physiology, evolutionary biology, ecology, comparative genomics and phylogenomics.

Database URL:asgard.rc.fas.harvard.edu

The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus

BACKGROUND

Most evolutionary developmental biology ("evo-devo") studies of emerging model organisms focus on small numbers of candidate genes cloned individually using degenerate PCR. However, newly available sequencing technologies such as 454 pyrosequencing have recently begun to allow for massive gene discovery in animals without sequenced genomes. Within insects, although large volumes of sequence data are available for holometabolous insects, developmental studies of basally branching hemimetabolous insects typically suffer from low rates of gene discovery.

RESULTS

We used 454 pyrosequencing to sequence over 500 million bases of cDNA from the ovaries and embryos of the milkweed bug Oncopeltus fasciatus, which lacks a sequenced genome. This indirectly developing insect occupies an important phylogenetic position, branching basal to Diptera (including fruit flies) and Hymenoptera (including honeybees), and is an experimentally tractable model for short-germ development. 2,087,410 reads from both normalized and non-normalized cDNA assembled into 21,097 sequences (isotigs) and 112,531 singletons. The assembled sequences fell into 16,617 unique gene models, and included predictions of splicing isoforms, which we examined experimentally. Discovery of new genes plateaued after assembly of ~1.5 million reads, suggesting that we have sequenced nearly all transcripts present in the cDNA sampled. Many transcripts have been assembled at close to full length, and there is a net gain of sequence data for over half of the pre-existing O. fasciatus accessions for developmental genes in GenBank. We identified 10,775 unique genes, including members of all major conserved metazoan signaling pathways and genes involved in several major categories of early developmental processes. We also specifically address the effects of cDNA normalization on gene discovery in de novo transcriptome analyses.

CONCLUSIONS

Our sequencing, assembly and annotation framework provide a simple and effective way to achieve high-throughput gene discovery for organisms lacking a sequenced genome. These data will have applications to the study of the evolution of arthropod genes and genetic pathways, and to the wider evolution, development and genomics communities working with emerging model organisms.[The sequence data from this study have been submitted to GenBank under study accession number SRP002610 (http://www.ncbi.nlm.nih.gov/sra?term=SRP002610). Custom scripts generated are available at http://www.extavourlab.com/protocols/index.html. Seven Additional files are available.].