Identification of early zygotic genes in the yellow fever mosquito Aedes aegypti and discovery of a motif involved in early zygotic genome activation

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.

Fertility decline in female mosquitoes is regulated by the orco olfactory co-receptor

Female Aedes aegypti mosquitoes undergo multiple rounds of reproduction, known as gonotrophic cycles. These cycles span the period from blood meal intake to oviposition. Understanding how reproductive success is maintained across gonotrophic cycles allows for the identification of molecular targets to reduce mosquito population growth. Odorant receptor co-receptor (orco) encodes a conserved insect-specific transmembrane ion channel that complexes with tuning odorant receptors (ORs) to form a functional olfactory receptor. orco expression has been identified in the male and female mosquito germline, but its role is unclear. We report an orco-dependent, maternal effect reduction in fertility after the first gonotrophic cycle. This phenotype was removed by CRISPR-Cas9 reversion of the orco mutant locus. Eggs deposited by orco mutant females are fertilized but the embryos reveal developmental defects, reduced hatching, and changes in ion channel signaling gene transcription. We present an unexpected role for an olfactory receptor pathway in mosquito reproduction.

Arbovirus-vector protein interactomics identifies Loquacious as a co-factor for dengue virus replication in Aedes mosquitoes

Efficient virus replication in Aedes vector mosquitoes is essential for the transmission of arboviral diseases such as dengue virus (DENV) in human populations. Like in vertebrates, virus-host protein-protein interactions are essential for viral replication and immune evasion in the mosquito vector. Here, 79 mosquito host proteins interacting with DENV non-structural proteins NS1 and NS5 were identified by label-free mass spectrometry, followed by a functional screening. We confirmed interactions with host factors previously observed in mammals, such as the oligosaccharyltransferase complex, and we identified protein-protein interactions that seem to be specific for mosquitoes. Among the interactors, the double-stranded RNA (dsRNA) binding protein Loquacious (Loqs), an RNA interference (RNAi) cofactor, was found to be essential for efficient replication of DENV and Zika virus (ZIKV) in mosquito cells. Loqs did not affect viral RNA stability or translation of a DENV replicon and its proviral activity was independent of its RNAi regulatory activity. Interestingly, Loqs colocalized with DENV dsRNA replication intermediates in infected cells and directly interacted with high affinity with DENV RNA in the 3’ untranslated region in vitro (K_D = 48–62 nM). Our study provides an interactome for DENV NS1 and NS5 and identifies Loqs as a key proviral host factor in mosquitoes. We propose that DENV hijacks a factor of the RNAi mechanism for replication of its own RNA.

Dengue virus is a mosquito-transmitted virus endemic to the tropics and subtropics, affecting an estimated 390 million people yearly. While the mechanisms of infection, pathogenesis and immune evasion have been extensively studied in humans, replication in Aedes mosquitoes has received much less attention, despite being a critical step in the arbovirus transmission cycle. Here, we used a proteomic approach to identify Aedes mosquito proteins recruited by dengue virus non-structural proteins NS1 and NS5. In addition to previously established host proteins that interact with DENV in mammals, we identified Loquacious, a double-stranded RNA binding protein involved in the RNAi-based antiviral immune response of mosquitoes. Unexpectedly, our data showed that Loquacious functions as a proviral factor that is recruited to replication organelles to facilitate viral replication. We propose that DENV exploits host immune components, such as Loquacious, for its own benefit.

Uncovering the Worldwide Diversity and Evolution of the Virome of the Mosquitoes Aedes aegypti and Aedes albopictus

Aedes aegypti, the yellow fever mosquito, and Aedes albopictus, the Asian tiger mosquito, are the most significant vectors of dengue, Zika, and Chikungunya viruses globally. Studies examining host factors that control arbovirus transmission demonstrate that insect-specific viruses (ISVs) can modulate mosquitoes’ susceptibility to arbovirus infection in both in vivo and in vitro co-infection models. While research is ongoing to implicate individual ISVs as proviral or antiviral factors, we have a limited understanding of the composition and diversity of the Aedes virome. To address this gap, we used a meta-analysis approach to uncover virome diversity by analysing ~3000 available RNA sequencing libraries representing a worldwide geographic range for both mosquitoes. We identified ten novel viruses and previously characterised viruses, including mononegaviruses, orthomyxoviruses, negeviruses, and a novel bi-segmented negev-like group. Phylogenetic analysis suggests close relatedness to mosquito viruses implying likely insect host range except for one arbovirus, the multi-segmented Jingmen tick virus (Flaviviridae) in an Italian colony of Ae. albopictus. Individual mosquito transcriptomes revealed remarkable inter-host variation of ISVs within individuals from the same colony and heterogeneity between different laboratory strains. Additionally, we identified striking virus diversity in Wolbachia infected Aedes cell lines. This study expands our understanding of the virome of these important vectors. It provides a resource for further assessing the ecology, evolution, and interaction of ISVs with their mosquito hosts and the arboviruses they transmit.

Nanopore long-read RNA-seq and absolute quantification delineate transcription dynamics in early embryo development of an insect pest

The olive fruit fly, Bactrocera oleae, is the most important pest for the olive fruit but lacks adequate transcriptomic characterization that could aid in molecular control approaches. We apply nanopore long-read RNA-seq with internal RNA standards allowing absolute transcript quantification to analyze transcription dynamics during early embryo development for the first time in this organism. Sequencing on the MinION platform generated over 31 million reads. Over 50% of the expressed genes had at least one read covering its entire length validating our full-length approach. We generated a de novo transcriptome assembly and identified 1768 new genes and a total of 79,810 isoforms; a fourfold increase in transcriptome diversity compared to the current NCBI predicted transcriptome. Absolute transcript quantification per embryo allowed an insight into the dramatic re-organization of maternal transcripts. We further identified Zelda as a possible regulator of early zygotic genome activation in B. oleae and provide further insights into the maternal-to-zygotic transition. These data show the utility of long-read RNA in improving characterization of non-model organisms that lack a fully annotated genome, provide potential targets for sterile insect technic approaches, and provide the first insight into the transcriptome landscape of the developing olive fruit fly embryo.

