The regulatory genome of the malaria vector Anopheles gambiae: integrating chromatin accessibility and gene expression

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

Mosquito salivary apyrase regulates blood meal hemostasis and facilitates malaria parasite transmission

The evolution of hematophagy involves a series of adaptations that allow blood-feeding insects to access and consume blood efficiently while managing and circumventing the host's hemostatic and immune responses. Mosquito, and other insects, utilize salivary proteins to regulate these responses at the bite site during and after blood feeding. We investigated the function of Anopheles gambiae salivary apyrase (AgApyrase) in regulating hemostasis in the mosquito blood meal and in Plasmodium transmission. Our results demonstrate that salivary apyrase, a known inhibitor of platelet aggregation, interacts with and activates tissue plasminogen activator, facilitating the conversion of plasminogen to plasmin, a human protease that degrades fibrin and facilitates Plasmodium transmission. We show that mosquitoes ingest a substantial amount of apyrase during blood feeding, which reduces coagulation in the blood meal by enhancing fibrin degradation and inhibiting platelet aggregation. AgApyrase significantly enhanced Plasmodium infection in the mosquito midgut, whereas AgApyrase immunization inhibited Plasmodium mosquito infection and sporozoite transmission. This study highlights a pivotal role for mosquito salivary apyrase for regulation of hemostasis in the mosquito blood meal and for Plasmodium transmission to mosquitoes and to the mammalian host, underscoring the potential for strategies to prevent malaria transmission.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele-specific expression

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb-resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and CRM sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele-specific expression (ASE) in hybrids of insecticide susceptible and resistant strains, suggesting cis-regulation is an important mechanism of gene expression regulation in A. gambiae. The genes showing ASE included a higher proportion of Anopheles-specific genes on average younger than genes with balanced allelic expression.

Peptides with Antimicrobial Activity in the Saliva of the Malaria Vector Anopheles coluzzii

Mosquito saliva plays a crucial physiological role in both sugar and blood feeding by helping sugar digestion and exerting antihemostatic functions. During meal acquisition, mosquitoes are exposed to the internalization of external microbes. Since mosquitoes reingest significant amounts of saliva during feeding, we hypothesized that salivary antimicrobial components may participate in the protection of mouthparts, the crop, and the gut by inhibiting bacterial growth. To identify novel potential antimicrobials from mosquito saliva, we selected 11 candidates from Anopheles coluzzii salivary transcriptomic datasets and obtained them either using a cell-free transcription/translation expression system or, when feasible, via chemical synthesis. Hyp6.2 and hyp13, which were predicted to be produced as propeptides and cleaved in shorter mature forms, showed the most interesting results in bacterial growth inhibition assays. Hyp6.2 (putative mature form, 35 amino acid residues) significantly inhibited the growth of Gram-positive (Staphylococcus aureus) and Gram-negative (Escherichia coli and Serratia marcescens) bacteria. Hyp13 (short form, 19 amino acid residues) dose-dependently inhibited E. coli and S. marcescens growth, inducing membrane disruption in both Gram-positive and Gram-negative bacteria as indicated with scanning electron microscopy. In conclusion, we identified two A. coluzzii salivary peptides inhibiting Gram-positive and Gram-negative bacteria growth and possibly contributing to the protection of mosquito mouthparts and digestive tracts from microbial infection during and/or after feeding.

Mechanisms of transcriptional regulation in Anopheles gambiae revealed by allele specific expression

Malaria control relies on insecticides targeting the mosquito vector, but this is increasingly compromised by insecticide resistance, which can be achieved by elevated expression of detoxifying enzymes that metabolize the insecticide. In diploid organisms, gene expression is regulated both in cis, by regulatory sequences on the same chromosome, and by trans acting factors, affecting both alleles equally. Differing levels of transcription can be caused by mutations in cis-regulatory modules (CRM), but few of these have been identified in mosquitoes. We crossed bendiocarb resistant and susceptible Anopheles gambiae strains to identify cis-regulated genes that might be responsible for the resistant phenotype using RNAseq, and cis-regulatory module sequences controlling gene expression in insecticide resistance relevant tissues were predicted using machine learning. We found 115 genes showing allele specific expression in hybrids of insecticide susceptible and resistant strains, suggesting cis regulation is an important mechanism of gene expression regulation in Anopheles gambiae. The genes showing allele specific expression included a higher proportion of Anopheles specific genes on average younger than genes those with balanced allelic expression.

Anopheles salivary apyrase regulates blood meal hemostasis and drives malaria parasite transmission

