Embryo polarity in moth flies and mosquitoes relies on distinct old genes with localized transcript isoforms

Temporal control of RNAi reveals both robust and labile feedback loops in the segmentation clock of the red flour beetle

Animals from all major clades have evolved a segmented trunk, reflected in the human spine or the insect segments. These units emerge during embryogenesis from a posterior segment addition zone (SAZ), where repetitive gene activity is regulated by a mechanism described by the clock and wavefront/speed gradient model. In the red flour beetle Tribolium castaneum, RNA interference (RNAi) has been used to continuously knock down the function of primary pair-rule genes (pPRGs), caudal or Wnt pathway components, which has led to the complete breakdown of segmentation. However, it has remained untested, if this breakdown was reversible by bringing the missing gene function back to the system. To fill this gap, we established a transgenic system in T. castaneum, which allows blocking an ongoing RNAi effect with temporal control by expressing a viral inhibitor of RNAi via heat shock. We show that the T. castaneum segmentation machinery was able to reestablish after RNAi targeting the pPRGs Tc-eve, Tc-odd, and Tc-runt was blocked. However, we observed no rescue after blocking RNAi targeting Wnt pathway components. We conclude that the insect segmentation system contains both robust feedback loops that can reestablish and labile feedback loops that break down irreversibly. This combination may reconcile conflicting needs of the system: Labile systems controlling initiation and maintenance of the SAZ ensure that only one SAZ is formed. Robust feedback loops confer developmental robustness toward external disturbances.

Improved methodology for fixation and preparation of Aedes aegypti embryos

The yellow fever mosquito Aedes aegypti is a major disease vector and an increasingly popular emerging model research organism. We present here an improved protocol for the collection, fixation, and preparation of A. aegypti embryos for immunohistochemical and in situ hybridization studies. The processing of A. aegypti embryos for such studies is complicated by the inability to easily remove the vitelline membrane, which prevents the reagents needed for staining from reaching their targets, and which furthermore obscures visualization of the embryo since the membrane is highly sclerotized. Previously described protocols for removal of the vitelline membrane are very low throughput, limiting the capacity of work that can be accomplished in a reasonable timeframe. Our adapted protocol increases the throughput capacity of embryos by an individual user, with experienced users able to prepare an average of 100-150 embryos per hour. The protocol provides high-quality intact embryos that can be used for morphological, immunohistochemical, and in situ hybridization studies. The protocol has been successfully tested on embryos of ages ranging from 14h after egg laying (AEL) at 27°C through to 55h AEL. Critical to the success of the optimized protocol is the selection, fabrication, and description of the tools required. To this end, a video-demonstrated protocol has been placed at protocols.io to clarify the protocol and provide easy access and training to anyone interested in the preparation of A. aegypti embryos for biological studies.

Neofunctionalization driven by positive selection led to the retention of the loqs2 gene encoding an Aedes specific dsRNA binding protein

BACKGROUND

Mosquito borne viruses, such as dengue, Zika, yellow fever and Chikungunya, cause millions of infections every year. These viruses are mostly transmitted by two urban-adapted mosquito species, Aedes aegypti and Aedes albopictus. Although mechanistic understanding remains largely unknown, Aedes mosquitoes may have unique adaptations that lower the impact of viral infection. Recently, we reported the identification of an Aedes specific double-stranded RNA binding protein (dsRBP), named Loqs2, that is involved in the control of infection by dengue and Zika viruses in mosquitoes. Preliminary analyses suggested that the loqs2 gene is a paralog of loquacious (loqs) and r2d2, two co-factors of the RNA interference (RNAi) pathway, a major antiviral mechanism in insects.

RESULTS

Here we analyzed the origin and evolution of loqs2. Our data suggest that loqs2 originated from two independent duplications of the first double-stranded RNA binding domain of loqs that occurred before the origin of the Aedes Stegomyia subgenus, around 31 million years ago. We show that the loqs2 gene is evolving under relaxed purifying selection at a faster pace than loqs, with evidence of neofunctionalization driven by positive selection. Accordingly, we observed that Loqs2 is localized mainly in the nucleus, different from R2D2 and both isoforms of Loqs that are cytoplasmic. In contrast to r2d2 and loqs, loqs2 expression is stage- and tissue-specific, restricted mostly to reproductive tissues in adult Ae. aegypti and Ae. albopictus. Transgenic mosquitoes engineered to express loqs2 ubiquitously undergo developmental arrest at larval stages that correlates with massive dysregulation of gene expression without major effects on microRNAs or other endogenous small RNAs, classically associated with RNA interference.

CONCLUSIONS

Our results uncover the peculiar origin and neofunctionalization of loqs2 driven by positive selection. This study shows an example of unique adaptations in Aedes mosquitoes that could ultimately help explain their effectiveness as virus vectors.

The chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus uncovers patterns of genome evolution in mosquitoes

BACKGROUND

Understanding genome organization and evolution is important for species involved in transmission of human diseases, such as mosquitoes. Anophelinae and Culicinae subfamilies of mosquitoes show striking differences in genome sizes, sex chromosome arrangements, behavior, and ability to transmit pathogens. However, the genomic basis of these differences is not fully understood.

METHODS

In this study, we used a combination of advanced genome technologies such as Oxford Nanopore Technology sequencing, Hi-C scaffolding, Bionano, and cytogenetic mapping to develop an improved chromosome-scale genome assembly for the West Nile vector Culex quinquefasciatus.

RESULTS

We then used this assembly to annotate odorant receptors, odorant binding proteins, and transposable elements. A genomic region containing male-specific sequences on chromosome 1 and a polymorphic inversion on chromosome 3 were identified in the Cx. quinquefasciatus genome. In addition, the genome of Cx. quinquefasciatus was compared with the genomes of other mosquitoes such as malaria vectors An. coluzzi and An. albimanus, and the vector of arboviruses Ae. aegypti. Our work confirms significant expansion of the two chemosensory gene families in Cx. quinquefasciatus, as well as a significant increase and relocation of the transposable elements in both Cx. quinquefasciatus and Ae. aegypti relative to the Anophelines. Phylogenetic analysis clarifies the divergence time between the mosquito species. Our study provides new insights into chromosomal evolution in mosquitoes and finds that the X chromosome of Anophelinae and the sex-determining chromosome 1 of Culicinae have a significantly higher rate of evolution than autosomes.

CONCLUSION

The improved Cx. quinquefasciatus genome assembly uncovered new details of mosquito genome evolution and has the potential to speed up the development of novel vector control strategies.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, limitations, and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for continuous extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies, while the containers and reproducible workflows could easily be deployed and extended to evaluate new methods or data sets.

Extensible benchmarking of methods that identify and quantify polyadenylation sites from RNA-seq data

The tremendous rate with which data is generated and analysis methods emerge makes it increasingly difficult to keep track of their domain of applicability, assumptions, and limitations and consequently, of the efficacy and precision with which they solve specific tasks. Therefore, there is an increasing need for benchmarks, and for the provision of infrastructure for continuous method evaluation. APAeval is an international community effort, organized by the RNA Society in 2021, to benchmark tools for the identification and quantification of the usage of alternative polyadenylation (APA) sites from short-read, bulk RNA-sequencing (RNA-seq) data. Here, we reviewed 17 tools and benchmarked eight on their ability to perform APA identification and quantification, using a comprehensive set of RNA-seq experiments comprising real, synthetic, and matched 3'-end sequencing data. To support continuous benchmarking, we have incorporated the results into the OpenEBench online platform, which allows for seamless extension of the set of methods, metrics, and challenges. We envisage that our analyses will assist researchers in selecting the appropriate tools for their studies. Furthermore, the containers and reproducible workflows generated in the course of this project can be seamlessly deployed and extended in the future to evaluate new methods or datasets.

Twenty-seven ZAD-ZNF genes of Drosophila melanogaster are orthologous to the embryo polarity determining mosquito gene cucoid

The C2H2 zinc finger gene cucoid establishes anterior-posterior (AP) polarity in the early embryo of culicine mosquitoes. This gene is unrelated to genes that establish embryo polarity in other fly species (Diptera), such as the homeobox gene bicoid, which serves this function in the traditional model organism Drosophila melanogaster. The cucoid gene is a conserved single copy gene across lower dipterans but nothing is known about its function in other species, and its evolution in higher dipterans, including Drosophila, is unresolved. We found that cucoid is a member of the ZAD-containing C2H2 zinc finger (ZAD-ZNF) gene family and is orthologous to 27 of the 91 members of this family in D. melanogaster, including M1BP, ranshi, ouib, nom, zaf1, odj, Nnk, trem, Zif, and eighteen uncharacterized genes. Available knowledge of the functions of cucoid orthologs in Drosophila melanogaster suggest that the progenitor of this lineage specific expansion may have played a role in regulating chromatin. We also describe many aspects of the gene duplication history of cucoid in the brachyceran lineage of D. melanogaster, thereby providing a framework for predicting potential redundancies among these genes in D. melanogaster.

Evolution and loss of ß-catenin and TCF-dependent axis specification in insects

Mechanisms and evolution of primary axis specification in insects are discussed in the context of the roles of ß-catenin and TCF in polarizing metazoan embryos. Three hypotheses are presented. First, insects with sequential segmentation and posterior growth use cell-autonomous mechanisms for establishing embryo polarity via the nuclear ratio of ß-catenin and TCF. Second, TCF homologs establish competence for anterior specification. Third, the evolution of simultaneous segmentation mechanisms, also known as long-germ development, resulted in primary axis specification mechanisms that are independent of ß-catenin but reliant on TCF, a condition that preceded the frequent replacement of anterior determinants in long germ insects.