Characterization of a female germline and early zygote promoter from the transcription factor bZip1 in the dengue mosquito Aedes aegypti

Background

The wide distribution of Aedes aegypti, the main vector of dengue and yellow fever viruses, currently puts three billion people in the world at risk of infection with these viruses. Continuous transmission of these and other viruses despite aggressive efforts to prevent this emphasizes the need to develop new control strategies. Proposals to control disease transmission based on vector engineering, including both population suppression and population replacement, rely on the development of transgenes under the control of regulatory elements able to drive molecules in a specific tissue, time and strength.

Methods

Here we report the characterization of a promoter active in both the female germline and early zygote, derived from the transcription factor bZip1 in the mosquito Ae. aegypti, using transposon-based methods and RT-qPCR.

Results

We generated seven transgenic lines carrying AabZip1-reporter constructs and observed expression in both the ovary and early embryo. RT-qPCR analysis was performed to evaluate transcript expression patterns for each line, confirming that transgenic expression from the AabZip1 promoter largely recapitulated the endogenous expression pattern, albeit the strength of maternal expression appeared to be strongly influenced by chromosomal position.

Conclusions

This study provides a new regulatory sequence that can be useful for generating transgenic lines that can become a tool in vector control strategies.

Alternative polyadenylation coordinates embryonic development, sexual dimorphism and longitudinal growth in Xenopus tropicalis

RNA alternative polyadenylation contributes to the complexity of information transfer from genome to phenome, thus amplifying gene function. Here, we report the first X. tropicalis resource with 127,914 alternative polyadenylation (APA) sites derived from embryos and adults. Overall, APA networks play central roles in coordinating the maternal-zygotic transition (MZT) in embryos, sexual dimorphism in adults and longitudinal growth from embryos to adults. APA sites coordinate reprogramming in embryos before the MZT, but developmental events after the MZT due to zygotic genome activation. The APA transcriptomes of young adults are more variable than growing adults and male frog APA transcriptomes are more divergent than females. The APA profiles of young females were similar to embryos before the MZT. Enriched pathways in developing embryos were distinct across the MZT and noticeably segregated from adults. Briefly, our results suggest that the minimal functional units in genomes are alternative transcripts as opposed to genes.

Pure early zygotic genes in the Asian malaria mosquito Anopheles stephensi

BACKGROUND

The Asian malaria mosquito, Anopheles stephensi, is a major urban malaria vector in the Middle East and on the Indian subcontinent. Early zygotic transcription, which marks the maternal-to-zygotic transition, has not been systematically studied in An. stephensi or any other Anopheles mosquitoes. Improved understanding of early embryonic gene expression in An. stephensi will facilitate genetic and evolutionary studies and help with the development of novel control strategies for this important disease vector.

RESULTS

We obtained RNA-seq data in biological triplicates from four early An. stephensi embryonic time points. Using these data, we identified 70 and 153 pure early zygotic genes (pEZGs) under stringent and relaxed conditions, respectively. We show that these pEZGs are enriched in functional groups related to DNA-binding transcription regulators, cell cycle modulators, proteases, transport, and cellular metabolism. On average these pEZGs are shorter and have less introns than other An. stephensi genes. Some of the pEZGs may arise de novo while others have clear non-pEZG paralogs. There is no or very limited overlap between An. stephensi pEZGs and Drosophila melanogaster or Aedes aegypti pEZGs. Interestingly, the upstream region of An. stephensi pEZGs lack significant enrichment of a previously reported TAGteam/VBRGGTA motif found in the regulatory region of pEZGs in D. melanogaster and Ae. aegypti. However, a GT-rich motif was found in An. stephensi pEZGs instead.

CONCLUSIONS

We have identified a number of pEZGs whose predicted functions and structures are consistent with their collective roles in the degradation of maternally deposited components, activation of the zygotic genome, cell division, and metabolism. The pEZGs appear to rapidly turn over within the Dipteran order and even within the Culicidae family. These pEZGs, and the shared regulatory motif, could provide the promoter or regulatory sequences to drive gene expression in the syncytial or early cellular blastoderm, a period when the developing embryo is accessible to genetic manipulation. In addition, these molecular resources may be used to achieve sex separation of mosquitoes for sterile insect technique.

Functional analysis of the promoter of an early zygotic gene KLC2 in Aedes aegypti

Background

Aedes aegypti is an important mosquito vector that transmits arboviruses that cause devastating diseases including Zika, dengue fever, yellow fever and chikungunya. Improved understanding of gene regulation in the early development of Ae. aegypti will facilitate genetic studies and help the development of novel control strategies of this important disease vector.

Results

In this study, we demonstrated through transgenic assays that the promoter of an endogenous early zygotic gene KLC2 could drive gene expression in the syncytial blastoderm and early cellular blastoderm, which is a stage that the developing germline and the rest of embryo are accessible to genetic manipulation. An unexpected expression of the reporter gene in transgenic male testes was also observed. Further analysis confirmed the expression of the endogenous KLC2 in the testes, which was not detected in the previous RNA sequencing data.

Conclusions

Our finding provided a new promoter element that can be used in future genetic studies and applications in Ae. aegypti. Moreover, our transgenic reporter assays showed that cautions are needed when interpreting RNA sequencing data as transient or tissue-specific transcription may go undetected by RNAseq.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3210-0) contains supplementary material, which is available to authorized users.

Regulatory principles governing the maternal-to-zygotic transition: insights from Drosophila melanogaster

The onset of metazoan development requires that two terminally differentiated germ cells, a sperm and an oocyte, become reprogrammed to the totipotent embryo, which can subsequently give rise to all the cell types of the adult organism. In nearly all animals, maternal gene products regulate the initial events of embryogenesis while the zygotic genome remains transcriptionally silent. Developmental control is then passed from mother to zygote through a process known as the maternal-to-zygotic transition (MZT). The MZT comprises an intimately connected set of molecular events that mediate degradation of maternally deposited mRNAs and transcriptional activation of the zygotic genome. This essential developmental transition is conserved among metazoans but is perhaps best understood in the fruit fly, Drosophila melanogaster. In this article, we will review our understanding of the events that drive the MZT in Drosophila embryos and highlight parallel mechanisms driving this transition in other animals.