Mosquito salivary proteins play a crucial role in regulating hemostatic responses at the bite site during blood feeding. In this study, we investigate the function of Anopheles gambiae salivary apyrase (AgApyrase) in Plasmodium transmission. Our results demonstrate that salivary apyrase interacts with and activates tissue plasminogen activator, facilitating the conversion of plasminogen to plasmin, a human protein previously shown to be required for Plasmodium transmission. Microscopy imaging shows that mosquitoes ingest a substantial amount of apyrase during blood feeding which reduces coagulation in the blood meal by enhancing fibrin degradation and inhibiting platelet aggregation. Supplementation of Plasmodium infected blood with apyrase significantly enhanced Plasmodium infection in the mosquito midgut. In contrast, AgApyrase immunization inhibited Plasmodium mosquito infection and sporozoite transmission. This study highlights a pivotal role for mosquito salivary apyrase for regulation of hemostasis in the mosquito blood meal and for Plasmodium transmission to mosquitoes and to the mammal host, underscoring the potential for new strategies to prevent malaria transmission.

Use of Insect Promoters in Genetic Engineering to Control Mosquito-Borne Diseases

Mosquito transgenesis and gene-drive technologies provide the basis for developing promising new tools for vector-borne disease prevention by either suppressing wild mosquito populations or reducing their capacity from transmitting pathogens. Many studies of the regulatory DNA and promoters of genes with robust sex-, tissue- and stage-specific expression profiles have supported the development of new tools and strategies that could bring mosquito-borne diseases under control. Although the list of regulatory elements available is significant, only a limited set of those can reliably drive spatial-temporal expression. Here, we review the advances in our ability to express beneficial and other genes in mosquitoes, and highlight the information needed for the development of new mosquito-control and anti-disease strategies.

Chromatin accessibility-based characterisation of brain gene regulatory networks in three distinct honey bee polyphenisms

The honey bee genome has the capacity to produce three phenotypically distinct organisms (two diploid female castes: queen and worker, and a haploid male drone). Previous studies have implicated metabolic flux acting via epigenetic regulation in directing nutrition-driven phenotypic plasticity in the honey bee. However, the cis-acting DNA regulatory elements that establish tissue and polyphenism -specific epigenomes and gene expression programmes, remain unclear. Using a high resolution multiomic approach including assay for transposase-accessible chromatin by sequencing (ATAC-seq), RNA-seq and ChIP-seq, we produce the first genome-wide maps of the regulatory landscape across all three adult honey bee phenotypes identifying > 5000 regulatory regions in queen, 7500 in worker and 6500 in drone, with the vast majority of these sites located within intronic regions. These regions are defined by positive enrichment of H3K27ac and depletion of H3K4me3 and show a positive correlation with gene expression. Using ATAC-seq footprinting we determine queen, worker and drone -specific transcription factor occupancy and uncover novel phenotype-specific regulatory networks identifying two key nuclear receptors that have previously been implicated in caste-determination and adult behavioural maturation in honey bees; ecdysone receptor and ultraspiracle. Collectively, this study provides novel insights into key gene regulatory networks that are associated with these distinct polyphenisms in the honey bee.

The Landscape of Accessible Chromatin during Yak Adipocyte Differentiation

Although significant advancement has been made in the study of adipogenesis, knowledge about how chromatin accessibility regulates yak adipogenesis is lacking. We here described genome-wide dynamic chromatin accessibility in preadipocytes and adipocytes by using the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), and thus revealed the unique characteristics of open chromatin during yak adipocyte differentiation. The chromatin accessibility of preadipocytes and adipocytes exhibited a similar genomic distribution, displaying a preferential location within the intergenic region, intron, and promoter. The pathway enrichment analysis identified that genes with differential chromatin accessibility were involved in adipogenic metabolism regulation pathways, such as the peroxisome proliferator activated receptor-γ (PPAR) signaling pathway, wingless-type MMTV integration site (Wnt) signaling pathway, and extracellular matrix-receptor (ECM–receptor) interaction. Integration of ATAC-seq and mRNA-seq revealed that genes with a high expression were associated with high levels of chromatin accessibility, especially within 1 kb upstream and downstream of the transcription start site. In addition, we identified a series of transcription factors (TFs) related to adipogenesis and created the TF regulatory network, providing the possible interactions between TFs during yak adipogenesis. This study is crucial for advancing the understanding of transcriptional regulatory mechanisms of adipogenesis and provides valuable information for understanding the adaptation of plateau species to high-altitude environments by maintaining whole body homeostasis through fat metabolism.

Identification of a rapidly-spreading triple mutant for high-level metabolic insecticide resistance in Anopheles gambiae provides a real-time molecular diagnostic for antimalarial intervention deployment

