Chromatin changes in Anopheles gambiae induced by Plasmodium falciparum infection

Conserved transcription factors coordinate synaptic gene expression through repression

Chemical synapses are the primary sites of communication in the nervous system. Synapse formation is a complex process involving hundreds of proteins that must be expressed in two cells at the same time. We find that synaptic genes are broadly and specifically coordinated at the level of transcription across developing nervous systems. How this spatiotemporal coordination is achieved remains an open question. Through genomic and functional studies in Drosophila, we demonstrate corresponding coordination of chromatin accessibility and identify chromatin regulators DEAF1 and CLAMP as broad repressors of synaptic gene expression outside windows of peak synaptogenesis. Disruption of either factor temporally dysregulates synaptic gene expression across neuronal subtypes, leading to excess synapse formation. We further find that DEAF1, which is linked to syndromic intellectual disability, is both necessary and sufficient to constrain synapse formation. Our findings reveal the critical importance of broad temporally coordinated repression of synaptic gene expression in regulating neuronal connectivity and identify two key repressors.

Epigenetic Regulation in Insect-Microbe Interactions

Insects have evolved diverse interactions with a variety of microbes, such as pathogenic fungi, bacteria, and viruses. The immune responses of insect hosts, along with the dynamic infection process of microbes in response to the changing host environment and defenses, require rapid and fine-tuned regulation of gene expression programs. Epigenetic mechanisms, including DNA methylation, histone modifications, and noncoding RNA regulation, play important roles in regulating the expression of genes involved in insect immunity and microbial pathogenicity. This review highlights recent discoveries and insights into epigenetic regulatory mechanisms that modulate insect-microbe interactions. A deeper understanding of these regulatory mechanisms underlying insect-microbe interactions holds promise for the development of novel strategies for biological control of insect pests and mitigation of vector-borne diseases.

Comparison of the effect of bacterial stimulation on the global epigenetic landscape and transcription of immune genes in primarily zoophilic members of the Anopheles gambiae complex (Diptera: Culicidae)

Members of the Anopheles gambiae complex vary in their vector competence, and this is often attributed to behavioural differences. Similarly, there are differences in transmission capabilities of the zoophilic members of this complex despite exhibiting similar behaviours. Therefore, behavioural differences alone cannot fully explain vector competence variation within members of the An. gambiae complex. The immune system of mosquitoes plays a key role in determining susceptibility to parasite infection and consequently transmission capacity. This study aimed to examine variations in the immune response of An. arabiensis, An. merus and An. quadriannulatus, a major, minor, and non-vector respectively. The global epigenetic landscape was characterised and the expression of Defensin-1 and Gambicin was assessed in response to Gram-positive (Streptococcus pyogenes) and Gram-negative (Escherichia coli) bacterial infections. The effect of insecticide resistance in An. arabiensis on these aspects was also assessed. The immune system was stimulated by a blood-borne bacterial supplementation. The 5mC, 5hmC, m6A methylation levels and Histone Acetyl Transferase activity were assessed with commercial ELISA kits. The transcript levels of Defensin-1 and Gambicin were assessed by quantitative Real-Time Polymerase Chain Reaction. Species-specific differences in 5mC and m6A methylation existed both constitutively as well as post immune stimulation. The epigenetic patterns observed in the laboratory strains were largely conserved in F1 offspring of wild-caught adults. The methylation patterns in the major vector typically differed from that of the minor/non-vectors. The differences between insecticide susceptible and resistant An. arabiensis were more reflected in the expression of Defensin-1 and Gambicin. The expression of these peptides differed in the strains only after bacterial stimulation. Anopheles merus and An. quadriannulatus expressed significantly higher levels of antimicrobial peptides, both constitutively and after immune stimulation. These findings suggest molecular variations in the immune response of members of the An. gambiae complex.

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.

Lysine Methylation and Histone Modifications during Cold Stress of Insects: Freeze-Tolerant Eurosta solidaginis and Freeze-Avoiding Epiblema scudderiana