The sequence of a male-specific genome region containing the sex determination switch in Aedes aegypti

Background

Aedes aegypti is the principal vector of several important arboviruses. Among the methods of vector control to limit transmission of disease are genetic strategies that involve the release of sterile or genetically modified non-biting males, which has generated interest in manipulating mosquito sex ratios. Sex determination in Ae. aegypti is controlled by a non-recombining Y chromosome-like region called the M locus, yet characterisation of this locus has been thwarted by the repetitive nature of the genome. In 2015, an M locus gene named Nix was identified that displays the qualities of a sex determination switch.

Results

With the use of a whole-genome bacterial artificial chromosome (BAC) library, we amplified and sequenced a ~200 kb region containing the male-determining gene Nix. In this study, we show that Nix is comprised of two exons separated by a 99 kb intron primarily composed of repetitive DNA, especially transposable elements.

Conclusions

Nix, an unusually large and highly repetitive gene, exhibits features in common with Y chromosome genes in other organisms. We speculate that the lack of recombination at the M locus has allowed the expansion of repeats in a manner characteristic of a sex-limited chromosome, in accordance with proposed models of sex chromosome evolution in insects.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3090-3) contains supplementary material, which is available to authorized users.

Multiscale Approach to Deciphering the Molecular Mechanisms Involved in the Direct and Intergenerational Effect of Ibuprofen on Mosquito Aedes aegypti

The anti-inflammatory ibuprofen is a ubiquitous surface water contaminant. However, the chronic impact of this pharmaceutical on aquatic invertebrate populations remains poorly understood. In model insect Aedes aegypti, we investigated the intergenerational consequences of parental chronic exposure to an environmentally relevant concentration of ibuprofen. While exposed individuals did not show any phenotypic changes, their progeny showed accelerated development and an increased tolerance to starvation. In order to understand the mechanistic processes underpinning the direct and intergenerational impacts of ibuprofen, we combined transcriptomic, metabolomics, and hormone kinetics studies at several life stages in exposed individuals and their progeny. This integrative approach revealed moderate transcriptional changes in exposed larvae consistent with the pharmacological mode of action of ibuprofen. Parental exposure led to lower levels of several polar metabolites in progeny eggs and to major transcriptional changes in the following larval stage. These transcriptional changes, most likely driven by changes in the expression of numerous transcription factors and epigenetic regulators, led to ecdysone signaling and stress response potentiation. Overall, the present study illustrates the complexity of the molecular basis of the intergenerational pollutant response in insects and the importance of considering the entire life cycle of exposed organisms and of their progeny in order to fully understand the mode of action of pollutants and their impact on ecosystems.

A conserved maternal-specific repressive domain in Zelda revealed by Cas9-mediated mutagenesis in Drosophila melanogaster

In nearly all metazoans, the earliest stages of development are controlled by maternally deposited mRNAs and proteins. The zygotic genome becomes transcriptionally active hours after fertilization. Transcriptional activation during this maternal-to-zygotic transition (MZT) is tightly coordinated with the degradation of maternally provided mRNAs. In Drosophila melanogaster, the transcription factor Zelda plays an essential role in widespread activation of the zygotic genome. While Zelda expression is required both maternally and zygotically, the mechanisms by which it functions to remodel the embryonic genome and prepare the embryo for development remain unclear. Using Cas9-mediated genome editing to generate targeted mutations in the endogenous zelda locus, we determined the functional relevance of protein domains conserved amongst Zelda orthologs. We showed that neither a conserved N-terminal zinc finger nor an acidic patch were required for activity. Similarly, a previously identified splice isoform of zelda is dispensable for viability. By contrast, we identified a highly conserved zinc-finger domain that is essential for the maternal, but not zygotic functions of Zelda. Animals homozygous for mutations in this domain survived to adulthood, but embryos inheriting these loss-of-function alleles from their mothers died late in embryogenesis. These mutations did not interfere with the capacity of Zelda to activate transcription in cell culture. Unexpectedly, these mutations generated a hyperactive form of the protein and enhanced Zelda-dependent gene expression. These data have defined a protein domain critical for controlling Zelda activity during the MZT, but dispensable for its roles later in development, for the first time separating the maternal and zygotic requirements for Zelda. This demonstrates that highly regulated levels of Zelda activity are required for establishing the developmental program during the MZT. We propose that tightly regulated gene expression is essential to navigate the MZT and that failure to precisely execute this developmental program leads to embryonic lethality.

Following fertilization, the one-celled zygote must be rapidly reprogrammed to enable the development of a new, unique organism. During these initial stages of development there is little or no transcription of the zygotic genome, and maternally deposited products control this process. Among the essential maternal products are mRNAs that encode transcription factors required for preparing the zygotic genome for transcriptional activation. This ensures that there is a precisely coordinated hand-off from maternal to zygotic control. In Drosophila melanogaster, the transcription factor Zelda is essential for activating the zygotic genome and coupling this activation to the degradation of the maternally deposited products. Nonetheless, the mechanism by which Zelda functions remains unclear. Here we used Cas9-mediated genome engineering to determine the functional requirements for highly conserved domains within Zelda. We identified a domain required specifically for Zelda’s role in reprogramming the early embryonic genome, but not essential for its functions later in development. Surprisingly, this domain restricts the ability of Zelda to activate transcription. These data demonstrate that Zelda activity is tightly regulated, and we propose that precise regulation of both the timing and levels of genome activation is required for the embryo to successfully transition from maternal to zygotic control.

Evolution and multiple roles of the Pancrustacea specific transcription factor zelda in insects

