Mosquito-Borne Diseases and Omics: Tissue-Restricted Expression and Alternative Splicing Revealed by Transcriptome Profiling of Anopheles stephensi

Circadian rhythms mediate malaria transmission potential

Malaria transmission begins when infected female Anopheles mosquitos deposit Plasmodium parasites into the mammalian host's skin during a bloodmeal. The salivary gland-resident sporozoite parasites migrate to the bloodstream, subsequently invading and replicating within hepatocytes. As Anopheles mosquitos are more active at night, with a 24-hour rhythm, we investigated whether their salivary glands are under circadian control, anticipating bloodmeals and modulating sporozoite biology for host encounters. Here we show that approximately half of the mosquito salivary gland transcriptome, particularly genes essential for efficient bloodmeals such as anti-blood clotting factors, exhibits circadian rhythmic expression. Furthermore, we demonstrate that mosquitoes prefer to feed during nighttime, with the amount of blood ingested varying cyclically throughout the day. Notably, we show a substantial subset of the sporozoite transcriptome cycling throughout the day. These include genes involved in parasite motility, potentially modulating the ability to initiate infection at different times of day. Thus, although sporozoites are typically considered quiescent, our results demonstrate their transcriptional activity, revealing robust daily rhythms of gene expression. Our findings suggest a circadian evolutionary relationship between the vector, parasite and mammalian host that together modulate malaria transmission.

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.

The impact of transgenesis on mosquito fitness: A review

Transgenic mosquitoes developed by genetic manipulation, offer a promising strategy for the sustainable and effective control of mosquito-borne diseases. This strategy relies on the mass release of transgenic mosquitoes into the wild, where their transgene is expected to persist in the natural environment, either permanently or transiently, within the mosquito population. In such circumstances, the fitness of transgenic mosquitoes is an important factor in determining their survival in the wild. The impact of transgene expression, insertional mutagenesis, inbreeding depression related to laboratory adaptation, and the hitchhiking effect involved in developing homozygous mosquito lines can all have an effect on the fitness of transgenic mosquitoes. Therefore, real-time estimation of transgene-associated fitness cost is imperative for modeling and planning transgenic mosquito release programs. This can be achieved by directly comparing fitness parameters in individuals homozygous or hemizygous for the transgene and their wild-type counterparts, or by cage invasion experiments to monitor the frequency of the transgenic allele over multiple generations. Recent advancements such as site-specific integration systems and gene drives, provide platforms to address fitness issues in transgenic mosquitoes. More research on the fitness of transgenic individuals is required to develop transgenic mosquitoes with a low fitness cost.

A transcriptomic atlas of Aedes aegypti reveals detailed functional organization of major body parts and gut regional specializations in sugar-fed and blood-fed adult females