Overwintering survival by insects, whether of the freeze-tolerant or freeze-avoiding types, is typically associated with a strong suppression of metabolic rate (e.g., entry into diapause) that involves the differential expression of many genes with regulation at the transcriptional, translational or post-translational levels. Epigenetic modifications have been suggested to play a vital role in regulating cold responses of insects. However, knowledge of the roles of epigenetic mechanisms in modulating gene expression for winter survival of the larvae of two goldenrod gall formers, the freeze-tolerant dipteran Eurosta solidaginis and the freeze-avoiding lepidopteran Epiblema scudderiana, remain unknown. The current study evaluates the role of cold-induced lysine methylation and histone modifications, with enzymes of lysine methylation (SETD8, SETD7, SUV39H1, SMYD2 and ASH2L), as well as relative levels of histone H3 acetylation (H3K9ac, H3K18ac, H3K27ac, H3K56ac) and methylation (H3K4me1, H3K9me3, H3K36me2) examined in two insects. Significant (p < 0.05) reductions were observed in most of the targets of histone methylation/acetylation for decreasing temperatures of Ep. scudderiana larvae, whereas selected histone methylation/acetylation targets were conversely elevated (p < 0.05) in E. solidaginis, particularly under conditions of 5 °C for 4 h. Histone H3 expression was found to be variable without statistical differences in larval goldenrod gall moths and gall flies. These results provide basic information on the patterns of epigenetic regulation involved in insect cold hardiness.

Genome-wide profiles of H3K9me3, H3K27me3 modifications, and DNA methylation during diapause of Asian corn borer (Ostrinia furnacalis)

Diapause represents a crucial adaptive strategy used by insects to cope with changing environmental conditions. In North China, the Asian corn borer (Ostrinia furnacalis) enters a winter larval diapause stage. Although there is growing evidence implicating epigenetic mechanisms in diapause regulation, it remains unclear whether dynamic genome-wide profiles of epigenetic modifications exist during this process. By investigating multiple histone modifications, we have discovered the essential roles of H3K9me3 and H3K27me3 during diapause of the Asian corn borer. Building upon previous findings in vertebrates highlighting the connection between DNA methylation and repressive histone methylations, we have examined changes in the genome-wide profile of H3K9me3, H3K27me3, and DNA methylation at the nondiapause, prediapause, and diapause stages. Data analysis reveals significant alterations in these three modifications during diapause. Moreover, we observe a correlation between the H3K9me3 and H3K27me3 modification sites during diapause, whereas DNA modifications show little association with either H3K9me3 or H3K27me3. Integrative analysis of epigenome and expression data unveils the relationship between these epigenetic modifications and gene expression levels at corresponding diapause stages. Furthermore, by studying the function of histone modifications on genes known to be important in diapause, especially those involved in the juvenile pathway, we discover that the juvenile hormone pathway lies downstream from H3K9me3 and H3K27me3 histone modifications. Finally, the analysis of gene loci with modified modifications unreported in diapause uncovers novel pathways potentially crucial in diapause regulation. This study provides a valuable resource for future investigations aiming to elucidate the underlying mechanisms of diapause.

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.

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.

Genomic advances in the study of the mosquito vector during avian malaria infection

Invertebrate host–parasite associations are one of the keystones in order to understand vector-borne diseases. The study of these specific interactions provides information not only about how the vector is affected by the parasite at the gene-expression level, but might also reveal mosquito strategies for blocking the transmission of the parasites. A very well-known vector for human malaria is Anopheles gambiae. This mosquito species has been the main focus for genomics studies determining essential key genes and pathways over the course of a malaria infection. However, to-date there is an important knowledge gap concerning other non-mammophilic mosquito species, for example some species from the Culex genera which may transmit avian malaria but also zoonotic pathogens such as West Nile virus. From an evolutionary perspective, these 2 mosquito genera diverged 170 million years ago, hence allowing studies in both species determining evolutionary conserved genes essential during malaria infections, which in turn might help to find key genes for blocking malaria cycle inside the mosquito. Here, we extensively review the current knowledge on key genes and pathways expressed in Anopheles over the course of malaria infections and highlight the importance of conducting genomic investigations for detecting pathways in Culex mosquitoes linked to infection of avian malaria. By pooling this information, we underline the need to increase genomic studies in mosquito–parasite associations, such as the one in Culex–Plasmodium, that can provide a better understanding of the infection dynamics in wildlife and reduce the negative impact on ecosystems.

Notch Signaling in Insect Development: A Simple Pathway with Diverse Functions