Gene regulatory networks (GRNs) evolve as a result of the coevolutionary processes acting on transcription factors (TFs) and the cis-regulatory modules they bind. The zinc-finger TF zelda (zld) is essential for the maternal-to-zygotic transition (MZT) in Drosophila melanogaster, where it directly binds over thousand cis-regulatory modules to regulate chromatin accessibility. D. melanogaster displays a long germ type of embryonic development, where all segments are simultaneously generated along the whole egg. However, it remains unclear if zld is also involved in the MZT of short-germ insects (including those from basal lineages) or in other biological processes. Here we show that zld is an innovation of the Pancrustacea lineage, being absent in more distant arthropods (e.g. chelicerates) and other organisms. To better understand zld´s ancestral function, we thoroughly investigated its roles in a short-germ beetle, Tribolium castaneum, using molecular biology and computational approaches. Our results demonstrate roles for zld not only during the MZT, but also in posterior segmentation and patterning of imaginal disc derived structures. Further, we also demonstrate that zld is critical for posterior segmentation in the hemipteran Rhodnius prolixus, indicating this function predates the origin of holometabolous insects and was subsequently lost in long-germ insects. Our results unveil new roles of zld in different biological contexts and suggest that changes in expression of zld (and probably other major TFs) are critical in the evolution of insect GRNs.

Pioneer transcription factors (TFs) are considered the first regulators of chromatin accessibility in fruit flies and vertebrates, modulating the expression of a large number of target genes. In fruit flies, zelda resembles a pioneer TF, being essential during early embryogenesis. However, the evolutionary origins and ancestral functions of zelda remain largely unknown. Through a number of gene silencing, microscopy and evolutionary analysis, the present work shows that zelda is an innovation of the Pancrustacea lineage, governing not only the MZT in the short-germ insect Tribolium castaneum, but also posterior segmentation and post-embryonic patterning of imaginal disc derived structures such as wings, legs and antennae. Further, zelda regulation of posterior segmentation predates the origin of insects with complete metamorphosis (holometabolous), as supported by gene silencing experiments in the kissing bug Rhodnius prolixus. We hypothesize that the emergence of zelda contributed to the evolution of gene regulatory networks and new morphological structures of insects.

Transcriptome Analysis of Honeybee (Apis Mellifera) Haploid and Diploid Embryos Reveals Early Zygotic Transcription during Cleavage

In honeybees, the haplodiploid sex determination system promotes a unique embryogenesis process wherein females develop from fertilized eggs and males develop from unfertilized eggs. However, the developmental strategies of honeybees during early embryogenesis are virtually unknown. Similar to most animals, the honeybee oocytes are supplied with proteins and regulatory elements that support early embryogenesis. As the embryo develops, the zygotic genome is activated and zygotic products gradually replace the preloaded maternal material. The analysis of small RNA and mRNA libraries of mature oocytes and embryos originated from fertilized and unfertilized eggs has allowed us to explore the gene expression dynamics in the first steps of development and during the maternal-to-zygotic transition (MZT). We localized a short sequence motif identified as TAGteam motif and hypothesized to play a similar role in honeybees as in fruit flies, which includes the timing of early zygotic expression (MZT), a function sustained by the presence of the zelda ortholog, which is the main regulator of genome activation. Predicted microRNA (miRNA)-target interactions indicated that there were specific regulators of haploid and diploid embryonic development and an overlap of maternal and zygotic gene expression during the early steps of embryogenesis. Although a number of functions are highly conserved during the early steps of honeybee embryogenesis, the results showed that zygotic genome activation occurs earlier in honeybees than in Drosophila based on the presence of three primary miRNAs (pri-miRNAs) (ame-mir-375, ame-mir-34 and ame-mir-263b) during the cleavage stage in haploid and diploid embryonic development.

Increased production of piRNAs from euchromatic clusters and genes in Anopheles gambiae compared with Drosophila melanogaster

BACKGROUND

Specific genomic loci, termed Piwi-interacting RNA (piRNA) clusters, manufacture piRNAs that serve as guides for the inactivation of complementary transposable elements (TEs). The piRNA pathway has been accurately detailed in Drosophila melanogaster, while it remains poorly examined in other insects. This pathway is increasingly recognized as critical for germline development and reproduction. Understanding of the piRNA functions in mosquitoes could offer an opportunity for disease vector control by the reduction of their reproductive potential.

RESULTS

To analyze the similarities and differences in this pathway between Drosophila and mosquito, we performed an in-depth analysis of the genomic loci producing piRNAs and their targets in the African malaria vector Anopheles gambiae. We identified 187 piRNA clusters in the An. gambiae genome and 155 piRNA clusters in the D. melanogaster genome. We demonstrate that many more piRNA clusters in the mosquito compared with the fruit fly are uni-directionally transcribed and are located outside pericentromeric heterochromatin. About 11 % of the An. gambiae piRNA population map to gene transcripts. This is a noticeable increase compared with the ~6 % of the piRNA population mapped to genes in D. melanogaster. A subset of the piRNA-enriched genes in An. gambiae has functions related to reproduction and development. At least 24 and 65 % of the mapped piRNAs correspond to genomic TE sequences in An. gambiae and D. melanogaster, respectively. DNA transposons and non-LTR retrotransposons are more abundant in An. gambiae, while LTR retrotransposons are more abundant in D. melanogaster. Yet, piRNAs predominantly target LTR retrotransposons in both species, which may point to a distinct feature of these elements compared to the other classes of TEs concerning their silencing by the piRNA pathway.

CONCLUSIONS

Here, we demonstrate that piRNA-producing loci have more ubiquitous distribution in the An. gambiae genome than in the genome of D. melanogaster. Also, protein-coding genes have an increased role in production of piRNAs in the germline of this mosquito. Genes involved in germline and embryonic development of An. gambiae generate a substantial portion of piRNAs, suggesting a role of the piRNA pathway in the epigenetic regulation of the reproductive processes in the African malaria vector.

SEX DETERMINATION. A male-determining factor in the mosquito Aedes aegypti

Sex determination in the mosquito Aedes aegypti is governed by a dominant male-determining factor (M factor) located within a Y chromosome-like region called the M locus. Here, we show that an M-locus gene, Nix, functions as an M factor in A. aegypti. Nix exhibits persistent M linkage and early embryonic expression, two characteristics required of an M factor. Nix knockout with clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 resulted in largely feminized genetic males and the production of female isoforms of two key regulators of sexual differentiation: doublesex and fruitless. Ectopic expression of Nix resulted in genetic females with nearly complete male genitalia. Thus, Nix is both required and sufficient to initiate male development. This study provides a foundation for mosquito control strategies that convert female mosquitoes into harmless males.

