Cave beetle lineages gained genes before going down under: An example of repeated genomic exaptation?

The adaptation of animals to subterranean habitats like caves and aquifers stereotypically leads to dramatic trait-loss consequences like the lack of eyes and body pigmentation. These body plan regression trends are expected to be tied to gene loss as well. Indeed, previous studies documented the degeneration of vision genes in obligate cave dwellers. Contradicting this picture, the first broad-scale comparative transcriptome-wide study of gene content evolution in separate subterranean Australian and Mediterranean beetle clades unearthed evidence of global gene gain and retention. This suggests that the transition to cave life may be more contingent on gene repertoire expansion than contraction. Future studies, however, will need to examine how much the observed patterns of gene content evolution reflect subfunctionalization and fitness-securing genetic redundancy outcomes following gene duplication as opposed to adaptive trajectories.

Coming into clear sight at last: Ancestral and derived events during chelicerate visual system development.

Pioneering molecular work on chelicerate visual system development in the horseshoe crab Limulus polyphemus surprised with the possibility that this process may not depend on the deeply conserved retinal determination function of Pax6 transcription factors. Genomic, transcriptomic, and developmental studies in spiders now reveal that the arthropod Pax6 homologs eyeless and twin of eyeless act as ancestral determinants of the ocular head segment in chelicerates, which clarifies deep gene regulatory and structural homologies and recommends more unified terminologies in the comparison of arthropod visual systems. Following this phylotypic stage, chelicerate visual system development differs fundamentally from other arthropods during the compartmentalization of the ocular segment in that eye and optic neuropil primordia originate independently from each other. Comparative analyses of this phase identified further gene regulatory homologies but also major differences, most notably the possibly compensatory replacement of Pax6 by Pax2 in lateral eye specification. Also see the video abstract here: https://youtu.be/Hdfr3z5kEXU.

Close to complete conservation of the brachyceran opsin repertoire in the stalk-eyed fly Teleopsis dalmanni.

Due to the unique morphology of their adult visual system, stalk-eyed flies represent an important model of exaggerated trait evolution through sexual selection. Early physiological measurements indicated wavelength sensitivity peaks in the ultraviolet (360 nm), blue (450), blue-green (490 nm), and red (>550 nm) ranges in the compound eye retina of the stalk-eyed fly Teleopsis dalmanni, consistent with the trichromatic color and broad range motion detection vision system of brachyceran Diptera. A previous study of dipteran opsin gene diversification, however, detected only homologs of members of the long wavelength range sensitive opsin subfamilies Rh2 and Rh6 in T. dalmanni. Here, I report findings from analyzing the most recent T. dalmanni genome assembly, which revealed the conservation of most brachyceran opsin homologs except for the UV wavelength range-sensitive homolog Rh4. These results and other examples highlight the caution that needs to be applied to gene loss conclusions.

The brachyceran de novo gene PIP82, a phosphorylation target of aPKC, is essential for proper formation and maintenance of the rhabdomeric photoreceptor apical domain in Drosophila.

The Drosophila apical photoreceptor membrane is defined by the presence of two distinct morphological regions, the microvilli-based rhabdomere and the stalk membrane. The subdivision of the apical membrane contributes to the geometrical positioning and the stereotypical morphology of the rhabdomeres in compound eyes with open rhabdoms and neural superposition. Here we describe the characterization of the photoreceptor specific protein PIP82. We found that PIP82's subcellular localization demarcates the rhabdomeric portion of the apical membrane. We further demonstrate that PIP82 is a phosphorylation target of aPKC. PIP82 localization is modulated by phosphorylation, and in vivo, the loss of the aPKC/Crumbs complex results in an expansion of the PIP82 localization domain. The absence of PIP82 in photoreceptors leads to misshapped rhabdomeres as a result of misdirected cellular trafficking of rhabdomere proteins. Comparative analyses reveal that PIP82 originated de novo in the lineage leading to brachyceran Diptera, which is also characterized by the transition from fused to open rhabdoms. Taken together, these findings define a novel factor that delineates and maintains a specific apical membrane domain, and offers new insights into the functional organization and evolutionary history of the Drosophila retina.

Semper's cells in the insect compound eye: Insights into ocular form and function.