Mosquitoes transmit numerous pathogens, but large gaps remain in our understanding of their physiology. To facilitate explorations of mosquito biology, we have created Aegypti-Atlas (http://aegyptiatlas.buchonlab.com/), an online resource hosting RNAseq profiles of Ae. aegypti body parts (head, thorax, abdomen, gut, Malpighian tubules, ovaries), gut regions (crop, proventriculus, anterior and posterior midgut, hindgut), and a gut time course of blood meal digestion. Using Aegypti-Atlas, we provide insights into regionalization of gut function, blood feeding response, and immune defenses. We find that the anterior and posterior midgut possess digestive specializations which are preserved in the blood-fed state. Blood feeding initiates the sequential induction and repression/depletion of multiple cohorts of peptidases. With respect to defense, immune signaling components, but not recognition or effector molecules, show enrichment in ovaries. Basal expression of antimicrobial peptides is dominated by holotricin and gambicin, which are expressed in carcass and digestive tissues, respectively, in a mutually exclusive manner. In the midgut, gambicin and other effectors are almost exclusively expressed in the anterior regions, while the posterior midgut exhibits hallmarks of immune tolerance. Finally, in a cross-species comparison between Ae. aegypti and Anopheles gambiae midguts, we observe that regional digestive and immune specializations are conserved, indicating that our dataset may be broadly relevant to multiple mosquito species. We demonstrate that the expression of orthologous genes is highly correlated, with the exception of a ‘species signature’ comprising a few highly/disparately expressed genes. With this work, we show the potential of Aegypti-Atlas to unlock a more complete understanding of mosquito biology.

FINDER: an automated software package to annotate eukaryotic genes from RNA-Seq data and associated protein sequences

BACKGROUND

Gene annotation in eukaryotes is a non-trivial task that requires meticulous analysis of accumulated transcript data. Challenges include transcriptionally active regions of the genome that contain overlapping genes, genes that produce numerous transcripts, transposable elements and numerous diverse sequence repeats. Currently available gene annotation software applications depend on pre-constructed full-length gene sequence assemblies which are not guaranteed to be error-free. The origins of these sequences are often uncertain, making it difficult to identify and rectify errors in them. This hinders the creation of an accurate and holistic representation of the transcriptomic landscape across multiple tissue types and experimental conditions. Therefore, to gauge the extent of diversity in gene structures, a comprehensive analysis of genome-wide expression data is imperative.

RESULTS

We present FINDER, a fully automated computational tool that optimizes the entire process of annotating genes and transcript structures. Unlike current state-of-the-art pipelines, FINDER automates the RNA-Seq pre-processing step by working directly with raw sequence reads and optimizes gene prediction from BRAKER2 by supplementing these reads with associated proteins. The FINDER pipeline (1) reports transcripts and recognizes genes that are expressed under specific conditions, (2) generates all possible alternatively spliced transcripts from expressed RNA-Seq data, (3) analyzes read coverage patterns to modify existing transcript models and create new ones, and (4) scores genes as high- or low-confidence based on the available evidence across multiple datasets. We demonstrate the ability of FINDER to automatically annotate a diverse pool of genomes from eight species.

CONCLUSIONS

FINDER takes a completely automated approach to annotate genes directly from raw expression data. It is capable of processing eukaryotic genomes of all sizes and requires no manual supervision-ideal for bench researchers with limited experience in handling computational tools.

Genetic engineering and bacterial pathogenesis against the vectorial capacity of mosquitoes

Mosquitoes are the main vector of multiple diseases worldwide and transmit viral (malaria, chikungunya, encephalitis, yellow fever, as well as dengue fever), as well as bacterial diseases (tularemia). To manage the outbreak of mosquito populations, various management programs include the application of chemicals, followed by biological and genetic control. Here we aimed to focus on the role of bacterial pathogenesis and molecular tactics for the management of mosquitoes and their vectorial capacity. Bacterial pathogenesis and molecular manipulations have a substantial impact on the biology of mosquitoes, and both strategies change the gene expression and regulation of disease vectors. The strategy for genetic modification is also proved to be excellent for the management of mosquitoes, which halt the development of population via incompatibility of different sex. Therefore, the purpose of the present discussion is to illustrate the impact of both approaches against the vectorial capacity of mosquitoes. Moreover, it could be helpful to understand the relationship of insect-pathogen and to manage various insect vectors as well as diseases.

GFF Utilities: GffRead and GffCompare

GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations. Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

GFF Utilities: GffRead and GffCompare

Summary: GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations.

Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

GFF Utilities: GffRead and GffCompare

GTF (Gene Transfer Format) and GFF (General Feature Format) are popular file formats used by bioinformatics programs to represent and exchange information about various genomic features, such as gene and transcript locations and structure. GffRead and GffCompare are open source programs that provide extensive and efficient solutions to manipulate files in a GTF or GFF format. While GffRead can convert, sort, filter, transform, or cluster genomic features, GffCompare can be used to compare and merge different gene annotations. Availability and implementation: GFF utilities are implemented in C++ for Linux and OS X and released as open source under an MIT license  ( https://github.com/gpertea/gffread, https://github.com/gpertea/gffcompare).

Proteome data of Anopheles stephensi salivary glands using high-resolution mass spectrometry analysis

The data article reports data of the proteins expressed in female Anopheles stephensi salivary glands. Proteomic data were acquired using high-resolution mass spectrometers - Orbitrap-Velos and Orbitrap-Elite. Samples derived from adult female A. stephensi salivary glands led to the identification of 4390 proteins. Mass spectrometry data were analyzed on Proteome Discoverer (Version 2.1) platform with Sequest and Mascot search engines. The identified proteins were analyzed for their Gene Ontology annotation, interaction network and their possible roles in vector-parasite interaction. The data provided here are related to our published article "Integrating transcriptomics and proteomics data for accurate assembly and annotation of genomes" (Prasad et al., 2017) [1].

Dataset on fat body proteome of Anopheles stephensi Liston

Fat body from Anopheles stephensi female mosquitoes were dissected and processed for proteomic analysis. Both SDS-PAGE and basic Reverse Phase Liquid Chromatography-based fractionation strategies were used to achieve a broad coverage of protein identification. The fractionated peptides were then analyzed on a high-resolution mass spectrometer. Searching the raw data against the protein database of An. stephensi resulted in identification of 4535 proteins, which is, to our knowledge, the largest catalog of fat body proteome in any mosquito vector species reported so far. Bioinformatics analysis on these fat body proteins suggested the enrichment of biological processes including carbon and lipid metabolism, amino acid metabolism, signal peptide processing and oxidation-reduction. In addition, using proteogenomic approaches, 43 novel proteins were identified, which were not listed in the annotated gene annotations of An. stephensi. The data used in the analysis are related to the article ‘Integrating transcriptomic and proteomic data for accurate assembly and annotation of genomes’ (Prasad et al., 2017).

Proteome data of Anopheles stephensi ovary using high-resolution mass spectrometry

This article contains data on the proteins expressed in the ovaries of Anopheles stephensi, a major vector of malaria in India. Data acquisition was performed using a high-resolution Orbitrap-Velos mass spectrometer. The acquired MS/MS data was searched against An. stephensi protein database comprising of 11,789 sequences. Overall, 4407 proteins were identified, functional analysis was performed for the identified proteins and a protein-protein interaction map predicted. The data provided here is also related to a published article - “Integrating transcriptomics and proteomics data for accurate assembly and annotation of genomes” (Prasad et al., 2017) [1].

Integrated Multi-Omic Analysis of Mycobacterium tuberculosis H37Ra Redefines Virulence Attributes

H37Ra is a virulence attenuated strain of Mycobacterium tuberculosis widely employed as a model to investigate virulence mechanisms. Comparative high-throughput studies have earlier correlated its avirulence to the presence of specific mutations or absence of certain proteins. However, a recent sequencing study of H37Ra, has disproved several genomic differences earlier reported to be associated with virulence. This warrants further investigations on the H37Ra proteome as well. In this study, we carried out an integrated analysis of the genome, transcriptome, and proteome of H37Ra. In addition to confirming single nucleotide variations (SNVs) and insertion-deletions that were reported earlier, our study provides novel insights into the mutation spectrum in the promoter regions of 7 genes. We also provide transcriptional and proteomic evidence for 3,900 genes representing ~80% of the total predicted gene count including 408 proteins that have not been identified previously. We identified 9 genes whose coding potential was hitherto reported to be absent in H37Ra. These include 2 putative virulence factors belonging to ESAT-6 like family of proteins. Furthermore, proteogenomic analysis enabled us to identify 63 novel proteins coding genes and correct 25 existing gene models in H37Ra genome. A majority of these were found to be conserved in the virulent strain H37Rv as well as in other mycobacterial species suggesting that the differences in the virulent and avirulent strains of M. tuberculosis are not entirely dependent on the expression of certain proteins or their absence but may possibly be ascertained to functional changes.

Response to Blood Meal in the Fat Body of Anopheles stephensi Using Quantitative Proteomics: Toward New Vector Control Strategies Against Malaria

Malaria remains a grand challenge for disruptive innovation in global health therapeutics and diagnostics. Anopheles stephensi is one of the major vectors of malaria in Asia. Vector and transmission control are key focus areas in the fight against malaria, a field of postgenomics research where proteomics can play a substantive role. Moreover, to identify novel strategies to control the vector population, it is necessary to understand the vector life processes at a global and molecular scale. In this context, fat body is a vital organ required for vitellogenesis, vector immunity, vector physiology, and vector-parasite interaction. Given its central role in energy metabolism, vitellogenesis, and immune function, the proteome profile of the fat body and the impact of blood meal (BM) ingestion on the protein abundances of this vital organ have not been investigated so far. Therefore, using a proteomics approach, we identified the proteins expressed in the fat body of An. stephensi and their differential expression in response to BM ingestion. In all, we identified 3,218 proteins in the fat body using high-resolution mass spectrometry, of which 483 were found to be differentially expressed in response to the BM ingestion. Bioinformatics analysis of these proteins underscored their role in amino acid metabolism, vitellogenesis, lipid transport, signal peptide processing, mosquito immunity, and oxidation-reduction processes. Interestingly, we identified five novel genes, which were found to be differentially expressed upon BM ingestion. Proteins that exhibited altered expression in the present study are potential targets for vector control strategies and development of transmission blocking vaccines in the fight against malaria.