MicroRNAs of two medically important mosquito species: Aedes aegypti and Anopheles stephensi

MicroRNAs (miRNAs) are endogenous, single-stranded small RNAs that have important regulatory functions at the post-transcriptional level. In the present study, we characterize miRNAs in two divergent mosquito species, Aedes aegypti and Anopheles stephensi, through deep sequencing of small RNAs spanning all developmental stages. We discovered eight novel miRNAs in Ae. aegypti and 20 novel miRNAs in An. stephensi, which enabled the first systematic analysis of miRNA evolution in mosquitos. We traced the phylogenetic history of all miRNAs in both species and report a rate of 0.055-0.13 miRNA net gain per million years. Most novel miRNAs originate de novo. Duplications that produced miRNA clusters and families are more common in Ae. aegypti than in An. stephensi. We also identified arm-switch as a source of new miRNAs. Expression profile analysis identified mosquito-specific miRNAs that showed strong stage-specific expression in one or both lineages. For example, the aae-miR-2941/2946 family represents the most abundant maternally deposited and zygotically transcribed miRNAs in Ae. aegypti. miR-2943 is a highly expressed zygotic miRNA in both Ae. aegypti and An. stephensi. Such information provides the basis from which to study the function of these miRNAs in biology common to all mosquitos or unique to one particular lineage.

How gene expression in fast-proliferating cells keeps pace

The development of living organisms requires a precise coordination of all basic cellular processes, in space and time. Early embryogenesis of most species with externally deposited eggs starts with a series of extremely fast cleavage cycles. These divisions have a strong influence on gene expression as mitosis represses transcription and pre-mRNA processing. In this review, we will describe the distinct adaptations for efficient gene expression and discuss the emerging role of the multifunctional NineTeen Complex (NTC) in gene expression and genomic stability during fast proliferation.

The hub protein loquacious connects the microRNA and short interfering RNA pathways in mosquitoes

Aedes aegypti mosquitoes vector several arboviruses of global health significance, including dengue viruses and chikungunya virus. RNA interference (RNAi) plays an important role in antiviral immunity, gene regulation and protection from transposable elements. Double-stranded RNA binding proteins (dsRBPs) are important for efficient RNAi; in Drosophila functional specialization of the miRNA, endo-siRNA and exo-siRNA pathway is aided by the dsRBPs Loquacious (Loqs-PB, Loqs-PD) and R2D2, respectively. However, this functional specialization has not been investigated in other dipterans. We were unable to detect Loqs-PD in Ae. aegypti; analysis of other dipteran genomes demonstrated that this isoform is not conserved outside of Drosophila. Overexpression experiments and small RNA sequencing following depletion of each dsRBP revealed that R2D2 and Loqs-PA cooperate non-redundantly in siRNA production, and that these proteins exhibit an inhibitory effect on miRNA levels. Conversely, Loqs-PB alone interacted with mosquito dicer-1 and was essential for full miRNA production. Mosquito Loqs interacted with both argonaute 1 and 2 in a manner independent of its interactions with dicer. We conclude that the functional specialization of Loqs-PD in Drosophila is a recently derived trait, and that in other dipterans, including the medically important mosquitoes, Loqs-PA participates in both the miRNA and endo-siRNA based pathways.

Transcriptomic study of the red palm weevil Rhynchophorus ferrugineus embryogenesis

The red palm weevil (RPW), Rhynchophorus ferrugineus (Coleoptera: Curculionidae), is an invasive, concealed and destructive tissue borer, and it becomes a lethal pest of the palm family of plants and has been reported to attack 20 palm species around the globe. Here we report a systematic transcriptomic study on embryogenesis of RPW, where we analyze the transcriptomes across five developmental stages of RPW embryogenesis, involving four embryonic stages (E1, E2, E3 and E4) and one larval stage (L1). Using the RNA-seq and next-generation platforms, we generated 80 to 91 million reads for each library and assemble 22 532 genes that are expressed at different embryonic stages. Among the total transcripts from the five embryonic development stages, we found that 30.45 % are differentially expressed, 10.10 % show stage-specificity and even a larger fraction, 62.88 %, exhibit constitutive expression in all the stages. We also analyzes the expression dynamics of several conserved signaling pathways (such as Hedgehog, JAK-STAT, Notch, TGF-β, Ras/MAPK and Wnt), as well as key developmental genes, including those related to apoptosis, axis formation, Hox complex, neurogenesis and segmentation. The datasets provide an essential resource for gene annotation and RPW functional genomics, including studies by using tools and concepts from multiple disciplines, such as development, physiology, biochemistry, molecular biology and genetics.

Functional characterization of calliphorid cell death genes and cellularization gene promoters for controlling gene expression and cell viability in early embryos

The New World screwworm fly, Cochliomyia hominivorax, and the Australian sheep blow fly, Lucilia cuprina, are major pests of livestock. The sterile insect technique was used to eradicate C. hominivorax from North and Central America. This involved area-wide releases of male and female flies that had been sterilized by radiation. Genetic systems have been developed for making 'male-only' strains that would improve the efficiency of genetic control of insect pests. One system involves induction of female lethality in embryos through activation of a pro-apoptotic gene by the tetracycline-dependent transactivator. Sex-specific expression is achieved using an intron from the transformer gene, which we previously isolated from several calliphorids. In the present study, we report the isolation of the promoters from the C. hominivorax slam and Lucilia sericata bnk cellularization genes and show that these promoters can drive expression of a GFP reporter gene in early embryos of transgenic L. cuprina. Additionally, we report the isolation of the L. sericata pro-apoptotic hid and rpr genes, identify conserved motifs in the encoded proteins and determine the relative expression of these genes at different stages of development. We show that widespread expression of the L. sericata pro-apoptotic genes was lethal in Drosophila melanogaster. The isolated gene promoters and pro-apoptotic genes could potentially be used to build transgenic embryonic sexing strains of calliphorid livestock pests.

Transcriptional activation is a conserved feature of the early embryonic factor Zelda that requires a cluster of four zinc fingers for DNA binding and a low-complexity activation domain

