Insights into the epigenomic landscape of the human malaria vector Anopheles gambiae

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.

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.

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.

Genome-Wide Mapping of Histone H3 Lysine 4 Trimethylation (H3K4me3) and Its Involvement in Fatty Acid Biosynthesis in Sunflower Developing Seeds

Histone modifications are of paramount importance during plant development. Investigating chromatin remodeling in developing oilseeds sheds light on the molecular mechanisms controlling fatty acid metabolism and facilitates the identification of new functional regions in oil crop genomes. The present study characterizes the epigenetic modifications H3K4me3 in relationship with the expression of fatty acid-related genes and transcription factors in developing sunflower seeds. Two master transcriptional regulators identified in this analysis, VIV1 (homologous to Arabidopsis ABI3) and FUS3, cooperate in the regulation of WRINKLED 1, a transcriptional factor regulating glycolysis, and fatty acid synthesis in developing oilseeds.

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.

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.

bsAS, an antisense long non-coding RNA, essential for correct wing development through regulation of blistered/DSRF isoform usage

Natural Antisense Transcripts (NATs) are long non-coding RNAs (lncRNAs) that overlap coding genes in the opposite strand. NATs roles have been related to gene regulation through different mechanisms, including post-transcriptional RNA processing. With the aim to identify NATs with potential regulatory function during fly development, we generated RNA-Seq data in Drosophila developing tissues and found bsAS, one of the most highly expressed lncRNAs in the fly wing. bsAS is antisense to bs/DSRF, a gene involved in wing development and neural processes. bsAS plays a crucial role in the tissue specific regulation of the expression of the bs/DSRF isoforms. This regulation is essential for the correct determination of cell fate during Drosophila development, as bsAS knockouts show highly aberrant phenotypes. Regulation of bs isoform usage by bsAS is mediated by specific physical interactions between the promoters of these two genes, which suggests a regulatory mechanism involving the collision of RNA polymerases transcribing in opposite directions. Evolutionary analysis suggests that bsAS NAT emerged simultaneously to the long-short isoform structure of bs, preceding the emergence of wings in insects.

Conserved Small Nucleotidic Elements at the Origin of Concerted piRNA Biogenesis from Genes and lncRNAs

PIWI-interacting RNAs (piRNAs) target transcripts by sequence complementarity serving as guides for RNA slicing in animal germ cells. The piRNA pathway is increasingly recognized as critical for essential cellular functions such as germline development and reproduction. In the Anopheles gambiae ovary, as much as 11% of piRNAs map to protein-coding genes. Here, we show that ovarian mRNAs and long non-coding RNAs (lncRNAs) are processed into piRNAs that can direct other transcripts into the piRNA biogenesis pathway. Targeting piRNAs fuel transcripts either into the ping-pong cycle of piRNA amplification or into the machinery of phased piRNA biogenesis, thereby creating networks of inter-regulating transcripts. RNAs of the same network share related genomic repeats. These repeats give rise to piRNAs, which target other transcripts and lead to a cascade of concerted RNA slicing. While ping-pong networks are based on repeats of several hundred nucleotides, networks that rely on phased piRNA biogenesis operate through short ~40-nucleotides long repeats, which we named snetDNAs. Interestingly, snetDNAs are recurring in evolution from insects to mammals. Our study brings to light a new type of conserved regulatory pathway, the snetDNA-pathway, by which short sequences can include independent genes and lncRNAs in the same biological pathway.

Chromatin changes in Anopheles gambiae induced by Plasmodium falciparum infection

BACKGROUND

Infection by the human malaria parasite leads to important changes in mosquito phenotypic traits related to vector competence. However, we still lack a clear understanding of the underlying mechanisms and, in particular, of the epigenetic basis for these changes. We have examined genome-wide distribution maps of H3K27ac, H3K9ac, H3K9me3 and H3K4me3 by ChIP-seq and the transcriptome by RNA-seq, of midguts from Anopheles gambiae mosquitoes blood-fed uninfected and infected with natural isolates of the human malaria parasite Plasmodium falciparum in Burkina Faso.

RESULTS

We report 15,916 regions containing differential histone modification enrichment between infected and uninfected, of which 8339 locate at promoters and/or intersect with genes. The functional annotation of these regions allowed us to identify infection-responsive genes showing differential enrichment in various histone modifications, such as CLIP proteases, antimicrobial peptides-encoding genes, and genes related to melanization responses and the complement system. Further, the motif analysis of regions differentially enriched in various histone modifications predicts binding sites that might be involved in the cis-regulation of these regions, such as Deaf1, Pangolin and Dorsal transcription factors (TFs). Some of these TFs are known to regulate immunity gene expression in Drosophila and are involved in the Notch and JAK/STAT signaling pathways.

CONCLUSIONS

The analysis of malaria infection-induced chromatin changes in mosquitoes is important not only to identify regulatory elements and genes underlying mosquito responses to P. falciparum infection, but also for possible applications to the genetic manipulation of mosquitoes and to other mosquito-borne systems.

Epigenetic regulation of Plasmodium falciparum clonally variant gene expression during development in Anopheles gambiae