The arthropod compound eye represents one of two major eye types in the animal kingdom and has served as an essential experimental paradigm for defining fundamental mechanisms underlying sensory organ formation, function, and maintenance. One of the most distinguishing features of the compound eye is the highly regular array of lens facets that define individual eye (ommatidial) units. These lens facets are produced by a deeply conserved quartet of cuticle-secreting cells, called Semper cells (SCs). Also widely known as cone cells, SCs were originally identified for their secretion of the dioptric system, i.e. the corneal lens and underlying crystalline cones. Additionally, SCs are now known to execute a diversity of patterning and glial functions in compound eye development and maintenance. Here, we present an integrated account of our current knowledge of SC multifunctionality in the Drosophila compound eye, highlighting emerging gene regulatory modules that may drive the diverse roles for these cells. Drawing comparisons with other deeply conserved retinal glia in the vertebrate single lens eye, this discussion speaks to glial cell origins and opens new avenues for understanding sensory system support programs.

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.

Genome-enabled insights into the biology of thrips as crop pests.

BACKGROUND: The western flower thrips, Frankliniella occidentalis (Pergande), is a globally invasive pest and plant virus vector on a wide array of food, fiber, and ornamental crops. The underlying genetic mechanisms of the processes governing thrips pest and vector biology, feeding behaviors, ecology, and insecticide resistance are largely unknown. To address this gap, we present the F. occidentalis draft genome assembly and official gene set. RESULTS: We report on the first genome sequence for any member of the insect order Thysanoptera. Benchmarking Universal Single-Copy Ortholog (BUSCO) assessments of the genome assembly (size = 415.8 Mb, scaffold N50 = 948.9 kb) revealed a relatively complete and well-annotated assembly in comparison to other insect genomes. The genome is unusually GC-rich (50%) compared to other insect genomes to date. The official gene set (OGS v1.0) contains 16,859 genes, of which ~ 10% were manually verified and corrected by our consortium. We focused on manual annotation, phylogenetic, and expression evidence analyses for gene sets centered on primary themes in the life histories and activities of plant-colonizing insects. Highlights include the following: (1) divergent clades and large expansions in genes associated with environmental sensing (chemosensory receptors) and detoxification (CYP4, CYP6, and CCE enzymes) of substances encountered in agricultural environments; (2) a comprehensive set of salivary gland genes supported by enriched expression; (3) apparent absence of members of the IMD innate immune defense pathway; and (4) developmental- and sex-specific expression analyses of genes associated with progression from larvae to adulthood through neometaboly, a distinct form of maturation differing from either incomplete or complete metamorphosis in the Insecta. CONCLUSIONS: Analysis of the F. occidentalis genome offers insights into the polyphagous behavior of this insect pest that finds, colonizes, and survives on a widely diverse array of plants. The genomic resources presented here enable a more complete analysis of insect evolution and biology, providing a missing taxon for contemporary insect genomics-based analyses. Our study also offers a genomic benchmark for molecular and evolutionary investigations of other Thysanoptera species.

Correction to: Genome-enabled insights into the biology of thrips as crop pests.

An amendment to this paper has been published and can be accessed via the original article.

Enhanced genome assembly and a new official gene set for Tribolium castaneum.

BACKGROUND: The red flour beetle Tribolium castaneum has emerged as an important model organism for the study of gene function in development and physiology, for ecological and evolutionary genomics, for pest control and a plethora of other topics. RNA interference (RNAi), transgenesis and genome editing are well established and the resources for genome-wide RNAi screening have become available in this model. All these techniques depend on a high quality genome assembly and precise gene models. However, the first version of the genome assembly was generated by Sanger sequencing, and with a small set of RNA sequence data limiting annotation quality. RESULTS: Here, we present an improved genome assembly (Tcas5.2) and an enhanced genome annotation resulting in a new official gene set (OGS3) for Tribolium castaneum, which significantly increase the quality of the genomic resources. By adding large-distance jumping library DNA sequencing to join scaffolds and fill small gaps, the gaps in the genome assembly were reduced and the N50 increased to 4753kbp. The precision of the gene models was enhanced by the use of a large body of RNA-Seq reads of different life history stages and tissue types, leading to the discovery of 1452 novel gene sequences. We also added new features such as alternative splicing, well defined UTRs and microRNA target predictions. For quality control, 399 gene models were evaluated by manual inspection. The current gene set was submitted to Genbank and accepted as a RefSeq genome by NCBI. CONCLUSIONS: The new genome assembly (Tcas5.2) and the official gene set (OGS3) provide enhanced genomic resources for genetic work in Tribolium castaneum. The much improved information on transcription start sites supports transgenic and gene editing approaches. Further, novel types of information such as splice variants and microRNA target genes open additional possibilities for analysis.