Delayed transcriptional activation of the zygotic genome is a nearly universal phenomenon in metazoans. Immediately following fertilization, development is controlled by maternally deposited products, and it is not until later stages that widespread activation of the zygotic genome occurs. Although the mechanisms driving this genome activation are currently unknown, the transcriptional activator Zelda (ZLD) has been shown to be instrumental in driving this process in Drosophila melanogaster. Here we define functional domains of ZLD required for both DNA binding and transcriptional activation. We show that the C-terminal cluster of four zinc fingers mediates binding to TAGteam DNA elements in the promoters of early expressed genes. All four zinc fingers are required for this activity, and splice isoforms lacking three of the four zinc fingers fail to activate transcription. These truncated splice isoforms dominantly suppress activation by the full-length, embryonically expressed isoform. We map the transcriptional activation domain of ZLD to a central region characterized by low complexity. Despite relatively little sequence conservation within this domain, ZLD orthologs from Drosophila virilis, Anopheles gambiae, and Nasonia vitripennis activate transcription in D. melanogaster cells. Transcriptional activation by these ZLD orthologs suggests that ZLD functions through conserved interactions with a protein cofactor(s). We have identified distinct DNA-binding and activation domains within the critical transcription factor ZLD that controls the initial activation of the zygotic genome.

Maternal germline-specific genes in the Asian malaria mosquito Anopheles stephensi: characterization and application for disease control

Anopheles stephensi is a principal vector of urban malaria on the Indian subcontinent and an emerging model for molecular and genetic studies of mosquito biology. To enhance our understanding of female mosquito reproduction, and to develop new tools for basic research and for genetic strategies to control mosquito-borne infectious diseases, we identified 79 genes that displayed previtellogenic germline-specific expression based on RNA-Seq data generated from 11 life stage-specific and sex-specific samples. Analysis of this gene set provided insights into the biology and evolution of female reproduction. Promoters from two of these candidates, vitellogenin receptor and nanos, were used in independent transgenic cassettes for the expression of artificial microRNAs against suspected mosquito maternal-effect genes, discontinuous actin hexagon and myd88. We show these promoters have early germline-specific expression and demonstrate 73% and 42% knockdown of myd88 and discontinuous actin hexagon mRNA in ovaries 48 hr after blood meal, respectively. Additionally, we demonstrate maternal-specific delivery of mRNA and protein to progeny embryos. We discuss the application of this system of maternal delivery of mRNA/miRNA/protein in research on mosquito reproduction and embryonic development, and for the development of a gene drive system based on maternal-effect dominant embryonic arrest.

Introns and gene expression: cellular constraints, transcriptional regulation, and evolutionary consequences

A gene's “expression profile” denotes the number of transcripts present relative to all other transcripts. The overall rate of transcript production is determined by transcription and RNA processing rates. While the speed of elongating RNA polymerase II has been characterized for many different genes and organisms, gene-architectural features – primarily the number and length of exons and introns – have recently emerged as important regulatory players. Several new studies indicate that rapidly cycling cells constrain gene-architecture toward short genes with a few introns, allowing efficient expression during short cell cycles. In contrast, longer genes with long introns exhibit delayed expression, which can serve as timing mechanisms for patterning processes. These findings indicate that cell cycle constraints drive the evolution of gene-architecture and shape the transcriptome of a given cell type. Furthermore, a tendency for short genes to be evolutionarily young hints at links between cellular constraints and the evolution of animal ontogeny.

Insights into the preservation of the homomorphic sex-determining chromosome of Aedes aegypti from the discovery of a male-biased gene tightly linked to the M-locus

The preservation of a homomorphic sex-determining chromosome in some organisms without transformation into a heteromorphic sex chromosome is a long-standing enigma in evolutionary biology. A dominant sex-determining locus (or M-locus) in an undifferentiated homomorphic chromosome confers the male phenotype in the yellow fever mosquito Aedes aegypti. Genetic evidence suggests that the M-locus is in a nonrecombining region. However, the molecular nature of the M-locus has not been characterized. Using a recently developed approach based on Illumina sequencing of male and female genomic DNA, we identified a novel gene, myo-sex, that is present almost exclusively in the male genome but can sporadically be found in the female genome due to recombination. For simplicity, we define sequences that are primarily found in the male genome as male-biased. Fluorescence in situ hybridization (FISH) on A. aegypti chromosomes demonstrated that the myo-sex probe localized to region 1q21, the established location of the M-locus. Myo-sex is a duplicated myosin heavy chain gene that is highly expressed in the pupa and adult male. Myo-sex shares 83% nucleotide identity and 97% amino acid identity with its closest autosomal paralog, consistent with ancient duplication followed by strong purifying selection. Compared with males, myo-sex is expressed at very low levels in the females that acquired it, indicating that myo-sex may be sexually antagonistic. This study establishes a framework to discover male-biased sequences within a homomorphic sex-determining chromosome and offers new insights into the evolutionary forces that have impeded the expansion of the nonrecombining M-locus in A. aegypti.

Zygotic genome activation during the maternal-to-zygotic transition

Embryogenesis depends on a highly coordinated cascade of genetically encoded events. In animals, maternal factors contributed by the egg cytoplasm initially control development, whereas the zygotic nuclear genome is quiescent. Subsequently, the genome is activated, embryonic gene products are mobilized, and maternal factors are cleared. This transfer of developmental control is called the maternal-to-zygotic transition (MZT). In this review, we discuss recent advances toward understanding the scope, timing, and mechanisms that underlie zygotic genome activation at the MZT in animals. We describe high-throughput techniques to measure the embryonic transcriptome and explore how regulation of the cell cycle, chromatin, and transcription factors together elicits specific patterns of embryonic gene expression. Finally, we illustrate the interplay between zygotic transcription and maternal clearance and show how these two activities combine to reprogram two terminally differentiated gametes into a totipotent embryo.

Requirement for highly efficient pre-mRNA splicing during Drosophila early embryonic development