Screens in fly and beetle reveal vastly divergent gene sets required for developmental processes

Background

Most of the known genes required for developmental processes have been identified by genetic screens in a few well-studied model organisms, which have been considered representative of related species, and informative—to some degree—for human biology. The fruit fly Drosophila melanogaster is a prime model for insect genetics, and while conservation of many gene functions has been observed among bilaterian animals, a plethora of data show evolutionary divergence of gene function among more closely-related groups, such as within the insects. A quantification of conservation versus divergence of gene functions has been missing, without which it is unclear how representative data from model systems actually are.

Results

Here, we systematically compare the gene sets required for a number of homologous but divergent developmental processes between fly and beetle in order to quantify the difference of the gene sets. To that end, we expanded our RNAi screen in the red flour beetle Tribolium castaneum to cover more than half of the protein-coding genes. Then we compared the gene sets required for four different developmental processes between beetle and fly. We found that around 50% of the gene functions were identified in the screens of both species while for the rest, phenotypes were revealed only in fly (~ 10%) or beetle (~ 40%) reflecting both technical and biological differences. Accordingly, we were able to annotate novel developmental GO terms for 96 genes studied in this work. With this work, we publish the final dataset for the pupal injection screen of the iBeetle screen reaching a coverage of 87% (13,020 genes).

Conclusions

We conclude that the gene sets required for a homologous process diverge more than widely believed. Hence, the insights gained in flies may be less representative for insects or protostomes than previously thought, and work in complementary model systems is required to gain a comprehensive picture. The RNAi screening resources developed in this project, the expanding transgenic toolkit, and our large-scale functional data make T. castaneum an excellent model system in that endeavor.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01231-4.

Evolution of a cytoplasmic determinant: evidence for the biochemical basis of functional evolution of the novel germ line regulator oskar

Germ line specification is essential in sexually reproducing organisms. Despite their critical role, the evolutionary history of the genes that specify animal germ cells is heterogeneous and dynamic. In many insects, the gene oskar is required for the specification of the germ line. However, the germ line role of oskar is thought to be a derived role resulting from co-option from an ancestral somatic role. To address how evolutionary changes in protein sequence could have led to changes in the function of Oskar protein that enabled it to regulate germ line specification, we searched for oskar orthologs in 1,565 publicly available insect genomic and transcriptomic data sets. The earliest-diverging lineage in which we identified an oskar ortholog was the order Zygentoma (silverfish and firebrats), suggesting that oskar originated before the origin of winged insects. We noted some order-specific trends in oskar sequence evolution, including whole gene duplications, clade-specific losses, and rapid divergence. An alignment of all known 379 Oskar sequences revealed new highly conserved residues as candidates that promote dimerization of the LOTUS domain. Moreover, we identified regions of the OSK domain with conserved predicted RNA binding potential. Furthermore, we show that despite a low overall amino acid conservation, the LOTUS domain shows higher conservation of predicted secondary structure than the OSK domain. Finally, we suggest new key amino acids in the LOTUS domain that may be involved in the previously reported Oskar−Vasa physical interaction that is required for its germ line role.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Regional patterning and regulation of melanin pigmentation in insects

Insects display an immense diversity in melanin pigmentation, which is generated by the interplay between the regulatory genes (that provide general patterning information) and effector genes (that provide coloration of the pattern). However, recent studies encompassing several different orders (Hemiptera, Blattodea, Coleoptera, and Lepidoptera) have shown that knockdowns of melanin producing genes alone can generate distinct region-specific patterns. This review surveys the most recent studies to further document the regional patterning of effector genes, and highlights the new advances and their implications for future research.

PAS-seq 2: A fast and sensitive method for global profiling of polyadenylated RNAs

Alternative polyadenylation (APA) is a widespread phenomenon in eukaryotes that contributes to regulating gene expression and generating proteomic diversity. APA plays critical roles in development and its mis-regulation has been implicated in a wide variety of human diseases, including cancer. To study APA on the transcriptome-wide level, numerous deep sequencing methods that capture 3' end of mRNAs have been developed in the past decade, but they generally require a large amount of hands-on time and/or high RNA input. Here, we introduce PAS-seq 2, a fast and sensitive method for global and quantitative profiling of polyadenylated RNAs. Compared to our original PAS-seq, this method takes less time and requires much lower total RNA input due to improvement in the reverse transcription process. PAS-seq 2 can be applied to both APA and differential gene expression analyses.

Juxtapose: a gene-embedding approach for comparing co-expression networks

BACKGROUND