Brown marmorated stink bug, Halyomorpha halys (Stål), genome: putative underpinnings of polyphagy, insecticide resistance potential and biology of a top worldwide pest.

BACKGROUND: Halyomorpha halys (Stål), the brown marmorated stink bug, is a highly invasive insect species due in part to its exceptionally high levels of polyphagy. This species is also a nuisance due to overwintering in human-made structures. It has caused significant agricultural losses in recent years along the Atlantic seaboard of North America and in continental Europe. Genomic resources will assist with determining the molecular basis for this species' feeding and habitat traits, defining potential targets for pest management strategies. RESULTS: Analysis of the 1.15-Gb draft genome assembly has identified a wide variety of genetic elements underpinning the biological characteristics of this formidable pest species, encompassing the roles of sensory functions, digestion, immunity, detoxification and development, all of which likely support H. halys' capacity for invasiveness. Many of the genes identified herein have potential for biomolecular pesticide applications. CONCLUSIONS: Availability of the H. halys genome sequence will be useful for the development of environmentally friendly biomolecular pesticides to be applied in concert with more traditional, synthetic chemical-based controls.

Evolutionary conservation of opsin gene expression patterns in the compound eyes of darkling beetles.

Recent large-scale studies of opsin gene contents in representatives of the largest order of insects, the Coleoptera (beetles), revealed that the blue wavelength-sensitive (B) opsin subfamily is absent in this clade, while the ultraviolet- (UV) and long wavelength-sensitive (LW) opsin subfamilies are broadly conserved with gene duplications possibly reintroducing blue sensitivity in select subclades. Little is known yet, however, how opsin genes are expressed in the compound eyes of beetles. In a previous study, we analyzed opsin gene expression in the red flour beetle Tribolium castaneum, a member of the family of darkling beetles (Tenebrionidae), and found that a singleton LW opsin homolog is homogeneously expressed in all photoreceptors of the compound eye retina with a singleton UV opsin homolog being co-expressed in the R7 subtype photoreceptors. To probe for the evolutionary conservation of these expression patterns, we isolated complete opsin transcript sequences from three additional species in the subfamily Tenebrionidae (Tribolium confusum, Tenebrio molitor, Zophobas morio) and studied their expression via whole mount in situ hybridization in the pupal retina. These experiments revealed very similar, if not identical, photoreceptor subtype-specific expression patterns in all three species compared with T. castaneum. Documenting a deep conservation of photoreceptor subtype-specific opsin gene expression in this range of darkling beetles, our study provides a first point of reference for broader comparative studies of retinal organization in the Coleoptera.

Towards reconstructing the dipteran demise of an ancient essential gene: E3 ubiquitin ligase Murine double minute.

Genome studies have uncovered many examples of essential gene loss, raising the question of how ancient genes transition from essentiality to dispensability. We explored this process for the deeply conserved E3 ubiquitin ligase Murine double minute (Mdm), which is lacking in Drosophila despite the conservation of its main regulatory target, the cellular stress response gene p53. Conducting gene expression and knockdown experiments in the red flour beetle Tribolium castaneum, we found evidence that Mdm has remained essential in insects where it is present. Using bioinformatics approaches, we confirm the absence of the Mdm gene family in Drosophila, mapping its loss to the stem lineage of schizophoran Diptera and Pipunculidae (big-headed flies), about 95-85 million years ago. Intriguingly, this gene loss event was preceded by the de novo origin of the gene Companion of reaper (Corp), a novel p53 regulatory factor that is characterized by functional similarities to vertebrate Mdm2 despite lacking E3 ubiquitin ligase protein domains. Speaking against a 1:1 compensatory gene gain/loss scenario, however, we found that hoverflies (Syrphidae) and pointed-wing flies (Lonchopteridae) possess both Mdm and Corp. This implies that the two p53 regulators have been coexisting for ~ 150 million years in select dipteran clades and for at least 50 million years in the lineage to Schizophora and Pipunculidae. Given these extensive time spans of Mdm/Corp coexistence, we speculate that the loss of Mdm in the lineage to Drosophila involved further acquisitions of compensatory gene activities besides the emergence of Corp. Combined with the previously noted reduction of an ancestral P53 contact domain in the Mdm homologs of crustaceans and insects, we conclude that the loss of the ancient Mdm gene family in flies was the outcome of incremental functional regression over long macroevolutionary time scales.

Ancient genetic redundancy of eyeless and twin of eyeless in the arthropod ocular segment.