Drosophila syncytial nuclear divisions limit transcription unit size of early zygotic genes. As mitosis inhibits not only transcription, but also pre-mRNA splicing, we reasoned that constraints on splicing were likely to exist in the early embryo, being splicing avoidance a possible explanation why most early zygotic genes are intronless. We isolated two mutant alleles for a subunit of the NTC/Prp19 complexes, which specifically impaired pre-mRNA splicing of early zygotic but not maternally encoded transcripts. We hypothesized that the requirements for pre-mRNA splicing efficiency were likely to vary during development. Ectopic maternal expression of an early zygotic pre-mRNA was sufficient to suppress its splicing defects in the mutant background. Furthermore, a small early zygotic transcript with multiple introns was poorly spliced in wild-type embryos. Our findings demonstrate for the first time the existence of a developmental pre-requisite for highly efficient splicing during Drosophila early embryonic development and suggest in highly proliferative tissues a need for coordination between cell cycle and gene architecture to ensure correct gene expression and avoid abnormally processed transcripts. DOI: http://dx.doi.org/10.7554/eLife.02181.001.

A selfish gene chastened: Tribolium castaneum Medea M4 is silenced by a complementary gene

Maternal-effect dominant embryonic arrest (Medea) of Tribolium castaneum are autosomal factors that act maternally to cause the death of any progeny that do not inherit them. This selfish behavior is thought to result from a maternally expressed poison and zygotically expressed antidote. Medea factors and the hybrid incompatibility factor, H, have a negative interaction consistent with complementary genes of the Dobzhansky-Muller model for post-zygotic isolation. This negative interaction may result from H suppression of Medea zygotic antidote, leaving zygotes incompletely protected from maternal poison. I report here a test of the hypothesis that H also suppresses the Medea maternal poison. Viable F1 females were generated from a cross of Medea M4 strain males to H strain females. These females, heterozygous for both M4 and H, failed to express M4 maternal lethal activity when crossed to their male sibs. Transmission of non-M4 homologues from these females was confirmed using a dominant transgenic enhanced green fluorescent protein eye color marker, tightly linked in cis to M4. M4 beetles, lacking H, were selected from the F2 population. Female descendants of these clearly expressed M4 maternal lethal activity, indicating restoration of this activity after H was segregated away. I conclude that H, or a factor tightly linked to H, suppresses Medea M4 maternal poison.

Integrated proteomic and transcriptomic analysis of the Aedes aegypti eggshell

Background

Mosquito eggshells show remarkable diversity in physical properties and structure consistent with adaptations to the wide variety of environments exploited by these insects. We applied proteomic, transcriptomic, and hybridization in situ techniques to identify gene products and pathways that participate in the assembly of the Aedes aegypti eggshell. Aedes aegypti population density is low during cold and dry seasons and increases immediately after rainfall. The survival of embryos through unfavorable periods is a key factor in the persistence of their populations. The work described here supports integrated vector control approaches that target eggshell formation and result in Ae. aegypti drought-intolerant phenotypes for public health initiatives directed to reduce mosquito-borne diseases.

Results

A total of 130 proteins were identified from the combined mass spectrometric analyses of eggshell preparations.

Conclusions

Classification of proteins according to their known and putative functions revealed the complexity of the eggshell structure. Three novel Ae. aegypti vitelline membrane proteins were discovered. Odorant-binding and cysteine-rich proteins that may be structural components of the eggshell were identified. Enzymes with peroxidase, laccase and phenoloxidase activities also were identified, and their likely involvements in cross-linking reactions that stabilize the eggshell structure are discussed.

Targeted genome editing in Aedes aegypti using TALENs

The Culicine mosquito, Aedes aegypti, is both a major vector of arthropod-borne viruses (arboviruses) and a genetic model organism for arbovirus transmission. TALE nucleases (TALENs), a group of artificial enzymes capable of generating site-specific DNA lesions, consist of a non-specific FokI endonuclease cleavage domain fused to an engineered DNA binding domain specific to a target site. While TALENs have become an important tool for targeted gene disruption in a variety of organisms, application to the mosquito genome is a new approach. We recently described the use of TALENs to perform heritable genetic disruptions in A. aegypti. Here, we provide detailed methods that will allow other research laboratories to capitalize on the potential of this technology for understanding mosquito gene function. We describe target site selection, transient embryo-based assays to rapidly assess TALEN activity, embryonic microinjection and downstream screening steps to identify target site mutations.

The earliest transcribed zygotic genes are short, newly evolved, and different across species

The transition from maternal to zygotic control is fundamental to the life cycle of all multicellular organisms. It is widely believed that genomes are transcriptionally inactive from fertilization until zygotic genome activation (ZGA). Thus, the earliest genes expressed probably support the rapid cell divisions that precede morphogenesis and, if so, might be evolutionarily conserved. Here, we identify the earliest zygotic transcripts in the zebrafish, Danio rerio, through metabolic labeling and purification of RNA from staged embryos. Surprisingly, the mitochondrial genome was highly active from the one-cell stage onwards, showing that significant transcriptional activity exists at fertilization. We show that 592 nuclear genes become active when cell cycles are still only 15 min long, confining expression to relatively short genes. Furthermore, these zygotic genes are evolutionarily younger than those expressed at other developmental stages. Comparison of fish, fly, and mouse data revealed different sets of genes expressed at ZGA. This species specificity uncovers an evolutionary plasticity in early embryogenesis that probably confers substantial adaptive potential.

Transcriptome analysis of Anopheles stephensi embryo using expressed sequence tags

Germ band retraction (GBR) stage is one of the important stages during insect development. It is associated with an extensive epithelial morphogenesis and may also be pivotal in generation of morphological diversity in insects. Despite its importance, only a handful of studies report the transcriptome repertoire of this stage in insects. Here, we report generation, annotation and analysis of ESTs from the embryonic stage (16-22 h post fertilization) of laboratoryreared Anopheles stephensi mosquitoes. A total of 1002 contigs were obtained upon clustering of 1140 high-quality ESTs, which demonstrates an astonishingly low transcript redundancy (12.1 percent). Putative functions were assigned only to 213 contigs (21 percent), comprising mainly of transcripts encoding protein synthesis machinery. Approximately 78 percent of the transcripts remain uncharacterized, illustrating a lack of sequence information about the genes expressed in the embryonic stages of mosquitoes. This study highlights several novel transcripts, which apart from insect development, may significantly contribute to the essential biological complexity underlying insect viability in adverse environments. Nonetheless, the generated sequence information from this work provides a comprehensive resource for genome annotation, microarray development, phylogenetic analysis and other molecular biology applications in entomology.