P. falciparum phenotypic plasticity is linked to the variant expression of clonal multigene families such as the var genes. We have examined changes in transcription and histone modifications that occur during sporogonic development of P. falciparum in the mosquito host. All var genes are silenced or transcribed at low levels in blood stages (gametocyte/ring) of the parasite in the human host. After infection of mosquitoes, a single var gene is selected for expression in the oocyst, and transcription of this gene increases dramatically in the sporozoite. The same PF3D7_1255200 var gene was activated in 4 different experimental infections. Transcription of this var gene during parasite development in the mosquito correlates with the presence of low levels of H3K9me3 at the binding site for the PF3D7_1466400 AP2 transcription factor. This chromatin state in the sporozoite also correlates with the expression of an antisense long non-coding RNA (lncRNA) that has previously been shown to promote var gene transcription during the intraerythrocytic cycle in vitro. Expression of both the sense protein-coding transcript and the antisense lncRNA increase dramatically in sporozoites. The findings suggest a complex process for the activation of a single particular var gene that involves AP2 transcription factors and lncRNAs.

Unusual chromatin state in Rhynchosciara americana (Diptera: Sciaridae)

Euchromatin and heterochromatin are usually defined by the degree of DNA compaction, gene content and combinations of histone and non-histone proteins. More recent studies on protein location have been able to specify a variety of chromatin types thus adding chromatin configurations other than the two basic reference states. Chromatin research exploiting non-model organisms has the potential to provide novel information related to epigenetic modifications and their impact on chromosome structure and function. Polytene chromosomes of Rhynchosciara americana display a particular region within the A9 sub-section characterised by lack of DNA compaction as well as an usual polytene banding pattern. DNA content in the sub-section seems to be low as deduced by DAPI staining. Antibodies to H3K4me, a conserved epigenetic transcription marker,labelled the A9 sub-section strongly. In contrast,transcriptional activity in the region, if any, seems to be low as inferred by detection of RNA polymerase II and RNA. Histone markers related to heterochromatin formation such as H3K9me and H3K27me are underrepresented in the A9 sub-section. However, a chromodomain-containing sciarid protein was detected in the region, displaying levels of fluorescence very close to those observed in pericentric heterochromatin.A plasmid micro-library constructed with microdissected DNA from the A9 sub-section was screened for repetitive DNA. The proportion of inserts containing repeats was found to be similar to that contained in another micro-library made with DNA from a single chromosome end of this species. The data suggest an unusual "chromatin colour" indicating that high levels of histone markers related to transcription coexist with a significant presence of chromodomain-containing proteins and the virtual absence of histone modifications observed in heterochromatin formation.

"Salivary gland cellular architecture in the Asian malaria vector mosquito Anopheles stephensi"

Background

Anopheles mosquitoes are vectors for malaria, a disease with continued grave outcomes for human health. Transmission of malaria from mosquitoes to humans occurs by parasite passage through the salivary glands (SGs). Previous studies of mosquito SG architecture have been limited in scope and detail.

Methods

We developed a simple, optimized protocol for fluorescence staining using dyes and/or antibodies to interrogate cellular architecture in Anopheles stephensi adult SGs. We used common biological dyes, antibodies to well-conserved structural and organellar markers, and antibodies against Anopheles salivary proteins to visualize many individual SGs at high resolution by confocal microscopy.

Results

These analyses confirmed morphological features previously described using electron microscopy and uncovered a high degree of individual variation in SG structure. Our studies provide evidence for two alternative models for the origin of the salivary duct, the structure facilitating parasite transport out of SGs. We compare SG cellular architecture in An. stephensi and Drosophila melanogaster, a fellow Dipteran whose adult SGs are nearly completely unstudied, and find many conserved features despite divergence in overall form and function. Anopheles salivary proteins previously observed at the basement membrane were localized either in SG cells, secretory cavities, or the SG lumen. Our studies also revealed a population of cells with characteristics consistent with regenerative cells, similar to muscle satellite cells or midgut regenerative cells.

Conclusions

This work serves as a foundation for linking Anopheles stephensi SG cellular architecture to function and as a basis for generating and evaluating tools aimed at preventing malaria transmission at the level of mosquito SGs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1229-z) contains supplementary material, which is available to authorized users.

Unraveling dual feeding associated molecular complexity of salivary glands in the mosquito Anopheles culicifacies

Mosquito salivary glands are well known to facilitate meal acquisition, however the fundamental question on how adult female salivary gland manages molecular responses during sugar versus blood meal uptake remains unanswered. To investigate these responses, we analyzed a total of 58.5 million raw reads generated from two independent RNAseq libraries of the salivary glands collected from 3–4 day-old sugar and blood fed Anopheles culicifacies mosquitoes. Comprehensive functional annotation analysis of 10,931 contigs unraveled that salivary glands may encode diverse nature of proteins in response to distinct physiological feeding status. Digital gene expression analysis and PCR validation indicated that first blood meal significantly alters the molecular architecture of the salivary glands. Comparative microscopic analysis also revealed that first blood meal uptake not only causes an alteration of at least 12–22% of morphological features of the salivary glands but also results in cellular changes e.g. apoptosis, confirming together that adult female salivary glands are specialized organs to manage meal specific responses. Unraveling the underlying mechanism of mosquito salivary gene expression, controlling dual feeding associated responses may provide a new opportunity to control vector borne diseases.