Studies of insecticide resistance provide insights into the capacity of populations to show rapid evolutionary responses to contemporary selection. Malaria control remains heavily dependent on pyrethroid insecticides, primarily in long lasting insecticidal nets (LLINs). Resistance in the major malaria vectors has increased in concert with the expansion of LLIN distributions. Identifying genetic mechanisms underlying high-level resistance is crucial for the development and deployment of resistance-breaking tools. Using the Anopheles gambiae 1000 genomes (Ag1000g) data we identified a very recent selective sweep in mosquitoes from Uganda which localized to a cluster of cytochrome P450 genes. Further interrogation revealed a haplotype involving a trio of mutations, a nonsynonymous point mutation in Cyp6p4 (I236M), an upstream insertion of a partial Zanzibar-like transposable element (TE) and a duplication of the Cyp6aa1 gene. The mutations appear to have originated recently in An. gambiae from the Kenya-Uganda border, with stepwise replacement of the double-mutant (Zanzibar-like TE and Cyp6p4-236 M) with the triple-mutant haplotype (including Cyp6aa1 duplication), which has spread into the Democratic Republic of Congo and Tanzania. The triple-mutant haplotype is strongly associated with increased expression of genes able to metabolize pyrethroids and is strongly predictive of resistance to pyrethroids most notably deltamethrin. Importantly, there was increased mortality in mosquitoes carrying the triple-mutation when exposed to nets cotreated with the synergist piperonyl butoxide (PBO). Frequencies of the triple-mutant haplotype remain spatially variable within countries, suggesting an effective marker system to guide deployment decisions for limited supplies of PBO-pyrethroid cotreated LLINs across African countries.

Comprehensive Genomic Discovery of Non-Coding Transcriptional Enhancers in the African Malaria Vector Anopheles coluzzii

Almost all regulation of gene expression in eukaryotic genomes is mediated by the action of distant non-coding transcriptional enhancers upon proximal gene promoters. Enhancer locations cannot be accurately predicted bioinformatically because of the absence of a defined sequence code, and thus functional assays are required for their direct detection. Here we used a massively parallel reporter assay, Self-Transcribing Active Regulatory Region sequencing (STARR-seq), to generate the first comprehensive genome-wide map of enhancers in Anopheles coluzzii, a major African malaria vector in the Gambiae species complex. The screen was carried out by transfecting reporter libraries created from the genomic DNA of 60 wild A. coluzzii from Burkina Faso into A. coluzzii 4a3A cells, in order to functionally query enhancer activity of the natural population within the homologous cellular context. We report a catalog of 3,288 active genomic enhancers that were significant across three biological replicates, 74% of them located in intergenic and intronic regions. The STARR-seq enhancer screen is chromatin-free and thus detects inherent activity of a comprehensive catalog of enhancers that may be restricted in vivo to specific cell types or developmental stages. Testing of a validation panel of enhancer candidates using manual luciferase assays confirmed enhancer function in 26 of 28 (93%) of the candidates over a wide dynamic range of activity from two to at least 16-fold activity above baseline. The enhancers occupy only 0.7% of the genome, and display distinct composition features. The enhancer compartment is significantly enriched for 15 transcription factor binding site signatures, and displays divergence for specific dinucleotide repeats, as compared to matched non-enhancer genomic controls. The genome-wide catalog of A. coluzzii enhancers is publicly available in a simple searchable graphic format. This enhancer catalogue will be valuable in linking genetic and phenotypic variation, in identifying regulatory elements that could be employed in vector manipulation, and in better targeting of chromosome editing to minimize extraneous regulation influences on the introduced sequences.

Importance: Understanding the role of the non-coding regulatory genome in complex disease phenotypes is essential, but even in well-characterized model organisms, identification of regulatory regions within the vast non-coding genome remains a challenge. We used a large-scale assay to generate a genome wide map of transcriptional enhancers. Such a catalogue for the important malaria vector, Anopheles coluzzii, will be an important research tool as the role of non-coding regulatory variation in differential susceptibility to malaria infection is explored and as a public resource for research on this important insect vector of disease.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Distinct mechanisms mediate X chromosome dosage compensation in Anopheles and Drosophila

CRISPR knockout of msl-2 and epigenome analyses in Anopheles reveal that X chromosome dosage compensation in mosquitos and Drosophila is achieved by two different molecular mechanisms.

Sex chromosomes induce potentially deleterious gene expression imbalances that are frequently corrected by dosage compensation (DC). Three distinct molecular strategies to achieve DC have been previously described in nematodes, fruit flies, and mammals. Is this a consequence of distinct genomes, functional or ecological constraints, or random initial commitment to an evolutionary trajectory? Here, we study DC in the malaria mosquito Anopheles gambiae. The Anopheles and Drosophila X chromosomes evolved independently but share a high degree of homology. We find that Anopheles achieves DC by a mechanism distinct from the Drosophila MSL complex–histone H4 lysine 16 acetylation pathway. CRISPR knockout of Anopheles msl-2 leads to embryonic lethality in both sexes. Transcriptome analyses indicate that this phenotype is not a consequence of defective X chromosome DC. By immunofluorescence and ChIP, H4K16ac does not preferentially enrich on the male X. Instead, the mosquito MSL pathway regulates conserved developmental genes. We conclude that a novel mechanism confers X chromosome up-regulation in Anopheles. Our findings highlight the pluralism of gene-dosage buffering mechanisms even under similar genomic and functional constraints.