The developmental transcriptome of the mosquito Aedes aegypti, an invasive species and major arbovirus vector

Mosquitoes are vectors of a number of important human and animal diseases. The development of novel vector control strategies requires a thorough understanding of mosquito biology. To facilitate this, we used RNA-seq to identify novel genes and provide the first high-resolution view of the transcriptome throughout development and in response to blood feeding in a mosquito vector of human disease, Aedes aegypti, the primary vector for Dengue and yellow fever. We characterized mRNA expression at 34 distinct time points throughout Aedes development, including adult somatic and germline tissues, by using polyA+ RNA-seq. We identify a total of 14,238 novel new transcribed regions corresponding to 12,597 new loci, as well as many novel transcript isoforms of previously annotated genes. Altogether these results increase the annotated fraction of the transcribed genome into long polyA+ RNAs by more than twofold. We also identified a number of patterns of shared gene expression, as well as genes and/or exons expressed sex-specifically or sex-differentially. Expression profiles of small RNAs in ovaries, early embryos, testes, and adult male and female somatic tissues also were determined, resulting in the identification of 38 new Aedes-specific miRNAs, and ~291,000 small RNA new transcribed regions, many of which are likely to be endogenous small-interfering RNAs and Piwi-interacting RNAs. Genes of potential interest for transgene-based vector control strategies also are highlighted. Our data have been incorporated into a user-friendly genome browser located at www.Aedes.caltech.edu, with relevant links to Vectorbase (www.vectorbase.org)

A unique Y gene in the Asian malaria mosquito Anopheles stephensi encodes a small lysine-rich protein and is transcribed at the onset of embryonic development

In many organisms the Y chromosome initiates sex determination and regulates male fertility and mating behaviour. However, molecular characterization of Y genes is rare outside of a few model species because it is difficult to clone and analyse repeat-rich heterochromatic Y sequences. In insects, Y genes are only well characterized in a small number of Drosophila species. Here we report the discovery of GUY1 (gene unique to the Y), a gene unique to the Y chromosome in the Asian malaria mosquito, Anopheles stephensi, using an approach that compares Illumina sequences separately obtained from male and female genomic DNA. Experimental evidence confirmed that GUY1 is a single copy gene found only on the Y chromosome. GUY1 is transcribed at the very onset of zygotic transcription and encodes a small lysine-rich protein that forms two alpha helices and shows DNA-binding properties. Interestingly, three helix-loop-helix proteins are key factors that determine sex in the early embryo in Drosophila melanogaster. Single embryo analysis indicated that GUY1 is only transcribed in male embryos and that the GUY1 promoter is functional in the early embryos. GUY1 may be used as a paternally inherited molecular marker. Further investigation of GUY1 will contribute to the genetic approaches to control mosquito-borne diseases.

An essential role for zygotic expression in the pre-cellular Drosophila embryo

The Drosophila embryo proceeds through thirteen mitotic divisions as a syncytium. Its nuclei distribute in the embryo's interior during the first six divisions, dividing synchronously with a cycle time of less than ten minutes. After seven divisions (nuclear cycle 8), the syncytial blastoderm forms as the nuclei approach the embryo surface and slow their cycle time; subsequent divisions proceed in waves that initiate at the poles. Because genetic studies have not identified zygotic mutants that affect the early divisions and because transcription has not been detected before cycle 8, the early, pre-blastoderm embryo has been considered to rely entirely on maternal contributions and to be transcriptionally silent. Our studies identified several abnormal phenotypes in live engrailed (en) mutant embryos prior to cycle 8, as well as a small group of genes that are transcribed in embryos prior to cycle 7. Nuclei in en embryos divide asynchronously, an abnormality that was detected as early as nuclear cycle 2–3. Anti-En antibody detected nuclear En protein in embryos at cycle 2, and expression of an En:GFP fusion protein encoded in the paternal genome was also detected in cycle 2 nuclei. These findings demonstrate that the Drosophila embryo is functionally competent for gene expression prior to the onset of its rapid nuclear divisions and that the embryo requires functions that are expressed in the zygote in order to faithfully prosecute its early, pre-cellularization mitotic cycles.

Genetic studies identified many genes that are required during Drosophila oogenesis to endow the embryo with structures and components it will need to develop; they have also identified many genes that the embryo must express. However, measures of transcription have detected zygotic transcripts only after seven nuclear divisions, and many studies have concluded that zygotic mutants do not affect embryos prior to cellularization. The model that has emerged is that the earliest stages of embryogenesis rely solely on maternal stores and do not receive input from the zygotic genome. The fact that the embryo's nuclei divide rapidly with a cycling time of less than ten minutes has been interpreted to support this model, because it has been assumed that the nuclear cycle is too short for productive gene expression. Using sensitive measures of transcription and histological procedures that detect subtle differences, we found evidence for expression as early as nuclear cycle 2, and we identified a requirement for zygotic gene expression in embryos with just 2–4 nuclei. These findings challenge the idea that the Drosophila embryo is entirely pre-programmed and that its early development is under exclusive maternal control.

Chromatin-linked determinants of zygotic genome activation

The merging of the maternal and paternal genomes into a single pronucleus after fertilization is accompanied by a remarkable reconfiguration of chromatin in the newly formed zygote. The first stages of embryonic chromatin remodeling take place in the absence of ongoing transcription, during a species-specific developmental time-frame. Once post-fertilization chromatin states are organized, zygotic genome activation (ZGA) is initiated, and embryonic transcripts gradually take control of development. We review here transitions in chromatin modifications associated with the onset of ZGA, and the role of transcription factors and DNA motifs in the regulation of ZGA. We propose a model of sequential chromatin remodeling events preceding ZGA, leading to the onset of embryonic transcription.