Pax6 transcription factors are essential upstream regulators in the developing anterior brain and peripheral visual system of most bilaterian animals. While a single homolog is in charge of these functions in vertebrates, two Pax6 genes are in Drosophila: eyeless (ey) and twin of eyeless (toy). At first glance, their co-existence seems sufficiently explained by their differential involvement in the specification of two types of insect visual organs: the lateral compound eyes (ey) and the dorsal ocelli (toy). Less straightforward to understand, however, is their genetic redundancy in promoting defined early and late growth phases of the precursor tissue to these organs: the eye-antennal imaginal disc. Drawing on comparative sequence, expression, and gene function evidence, I here conclude that this gene regulatory network module dates back to the dawn of arthropod evolution, securing the embryonic development of the ocular head segment. Thus, ey and toy constitute a paradigm to explore the organization and functional significance of longterm conserved genetic redundancy of duplicated genes. Indeed, as first steps in this direction, recent studies uncovered the shared use of binding sites in shared enhancers of target genes that are under redundant (string) and, strikingly, even subfunctionalized control by ey and toy (atonal). Equally significant, the evolutionarily recent and paralog-specific function of ey to repress the transcription of the antenna fate regulator Distal-less offers a functionally and phylogenetically well-defined opportunity to study the reconciliation of shared, partitioned, and newly acquired functions in a duplicated developmental gene pair.

Distinct regulation of atonal in a visual organ of Drosophila: Organ-specific enhancer and lack of autoregulation in the larval eye.

Drosophila has three types of visual organs, the larval eyes or Bolwig's organs (BO), the ocelli (OC) and the compound eyes (CE). In all, the bHLH protein Atonal (Ato) functions as the proneural factor for photoreceptors and effects the transition from progenitor cells to differentiating neurons. In this work, we investigate the regulation of ato expression in the BO primordium (BOP). Surprisingly, we find that ato transcription in the BOP is entirely independent of the shared regulatory DNA for the developing CE and OC. The core enhancer for BOP expression, atoBO, lies ~6kb upstream of the ato gene, in contrast to the downstream location of CE and OC regulatory elements. Moreover, maintenance of ato expression in the neuronal precursors through autoregulation-a common and ancient feature of ato expression that is well-documented in eyes, ocelli and chordotonal organs-does not occur in the BO. We also show that the atoBO enhancer contains two binding sites for the transcription factor Sine oculis (So), a core component of the progenitor specification network in all three visual organs. These binding sites function in vivo and are specifically bound by So in vitro. Taken together, our findings reveal that the control of ato transcription in the evolutionarily derived BO has diverged considerably from ato regulation in the more ancestral compound eyes and ocelli, to the extent of acquiring what appears to be a distinct and evolutionarily novel cis-regulatory module.

The genome of the water strider Gerris buenoi reveals expansions of gene repertoires associated with adaptations to life on the water.

BACKGROUND: Having conquered water surfaces worldwide, the semi-aquatic bugs occupy ponds, streams, lakes, mangroves, and even open oceans. The diversity of this group has inspired a range of scientific studies from ecology and evolution to developmental genetics and hydrodynamics of fluid locomotion. However, the lack of a representative water strider genome hinders our ability to more thoroughly investigate the molecular mechanisms underlying the processes of adaptation and diversification within this group. RESULTS: Here we report the sequencing and manual annotation of the Gerris buenoi (G. buenoi) genome; the first water strider genome to be sequenced thus far. The size of the G. buenoi genome is approximately 1,000 Mb, and this sequencing effort has recovered 20,949 predicted protein-coding genes. Manual annotation uncovered a number of local (tandem and proximal) gene duplications and expansions of gene families known for their importance in a variety of processes associated with morphological and physiological adaptations to a water surface lifestyle. These expansions may affect key processes associated with growth, vision, desiccation resistance, detoxification, olfaction and epigenetic regulation. Strikingly, the G. buenoi genome contains three insulin receptors, suggesting key changes in the rewiring and function of the insulin pathway. Other genomic changes affecting with opsin genes may be associated with wavelength sensitivity shifts in opsins, which is likely to be key in facilitating specific adaptations in vision for diverse water habitats. CONCLUSIONS: Our findings suggest that local gene duplications might have played an important role during the evolution of water striders. Along with these findings, the sequencing of the G. buenoi genome now provides us the opportunity to pursue exciting research opportunities to further understand the genomic underpinnings of traits associated with the extreme body plan and life history of water striders.

Massive parallel expansions of Methuselah/Methuselah-like receptors in schizophoran Diptera.