Gene co-expression networks (GCNs) are not easily comparable due to their complex structure. In this paper, we propose a tool, Juxtapose, together with similarity measures that can be utilized for comparative transcriptomics between a set of organisms. While we focus on its application to comparing co-expression networks across species in evolutionary studies, Juxtapose is also generalizable to co-expression network comparisons across tissues or conditions within the same species.

METHODS

A word embedding strategy commonly used in natural language processing was utilized in order to generate gene embeddings based on walks made throughout the GCNs. Juxtapose was evaluated based on its ability to embed the nodes of synthetic structures in the networks consistently while also generating biologically informative results. Evaluation of the techniques proposed in this research utilized RNA-seq datasets from GTEx, a multi-species experiment of prefrontal cortex samples from the Gene Expression Omnibus, as well as synthesized datasets. Biological evaluation was performed using gene set enrichment analysis and known gene relationships in literature.

RESULTS

We show that Juxtapose is capable of globally aligning synthesized networks as well as identifying areas that are conserved in real gene co-expression networks without reliance on external biological information. Furthermore, output from a matching algorithm that uses cosine distance between GCN embeddings is shown to be an informative measure of similarity that reflects the amount of topological similarity between networks.

CONCLUSIONS

Juxtapose can be used to align GCNs without relying on known biological similarities and enables post-hoc analyses using biological parameters, such as orthology of genes, or conserved or variable pathways.

AVAILABILITY

A development version of the software used in this paper is available at https://github.com/klovens/juxtapose.

Zygotic pioneer factor activity of Odd-paired/Zic is necessary for late function of the Drosophila segmentation network

Because chromatin determines whether information encoded in DNA is accessible to transcription factors, dynamic chromatin states in development may constrain how gene regulatory networks impart embryonic pattern. To determine the interplay between chromatin states and regulatory network function, we performed ATAC-seq on Drosophila embryos during the establishment of the segmentation network, comparing wild-type and mutant embryos in which all graded maternal patterning inputs are eliminated. While during the period between zygotic genome activation and gastrulation many regions maintain stable accessibility, cis-regulatory modules (CRMs) within the network undergo extensive patterning-dependent changes in accessibility. A component of the network, Odd-paired (opa), is necessary for pioneering accessibility of late segmentation network CRMs. opa-driven changes in accessibility are accompanied by equivalent changes in gene expression. Interfering with the timing of opa activity impacts the proper patterning of expression. These results indicate that dynamic systems for chromatin regulation directly impact the reading of embryonic patterning information.

Oncopeltus-like gene expression patterns in Murgantia histrionica, a new hemipteran model system, suggest ancient regulatory network divergence

Background

Much has been learned about basic biology from studies of insect model systems. The pre-eminent insect model system, Drosophila melanogaster, is a holometabolous insect with a derived mode of segment formation. While additional insect models have been pioneered in recent years, most of these fall within holometabolous lineages. In contrast, hemimetabolous insects have garnered less attention, although they include agricultural pests, vectors of human disease, and present numerous evolutionary novelties in form and function. The milkweed bug, Oncopeltus fasciatus (order: Hemiptera)—close outgroup to holometabolous insects—is an emerging model system. However, comparative studies within this order are limited as many phytophagous hemipterans are difficult to stably maintain in the lab due to their reliance on fresh plants, deposition of eggs within plant material, and long development time from embryo to adult.

Results

Here we present the harlequin bug, Murgantia histrionica, as a new hemipteran model species. Murgantia—a member of the stink bug family Pentatomidae which shares a common ancestor with Oncopeltus ~ 200 mya—is easy to rear in the lab, produces a large number of eggs, and is amenable to molecular genetic techniques. We use Murgantia to ask whether Pair-Rule Genes (PRGs) are deployed in ways similar to holometabolous insects or to Oncopeltus. Specifically, PRGs even-skipped, odd-skipped, paired and sloppy-paired are initially expressed in PR-stripes in Drosophila and a number of holometabolous insects but in segmental-stripes in Oncopeltus. We found that these genes are likewise expressed in segmental-stripes in Murgantia, while runt displays partial PR-character in both species. Also like Oncopeltus, E75A is expressed in a clear PR-pattern in blastoderm- and germband-stage Murgantia embryos, although it plays no role in segmentation in Drosophila. Thus, genes diagnostic of the split between holometabolous insects and Oncopeltus are expressed in an Oncopeltus-like fashion during Murgantia development.

Conclusions

The similarity in gene expression between Murgantia and Oncopeltus suggests that Oncopeltus is not a sole outlier species in failing to utilize orthologs of Drosophila PRGs for PR-patterning. Rather, strategies deployed for PR-patterning, including the use of E75A in the PRG-network, are likely conserved within Hemiptera, and possibly more broadly among hemimetabolous insects.