Notch signaling is an evolutionarily conserved pathway which functions between adjacent cells to establish their distinct identities. Despite operating in a simple mechanism, Notch signaling plays remarkably diverse roles in development to regulate cell fate determination, organ growth and tissue patterning. While initially discovered and characterized in the model insect Drosophila melanogaster, recent studies across various insect species have revealed the broad involvement of Notch signaling in shaping insect tissues. This review focuses on providing a comprehensive picture regarding the roles of the Notch pathway in insect development. The roles of Notch in the formation and patterning of the insect embryo, wing, leg, ovary and several specific structures, as well as in physiological responses, are summarized. These results are discussed within the developmental context, aiming to deepen our understanding of the diversified functions of the Notch signaling pathway in different insect species.

Single nucleotide polymorphism (SNP) in the doublesex ( dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex ( dsx) gene, involved in mosquitos' somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene ( Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

ChIP-seq profiling of H3K4me3 and H3K27me3 in an invasive insect, Bactrocera dorsalis

Introduction: While it has been suggested that histone modifications can facilitate animal responses to rapidly changing environments, few studies have profiled whole-genome histone modification patterns in invasive species, leaving the regulatory landscape of histone modifications in invasive species unclear. Methods: Here, we screen genome-wide patterns of two important histone modifications, trimethylated Histone H3 Lysine 4 (H3K4me3) and trimethylated Histone H3 Lysine 27 (H3K27me3), in adult thorax muscles of a notorious invasive pest, the Oriental fruit fly Bactrocera dorsalis (Hendel) (Diptera: Tephritidae), using Chromatin Immunoprecipitation with high-throughput sequencing (ChIP-seq). Results: We identified promoters featured by the occupancy of H3K4me3, H3K27me3 or bivalent histone modifications that were respectively annotated with unique genes key to muscle development and structure maintenance. In addition, we found H3K27me3 occupied the entire body of genes, where the average enrichment was almost constant. Transcriptomic analysis indicated that H3K4me3 is associated with active gene transcription, and H3K27me3 is mostly associated with transcriptional repression. Importantly, we identified genes and putative motifs modified by distinct histone modification patterns that may possibly regulate flight activity. Discussion: These findings provide the first evidence of histone modification signature in B. dorsalis, and will be useful for future studies of epigenetic signature in other invasive insect species.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: Malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, involved in mosquitos’ somatic sexually dimorphic traits determination. However, the pathways that trigger the signal of dsx gene exon skipping alternative splicing mechanism in anopheline mosquitoes are not well characterized. This study aims to screen the An. gambiae dsx gene splice site sequences for single-nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

Zika virus infection drives epigenetic modulation of immunity by the histone acetyltransferase CBP of Aedes aegypti

Epigenetic mechanisms are responsible for a wide range of biological phenomena in insects, controlling embryonic development, growth, aging and nutrition. Despite this, the role of epigenetics in shaping insect-pathogen interactions has received little attention. Gene expression in eukaryotes is regulated by histone acetylation/deacetylation, an epigenetic process mediated by histone acetyltransferases (HATs) and histone deacetylases (HDACs). In this study, we explored the role of the Aedes aegypti histone acetyltransferase CBP (AaCBP) after infection with Zika virus (ZIKV), focusing on the two main immune tissues, the midgut and fat body. We showed that the expression and activity of AaCBP could be positively modulated by blood meal and ZIKV infection. Nevertheless, Zika-infected mosquitoes that were silenced for AaCBP revealed a significant reduction in the acetylation of H3K27 (CBP target marker), followed by downmodulation of the expression of immune genes, higher titers of ZIKV and lower survival rates. Importantly, in Zika-infected mosquitoes that were treated with sodium butyrate, a histone deacetylase inhibitor, their capacity to fight virus infection was rescued. Our data point to a direct correlation among histone hyperacetylation by AaCBP, upregulation of antimicrobial peptide genes and increased survival of Zika-infected-A. aegypti.

Characterisation of the epigenetic architecture of the major malaria vector Anopheles arabiensis (Diptera: Culicidae) after treatment with epigenetic modulators and heavy metals

Anopheles arabiensis (a member of the An. gambiae species complex) is a major vector of malaria in sub-Saharan Africa. Despite its disease vector status, there is currently a paucity of epigenetic information for this species. The aim this study was therefore to analyse global epigenetic markers and their response to metal exposure in insecticide susceptible and resistant laboratory strains of An. arabiensis. This was done using commercially available epigenetic marker quantification kits. In order to validate the efficacy of the kits, several kits were assessed to determine whether changes induced by known epigenetic modulators were detectable using these platforms. The efficacy of the dosages used were determined by examining the effect of the dosages used on insecticide resistant phenotypes. Upon confirmation that the dosages used were sufficient to induce a phenotypic change, the effect on epigenetic markers was assessed. Commercial kits were used to quantify 5-methylcysteine (5-mC) and 5-hydroxymethylcysteine (5-hmC) methylation in DNA, m6A methylation in mRNA as well as Histone Acetyl Transferase (HAT) activity. There was a marked difference in the phenotypic response in adult mosquitoes of the insecticide susceptible strain compared to that of its' resistant counterpart. For males and females of the resistant strain, exposure to nucleic acid modifying drugs typically increased their tolerance to insecticides. The patterns of changes in 5-mC methylation by epigenetic modulators was congruent with previous studies which quantified by mass spectrometry. The two strains differed in methylation patterns under control conditions and responded differentially to larval metal exposure. In the resistant strain, which previously was demonstrated to show increased detoxification enzyme activity and insecticide tolerance after the same treatment, the potential increase in transcriptional activity appeared to be modulated by reduced methylation and increased HAT activity. This study suggests that the commercial epigenetic quantification kits can be used to characterise phenotypic changes in An. arabiensis, and also shows that epigenetic regulation of the response to metal exposure is regulated at the DNA as opposed to the RNA level.

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.

Single nucleotide polymorphism (SNP) in the doublesex (dsx) gene splice sites and relevance for its alternative splicing in the malaria vector Anopheles gambiae

Background: The malaria burden continues to be significant in tropical regions, and conventional vector control methods are faced with challenges such as insecticide resistance. To overcome these challenges, additional vector control interventions are vital and include modern genetic approaches as well as classical methods like the sterile insect technique (SIT). In the major human malaria vector Anopheles gambiae, a candidate gene favourable for sterility induction is the doublesex (dsx) gene, encoding somatic sexually dimorphic traits in mosquitoes. However, the mechanism that regulates the expression of this gene in anopheline mosquitoes is poorly understood. This study aimed to screen the An. gambiae dsx gene splice site sequences for single nucleotide polymorphisms (SNPs) that could be critical to its alternative splicing. Methods: Variant annotation data from Ag1000G project phase 2 was analysed, in order to identify splice-relevant SNPs within acceptor and donor splice sites of the An. gambiae dsx gene (Agdsx). Results: SNPs were found in both donor and acceptor sites of the Agdsx. No splice-relevant SNPs were identified in the female-specific intron 4 acceptor site and the corresponding region in males. Two SNPs (rs48712947, rs48712962) were found in the female-specific donor site of exon 5. They were not specific to either males or females as the rs48712947 was found in female mosquitoes from Cameroon, and in both males and females from Burkina Faso. In the other splice sites, the intron 3 acceptor site carried the greatest abundance of SNPs.   Conclusions: There were no gender association between the identified SNPs and the random distribution of these SNPs in mosquito populations. The SNPs in Agdsx splice sites are not critical for the alternative splicing. Other molecular mechanisms should be considered and investigated.

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.

Factors Determining Transmission of Persistent Viruses by Bemisia tabaci and Emergence of New Virus-Vector Relationships

Many plant viruses depend on insect vectors for their transmission and dissemination. The whitefly Bemisia tabaci (Hemiptera: Aleyrodidae) is one of the most important virus vectors, transmitting more than four hundred virus species, the majority belonging to begomoviruses (Geminiviridae), with their ssDNA genomes. Begomoviruses are transmitted by B. tabaci in a persistent, circulative manner, during which the virus breaches barriers in the digestive, hemolymph, and salivary systems, and interacts with insect proteins along the transmission pathway. These interactions and the tissue tropism in the vector body determine the efficiency and specificity of the transmission. This review describes the mechanisms involved in circulative begomovirus transmission by B. tabaci, focusing on the most studied virus in this regard, namely the tomato yellow leaf curl virus (TYLCV) and its closely related isolates. Additionally, the review aims at drawing attention to the recent knowhow of unorthodox virus—B. tabaci interactions. The recent knowledge of whitefly-mediated transmission of two recombinant poleroviruses (Luteoviridae), a virus group with an ssRNA genome and known to be strictly transmitted with aphids, is discussed with its broader context in the emergence of new whitefly-driven virus diseases.

Effects of insemination and blood-feeding on locomotor activity of wild-derived females of the malaria mosquito Anopheles coluzzii

BACKGROUND

Behavioural shifts in the canonical location and timing of biting have been reported in natural populations of anopheline malaria vectors following the implementation of insecticide-based indoor vector control interventions. These modifications increase the likelihood of human-vector contact and allow mosquitoes to avoid insecticides, both conditions being favourable to residual transmission of the malarial parasites. The biting behaviour of mosquitoes follows rhythms that are under the control of biological clocks and environmental conditions, modulated by physiological states. In this work we explore modifications of spontaneous locomotor activity expressed by mosquitoes in different physiological states to highlight phenotypic variability associated to circadian control that may contribute to explain residual transmission in the field.

METHODS

The F10 generation progeny of field-collected Anopheles coluzzii from southwestern Burkina Faso was tested using an automated recording apparatus (Locomotor Activity Monitor, TriKinetics Inc.) under LD 12:12 or DD light regimens in laboratory-controlled conditions. Activity recordings of each test were carried out for a week with 6-day-old females belonging to four experimental treatments, representing factorial combinations of two physiological variables: insemination status (virgin vs inseminated) and gonotrophic status (glucose fed vs blood fed). Chronobiological features of rhythmicity in locomotor activity were explored using periodograms, diversity indices, and generalized linear mixed modelling.

RESULTS

The average strength of activity, onset of activity, and acrophase were modulated by both nutritional and insemination status as well as by the light regimen. Inseminated females showed a significant excess of arrhythmic activity under DD. When rhythmicity was observed in DD, females displayed sustained activity also during the subjective day.

CONCLUSIONS

Insemination and gonotrophic status influence the underlying light and circadian control of chronobiological features of locomotor activity. Overrepresentation of arrhythmic chronotypes as well as the sustained activity of inseminated females during the subjective day under DD conditions suggests potential activity of natural populations of A. coluzzii during daytime under dim conditions, with implications for residual transmission of malarial parasites.

From Global to Local-New Insights into Features of Pyrethroid Detoxification in Vector Mosquitoes

The threat of mosquito-borne diseases continues to be a problem for public health in subtropical and tropical regions of the world; in response, there has been increased use of adulticidal insecticides, such as pyrethroids, in human habitation areas over the last thirty years. As a result, the prevalence of pyrethroid-resistant genetic markers in natural mosquito populations has increased at an alarming rate. This review details recent advances in the understanding of specific mechanisms associated with pyrethroid resistance, with emphasis on features of insecticide detoxification and the interdependence of multiple cellular pathways. Together, these advances add important context to the understanding of the processes that are selected in resistant mosquitoes. Specifically, before pyrethroids bind to their targets on motoneurons, they must first permeate the outer cuticle and diffuse to inner tissues. Resistant mosquitoes have evolved detoxification mechanisms that rely on cytochrome P450s (CYP), esterases, carboxyesterases, and other oxidation/reduction (redox) components to effectively detoxify pyrethroids to nontoxic breakdown products that are then excreted. Enhanced resistance mechanisms have evolved to include alteration of gene copy number, transcriptional and post-transcriptional regulation of gene expression, as well as changes to cellular signaling mechanisms. Here, we outline the variety of ways in which detoxification has been selected in various mosquito populations, as well as key gene categories involved. Pathways associated with potential new genes of interest are proposed. Consideration of multiple cellular pathways could provide opportunities for development of new insecticides.

Filtering the Junk: Assigning Function to the Mosquito Non-Coding Genome

The portion of the mosquito genome that does not code for proteins contains regulatory elements that likely underlie variation for important phenotypes including resistance and susceptibility to infection with arboviruses and Apicomplexan parasites. Filtering the non-coding genome to uncover these functional elements is an expanding area of research, though identification of non-coding regulatory elements is challenging due to the lack of an amino acid-like code for the non-coding genome and a lack of sequence conservation across species. This review focuses on three types of non-coding regulatory elements: (1) microRNAs (miRNAs), (2) long non-coding RNAs (lncRNAs), and (3) enhancers, and summarizes current advances in technical and analytical approaches for measurement of each of these elements on a genome-wide scale. The review also summarizes and highlights novel findings following application of these techniques in mosquito-borne disease research. Looking beyond the protein-coding genome is essential for understanding the complexities that underlie differential gene expression in response to arboviral or parasite infection in mosquito disease vectors. A comprehensive understanding of the regulation of gene and protein expression will inform transgenic and other vector control methods rooted in naturally segregating genetic variation.

PeakMatcher facilitates updated Aedes aegypti embryonic cis-regulatory element map

BACKGROUND

The Aedes aegypti mosquito is a threat to human health across the globe. The A. aegypti genome was recently re-sequenced and re-assembled. Due to a combination of long-read PacBio and Hi-C sequencing, the AaegL5 assembly is chromosome complete and significantly improves the assembly in key areas such as the M/m sex-determining locus. Release of the updated genome assembly has precipitated the need to reprocess historical functional genomic data sets, including cis-regulatory element (CRE) maps that had previously been generated for A. aegypti.

RESULTS

We re-processed and re-analyzed the A. aegypti whole embryo FAIRE seq data to create an updated embryonic CRE map for the AaegL5 genome. We validated that the new CRE map recapitulates key features of the original AaegL3 CRE map. Further, we built on the improved assembly in the M/m locus to analyze overlaps of open chromatin regions with genes. To support the validation, we created a new method (PeakMatcher) for matching peaks from the same experimental data set across genome assemblies.

CONCLUSION

Use of PeakMatcher software, which is available publicly under an open-source license, facilitated the release of an updated and validated CRE map, which is available through the NIH GEO. These findings demonstrate that PeakMatcher software will be a useful resource for validation and transferring of previous annotations to updated genome assemblies.

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

Anopheles gambiae mosquitoes are primary human malaria vectors, but we know very little about their mechanisms of transcriptional regulation. We profiled chromatin accessibility by the assay for transposase-accessible chromatin by sequencing (ATAC-seq) in laboratory-reared A. gambiae mosquitoes experimentally infected with the human malaria parasite Plasmodium falciparum. By integrating ATAC-seq, RNA-seq and ChIP-seq data, we showed a positive correlation between accessibility at promoters and introns, gene expression and active histone marks. By comparing expression and chromatin structure patterns in different tissues, we were able to infer cis-regulatory elements controlling tissue-specific gene expression and to predict the in vivo binding sites of relevant transcription factors. The ATAC-seq assay also allowed the precise mapping of active regulatory regions, including novel transcription start sites and enhancers that were annotated to mosquito immune-related genes. Not only is this study important for advancing our understanding of mechanisms of transcriptional regulation in the mosquito vector of human malaria, but the information we produced also has great potential for developing new mosquito-control and anti-malaria strategies.

Capturing the transcription factor interactome in response to sub-lethal insecticide exposure

The increasing levels of pesticide resistance in agricultural pests and disease vectors represents a threat to both food security and global health. As insecticide resistance intensity strengthens and spreads, the likelihood of a pest encountering a sub-lethal dose of pesticide dramatically increases. Here, we apply dynamic Bayesian networks to a transcriptome time-course generated using sub-lethal pyrethroid exposure on a highly resistant Anopheles coluzzii population. The model accounts for circadian rhythm and ageing effects allowing high confidence identification of transcription factors with key roles in pesticide response. The associations generated by this model show high concordance with lab-based validation and identifies 44 transcription factors putatively regulating insecticide-responsive transcripts. We identify six key regulators, with each displaying differing enrichment terms, demonstrating the complexity of pesticide response. The considerable overlap of resistance mechanisms in agricultural pests and disease vectors strongly suggests that these findings are relevant in a wide variety of pest species.

Chromatin Structure and Function in Mosquitoes

The principles and function of chromatin and nuclear architecture have been extensively studied in model organisms, such as Drosophila melanogaster. However, little is known about the role of these epigenetic processes in transcriptional regulation in other insects including mosquitoes, which are major disease vectors and a worldwide threat for human health. Some of these life-threatening diseases are malaria, which is caused by protozoan parasites of the genus Plasmodium and transmitted by Anopheles mosquitoes; dengue fever, which is caused by an arbovirus mainly transmitted by Aedes aegypti; and West Nile fever, which is caused by an arbovirus transmitted by Culex spp. In this contribution, we review what is known about chromatin-associated mechanisms and the 3D genome structure in various mosquito vectors, including Anopheles, Aedes, and Culex spp. We also discuss the similarities between epigenetic mechanisms in mosquitoes and the model organism Drosophila melanogaster, and advocate that the field could benefit from the cross-application of state-of-the-art functional genomic technologies that are well-developed in the fruit fly. Uncovering the mosquito regulatory genome can lead to the discovery of unique regulatory networks associated with the parasitic life-style of these insects. It is also critical to understand the molecular interactions between the vectors and the pathogens that they transmit, which could hold the key to major breakthroughs on the fight against mosquito-borne diseases. Finally, it is clear that epigenetic mechanisms controlling mosquito environmental plasticity and evolvability are also of utmost importance, particularly in the current context of globalization and climate change.

Recent advances and future perspectives in vector-omics

We have reviewed recent progress and the remaining challenges in vector-omics. We have highlighted several technologies and applications that facilitate novel biological insights beyond achieving a reference-quality genome assembly. Among other topics, we have discussed the applications of chromatin conformation capture, chromatin accessibility assays, optical mapping, full-length RNA sequencing, single cell RNA analysis, proteomics, and population genomics. We anticipate that we will witness a great expansion in vector-omics research not only in its application in a broad range of species, but also its ability to uncover novel genetic elements and tackle previously inaccessible regions of the genome. It is our hope that the continued innovation in device portability, cost reduction, and informatics support will in the foreseeable future bring vector-omics to every vector laboratory and field station in the world, which will have an unparalleled impact on basic research and the control of vector-borne infectious diseases.

DNA Methylation in Anopheles albimanus Modulates the Midgut Immune Response Against Plasmodium berghei

Epigenetic mechanisms such as DNA methylation and histone post-translational modifications are fundamental for the phenotypic plasticity of insects during their interaction with the environment. In response to environmental cues, the methylation pattern in DNA is dynamically remodeled to achieve an epigenetic control of gene expression. DNA methylation is the focus of study in insects for its evolutionarily conserved character; however, there is scant knowledge about the epigenetic regulation in vector mosquitoes, especially during their infection by parasites. The aim of the present study was to evaluate the participation of DNA methylation in the immune response of Anopheles albimanus to a Plasmodium infection. For this, we first investigated the presence of a fully functional DNA methylation system in A. albimanus by assessing its potential role in larval development. Subsequently, we evaluated the transcriptional response to Plasmodium berghei of two mosquito phenotypes with different degrees of susceptibility to the parasite, in a scenario where their global DNA methylation had been pharmacologically inhibited. Our study revealed that A. albimanus has a functional DNA methylation system that is essential to larval viability, and that is also responsive to feeding and parasite challenges. The pharmacological erasure of the methylome with azacytidine or decitabine abolished the divergent responses of both mosquito phenotypes, leading to a transcriptionally similar response upon parasite challenge. This response was more specific, and the infection load in both phenotypes was lowered. Our findings suggest that DNA methylation may constitute a key factor in vector competence, and a promising target for preventing malaria transmission.

Influence of genetic polymorphism on transcriptional enhancer activity in the malaria vector Anopheles coluzzii

Enhancers are cis-regulatory elements that control most of the developmental and spatial gene expression in eukaryotes. Genetic variation of enhancer sequences is known to influence phenotypes, but the effect of enhancer variation upon enhancer functional activity and downstream phenotypes has barely been examined in any species. In the African malaria vector, Anopheles coluzzii, we identified candidate enhancers in the proximity of genes relevant for immunity, insecticide resistance, and development. The candidate enhancers were functionally validated using luciferase reporter assays, and their activity was found to be essentially independent of their physical orientation, a typical property of enhancers. All of the enhancers segregated genetically polymorphic alleles, which displayed significantly different levels of functional activity. Deletion mutagenesis and functional testing revealed a fine structure of positive and negative regulatory elements that modulate activity of the enhancer core. Enhancer polymorphisms segregate in wild A. coluzzii populations in West Africa. Thus, enhancer variants that modify target gene expression leading to likely phenotypic consequences are frequent in nature. These results demonstrate the existence of naturally polymorphic A. coluzzii enhancers, which may help explain important differences between individuals or populations for malaria transmission efficiency and vector adaptation to the environment.