Skip to main content

Advertisement

The transcriptome of the mosquito Aedes fluviatilis (Diptera: Culicidae), and transcriptional changes associated with its native Wolbachia infection

BMC Genomics volume 18, Article number: 6 (2017) | Download Citation

Article metrics

Abstract

Background

Wolbachia is a bacterial endosymbiont that naturally infects a wide range of insect species, and causes drastic changes to host biology. Stable infections of Wolbachia in mosquitoes can inhibit infection with medically important pathogens such as dengue virus and malaria-causing Plasmodium parasites. However, some native Wolbachia strains can enhance infection with certain pathogens, as is the case for the mosquito Aedes fluviatilis, where infection with Plasmodium gallinaceum is enhanced by the native wFlu Wolbachia strain. To better understand the biological interactions between mosquitoes and native Wolbachia infections, and to investigate the process of pathogen enhancement, we used RNA-Seq to generate the transcriptome of Ae. fluviatilis with and without Wolbachia infection.

Results

In total, we generated 22,280,160 Illumina paired-end reads from Wolbachia-infected and uninfected mosquitoes, and used these to make a de novo transcriptome assembly, resulting in 58,013 contigs with a median sequence length of 443 bp and an N50 of 2454 bp. Contigs were annotated through local alignments using BlastX, and associated with both gene ontology and KEGG orthology terms. Through baySeq, we identified 159 contigs that were significantly upregulated due to Wolbachia infection, and 98 that were downregulated. Critically, we saw no changes to Toll or IMD immune gene transcription, but did see evidence that wFlu infection altered the expression of several bacterial recognition genes, and immune-related genes that could influence Plasmodium infection. wFlu infection also had a widespread effect on a number of host physiological processes including protein, carbohydrate and lipid metabolism, and oxidative stress. We then compared our data set with transcriptomic data for other Wolbachia infections in Aedes aegypti, and identified a core set of 15 gene groups associated with Wolbachia infection in mosquitoes.

Conclusions

While the scale of transcriptional changes associated with wFlu infection might be small, the scope is rather large, which confirms that native Wolbachia infections maintain intricate molecular relationships with their mosquito hosts even after lengthy periods of co-evolution. We have also identified several potential means through which wFlu infection might influence Plasmodium infection in Ae. fluviatilis, and these genes should form the basis of future investigation into the enhancement of Plasmodium by Wolbachia.

Background

Wolbachia pipientis is a maternally inherited, bacterial endosymbiont that naturally infects an estimated 40% of terrestrial insect species [1]. Wolbachia induce a wide range of physiological manipulations in different insect hosts, with manipulation of reproductive biology promoting maternal transmission and thus bacterial propagation [2, 3]. It is through this ability to alter host biology that Wolbachia have gained interest as a form of biological control for the mosquito-transmitted pathogens that are responsible for diseases such as malaria, dengue fever, chikungunya and Zika fever, which represent a serious threat to human health across the globe [4, 5].

Many Wolbachia strains induce the reproductive manipulation cytoplasmic incompatibility (CI) in their hosts. This occurs when Wolbachia-infected male insects mate with uninfected females, which then produce unviable eggs. In contrast, Wolbachia-infected females successfully produce viable progeny after mating with either infected or uninfected males [3]. CI increases the proportion of Wolbachia-infected insects over subsequent generations, and serves to replace Wolbachia-uninfected individuals in population a population with those infected by the bacterium [6, 7]. CI-causing strains can be used to suppress mosquito populations that are uninfected by Wolbachia through the release of infected males, similar to the sterile insect technique, or to control Wolbachia-infected populations by releasing mosquitoes infected with a different strain, as this also crashes the population [8].

Several Wolbachia strains also produce anti-pathogenic effects in their hosts through the pathogen interference phenotype. The mechanics, scope and effectiveness of pathogen interference vary significantly between Wolbachia strains and insect hosts [9–12]. More effective pathogen interference severely inhibits pathogen development and transmission within the host. In Aedes aegypti, a prominent mosquito vector of human diseases, artificial Wolbachia infection interferes with the dengue, Zika, chikungunya, yellow fever and West Nile viruses, and other pathogens including the avian malaria Plasmodium gallinaceum, pathogenic bacteria, and the filarial nematode Brugia pahangi [13–19]. Recent work has also demonstrated that artificial Wolbachia infection in the malaria vector Anopheles stephensi can interfere with infection by the human malaria parasite Plasmodium falciparum, indicating that pathogen interference has broad applicability against human pathogens transmitted by mosquitoes [20].