The Methuselah/Methuselah-like (Mth/Mthl) family of G-protein coupled receptors (GPCRs) is represented by 16 homologs in the fruit fly Drosophila melanogaster. Three of them have thus far been functionally characterized and found to play critical roles in cell adhesion, immunity, lifespan, and oxidative stress regulation. Evolutionary studies have shown that the large number of D. melanogaster Mth/Mthl GPCRs arose by at least two rounds of gene duplications. The first produced the "mth superclade" subfamily and was followed by the expansion of the "melanogaster subgroup" cluster within the "mth superclade" of Mth/Mthl GPCRs. The adaptive significance of the Mth/Mthl receptor repertoire expansion in Drosophila remains elusive. Studying the Mth/Mthl gene family content in newly available dipteran genomes, we find that the first expansion of the mthl superclade predates the diversification of schizophoran Diptera approximately 65 million years ago. Unexpectedly, we further find that the subsequent expansion of the melanogaster subgroup cluster was paralleled by independent mth superclade Mth/Mthl GPCR expansions in other schizophoran clades (Muscidae and Tephritidae). Our study thus reveals an even more dynamic diversification of mth superclade GPCRs than previously appreciated and linked to the emergence of schizophoran flies, the most dramatic radiation in the dipteran tree of life.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology.

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

Shared and distinct mechanisms of atonal regulation in Drosophila ocelli and compound eyes.

The fruit fly Drosophila melanogaster has two types of external visual organs, a pair of compound eyes and a group of three ocelli. At the time of neurogenesis, the proneural transcription factor Atonal mediates the transition from progenitor cells to differentiating photoreceptor neurons in both organs. In the developing compound eye, atonal (ato) expression is directly induced by transcriptional regulators that confer retinal identity, the Retinal Determination (RD) factors. Little is known, however, about control of ato transcription in the ocelli. Here we show that a 2kb genomic DNA fragment contains distinct and common regulatory elements necessary for ato induction in compound eyes and ocelli. The three binding sites that mediate direct regulation by the RD factors Sine oculis and Eyeless in the compound eye are also required in the ocelli. However, in the latter, these sites mediate control by Sine oculis and the other Pax6 factor of Drosophila, Twin of eyeless, which can bind the Pax6 sites in vitro. Moreover, the three sites are differentially utilized in the ocelli: all three are similarly essential for atonal induction in the posterior ocelli, but show considerable redundancy in the anterior ocellus. Strikingly, this difference parallels the distinct control of ato transcription in the posterior and anterior progenitors of the developing compound eyes. From a comparative perspective, our findings suggest that the ocelli of arthropods may have originated through spatial partitioning from the dorsal edge of an ancestral compound eye.

Common transcriptional mechanisms for visual photoreceptor cell differentiation among Pancrustaceans.

A hallmark of visual rhabdomeric photoreceptors is the expression of a rhabdomeric opsin and uniquely associated phototransduction molecules, which are incorporated into a specialized expanded apical membrane, the rhabdomere. Given the extensive utilization of rhabdomeric photoreceptors in the eyes of protostomes, here we address whether a common transcriptional mechanism exists for the differentiation of rhabdomeric photoreceptors. In Drosophila, the transcription factors Pph13 and Orthodenticle (Otd) direct both aspects of differentiation: rhabdomeric opsin transcription and rhabdomere morphogenesis. We demonstrate that the orthologs of both proteins are expressed in the visual systems of the distantly related arthropod species Tribolium castaneum and Daphnia magna and that their functional roles are similar in these species. In particular, we establish that the Pph13 homologs have the ability to bind a subset of Rhodopsin core sequence I sites and that these sites are present in key phototransduction genes of both Tribolium and Daphnia. Furthermore, Pph13 and Otd orthologs are capable of executing deeply conserved functions of photoreceptor differentiation as evidenced by the ability to rescue their respective Drosophila mutant phenotypes. Pph13 homologs are equivalent in their ability to direct both rhabdomere morphogenesis and opsin expression within Drosophila, whereas Otd paralogs demonstrate differential abilities to regulate photoreceptor differentiation. Finally, loss-of-function analyses in Tribolium confirm the conserved requirement of Pph13 and Otd in regulating both rhabdomeric opsin transcription and rhabdomere morphogenesis. Taken together, our data identify components of a regulatory framework for rhabdomeric photoreceptor differentiation in Pancrustaceans, providing a foundation for defining ancestral regulatory modules of rhabdomeric photoreceptor differentiation.