CI and pathogen interference are the basis for the population replacement form of mosquito control utilised by the Eliminate Dengue Project (http://www.eliminatedengue.com/program). This strategy involves the release of Wolbachia-infected mosquitoes; CI allows the bacterium to spread and become stable within the target, wild population, while pathogen interference makes these mosquitoes less likely to transmit important viruses [17, 21]. Wolbachia has been successfully spread into a wild Ae. aegypti population [6], with the infection and strong pathogen interference against dengue virus persisting after several years of co-evolution [22, 23]

Neither Ae. aegypti nor An. stephensi are known to be naturally infected by Wolbachia. The infections of these mosquitoes described above were generated through transinfection, where Wolbachia is taken from a donor species and then injected into the eggs of the target species to create a stable germline infection transmitted to offspring [20, 24, 25]. In comparison to natural Wolbachia infections, such transinfections typically have a higher bacterial density, and infect a wider range of host tissues, which makes them far more likely to produce pathogen interference, and other extreme manipulations of host physiology [10, 20, 24]. Pathogen interference could potentially be lost from transinfected mosquitoes due to co-evolution between Wolbachia and host, or adaptation on the part of the pathogen [26]. Native Wolbachia infections typically produce minimal pathogen interference, and have little apparent utility to mosquito control strategies that require that trait. However, their low bacterial density, and presumed lower levels of virulence may be reflective the future biological state of transinfected mosquitoes after a long period of adaptation between host and symbiont.

Other native Wolbachia associations can enhance pathogen infection, as is the case for wPip in Culex pipiens when challenged by Plasmodium relictum [27]. Enhancement is commonly associated with artificial transient somatic Wolbachia infections, and has not been seen with stable germline transinfections [2]. Its mechanism is unknown, but may involve changes to host immunity, metabolism or transcription [27–30]. Needless to say that both loss of pathogen interference, and the development of enhancement would be undesirable consequences if they were to occur in Wolbachia used for vector control. To that end, understanding how native Wolbachia strains influence host physiology at the molecular level will provide useful information about how these strains influence response to pathogens, and potentially highlight a mechanism for enhancement.

The mosquito Aedes fluviatilis inhabits non-urbanized regions throughout Latin America. It is not regarded as a vector for human pathogens in the field, however it is a good laboratory model for P. gallinaceum infection [31]. It is naturally infected by the Wolbachia strain wFlu, which grows only to low density, causes CI, and does not induce noticeable fitness costs [15, 32]. The effect of wFlu on dengue virus has not been investigated, however wFlu was shown to enhance P. gallinaceum oocyst numbers during some experimental infections [32], making it an interesting model to understand both native Wolbachia infections and pathogen enhancement. To determine whether there was a transcriptional basis for this enhancement and to further the understanding of native Wolbachia strains, and the extent to which they impact host biology, we used RNA-Seq to generate the transcriptome of Ae. fluviatilis mosquitoes both with and without their native Wolbachia infection.

Results and Discussion

RNA sequencing and de novo transcriptome assembly

We generated a total of 22,280,160 Illumina paired-end reads across 6 Ae. fluviatilis libraries - 3 with Wolbachia infection (wFlu), and 3 where the native infection had been cleared by treatment with tetracycline (Tet). Each library was sequenced from a pool of 16 whole adult female mosquitoes, collected 6 days after eclosion. After the trimming of adaptors and filtering for low quality reads we were left with 19,919,299 high quality paired-end reads (Q30% = 91), across all six libraries. As there was no published genome for Ae. fluviatilis, we attempted to map against the Ae. aegypti predicted transcriptome with Bowtie2, but less than 25% could be successfully mapped, which was unsuitable for further analysis. We then used the complete set of reads to make a de novo transcriptome assembly with Trinity (see methods). A total of 58,013 contigs were assembled encompassing 64 million base pairs (bp), with a median sequence length of 443 bp and an N50 of 2454 bp (Table 1). Over 19,000 contigs were larger than 1Kb in size.

Table 1 Summary of sequencing results
Full size table

Functional annotation

Annotation was performed through local alignments using BlastX (with an e-value of 1e-10) against the NCBI non-redundant database (NR), which returned annotation for 26,066 contigs. Of these, 14,209 contained unique hits (Table 1). As previous molecular characterization of Ae. fluviatilis has been limited to a few microsatellites [33], our analysis relied entirely on finding similarity to other genomes. Consequently, 31,947 contigs (55% of the entire data set) did not get a significant BlastX hit, which could be explained by misassembly or from lack of representation in the NR database. 920 of these were considered to be protein-coding sequences based on predictions made by Transdecoder, but no further characterization was possible for the remaining 31,027 contigs. For the purposes of our analysis we have assumed that gene function is identical to what was described in the blast hits, but this may not be the case, as genetic homology does not always imply functional homology. However, we are fairly confident given that majority of annotations were previously described in closely related mosquitoes.

The assembled transcriptome of Ae. fluviatilis showed a high degree of similarity to Ae. aegypti, with 18,082 (69.38%) of the annotated contigs most closely matched to that species (Fig. 1). This was not surprising as Ae. aegypti has the most thoroughly annotated genome of the mosquitoes in the Aedes genus, meaning that this observed similarity may have more to do with the prevalence of Ae. aegypti sequences in public databases than the closeness of the relationship between the two species [34]. A further 5,616 (21.55%) of the annotated contigs were homologous to genes from other mosquito species, and overall 94.9% of the matches from BlastX were arthropod in origin. We conducted divergence analyses to better clarify the relationship between Ae. fluviatilis and other mosquito species (Additional file 1), and determined that Ae. fluviatilis diverged from these species approximately 98 million years ago (95% highest posterior density interval: 64.1 to 133.5 million years ago). A further 383 contigs matched to non-arthropod animals.

Fig. 1
figure1

Breakdown of contig annotations by organism of origin. After de novo assembly, contigs were annotated by local alignment with the BlastX database. The chart depicts the percentage of contigs where the most significant BlastX hit matched to a particular species of clade. 90.93% of the hits were to a mosquito species. Contigs without a BlastX hit (31,951) are not shown

Full size image

Across both the Tet and wFlu libraries we identified 751 contigs of bacterial origin (2.86%), and 112 of fungal origin (0.43%), which may potentially represent part of the Ae. fluviatilis intestinal microbiota, or could be the result of environmental contamination during sample collection. The bacterial sequences represented 250 distinct taxa, with the majority associated with a single contig. The diversity of sequence origins may indicate that the microbiota comprises many species, or that the majority of gut flora belonged to species that were not represented in the NR database. It is also possible that some of these contigs actually came from Wolbachia, but were different in composition from previously sequenced Wolbachia genes. Given the relatively low number of contigs of bacterial origin, the low copy numbers of these contigs and the low number of hits for each bacterial taxon in our data set, it was not possible to gain significant insight into the potential influence of wFlu on the composition or role of the microbiota.

Only 42 contigs of Wolbachia origin were identified in our data (Additional file 2). These contigs came from a variety of strains including wAna, wBm, wMelPop, and wPip. 38 were expressed in the wFlu libraries, and 4 in the Tet libraries - potentially via integration into the host genome. Wolbachia genomes typically contain between 900–1400 coding sequences [35]. Of the 42 identified Wolbachia contigs, 22 were hypothetical proteins with unknown function. Interestingly, one contig was a heme biosynthesis protein, which suggests that wFlu promotes heme production, and provides further evidence to support the theory that Wolbachia alters host iron and redox homeostasis, even amongst native associations [36, 37]. We also identified a glycosidase hydrolase, which plays a role in sugar metabolism. The remainder of the Wolbachia contigs fell into two broad functional categories. The first included genes involved in DNA/RNA processing, DNA repair and RNA synthesis, and are likely part of normal Wolbachia replication and transcription processes, as similar genes have been identified in the wMel genome [38]. A further group of contigs included Wolbachia membrane proteins, and ankyrin genes, which are used to attach the bacterial membrane to the host cytoskeleton. These genes likely represent part of the machinery used to mediate Wolbachia-host interactions [39], are typically present in large numbers in Wolbachia genomes [35, 38, 40–42], and may also be involved in CI [43].

The genomes of both Wolbachia strains and their insect hosts typically contain large quantities of mobile elements including transposons [38]. These can alter or disrupt the expression of genes, depending on their point of insertion into the host genome [34, 44]. We found 61 contigs related to transposable elements, none of which were matched to a Wolbachia genome. We also found a further 28 of viral origin, which could represent sequences from past or present members of the Ae. fluviatilis viral flora, which have not been well characterised.

As part of the contig annotation process, gene ontology (GO) and KEGG orthology (KO) terms for each contig were identified using Blast2GO and the KEGG database, respectively (Fig. 2). We identified genes associated with a wide range of biological processes, functions and structures. From the whole set of contigs, 23,035 were assigned to at least one GO term. The majority of contigs in the biological process category were associated with cellular or metabolic process. Cell, and cell part were the cellular component terms with the greatest frequency, while binding and catalytic activity were the most common terms associated with the molecular function category. The most common KO terms were related to diseases, and molecular information processing.

Fig. 2
figure2

GO and KO terms associated with the Aedes fluviatilis transcriptome. a The 40 most enriched gene ontology terms in the biological process (blue), cell component (red) and molecular function (yellow) categories at level 2. b First level (orange) and second level (purple) KEGG orthology functional category terms associated with Ae. fluviatilis. GO and KO term lists were generated using combined data for both the Wolbachia-infected and -uninfected libraries

Full size image

Differentially expressed contigs

Each RNA-Seq library was independently mapped to the assembled transcriptome using Bowtie. An average of 95% of the reads from each library were successfully mapped, and these data were used to generate counts for each contig. A list of differentially expressed contigs was generated using the baySeq package from Bioconductor, with 66% of these contigs (data not shown) also determined to be differentially expressed via analysis with the DESeq2 package, again from Bioconductor. Given this high level of concordance between the lists, we chose to proceed with further analysis of the baySeq list, as that method of analysis is known to be more sensitive [45]. Through baySeq we calculated the FPKM (Fragments Per Kilobase of transcript per Million mapped reads) for each contig, which considers the size of the contig in base pairs, and the overall data set size, as a measure of expression.

A total of 257 differentially expressed contigs were identified using baySeq, 159 were associated with the wFlu libraries, and 98 with the Tet libraries. Of these, 50 (19.4%) from the wFlu libraries (upregulated by wFlu) and 32 (32.7%) from the Tet libraries (downregulated by wFlu) had no matches after the BlastX search. This left 109 and 66 possible orthologs to known genes with increased and decreased expression due to Wolbachia infection, respectively (Additional file 3). Excluding multiple matches to genes with the same accession number, there were 95 unique contigs that were upregulated by Wolbachia, and 59 unique contigs that were downregulated, on the whole mosquito level. It is possible that analysis of whole mosquitoes precluded detection of genes that experienced a pattern of tissue-specific differential expression.

To confirm the accuracy of the sequencing data, the expression of six contigs was examined using RT-qPCR (Fig. 3). Two contigs, comp10645_c1_seq2 and comp10453_c0_seq5 were predicted to be upregulated by wFlu, and further analysis confirmed that expression was significantly higher in wFlu mosquitoes (Mann Whitney U tests: 10645 - P < 0.0001; 10453 - P < 0.0001). A further 2 contigs were predicted to be downregulated by wFlu, and the first of these, comp15178_c0_seq1, demonstrated that pattern (Mann Whitney U test: P < 0.0001). However for the second, comp14155_c0_seq1, expression levels were not significantly different due to the presence of Wolbachia (Unpaired t test: t = 0.1196 P = 0.9057). All four of those contigs were predicted to display differential expression in both the bayseq and DESeq2 analyses. Only low levels of comp14155_c0_seq1 were detected during sequencing, which could explain the lack of differential expression observed during RT-qPCR. The final 2 contigs that were examined were predicted to have equivalent expression levels between treatments, and both fit that pattern (Unpaired t tests: comp2025_c0_seq1 - t = 1.500, P = 0.1449; comp2041_c0_seq1 - t = 0.7994, P = 0.4308). Further validation of the differential expression of these genes using different techniques will likely prove valuable for future research.

Fig. 3
figure3

Confirmation of contig differential expression using RT-qPCR. 4 differentially expressed contigs (a-d), and 2 non-differentially expressed contigs (e-f) were selected at random, and expression levels were quantified via RT-qPCR. 5/6 contigs performed as expected, while one that was expected to show higher expression in Wolbachia-infected mosquitoes did not

Full size image

GO information for the annotated contigs was retrieved from Blast2GO, or repositories of transcriptional data such as FlyBase and VectorBase for some genes where Blast2GO produced no information (see methods). GO information could not be found for 28 upregulated contigs, and 17 downregulated contigs. The remaining upregulated contigs were associated with 286 GO terms, of which 93 terms had multiple hits. While downregulated contigs were associated with 190 GO terms, and 62 GO had multiple hits (Additional file 4). These GO terms and GenBank annotations were used to group the differentially expressed contigs based on their putative functions, and this information was used to develop profiles of the transcriptional changes that occurred both with (Table 2) and without (Table 3) wFlu infection. Some contigs had more than one annotated function, and are listed in multiple categories.

Table 2 Significantly upregulated contigs
Full size table
Table 3 Significantly downregulated contigs
Full size table

Immune stimulation and suppression

Transinfection with Wolbachia in mosquito species has typically led to widespread increases in the expression of immune genes, including those involved in the Toll and IMD immune pathways, and a large number of antimicrobial peptides [10, 11, 14, 20]. This immune activation has been linked to interference of dengue viruses in mosquitoes [10, 14, 15], and immune suppression has been hypothesized as a potential cause of the enhancement of Plasmodium infection seen in some native and transient Wolbachia infections [28, 29]. We observed no changes in the transcription of genes in the Toll or IMD pathways, including genes such as rel1, or rel2, which might have been indicative of a change in pathway regulation. We also saw no changes in the expression of any of the antimicrobial peptides commonly affected by Wolbachia. The absence of systemic immune activation is common amongst native Wolbachia infections, and may be symptomatic of increased tolerance on the part of the host, and reduced pathogenicity on the part of the symbiont [46, 47].

We did observe differential expression of 6 contigs directly involved in mosquito immunity, 4 upregulated (2 cell wall hydrolases, a galactose specific c-type lectin, and a gram negative bacteria binding protein) and 2 downregulated (a mucin like protein, and a galactose specific c-type lectin). These genes are typically associated with bacterial binding and degradation, but many have also been linked to Plasmodium infection and could have contributed to the enhancement of P. gallinaceum infection in Ae. fluviatilis [32]. Gram negative bacteria binding proteins, for instance can have a broader role in immune stimulation, and are involved in the response to Plasmodium infection in An. gambiae [48]. Some galactose-specific c-type lectins are highly differentially expressed by Wolbachia [11], their expression is stimulated by Plasmodium infection, and they have been shown to protect Plasmodium against melanisation by the host immune system [49]. There are also examples of mucin like genes and cell wall hydrolases that are critical to Plasmodium development in mosquitoes [50, 51].

We also identified other differentially expressed genes have also been linked with insect immunity or Plasmodium infection. Plasmodium infection is dependent on glycolysis, and the metabolism of amino acids and lipids [52], all affected by wFlu infection (see below). There was evidence of altered sugar metabolism and transport, which could promote the development of Plasmodium [53]. Additionally, there was increased expression of a β-hexosaminidase, and these genes have a peptidoglycan hydrolase effect, can be bactericidal, and have a broader immune function that could facilitate interaction with Plasmodium [54]. Finally, we also saw upregulation and downregulation of serine proteases, which are involved in both digestion and immunity. Serine proteases have been linked to Plasmodium growth and development, invasion of host cells, and can protect the parasite against the host oxidative stress response [55, 56]. Any of these transcriptional changes could be linked to the enhancement of P. gallinaceum, and therefore merit further investigation.

Redox process

One interesting characteristic of Wolbachia infection is the breakdown of redox homeostasis, which has been seen across multiple infected hosts [57, 58]. This effect often manifests through the induction of reactive oxygen species (ROS), and altered expression of genes involved in oxidative stress response, which occur with both native Wolbachia associations and transinfections [10, 12, 58], but it is unclear if they represent an immune response to Wolbachia, or normalization of the redox processes altered by infection [37]. Increased oxidative stress has been linked to pathogen interference against viruses and Plasmodium [10, 11, 20], and may be part of Wolbachia’s immune evasion strategy [37]. Similarly, higher oxidative stress levels are an important part of the anti-Plasmodium immune response [59].

In Ae. fluviatilis we saw evidence of an altered redox response due to wFlu infection, in the form of higher levels of a short-chain dehydrogenase, an oxidoreductase, which can induce ROS, and 3 cytochromes p450, which act as oxidases [60, 61]. Wolbachia have also been linked with iron metabolism and storage, and this can influence key physiological traits such as fecundity [62–64]. wFlu infection induced higher levels of neuferricin, a protein that binds iron-rich heme, which is interesting given that Wolbachia produce enzymes involved in heme synthesis [38], and may utilise it as an energy source [36]. Likewise, heme is essential to the development of Plasmodium in mosquitoes [65]. Given that wFlu appears to alter redox homeostasis, it would be interesting to see if wFlu infection induces higher ROS levels. It should be noted that tetracycline treatment impacts mitochondrial function, which can lead to changes in insect oxidative stress response [66]. While we did use tetracycline to clear the wFlu infection, our experiments were performed more than 2 years (approximately 30 generations) after antibiotic treatment.

Metabolism

We observed that 33 genes associated with metabolism and digestion were differentially expressed as a result of wFlu infection. Wolbachia demonstrate many nutrition-based physiological changes in different insect hosts, including mutualism through nutritional provision [40, 67], better performance under nutritional stress [62], and competition for nutrients leading to effects on host fecundity and fertility, and on pathogen interference [68, 69]. These processes are linked to a variety of nutrients and micronutrients including amino acids, iron, and flavin adenine dinucleotide, which indicates that Wolbachia interact with a wide range of host metabolic pathways.

Wolbachia infections in Ae. aegypti show differential expression of large numbers of digestive enzymes, particularly serine proteases and trypsins, which may indicate cannibalisation of host resources to promote Wolbachia propagation. Infection with wFlu increased the expression of 3 serine proteases, and 2 trypsins, and decreased levels of two other serine proteases and 2 zinc metallopeptidases. While the number of affected genes were fewer than seen for the Ae. aegypti transinfections it is sufficient to suggest that a similar process operates in Ae. fluviatilis. Likewise, wFlu infection elevated levels of brain chitinase and chia, an enzyme involved in breakdown of the peritrophic matrix, which may be further evidence of an altered digestive process, and could compromise the integrity of the peritrophic matrix and potentially facilitate Plasmodium invasion.

We observed that wFlu infection had a broad effect on several aspects of host metabolism, including carbohydrate and lipid metabolism, both areas where Wolbachia is lacking key biosynthesis genes, and where Wolbachia transinfection alters transcription in Ae. aegypti [11, 38, 70]. Previous work indicates that wFlu infection leads to elevated levels of glycogen [71], a major carbohydrate reserve, in developing embryos. While it is unclear if this effect occurs in adults, we did see elevated levels of β-glucosidase, which encodes an enzyme involved in the breakdown of complex carbohydrates including glycogen. Similarly, β-hexosaminidase b levels were also higher, with this enzyme involved in the hydrolysis of amino-sugars. Wolbachia heavily utilise host lipids and sterols, and alter the cellular lipid profile [68, 72, 73], potentially to serve bacterial propagation [74], with similar processes underlying Plasmodium development [75]. Lipid metabolism genes affected by wFlu had a diverse range of functions and include glucosyl glucuronosyl transferases, enzymes required for glucuronidation, which has a role in the metabolism of fatty acids and other compounds, and also pancreatic triacylglycerol lipases, which are used to digest triacylglycerides and process fats. Interestingly, wFlu also increased levels of phosphatidylcholine-sterol acyltransferase, which is involved in sterol and phospholipid metabolism, and may be utilised by Wolbachia to fabricate its membranous vacuoles. wFlu infection also decreased levels of diacylglyceride (DAG) synthesis protein, suggesting that the infection likely interferes with DAG synthesis, and therefore lipid transport as DAGs are prominent transport lipids. Likewise, changes in the expression of these genes could potentially promote P. gallinaceum development.

Infection also induced expression of chorion peroxidase, an enzyme involved in ovarian follicle maturation, and decreased the levels of a vitellogenic carboxypeptidase, which is a yolk protein produced in the fat body [76]. The effect on these genes likely arises due to the strong presence of wFlu in the ovaries, and may contribute to further differences in egg nutritional content, as is seen with glycogen [71].

Membranes

We observed that infection with wFlu also affected 11 transcripts involved in membrane structure, or transmembrane transport. This included an ompa-like protein that is part of the Wolbachia membrane, which was indicative of Wolbachia replication. We also saw the upregulation of membrane proteins involved in the transport of glucose and cholesterol through niemann-pick type c2 protein, and these could potentially be involved in altering Plasmodium infection dynamics, given the varied metabolic requirements of the parasite [77]. This is further evidence of a co-opting of metabolic machinery, as Wolbachia has been strongly linked to carbohydrate metabolism as an energy source [38], and cholesterol, as a critical component of the vacuoles that surround the bacteria [78].

Infection also caused the upregulation of 3 membrane-bound chemosensory receptors potentially linked to Ae. fluviatilis behaviour. These three include sulfakinin, a digestive regulator in insects [79], and 2 taste receptors, including an ionotropic receptor and a gustatory receptor linked to sugar feeding [80]. Wolbachia infection can influence host olfaction at the behavioural [81, 82] and molecular level [11], and in the case of wMelPop, also dramatically alters feeding behaviour [83]. These effects could possibly contribute to a minor change in olfactory response or potentially behaviour in response to wFlu infection. In contrast, only 3 membrane-related genes were differentially expressed in the Tet dataset, including an ATP-dependent transporter molecule of bacterial origin, and two membrane receptors with putative roles in immune signalling, and membrane-protein interactions.

Physiology

We saw that wFlu upregulated several genes with a potential role in neurological function and development. These included calcium and sodium channels, a major allergen, a fasciculation and elongation protein, which is involved in nerve signal transduction, and neuroblast formation protein, and BMP-induced dendrite morphogenesis factor, which are both involved in neural development. While wFlu does not heavily infect many non-reproductive tissues, it is present in the head at low density [15], and so could affect neurological function. Both wMel and wMelPop infect neural tissues [15, 17], and wMelPop has a pronounced effect, causing neurological degradation in Drosophila melanogaster [84]. Wolbachia also affect hormone levels and contribute to host behavioural changes [85–87]. The overall physiological effect of wFlu infection is unclear, however given that genes involved in neurological development, signalling and neurotransmitter trafficking and release were all upregulated by wFlu, there could be critical effects.

While wFlu does infect Ae. fluviatilis heads, it is not found in the omatidia cells of the eye, unlike wMelPop [15], yet wFlu infection decreased the expression of 4 photoreceptor proteins involved in phototransduction, and 3 visual receptors, indicating that wFlu may influence host visual perception. There have been no categorised effects of Wolbachia on host visual process, and no similar genes were affected by wMel or wMelPop [11]. Potential physiological consequences of these changes could be a decreased sensitivity to light, which would be disadvantageous to the mosquito, or may indicate decreased activity during low level light conditions [88–90].

Cellular processes

There was also evidence that wFlu infection altered common cellular processes by influencing the expression of genes involved in DNA repair and replication, and DNA packaging (Tables 2 and 3), with similar genes affected by other Wolbachia infection in mosquitoes - particularly wMelPop [11]. Other contigs were also linked to the processes of transcription and translation, and more specifically to mRNA processing, snRNA and rRNA processing, and protein folding, which is not unexpected given that Wolbachia produces and affects the production of small RNAs [91, 92].

We observed that wFlu altered the expression of chromobox protein homolog 1, a heterochromatin protein that could potentially be involved in epigenetic silencing of gene expression, or in chromosome integrity. Likewise, we observed altered levels of histone 2b, which is involved in DNA packaging as a key component of the nucleosome, and in gene silencing. Wolbachia infection in Ae. aegypti mosquitoes is known to affect cytosine methylation and gene silencing [93], and our results could indicate that there could potentially be a similar influence on silencing through histone modification. In that study, Wolbachia strongly altered the methylation profile of genes involved in membrane function including transport, and here we observed that a similar group of genes were affected during wFlu infection in Ae. fluviatilis.

Reproductive manipulation in Wolbachia-infected insects has also been associated with histones and chromatin. For instance, CI has been linked to delays in the deposition of histones in the male pronucleus [94], while male-killing has been linked to defective chromatin packaging, and altered chromosome behaviour [95]. Critically, the molecular effects of these processes are associated with adult males (CI), or in early stage embryos, and while our data were generated from adult females, these results may suggest that there are broader effects of Wolbachia infection on cellular processes related to DNA packaging, chromatin formation, and the regulation of gene expression in mosquitoes.

Comparison to other Wolbachia infections in mosquitoes

As wFlu is a native Wolbachia infection, and has likely undergone a lengthy period of co-evolution with Ae. fluviatilis, we were quite surprised to discover that infection affected the expression of genes involved in a wide range of host biological processes. The scope of the transcriptional changes associated with infection (257 genes), appears to be similar to that of wMel in Ae. aegypti (327 genes), although that study utilized a microarray rather than RNA-Seq, and the difference in sensitivity of the techniques may be a factor [11]. This similarity is interesting given that the latter association is both relatively novel, and is a transinfection that produces a broader range of physiological effects, and has greater bacterial density [15, 17, 32]. In contrast, transinfection of Ae. aegypti with the more virulent wMelPop strain has been demonstrated to alter the expression of 244 genes in one study [14], and 2723 in another [11], with the studies utilizing different methods of analysis. A larger set of differentially expressed genes would likely reflect the pathogenicity associated with wMelPop infection.

While these wMel and wMelPop data come from transinfections in a different species, we sought to determine if there was a core set of genes or functional groups of genes that were differentially expressed in all three associations, as this might explain something fundamental about the nature of Wolbachia infections in mosquitoes. 75 of the upregulated contigs, and 54 of the downregulated contigs were homologous to previously described Ae. aegypti transcripts in VectorBase. However, only 42 of these were also differentially expressed during either wMel or wMelPop infection. Interestingly, the majority of these were upregulated by wMel and wMelPop, even if they were downregulated by wFlu (Additional file 5), with this difference likely due to the relative novelty of the former transinfections.

Taking a broader approach, we then compared types of genes affected by all three strains, for example looking at all serine proteases, rather than a specific serine protease. Forty-six of the same types of genes were affected by both wMelPop and wFlu, with 33 upregulated and 13 downregulated (Additional file 5). For wMel and wFlu, 15 genes of the same type were affected, with 12 upregulated and 2 downregulated, all of which were also affected by wMelPop (Table 4). Genes that are upregulated in the presence of Wolbachia reflect 5 key areas, protein and fat metabolism, redox process, membrane transport, DNA/RNA processing, and bacterial recognition, all of which have been previously characterised in Wolbachia infections [10, 11, 38, 72, 96]. Genes that were downregulated include cuticle proteins and carboxypeptidases, which are involved in protein digestion. Cuticle proteins can be downregulated in response to tetracycline treatment in Wolbachia-infected Brugia malayi worms [97]. They can also be downregulated in response to viral infections, and potentially play a role in host resistance to infection [98]. These processes likely contribute to making the mosquito host environment more favourable for Wolbachia propagation, and thus represent areas of interest for further study.

Table 4 Core genes affected by Wolbachia infection in mosquitoes
Full size table

Conclusions

Here, we present the transcriptome of the mosquito Ae. fluviatilis, and consider the transcriptomic effects of its native Wolbachia strain, wFlu. Previous results suggest that wFlu infects host tissues at relatively low densities, causes incomplete CI and has no observable fitness cost [15, 32], in accordance with theories that suggest native Wolbachia strains have lost bacterial density and pathogenicity during long periods of co-evolution, and the development of tolerance on the part of their hosts [26]. Our data indicated that wFlu infection led to the differential expression of 257 genes, and while the scale of these changes was not as extreme as what is sometimes seen with Wolbachia transinfections in Ae. aegypti [10, 11, 14], the effect was still broad in scope and encompassed a wide range of biological processes, many of which are held in common with Wolbachia infections in other mosquitoes. Metabolic effects of wFlu infection appear to be particularly prominent [71], especially those associated with protein and lipid metabolism, and it is possible that the strain maintains the characteristic of a Jekyll and Hyde infection by both parasitising and providing key metabolites [26, 99]. And our results suggest that native strains such as wFlu likely have a greater impact on mosquito host biology than previously thought.

Critically, we did not see evidence that wFlu infection activated or suppressed the immune pathways typically associated with Wolbachia or Plasmodium infection [100, 101]. However, we did observe changes to a range of genes involved in immunity, oxidative stress and metabolism that have previously been associated with Plasmodium infection [59, 102], and could feasibly play a role in the enhancement of P. gallinaceum. Further molecular changes contributing to enhancement could be induced by infection or blood feeding, or under certain conditions associated with infections in the vertebrate host, and these could explain why enhancement of P. gallinaceum does not occur consistently [32]. Plasmodium enhancement is more prominent amongst transient artificial Wolbachia infections, where there is typically more extreme manipulation of host immunity and bacterial density that could contribute to the phenotype [103], and where temperature appears to be a key determinant [30]. While transient infections will likely never be utilised outside of the laboratory, determining the prevalence of enhancement amongst mosquito vectors naturally infected by Wolbachia, and identifying the causal mechanism of enhancement remain important issues going forward.

Methods

Mosquito material

Two mosquito lines were used in these experiments. The original, Wolbachia-infected Aedes fluviatilis (wFlu) line was derived from mosquitoes originally captured in Belo Horizonte in 1975 [32]. The mosquito colony was maintained in the laboratory until 2013 when a subset (Tet) was treated with tetracycline hydrochloride to remove the native Wolbachia infection, and then had their gut microbiota recolonized, as previously described [32]. Colony larvae were reared at low density in dechlorinated water, and were fed with fish food (Goldfish Colour, Alcon, Camboriú, Santa Catarina, Cat. No. 0504-2). Adults were maintained in low-density cages in a climate-controlled insectary (temperature: 27 ± 1 °C, RH: 70 ± 10%, photoperiod: 12 h light: 12 h dark), and provided 10% sucrose solution ad libitum. Mosquitoes used in experiments were maintained in small cylindrical cages (diameter – 16 cm, height – 18 cm) of approximately 80–90 individuals. Experiments were conducted more than 2 years after tetracycline treatment. In all experiments, the Tet line served as a Wolbachia-uninfected control line for the wFlu line, to study the effects of Wolbachia infection.

RNA extractions & library preparation

Six groups of samples were prepared for sequencing, with each group going on to form an independent library. Three groups each of 16 6-day old whole adult females from the Flu and Tet lines were collected and total RNA extracted using the Trizol® protocol according to manufacturer’s instructions (Invitrogen), for a total of 6 independent samples, with 3 biological replicates per treatment. Mosquitoes were fed only 10% sucrose prior to collection. RNA levels in each sample were quantified using a NanoDrop ND1000 (ThermoFisher Scientific). Sample degradation levels were checked by running a portion of the samples on a standard non-denaturing agarose gel containing bleach [104].

cDNA libraries were constructed using Illumina Truseq RNA Sample Preparation Kit (Illumina Inc.) according to the manufacturer instructions starting with 4 μg of total RNA. The library products were then sequenced using an Illumina MiSeq platform on a paired-end 300 bp run. After cleaning reads from adaptor sequences, the quality of the reads was assessed using the FastQC program (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Cleaned reads are available for download at the National Center for Biotechnology Information - Sequence Read Archive under the BioProject ID PRJNA320882.

De novo transcriptome assembly and contig annotation

Since no reference genome was available for Aedes fluviatilis, a de novo transcriptome assembly was built with Trinity (https://github.com/trinityrnaseq/trinityrnaseq/wiki) using default parameters. All six Illumina RNAseq datasets were combined in order to assemble a more reliable transcriptome, with a total of 19,919,299 paired-end, high quality reads. Contig sequences were searched for candidate proteins with TransDecoder (https://transdecoder.github.io/), again using standard parameters. The assembled contigs were annotated through local alignments with BlastX (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to the NCBI non-redundant (NR) and KEGG databases. BlastX parameters were set with an -e value of 1e-10. Blast2GO (https://www.blast2go.com/) was used to retrieve Gene Ontologies to annotated transcripts. Phylogenetic and divergence analyses were conducted using sequence data obtained during this study, or from UniProt (www.uniprot.org) or VectorBase (https://www.vectorbase.org). Methods and references for these analyses are described in Additional file 1.

Read mapping and differential gene expression

All Illumina paired-end reads libraries were mapped separately against the Ae. aegypti predicted transcriptome, available at VectorBase, and the Trinity assembled contigs, both with Bowtie2 (http://bowtie-bio.sourceforge.net/bowtie2/index.shtml) using the default parameters while configuring fragment length. The Integrative Genomics Viewer (http://www.broadinstitute.org/igv/) was used to visualize the reads that were mapped back to the assembled transcriptome. Read counts mapped to each transcript were acquired with a custom Pearl script (available upon request).

Two R packages from Bioconductor (www.bioconductor.org), baySeq and DESeq2, were selected in order to identify the contigs that were significantly differentially expressed due to the presence of Wolbachia. The baySeq method is more sensitive, but also carries a greater false positive call rate, hence differential expression was confirmed using DESeq2 [45, 105]. The baySeq FDR (false discovery rate) P value for multiple tests was set to 0.05. The DESeq2 adjusted P value was set to 0.01, and transcripts were filtered according to their log2FoldChange (higher than +1 or lower than -1). We then obtained an estimate of the average fold change for each of the differentially expressed contigs using FPKM (fragments per kilobase of exon per million fragments mapped). Contigs that were differentially expressed but had no BlastX hit were excluded from analysis.

GO terms for the remaining contigs were generated through Blast2GO (https://www.blast2go.com) and VectorBase, or, when data could not be discovered from these sources, through FlyBase (http://flybase.org) or InterPro (http://www.ebi.ac.uk/interpro/), based on the analysis of homologous genes. GO information from contigs that matched to the same BlastX hit was considered only once, to avoid double counting. These GO lists were compiled and used to create two lists, one for the Flu libraries with the native wFlu Wolbachia infection, and the other for the Tet libraries, where the Wolbachia had been removed. Contigs were grouped according to putative function, and these lists were then compared to develop a profile of the transcriptomic effects of wFlu on its mosquito host.

Differentially expressed contig functions and GOs were then compared against the transcriptomic data from wMel and wMelPop-infected Aedes aegypti arrays [11, 14], in order to develop a profile of types of transcriptional changes held in common between the three Wolbachia strains and their mosquito hosts. Contigs that had BlastX matches against Ae. aegypti were compared directly with Ae. aegypti expression data. Those contigs that matched to another species were compared directly against the Ae. aegypti genome in NCBI using BlastN (http://blast.ncbi.nlm.nih.gov/Blast.cgi). Any hits with a substantial match percentage (>80%) and a significant e-value were used for further analysis, using the first hit in a comparison in VectorBase while those that did not were not considered. This information was used to determine if any of the specific genes affected by wFlu infection were also affected by wMel or wMelPop infection.

Confirmation of differential expression

To assess the accuracy of the transcriptomic data set, six contigs were selected at random for expression analysis with RT-qPCR. Two of these contigs were indicated to have higher expression for Flu mosquitoes (AF10645; comp10645_c1_seq2 and AF10453; comp10453_c0_seq5), two had higher expression in Tet mosquitoes (AF15178; comp15178_c0_seq1 and AF14155; comp14155_c0_seq1), while two had equivalent expression between Flu and Tet mosquitoes (AF2025; comp2025_c0_seq1 and AF2041; comp2041_c0_seq1). The first four of these contigs were predicted to be differentially expressed through both the bayseq and DESeq2 analyses. Primers for these contigs were designed from the sequences generated during sequencing, which meant that they were suitable for cDNA. Primer sequences were designed using Primer 3 V0.4.0 (http://bioinfo.ut.ee/primer3-0.4.0/) to have a Tm of 55–60 °C and a product size range of 80–120 bp (Additional file 6).

Mosquitoes from both the Flu and Tet lines were reared as described above, and then females were collected individually at 6 days post eclosion. Total RNA was extracted and quantified as described above. First strand cDNA synthesis was conducted with 1 μg of total RNA from each sample using the M-MLV reverse transcriptase assay according to manufacturer’s instructions (Promega cat: C118A). cDNAs were then diluted 1:10 in nuclease-free water and stored at -30 °C. SYBR-based PCR was used to confirm the expression of each of the test genes, with 15 samples tested per mosquito line. Each gene was quantified in duplicate for all samples using the following mix: SYBR -5 μL, forward and reverse primers (10 μM) - 0.5 μL each, sterile RNase free water - 2 μL, sample 2 μL). RT-qPCR for samples was run on a LightCycler® 96 System (Roche) using the following profile: 10 min pre-incubation at 95 °C, 40 cycles of 15 s at 95 °C, 60 s at 60 °C, melt curve - 95 °C for 15 s, ramp from 60 °C to 95 °C at 1.6 °C/s. Expression values for each gene were normalised against actin1 expression, which had previously been demonstrated to be a good control gene for Ae. fluviatilis [32]. Mean normalised expression values for each gene were calculated using Q-Gene [106], and were compared statistically between Flu and Tet mosquitoes using Mann Whiney U tests (Prism v 6.0 g, Graphpad).

Abbreviations

ATP:

Adenine triphosphate

bp:

Base pairs

CI:

Cytoplasmic incompatibility

DAG:

Diacylglyceride

FPKM:

Fragments per kilobase of transcript per million mapped reads

GO:

Gene ontology

Kb:

Kilobase

KEGG:

Kyoto encyclopedia of genes and genomes

KO:

Kegg orthology

NR:

NCBI non-redundant database

RNA-Seq:

RNA sequencing (whole transcriptome shotgun sequencing)

ROS:

Reactive oxygen species

rRNA:

Ribosomal RNA

RT-qPCR:

Quantitative reverse transcriptase PCR

snRNA:

Small nuclear RNA

References

  1. 1.

    Zug R, Hammerstein P. Still a host of hosts for Wolbachia: analysis of recent data suggests that 40% of terrestrial arthropod species are infected. PLoS One. 2012;7(6), e38544.

  2. 2.

    Caragata EP, Dutra HL, Moreira LA: Exploiting intimate relationships: Controlling mosquito-transmitted disease with Wolbachia. Trends Parasitol. 2016

  3. 3.

    Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6(10):741–51.

  4. 4.

    McGraw EA, O’Neill SL. Beyond insecticides: new thinking on an ancient problem. Nat Rev Microbiol. 2013;11(3):181–93.

  5. 5.

    Vector-borne diseases. WHO Fact sheets. World Health Organization; 2016.

  6. 6.

    Hoffmann AA, Montgomery BL, Popovici J, Iturbe-Ormaetxe I, Johnson PH, Muzzi F, Greenfield M, Durkan M, Leong YS, Dong Y, et al. Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission. Nature. 2011;476(7361):454–U107.

  7. 7.

    Turelli M, Hoffmann AA. Rapid spread of an inherited incompatibility factor in California Drosophila. Nature. 1991;353(6343):440–2.

  8. 8.

    Zhang D, Lees RS, Xi Z, Gilles JR, Bourtzis K. Combining the Sterile Insect Technique with Wolbachia-Based Approaches: II--A Safer Approach to Aedes albopictus Population Suppression Programmes, Designed to Minimize the Consequences of Inadvertent Female Release. PLoS One. 2015;10(8), e0135194.

  9. 9.

    Molloy JC, Sinkins SP. Wolbachia Do Not Induce Reactive Oxygen Species-Dependent Immune Pathway Activation in Aedes albopictus. Viruses. 2015;7(8):4624–39.

  10. 10.

    Pan X, Zhou G, Wu J, Bian G, Lu P, Raikhel AS, Xi Z. Wolbachia induces reactive oxygen species (ROS)-dependent activation of the Toll pathway to control dengue virus in the mosquito Aedes aegypti. Proc Natl Acad Sci U S A. 2012;109(1):E23–31.

  11. 11.

    Rancès E, Ye YH, Woolfit M, McGraw EA, O’Neill SL. The relative importance of innate immune priming in Wolbachia-mediated Dengue interference. Plos Pathogens. 2012;8(2):e1002548.

  12. 12.

    Wong ZS, Brownlie JC, Johnson KN. Oxidative stress correlates with Wolbachia-mediated antiviral protection in Wolbachia-Drosophila associations. Appl Environ Microbiol. 2015;81(9):3001–5.

  13. 13.

    Hussain M, Lu G, Torres S, Edmonds JH, Kay BH, Khromykh AA, Asgari S. Effect of Wolbachia on replication of West Nile virus in a mosquito cell line and adult mosquitoes. J Virol. 2013;87(2):851–8.

  14. 14.

    Kambris Z, Cook PE, Phuc HK, Sinkins SP. Immune activation by life-shortening Wolbachia and reduced filarial competence in mosquitoes. Science. 2009;326(5949):134–6.

  15. 15.

    Moreira LA, Iturbe-Ormaetxe I, Jeffery JA, Lu GJ, Pyke AT, Hedges LM, Rocha BC, Hall-Mendelin S, Day A, Riegler M, et al. A Wolbachia symbiont in Aedes aegypti limits infection with Dengue, Chikungunya, and Plasmodium. Cell. 2009;139(7):1268–78.

  16. 16.

    van den Hurk AF, Hall-Mendelin S, Pyke AT, Frentiu FD, McElroy K, Day A, Higgs S, O’Neill SL. Impact of Wolbachia on infection with chikungunya and yellow fever viruses in the mosquito vector Aedes aegypti. PLoS Negl Trop Dis. 2012;6(11), e1892.

  17. 17.

    Walker T, Johnson PH, Moreira LA, Iturbe-Ormaetxe I, Frentiu FD, McMeniman CJ, Leong YS, Dong Y, Axford J, Kriesner P, et al. A non-virulent Wolbachia infection blocks dengue transmission and rapidly invades Aedes aegypti populations. Nature. 2011;476:450–5.

  18. 18.

    Ye YH, Woolfit M, Rances E, O’Neill SL, McGraw EA. Wolbachia-associated bacterial protection in the mosquito Aedes aegypti. PLoS Negl Trop Dis. 2013;7(8), e2362.

  19. 19.

    Dutra HL, Rocha MN, Dias FB, Mansur SB, Caragata EP, Moreira LA. Wolbachia Blocks Currently Circulating Zika Virus Isolates in Brazilian Aedes aegypti Mosquitoes. Cell Host Microbe. 2016;19(6):771–4.

  20. 20.

    Bian G, Joshi D, Dong Y, Lu P, Zhou G, Pan X, Xu Y, Dimopoulos G, Xi Z. Wolbachia invades Anopheles stephensi populations and induces refractoriness to Plasmodium infection. Science. 2013;340(6133):748–51.

  21. 21.

    Ferguson NM, Kien DT, Clapham H, Aguas R, Trung VT, Chau TN, Popovici J, Ryan PA, O’Neill SL, McGraw EA, et al. Modeling the impact on virus transmission of Wolbachia-mediated blocking of dengue virus infection of Aedes aegypti. Sci Transl Med. 2015;7(279):279ra237.

  22. 22.

    Frentiu FD, Zakir T, Walker T, Popovici J, Pyke AT, van den Hurk A, McGraw EA, O’Neill SL. Limited dengue virus replication in field-collected Aedes aegypti mosquitoes infected with Wolbachia. PLoS Negl Trop Dis. 2014;8(2), e2688.

  23. 23.

    Hoffmann AA, Iturbe-Ormaetxe I, Callahan AG, Phillips BL, Billington K, Axford JK, Montgomery B, Turley AP, O’Neill SL. Stability of the wMel Wolbachia Infection following invasion into Aedes aegypti populations. PLoS Negl Trop Dis. 2014;8(9), e3115.

  24. 24.

    McMeniman CJ, Lane RV, Cass BN, Fong AWC, Sidhu M, Wang YF, O’Neill SL. Stable introduction of a life-shortening Wolbachia infection into the mosquito Aedes aegypti. Science. 2009;323(5910):141–4.

  25. 25.

    Xi Z, Khoo CC, Dobson SL. Wolbachia establishment and invasion in an Aedes aegypti laboratory population. Science. 2005;310(5746):326–8.

  26. 26.

    Zug R, Hammerstein P. Bad guys turned nice? A critical assessment of Wolbachia mutualisms in arthropod hosts. Biol Rev Camb Philos Soc. 2014;90(1):89–111.

  27. 27.

    Zele F, Nicot A, Berthomieu A, Weill M, Duron O, Rivero A. Wolbachia increases susceptibility to Plasmodium infection in a natural system. Proc Biol Sci. 2014;281(1779):20132837.

  28. 28.

    Dodson BL, Hughes GL, Paul O, Matacchiero AC, Kramer LD, Rasgon JL. Wolbachia enhances West Nile virus (WNV) infection in the mosquito Culex tarsalis. PLoS Negl Trop Dis. 2014;8(7), e2965.

  29. 29.

    Hughes GL, Vega-Rodriguez J, Xue P, Rasgon JL. Wolbachia strain wAlbB enhances infection by the rodent malaria parasite Plasmodium berghei in Anopheles gambiae mosquitoes. Appl Environ Microbiol. 2012;78(5):1491–5.

  30. 30.

    Murdock CC, Blanford S, Hughes GL, Rasgon JL, Thomas MB. Temperature alters Plasmodium blocking by Wolbachia. Sci Rep. 2014;4:3932.

  31. 31.

    Tason de Camargo M, Krettli AU. Aedes fluviatilis (Lutz), a new experimental host for Plasmodium gallinaceum brumpt. J Parasitol. 1978;64(5):924–5.

  32. 32.

    Baton LA, Pacidonio EC, Goncalves DS, Moreira LA. wFlu: characterization and evaluation of a native Wolbachia from the mosquito Aedes fluviatilis as a potential vector control agent. PLoS One. 2013;8(3), e59619.

  33. 33.

    Multini LC, Marrelli MT, Wilke AB. Microsatellite loci cross-species transferability in Aedes fluviatilis (Diptera:Culicidae): a cost-effective approach for population genetics studies. Parasit Vectors. 2015;8:635.

  34. 34.

    Nene V, Wortman JR, Lawson D, Haas B, Kodira C, Tu ZJ, Loftus B, Xi Z, Megy K, Grabherr M, et al. Genome sequence of Aedes aegypti, a major arbovirus vector. Science. 2007;316(5832):1718–23.

  35. 35.

    Duplouy A, Iturbe-Ormaetxe I, Beatson SA, Szubert JM, Brownlie JC, McMeniman CJ, McGraw EA, Hurst GD, Charlat S, O’Neill SL, et al. Draft genome sequence of the male-killing Wolbachia strain wBol1 reveals recent horizontal gene transfers from diverse sources. BMC Genomics. 2013;14:20.

  36. 36.

    Gill AC, Darby AC, Makepeace BL. Iron necessity: the secret of Wolbachia’s success? PLoS Negl Trop Dis. 2014;8(10), e3224.

  37. 37.

    Zug R, Hammerstein P. Wolbachia and the insect immune system: what reactive oxygen species can tell us about the mechanisms of Wolbachia-host interactions. Front Microbiol. 2015;6:1201.

  38. 38.

    Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, Brownlie JC, McGraw EA, Martin W, Esser C, Ahmadinejad N, et al. Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: A streamlined genome overrun by mobile genetic elements. Plos Biol. 2004;2(3):327–41.

  39. 39.

    Iturbe-Ormaetxe I, Burke GR, Riegler M, O’Neill SL. Distribution, expression, and motif variability of ankyrin domain genes in Wolbachia pipientis. J Bacteriol. 2005;187(15):5136–45.

  40. 40.

    Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, Bhattacharyya A, Kapatral V, Kumar S, Posfai J, et al. The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. Plos Biol. 2005;3(4), e121.

  41. 41.

    Sutton ER, Harris SR, Parkhill J, Sinkins SP. Comparative genome analysis of Wolbachia strain wAu. BMC Genomics. 2014;15:928.

  42. 42.

    Woolfit M, Iturbe-Ormaetxe I, Brownlie JC, Walker T, Riegler M, Seleznev A, Popovici J, Rances E, Wee BA, Pavlides J, et al. Genomic evolution of the pathogenic Wolbachia strain, wMelPop. Genome Biol Evol. 2013;5(11):2189–204.

  43. 43.

    Duron O, Boureux A, Echaubard P, Berthomieu A, Berticat C, Fort P, Weill M. Variability and expression of ankyrin domain genes in Wolbachia variants infecting the mosquito Culex pipiens. J Bacteriol. 2007;189(12):4442–8.

  44. 44.

    Holt RA, Subramanian GM, Halpern A, Sutton GG, Charlab R, Nusskern DR, Wincker P, Clark AG, Ribeiro JM, Wides R, et al. The genome sequence of the malaria mosquito Anopheles gambiae. Science. 2002;298(5591):129–49.

  45. 45.

    Kvam VM, Liu P, Si Y. A comparison of statistical methods for detecting differentially expressed genes from RNA-seq data. Am J Bot. 2012;99(2):248–56.

  46. 46.

    Bourtzis K, Pettigrew MM, O’Neill SL. Wolbachia neither induces nor suppresses transcripts encoding antimicrobial peptides. Insect Mol Biol. 2000;9(6):635–9.

  47. 47.

    Chrostek E, Marialva MS, Yamada R, O’Neill SL, Teixeira L. High anti-viral protection without immune upregulation after interspecies Wolbachia transfer. PLoS One. 2014;9(6), e99025.

  48. 48.

    Warr E, Das S, Dong Y, Dimopoulos G. The Gram-negative bacteria-binding protein gene family: its role in the innate immune system of Anopheles gambiae and in anti-Plasmodium defence. Insect Mol Biol. 2008;17(1):39–51.

  49. 49.

    Osta MA, Christophides GK, Kafatos FC. Effects of mosquito genes on Plasmodium development. Science. 2004;303(5666):2030–2.

  50. 50.

    Berois M, Romero-Severson J, Severson DW. RNAi knock-downs support roles for the mucin-like (AeIMUC1) gene and short-chain dehydrogenase/reductase (SDR) gene in Aedes aegypti susceptibility to Plasmodium gallinaceum. Med Vet Entomol. 2012;26(1):112–5.

  51. 51.

    Groat-Carmona AM, Kain H, Brownell J, Douglass AN, Aly AS, Kappe SH. A Plasmodium α/β-hydrolase modulates the development of invasive stages. Cell Microbiol. 2015;17(12):1848–67.

  52. 52.

    Lamour SD, Straschil U, Saric J, Delves MJ. Changes in metabolic phenotypes of Plasmodium falciparum in vitro cultures during gametocyte development. Malar J. 2014;13:468.

  53. 53.

    Salcedo-Sora JE, Caamano-Gutierrez E, Ward SA, Biagini GA. The proliferating cell hypothesis: a metabolic framework for Plasmodium growth and development. Trends Parasitol. 2014;30(4):170–5.

  54. 54.

    Koo IC, Ohol YM, Wu P, Morisaki JH, Cox JS, Brown EJ. Role for lysosomal enzyme beta-hexosaminidase in the control of mycobacteria infection. Proc Natl Acad Sci U S A. 2008;105(2):710–5.

  55. 55.

    Rodrigues J, Oliveira GA, Kotsyfakis M, Dixit R, Molina-Cruz A, Jochim R, Barillas-Mury C. An epithelial serine protease, AgESP, is required for Plasmodium invasion in the mosquito Anopheles gambiae. PLoS One. 2012;7(4), e35210.

  56. 56.

    Sharma S, Jadli M, Singh A, Arora K, Malhotra P. A secretory multifunctional serine protease, DegP of Plasmodium falciparum, plays an important role in thermo-oxidative stress, parasite growth and development. FEBS J. 2014;281(6):1679–99.

  57. 57.

    Brennan LJ, Haukedal JA, Earle JC, Keddie B, Harris HL. Disruption of redox homeostasis leads to oxidative DNA damage in spermatocytes of Wolbachia-infected Drosophila simulans. Insect Mol Biol. 2012;21(5):510–20.

  58. 58.

    Zhang YK, Ding XL, Rong X, Hong XY. How do hosts react to endosymbionts? A new insight into the molecular mechanisms underlying the Wolbachia-host association. Insect Mol Biol. 2015;24(1):1–12.

  59. 59.

    Molina-Cruz A, DeJong RJ, Charles B, Gupta L, Kumar S, Jaramillo-Gutierrez G, Barillas-Mury C. Reactive oxygen species modulate Anopheles gambiae immunity against bacteria and Plasmodium. J Biol Chem. 2008;283(6):3217–23.

  60. 60.

    Esterhazy D, King MS, Yakovlev G, Hirst J. Production of reactive oxygen species by complex I (NADH:ubiquinone oxidoreductase) from Escherichia coli and comparison to the enzyme from mitochondria. Biochemistry. 2008;47(12):3964–71.

  61. 61.

    Raha S, Robinson BH. Mitochondria, oxygen free radicals, disease and ageing. Trends Biochem Sci. 2000;25(10):502–8.

  62. 62.

    Brownlie JC, Cass BN, Riegler M, Witsenburg JJ, Iturbe-Ormaetxe I, McGraw EA, O’Neill SL. Evidence for metabolic provisioning by a common invertebrate endosymbiont, Wolbachia pipientis, during periods of nutritional stress. PLoS Pathog. 2009;5(4), e1000368.

  63. 63.

    Kosmidis S, Missirlis F, Botella JA, Schneuwly S, Rouault TA, Skoulakis EM. Behavioral decline and premature lethality upon pan-neuronal ferritin overexpression in Drosophila infected with a virulent form of Wolbachia. Front Pharmacol. 2014;5:66.

  64. 64.

    Kremer N, Voronin DA, Charif D, Mavingui P, Mollereau B, Vavre F. Wolbachia interferes with ferritin expression and iron metabolism in insects. Plos Pathogens. 2009;5(10), e1000630.

  65. 65.

    Ke H, Sigala PA, Miura K, Morrisey JM, Mather MW, Crowley JR, Henderson JP, Goldberg DE, Long CA, Vaidya AB. The heme biosynthesis pathway is essential for Plasmodium falciparum development in mosquito stage but not in blood stages. J Biol Chem. 2014;289(50):34827–37.

  66. 66.

    Ballard JW, Melvin RG. Tetracycline treatment influences mitochondrial metabolism and mtDNA density two generations after treatment in Drosophila. Insect Mol Biol. 2007;16(6):799–802.

  67. 67.

    Moriyama M, Nikoh N, Hosokawa T, Fukatsu T. Riboflavin Provisioning Underlies Wolbachia’s Fitness Contribution to Its Insect Host. MBio. 2015;6(6):e01732–01715.

  68. 68.

    Caragata EP, Rances E, Hedges LM, Gofton AW, Johnson KN, O’Neill SL, McGraw EA. Dietary cholesterol modulates pathogen blocking by Wolbachia. Plos Pathogens. 2013;9(6), e1003459.

  69. 69.

    Caragata EP, Rances E, O’Neill SL, McGraw EA. Competition for amino acids between Wolbachia and the mosquito host, Aedes aegypti. Microb Ecol. 2014;67(1):205–18.

  70. 70.

    Atyame CM, Labbe P, Lebon C, Weill M, Moretti R, Marini F, Gouagna LC, Calvitti M, Tortosa P. Comparison of Irradiation and Wolbachia Based Approaches for Sterile-Male Strategies Targeting Aedes albopictus. PLoS One. 2016;11(1), e0146834.

  71. 71.

    da Rocha FM, Martins R, Pessoa Costa E, Pacidonio EC, Araujo de Abreu L, da Silva Vaz Jr I, Moreira LA, da Fonseca RN, Logullo C. The modulation of the symbiont/host interaction between Wolbachia pipientis and Aedes fluviatilis embryos by glycogen metabolism. PLoS One. 2014;9(6), e98966.

  72. 72.

    Molloy JC, Sommer U, Viant MR, Sinkins SP. Wolbachia modulates lipid metabolism in Aedes albopictus mosquito cells. Appl Environ Microbiol. 2016;82(10):3109–20.

  73. 73.

    Vollmer J, Schiefer A, Schneider T, Julicher K, Johnston KL, Taylor MJ, Sahl HG, Hoerauf A, Pfarr K. Requirement of lipid II biosynthesis for cell division in cell wall-less Wolbachia, endobacteria of arthropods and filarial nematodes. Int J Med Microbiol. 2013;303(3):140–9.

  74. 74.

    Herren JK, Paredes JC, Schupfer F, Arafah K, Bulet P, Lemaitre B. Insect endosymbiont proliferation is limited by lipid availability. Elife. 2014;3, e02964.

  75. 75.

    Atella GC, Bittencourt-Cunha PR, Nunes RD, Shahabuddin M, Silva-Neto MA. The major insect lipoprotein is a lipid source to mosquito stages of malaria parasite. Acta Trop. 2009;109(2):159–62.

  76. 76.

    Deitsch KW, Raikhel AS. Cloning and analysis of the locus for mosquito vitellogenic carboxypeptidase. Insect Mol Biol. 1993;2(4):205–13.

  77. 77.

    Labaied M, Jayabalasingham B, Bano N, Cha SJ, Sandoval J, Guan G, Coppens I. Plasmodium salvages cholesterol internalized by LDL and synthesized de novo in the liver. Cell Microbiol. 2011;13(4):569–86.

  78. 78.

    Cho KO, Kim GW, Lee OK. Wolbachia bacteria reside in host Golgi-related vesicles whose position is regulated by polarity proteins. PLoS One. 2011;6(7), e22703.

  79. 79.

    Zels S, Dillen S, Crabbe K, Spit J, Nachman RJ, Vanden Broeck J. Sulfakinin is an important regulator of digestive processes in the migratory locust, Locusta migratoria. Insect Biochem Mol Biol. 2015;61:8–16.

  80. 80.

    Arensburger P, Megy K, Waterhouse RM, Abrudan J, Amedeo P, Antelo B, Bartholomay L, Bidwell S, Caler E, Camara F, et al. Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics. Science. 2010;330(6000):86–8.

  81. 81.

    Caragata EP, Real KM, Zalucki MP, McGraw EA. Wolbachia infection increases recapture rate of field-released Drosophila melanogaster. Symbiosis. 2011;54(1):55–60.

  82. 82.

    Peng Y, Nielsen JE, Cunningham JP, McGraw EA. Wolbachia infection alters olfactory-cued locomotion in Drosophila spp. Appl Environ Microbiol. 2008;74(13):3943–8.

  83. 83.

    Turley AP, Moreira LA, O’Neill SL, McGraw EA. Wolbachia infection reduces blood-feeding success in the Dengue fever mosquito. Aedes aegypti Plos Neglect Trop D. 2009;3(9), e516.

  84. 84.

    Strunov A, Kiseleva E, Gottlieb Y. Spatial and temporal distribution of pathogenic Wolbachia strain wMelPop in Drosophila melanogaster central nervous system under different temperature conditions. J Invertebr Pathol. 2013;114(1):22–30.

  85. 85.

    Liu C, Wang JL, Zheng Y, Xiong EJ, Li JJ, Yuan LL, Yu XQ, Wang YF. Wolbachia-induced paternal defect in Drosophila is likely by interaction with the juvenile hormone pathway. Insect Biochem Mol Biol. 2014;49:49–58.

  86. 86.

    Moreira LA, Ye YH, Turner K, Eyles DW, McGraw EA, O’Neill SL. The wMelPop strain of Wolbachia interferes with dopamine levels in Aedes aegypti. Parasit Vectors. 2011;4:28.

  87. 87.

    Rohrscheib CE, Bondy E, Josh P, Riegler M, Eyles D, van Swinderen B, Weible 2nd MW, Brownlie JC. Wolbachia Influences the Production of Octopamine and Affects Drosophila Male Aggression. Appl Environ Microbiol. 2015;81(14):4573–80.

  88. 88.

    Jenkins AM, Muskavitch MA. Crepuscular Behavioral Variation and Profiling of Opsin Genes in Anopheles gambiae and Anopheles stephensi (Diptera: Culicidae). J Med Entomol. 2015;52(3):296–307.

  89. 89.

    Moon YM, Metoxen AJ, Leming MT, Whaley MA, O’Tousa JE. Rhodopsin management during the light-dark cycle of Anopheles gambiae mosquitoes. J Insect Physiol. 2014;70:88–93.

  90. 90.

    Rund SS, Gentile JE, Duffield GE. Extensive circadian and light regulation of the transcriptome in the malaria mosquito Anopheles gambiae. BMC Genomics. 2013;14:218.

  91. 91.

    Hussain M, O’Neill SL, Asgari S. Wolbachia interferes with the intracellular distribution of Argonaute 1 in the dengue vector Aedes aegypti by manipulating the host microRNAs. RNA Biol. 2013;10(12):1868–75.

  92. 92.

    Woolfit M, Algama M, Keith JM, McGraw EA, Popovici J. Discovery of putative small non-coding RNAs from the obligate intracellular bacterium Wolbachia pipientis. PLoS One. 2015;10(3), e0118595.

  93. 93.

    Ye YH, Woolfit M, Huttley GA, Rances E, Caragata EP, Popovici J, O’Neill SL, McGraw EA. Infection with a virulent strain of Wolbachia disrupts genome-wide patterns of cytosine methylation in the mosquito Aedes aegypti. PLoS One. 2013;8(6), e66482.

  94. 94.

    Landmann F, Orsi GA, Loppin B, Sullivan W. Wolbachia-mediated cytoplasmic incompatibility is associated with impaired histone deposition in the male pronucleus. PLoS Pathog. 2009;5(3), e1000343.

  95. 95.

    Riparbelli MG, Giordano R, Ueyama M, Callaini G. Wolbachia-mediated male killing is associated with defective chromatin remodeling. PLoS One. 2012;7(1), e30045.

  96. 96.

    Ikeya T, Broughton S, Alic N, Grandison R, Partridge L. The endosymbiont Wolbachia increases insulin/IGF-like signalling in Drosophila. Proc Biol Sci. 2009;276(1674):3799–807.

  97. 97.

    Ghedin E, Hailemariam T, DePasse JV, Zhang X, Oksov Y, Unnasch TR, Lustigman S. Brugia malayi gene expression in response to the targeting of the Wolbachia endosymbiont by tetracycline treatment. PLoS Negl Trop Dis. 2009;3(10), e525.

  98. 98.

    Colpitts TM, Cox J, Vanlandingham DL, Feitosa FM, Cheng G, Kurscheid S, Wang P, Krishnan MN, Higgs S, Fikrig E. Alterations in the Aedes aegypti transcriptome during infection with West Nile, dengue and yellow fever viruses. PLoS Pathog. 2011;7(9), e1002189.

  99. 99.

    Jiggins FM, Hurst GD. Rapid insect evolution by symbiont transfer. Science. 2011;332(6026):185–6.

  100. 100.

    Bahia AC, Kubota MS, Tempone AJ, Araujo HR, Guedes BA, Orfano AS, Tadei WP, Rios-Velasquez CM, Han YS, Secundino NF, et al. The JAK-STAT pathway controls Plasmodium vivax load in early stages of Anopheles aquasalis infection. PLoS Negl Trop Dis. 2011;5(11), e1317.

  101. 101.

    Garver LS, Dong Y, Dimopoulos G. Caspar controls resistance to Plasmodium falciparum in diverse anopheline species. PLoS Pathog. 2009;5(3), e1000335.

  102. 102.

    Liu K, Dong Y, Huang Y, Rasgon JL, Agre P. Impact of trehalose transporter knockdown on Anopheles gambiae stress adaptation and susceptibility to Plasmodium falciparum infection. Proc Natl Acad Sci U S A. 2013;110(43):17504–9.

  103. 103.

    Hughes GL, Koga R, Xue P, Fukatsu T, Rasgon JL. Wolbachia infections are virulent and inhibit the human malaria parasite Plasmodium falciparum in Anopheles gambiae. PLoS Pathog. 2011;7(5), e1002043.

  104. 104.

    Aranda PS, LaJoie DM, Jorcyk CL. Bleach gel: a simple agarose gel for analyzing RNA quality. Electrophoresis. 2012;33(2):366–9.

  105. 105.

    Xu X, Zhang Y, Williams J, Antoniou E, McCombie WR, Wu S, Zhu W, Davidson NO, Denoya P, Li E. Parallel comparison of Illumina RNA-Seq and Affymetrix microarray platforms on transcriptomic profiles generated from 5-aza-deoxy-cytidine treated HT-29 colon cancer cells and simulated datasets. BMC Bioinformatics. 2013;14 Suppl 9:S1.

  106. 106.

    Simon P. Q-Gene: processing quantitative real-time RT-PCR data. Bioinformatics. 2003;19(11):1439–40.

Download references

Acknowledgements

The authors wish to thank Dr Paulo Ribolla from UNESP for use of the Illumina MiSeq platform and the INCT-EM for the Real-Time PCR machine. We also wish to thank Alice Sabatino and Anna Paula Canuto for technical support. Finally, we thank Dr Victor Pylro and Dr Larissa Scholte from FIOCRUZ/Minas for assistance with the phylogenetic analyses.

Funding

This work was supported by FAPEMIG, CNPq, and CAPES. FSP received a CNPq post-doctoral fellowship (process number 168223/2014-7).

Availability of data and materials

All RNA sequencing data are available for download at the National Center for Biotechnology Information - Sequence Read Archive under the BioProject ID PRJNA320882.

Authors’ contributions

EPC participated in study design, performed the RT-qPCR validation assay, carried out comparative data analysis, and drafted the manuscript. FSP performed the sequence assembly, annotation, analysis of differential expression, phylogenetic analysis, and helped to draft the manuscript. LAB conceived and designed the study, prepared samples for sequencing, conducted the divergence analysis, and helped prepare the manuscript. JBLS prepared samples and performed the RT-qPCR validation assay. MHFS helped to design the study, and prepared RNA libraries for sequencing. LAM conceived and designed the study, drafted the manuscript and supervised the research. All authors contributed to and approved the final version of the manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Affiliations

  1. Grupo Mosquitos Vetores: Endossimbiontes e Interação Patógeno Vetor, Centro de Pesquisas René Rachou - Fiocruz, Belo Horizonte, Minas Gerais, Brazil

    • E. P. Caragata
    • , L. A. Baton
    • , J. B. L. Silva
    •  & L. A. Moreira

  2. Grupo de Informática de Biossistemas e Genômica, Centro de Pesquisas René Rachou - Fiocruz, Belo Horizonte, Minas Gerais, Brazil

    • F. S. Pais

  3. Instituto de Bioquímica Médica, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil

    • M. H. F. Sorgine

Authors

  1. Search for E. P. Caragata in:

  2. Search for F. S. Pais in:

  3. Search for L. A. Baton in:

  4. Search for J. B. L. Silva in:

  5. Search for M. H. F. Sorgine in:

  6. Search for L. A. Moreira in:

Corresponding author

Correspondence to L. A. Moreira.

Additional files

Additional file 1:

Divergence and phylogenetic analyses. Phylogenetic trees, and an estimate of the time of divergence for Aedes fluviatilis and other key mosquito species. (DOCX 1815 kb)

Additional file 2:

Contigs of Wolbachia origin. List of contigs identified as being from Wolbachia pipientis. (XLSX 59 kb)

Additional file 3:

Annotated, differentially expressed Ae. fluviatilis contigs. Table depicting all differentially expressed contigs, relevant annotations, and expression data. (XLSX 55 kb)

Additional file 4:

GO terms with multiple hits associated with wFlu infection. Table depicting GO terms that were associated with more than one differentially expressed contig, for Wolbachia-infected and -uninfected mosquitoes. (DOCX 124 kb)

Additional file 5:

Comparison of the transcriptomic effects of wFlu with wMel and wMelPop infection in Ae. aegypti. Two tables. The first depicts gene homologs that are differentially expressed in wFlu-infected Aedes fluviatilis and/or wMel-infected Aedes aegypti, and wMelPop-infected Aedes aegypti. The second depicts gene functions that are commonly differentially expressed between the three associations. (DOCX 154 kb)

Additional file 6:

List of primers used for RT-qPCR confirmation of differential expression. Oligonucleotide sequences for genes used in this work. (DOCX 63 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Aedes fluviatilis
  • Wolbachia
  • Transcriptome
  • RNA-Seq
  • Metabolism
  • Mosquito
  • Vector control
  • Symbiont
  • Oxidative stress

BMC Genomics

ISSN: 1471-2164

Contact us