Corrigendum: RNA-Seq Comparison of Larval and Adult Malpighian Tubules of the Yellow Fever Mosquito Aedes aegypti Reveals Life Stage-Specific Changes in Renal Function
- 1Department of Biology, New Mexico State University, Las Cruces, NM, USA
- 2Department of Computer Science, New Mexico State University, Las Cruces, NM, USA
- 3Department of Entomology, Ohio Agricultural Research and Development Center, The Ohio State University, Wooster, OH, USA
- 4National Center for Genome Resources, Santa Fe, NM, USA
- 5Institute of Applied Biosciences, New Mexico State University, Las Cruces, NM, USA
Introduction: The life history of Aedes aegypti presents diverse challenges to its diuretic system. During the larval and pupal life stages mosquitoes are aquatic. With the emergence of the adult they become terrestrial. This shifts the organism within minutes from an aquatic environment to a terrestrial environment where dehydration has to be avoided. In addition, female mosquitoes take large blood meals, which present an entirely new set of challenges to salt and water homeostasis.
Methods: To determine differences in gene expression associated with these different life stages, we performed an RNA-seq analysis of the main diuretic tissue in A. aegypti, the Malpighian tubules. We compared transcript abundance in 4th instar larvae to that of adult females and analyzed the data with a focus on transcripts that encode proteins potentially involved in diuresis, like water and solute channels as well as ion transporters. We compared our results against the model of potassium- and sodium chloride excretion in the Malpighian tubules proposed by Hine et al. (2014), which involves at least eight ion transporters and a proton-pump.
Results: We found 3,421 of a total number of 17,478 (19.6%) unique transcripts with a P < 0.05 and at least a 2.5 fold change in expression levels between the two groups. We identified two novel transporter genes that are highly expressed in the adult Malpighian tubules, which have not previously been part of the transport model in this species and may play important roles in diuresis. We also identified candidates for hypothesized sodium and chloride channels. Detoxification genes were generally higher expressed in larvae.
Significance: This study represents the first comparison of Malpighian tubule transcriptomes between larval and adult A. aegypti mosquitoes, highlighting key differences in their renal systems that arise as they transform from an aquatic filter-feeding larval stage to a terrestrial, blood-feeding adult stage.
Introduction
Aedes aegypti
Since the abandonment of the wide-spread vector control programs in the late 1960's, the yellow fever mosquito, Aedes aegypti, has rapidly re-emerged across the globe. The current range of this important disease vector puts ~2.5 billion people at risk for dengue fever, with an estimated 50–100 million cases per year. A. aegypti is also the primary vector of the viruses which cause yellow fever, chikungunya, and Zika (Nene et al., 2007; CDC, 2009; WHO, 2011; Rodriguez-Morales, 2015; Benelli and Mehlhorn, 2016; Lazear et al., 2016).
Osmoregulation in A. aegypti
Larval and adult mosquitoes inhabit very different environments with diverse requirements regarding osmoregulation and excretion. During the larval (and pupal) life stages, mosquitoes are aquatic organisms that are immersed in a hypotonic environment (Bradley, 1987; Marusalin et al., 2012). Post-eclosion, adult A. aegypti are terrestrial organisms that face the threat of rapid dehydration. Accordingly, the osmoregulatory system must change rapidly for the mosquito to acclimate to the new environment (Bradley, 1987; Piermarini, 2016). Adult females are presented with an additional challenge when feeding on vertebrate blood. Within minutes, a female A. aegypti can ingest more than her own body weight in blood which drastically reduces her mobility (Clements, 1992). In order to regain mobility and maintain osmotic balance, A. aegypti females excrete at least 40% of the volume contained in a blood meal within 1–2 h (Williams et al., 1983; Drake et al., 2010). Vertebrate blood is composed of ~85% water, and a mixture of ions, sugars, and proteins (Scanlon and Tina, 2007). Thus, ingesting large volumes of blood affects osmotic balance and leads to the release and production of toxic metabolic wastes (e.g., heme, NH3) during digestion.
Malpighian Tubules
The Malpighian tubules (MT) are the key excretory tissues in mosquitoes for osmotic balance and excretion of wastes, such as xenobiotics and excess nitrogen. In A. aegypti, there are five MT that attach to the alimentary canal at the midgut/hindgut junction. The tubules are comprised of two main cell types (principal cells and stellate cells) that form a single epithelial layer around the lumen (Beyenbach and Piermarini, 2011). Stellate cells intercalate between principal cells along the distal segment of the tubule. Principal and stellate cells are thought to work in conjunction to maintain proper osmotic balance within the hemolymph (Beyenbach et al., 2010; Beyenbach and Piermarini, 2011).
Mechanisms of Transepithelial Fluid Secretion
In the current model of transepithelial fluid secretion in the MT of A. aegypti, a proton motive force is created by the principal cells through a V-type H+ ATPase (V-ATPase) located in the apical membrane (Beyenbach et al., 2010; Beyenbach and Piermarini, 2011; Piermarini, 2016). The V-ATPase pumps protons into the lumen creating proton and voltage gradients to power other electrogenic exchange transporters and ion channels that move Na+ and K+ from the hemolymph into the lumen of the tubules. Only some of the molecular mechanisms have been identified and characterized so far. The stellate cells exchange intracellular HCO3 with Cl− from the hemolymph, then the Cl− is transported to the lumen, presumably by chloride channels in stellate cells and/or a K,Cl cotransporter (KCC) in principal cells (Beyenbach, 2003; Beyenbach et al., 2010; Beyenbach and Piermarini, 2011; Piermarini et al., 2011). It is unknown exactly how water moves from the hemolymph to the lumen, however in A. aegypti, several aquaporin (AQP) mRNAs are enriched in the MT (Pietrantonio et al., 2000; Drake et al., 2010, 2015). In MT of Anopheles gambiae, Prip (AQP2 in A. aegypti) immunoreactivity has been located on the basolateral membrane of stellate cells in the distal segment and principal cells of the proximal segment (Liu et al., 2011; Tsujimoto et al., 2013). RNAi-mediated knockdown of selected AQPs in adult female A. aegypti resulted in greatly reduced whole-mosquito diuresis capabilities (Drake et al., 2010), consistent with roles in transepithelial water transport in the MT. Mosquito diuresis is regulated by a number of neuropeptides including kinins, cardioacceleratory peptides (CAPs), adiokinetic hormones (AKHs), and corticotropin-releasing factor related (CRF related) hormones (Gade, 2004).
Excretion and Detoxification
In addition to maintaining osmotic balance, the MT have a role in detoxification and excretion of metabolic waste products. Uric acid is thought to be transported into portions of the MT where it accumulates into crystals and assists with transport of water into the tubules by contributing to the osmotic gradient (O'Donnell et al., 1983). Enrichment of cytochrome p450's and glutathione S-transferases (GST) in the MT of Drosophila, A. gambiae, and Aedes albopictus indicates that this organ may also play a role in detoxification (Beyenbach et al., 2010; Ingham et al., 2014; Esquivel et al., 2016).
This Paper
Previous studies have characterized the transcriptome of the MT in A. albopictus (before and after blood feeding) (Esquivel et al., 2014, 2016) and A. gambiae (larval, adult, before, and after blood feeding) (Overend et al., 2015), but transcriptomic studies in the MT of A. aegypti have not previously been performed. The goal of this study was to characterize the changes occurring in the MT of A. aegypti mosquitoes when switching from an aquatic to a terrestrial life. We conducted an in depth RNA-seq analysis of MT from 4th instar larvae and adult females 3 days post-eclosion.
We found a high number of transcripts differentially expressed in the MTs of adult mosquitoes and larvae. We also identified potential sodium and chloride channels that expand our current model of diuresis in mosquitoes.
Materials and Methods
Mosquito Culture
Mosquitoes from the A. aegypti Rockefeller (ROCK) strain were raised as previously described (Price et al., 2011) with the exception that they were fed as larvae solely on cat food (Special Kitty Original, Wal-Mart stores, Bentonville, AR).
Mosquito Dissection and RNA Isolation
MT were isolated from ~30 adult female or 30 larval (4th instar) A. aegypti in PBS and placed in Trizol® (Thermo Fisher Science). This was done in duplicate, creating two replicates for each group. Total RNA was then extracted from these tubule samples according to the manufacturers protocol (Chomczynski, 1993). A Nanodrop 1000 (Thermo scientific) was used to quantify total RNA concentration. RNA quality was assessed visually by an RNA gel.
Illumina Library Preparation
Four micrograms of total RNA (the recommended maximum) from each sample was used to prepare a cDNA library for each sample using the TruSeq RNA Sample Preparation Kit v2 (Illumina), according to the manufacturer's protocol for low-throughput sample preparation, with the following modifications (Tsujimoto et al., 2017). Differences in the protocol and our preparation procedure were: (1) using PCR strip tubes instead of PCR plates, (2) Elute, Prime, Fragment mix was thawed on ice and mixed into each well of the RBP plate on ice, and (3) Ligation mix was thawed on ice and mixed into each well of the ALP plate on ice.
The resulting libraries were quantified using a Nanodrop 1000 and Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA) at New Mexico State University, then sent for sequencing to the National Center for Genome Resources (Santa Fe, NM). The sequencing lab further analyzed the libraries and pooled them for sequencing on a single lane of Illumina HiSeq2000 1 × 100 bp reads.
Bioinformatics
Illumina reads were aligned using bowtie2 (v2.2.9) and Tophat2 (v2.1.1) to the A. aegypti reference transcripts available from Vectorbase (downloaded May. 5, 2016) (Lawson et al., 2009; Langmead and Salzberg, 2012). Each library was aligned separately. Alignments in each library for each transcript were tallied from bowtie2/Tophat2 results by htseq-count (Anders et al., 2015), then expression was compared between larvae and adults using the DESeq package in R (Anders and Huber, 2010). After alignment by bowtie2 and Tophat2, the FPKM (Fragments Per Kilobase of transcript per Million mapped reads) values and expression comparison of transcript isoforms were also calculated using Cuffdiff (Trapnell et al., 2010). A Pearson correlation test (Galton, 1877; Pearson, 1895) was used to compare DESeq results and Cuffdiff results with a correlation value of 0.98 and a P-value of lower than 2.2e-16.
Functional Clustering Analysis of Transcripts Up-/Down-Regulated
For the clustering analysis, we followed a similar protocol to Esquivel et al. (2014, 2016). In brief, transcripts were separated based on their expression in the MT of adults relative to MT of larvae (i.e., down-regulated or up-regulated in adult MT). Then, transcripts were submitted to a standalone BLASTn (version 2.2.31) analysis against the A. gambiae transcriptome (PEST strain transcript sequences, AgampP4.4 geneset, v1.00; https://www.vectorbase.org/). The “best hit ortholog” for each A. aegypti transcript was retrieved. A. gambiae orthologs were subjected to a Database for Annotation, Visualization, and Integrated Discovery (DAVID, version 6.7) functional clustering analysis (Niaid, 2006). Only those functional clusters with an enrichment value >1.3 (corresponding to P < 0.05) were retrieved and recorded for further analysis.
Results and Discussion
General Sequencing Results
Two cDNA libraries for each group were sequenced, resulting in 26.3 and 26.2 million reads respectively for the adult MT libraries. All raw sequence data generated were submitted to the sequence read archive (Leinonen et al., 2011) and accepted under accession number SRX1884160 for larvae and SRX539939 for adult. Sequencing of larval libraries resulted in 31.0 and 31.2 million reads. Seventy-three to seventy-six percent of the reads in each library successfully aligned to the A. aegypti reference transcriptome. We calculated the Pearson correlation coefficient for our sample repeats. Coefficients between adult: adult and larval: larval libraries were 0.996 and 0.991, respectively. Our replicates exhibited a high degree of similarity, meaning that the transcriptome of A. aegypti Malpighian tubules was stable across biological repeats and not highly variable (Bonizzoni et al., 2011). The Pearson correlation between adult and larval libraries was 0.732, indicating potential changes in transcript expression between larval and adult MT. Figure 1 shows a volcano plot visualizing the differences in gene expression between the larval and adult transcriptomes. The overall sequence results are shown in Table S1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/243161/fphys-08-00283-HTML/image_m/fphys-08-00283-g001.jpg)
Figure 1. Volcano plot generated with larval vs. adult transcript data. Transcripts altered at a less than log 2 fold level are in red, transcripts altered at a greater than log 2 fold level, but with a P < 0.05 are in blue. Transcripts with an expression level altered more than log 2 fold and have a P > 0.05 are in black.
General Comparison of Larval and Adult Transcriptomes
Due to the vastly different osmoregulatory requirements of aquatic larvae and terrestrial adults, we expected to see significant differences between the transcript expression in their MT.
The pathway analysis performed using DAVID v6.7 (Huang et al., 2009) showed a total of 23 functional clusters among the transcripts up-regulated in adult MT, and 10 clusters among the transcripts down-regulated in adult MT (Table 1). Clusters were then manually categorized based on their putative biological role: transcription and translation (9 clusters), cell signaling and voltage-gated channels (10 clusters), protein sorting and trafficking (5 clusters), organelle sorting and trafficking (7 clusters), cell signaling (1 cluster), and redox and detoxification (1 cluster; Table 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/243161/fphys-08-00283-HTML/image_m/fphys-08-00283-t001.jpg)
Table 1. Categorization of the enriched functional clusters of transcripts in the MT of larval and adult mosquitoes.
Functional Clusters from Up-Regulated Transcripts in Adult MT
Inspection of the “transcription and translation” and “protein sorting and trafficking” categories suggest that MT of adult A. aegypti have enhanced activities related with (1) transcription of genes, (2) synthesis of proteins, (3) post-translational modification of proteins, and import/export of proteins among cellular organelles. In addition, the “cell signaling and voltage-gated channels” category gather a significant number of transcripts associated with active movement of ions (e.g., Na+, K+) through membranes for osmotic cellular regulation and/or cell signaling.
Functional Clusters from Down-Regulated Transcripts in Adult MT
The category “organelle sorting and trafficking” and “protein sorting and trafficking” contains several transcripts associated with movement/re-arrangement of organelles by proteins (e.g., tubulin and myosin), suggesting that MT of non-blood fed females have reduced internal trafficking of cellular components. Furthermore, MT of adults are also showing a down-regulation of transcripts associated with cell signaling and redox/detoxification processes. These changes might correspond to the switch from an aquatic stage living in a relatively closed environment (i.e., a container) where toxins and wastes can accumulate and cannot be avoided until after pupal metamorphosis to the non-blood fed terrestrial/aerial stage living in a relatively open environment, where potential toxins and wastes can be avoided.
Most Highly Expressed Transcripts in Both Samples
Tables 2, 3 show the top twenty most highly expressed genes in larval and adult MTs, respectively. Five transcripts are found in both tables:
• AAEL018662: predicted to encode cytochrome c oxidase subunit I
• AAEL004851: unknown protein
• AAEL017413: unknown protein
• AAEL018672: no protein (tRNA)
• AAEL017096: unknown protein.
In larvae, cysteine/serine peptidase domain, synaptic vesicle protein, cytochrome c oxidase subunit I, metalloproteinase, putative, and 16 hypothetical proteins round out the top 20.
In the adult MTs, ribosomal proteins make up 4 of the 20 most highly expressed transcripts, while in larvae there are no ribosomal proteins in the twenty mostly highly expressed transcripts.
Most Highly Changed Transcripts between Larva and Adult
Table 4 shows the top 10 most upregulated genes in larvae compared to adult while Table 5 shows the top 10 most upregulated genes in adults compared to larvae. Four of 10 genes that are highly upregulated in the larval stage are annotated as hypothetical proteins while seven of the 10 adult genes are annotated. In adult MT, ribosomal genes were expressed at very high levels and exhibited very high degrees of upregulation compared to larvae.
Based upon these results (data presented in Tables 1–5), the adult MTs appear to be more heavily geared toward protein synthesis than those of larvae, which is consistent with results from the MT of non-blood fed A. albopictus (Esquivel et al., 2016).
Ion and Water Transport Machinery
The DAVID functional cluster analysis identified four groups associated with ion transport: “voltage-gated channel activity,” “transmembrane,” “ligand-gated ion channel activity,” and “voltage-gated potassium channel activity” (Table 1). Other transcripts associated with epithelial ion transport are also found within the “cell signaling and voltage-gated channels” clusters (Table 1). In addition, a manual search was performed to find transcripts associated with the current model of transepithelial fluid secretion in MT (Hine et al., 2014; Piermarini, 2016). We were able to identify at least one member of each gene family represented in this model in both larval as well as adult MT transcriptomes. Figure 2 shows the relative expression levels of some model-associated transcripts. This model includes sodium/potassium ATPases (Na/K ATPase), a chloride/bicarbonate anion exchanger (Cl/HCO3), sodium/proton antiporters (NHA1, NHA2), a sodium/proton exchanger (NHE3), inward rectifier potassium channels (Kir1 & 3), a sodium channel (NaC), sodium/potassium/chloride cotransporters (NKCC), chloride channels (ClC), potassium/chloride cotransporters (KCC), and V-ATPase subunits. Our analysis resulted in the identification of several novel paralogous transcripts expressed in the MT and also allowed the identification of some genes that are dominantly expressed. We found that many transcripts show life stage-dependent expression. Most transcripts that we analyzed had a greater level of expression in the larval MTs with only a few exceptions. For example, the putative sodium channel (NaC) AAEL014228, the only one of three NaC genes that is highly expressed in MTs, is 3.1 fold increased in adult relative to larval MTs, suggesting that it may play a role in the post-prandial diuresis of females.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/243161/fphys-08-00283-HTML/image_m/fphys-08-00283-g002.jpg)
Figure 2. Model of NaCl and KCl excretion. Expression levels of transporters hypothesized to be involved in excretion of NaCl and KCl in Aedes aegypti. Shown is the log scale transcript expression of each transporter we found to be expressed, larvae in red, adults in blue. The cellular localization of several of these transporters has not been determined and is therefore hypothetical. PC, principal cell; SC, stellate cell; v, vesicle.
We found that Kir1 (AAEL008932) was up-regulated in adults while Kir3 (AAEL001646) was up-regulated in larvae which confirms qPCR results published earlier (Yang et al., 2017). In adults, Kir1 is expressed on the basolateral membrane of stellate cells where it contributes to most of the transepithelial secretion of K+ (Piermarini et al., 2015). On the other hand, Kir3 is expressed in intracellular compartments of principal and stellate cells in adult MT, but its physiological role is unknown (Piermarini et al., 2015).
NHA2 (AAEL000137) was highly up-regulated in larvae (6.46x), consistent with a putative role in the absorption of luminal Na+ as proposed by Xiang et al. (2012).
The NKCC transcripts AAEL009888 and AAEL006180 were not differentially expressed.
Catalytic subunit A of V-ATPase (AAEL008787) was 2.96x higher in larvae, while many other subunits show moderate up-regulation in larvae (see Table 6). This is in contrast to microarray data from A. gambiae which indicates that in this mosquito species V-ATPase subunits are generally higher expressed in adults compared to larvae (Overend et al., 2015).
These results suggest there may be potential differences in transport rates of certain ions between larval and adult MT, due not only to changes in expression levels but also of transporter types which may have different kinetics.
As mentioned earlier, we found a transcript encoding the alpha-subunit of the sodium/potassium ATPase (AAEL008787) to be 2.96x more abundant in MTs of larvae vs. adults (Figure 2). The Na/K-ATPase has recently been considered to play an important role in transepithelial fluid secretion by MT of A. aegypti (Hine et al., 2014). Immunoreactivity for the alpha subunit localizes to the basolateral membranes of stellate cells in the distal segment and principal cells in the proximal segment (Patrick et al., 2006), and ouabain significantly inhibits K+, Cl−, and fluid secretion in isolated MTs (Hine et al., 2014). Thus, our finding is consistent with a potentially important role of this ATPase in transepithelial fluid and cation secretion by adult MTs.
The larval challenge is to retain NaCl and KCl and to remain hypertonic by excreting a dilute urine, in the face of a constant hypotonic external force. In contrast, the adult female mosquito is challenged with the sudden need to excrete a large amount of fluid following a blood meal. The larvae need to get rid of water constantly and retain ions, whereas the adult female needs to retain water to prevent dehydration until the moment she takes a blood meal. At that point she needs to get rid of water and ions very quickly (Patrick et al., 2001; Beyenbach, 2003; Beyenbach and Piermarini, 2011). Thus, the transcriptome of adult female MTs is likely primed with molecular mechanisms to fulfill her diuretic needs in response to a blood meal.
The model hypothesizes the existence of a basolateral sodium channel (NaC), which allows the entry of Na+ into the principal cells from the hemolymph, and an apical chloride channel, which allows the movement of Cl− from the stellate cells to the tubule lumen (Figure 2). In our analysis, we found one NaC (AAEL014228), and two ClC's (AAEL005950, AAEL001752) to be highly expressed in adult tubules (Figure 2). These two ClCs are orthologous to ClC-a and ClC-b in Drosophila, respectively. ClC-a (2x increase) and b (1.1x increase) were found to be more highly expressed in adults. ClC-a is known in Drosophila to be exclusively expressed in stellate cells and responds to diuretic hormone-induced secretion, while ClC-b is a housekeeping channel (Cabrero et al., 2014). There are three ClC-type channels in Drosophila, whereas in A. aegypti there are seven. The large expansion of genes in A. aegypti is potentially due to differences in selective pressure as a result of different feeding habits of fruit flies vs. mosquitoes (Wang et al., 2004).
Notably, a potassium-dependent sodium-calcium exchanger (NCKX1) was very highly expressed in larval and adult MT (AAEL004805). An ortholog of this transcript was also highly abundant in the MT of adult female A. albopictus and A. gambiae (Overend et al., 2015; Esquivel et al., 2016). NCKX1 may be associated with the response of MTs to kinins. That is, after stimulation with kinin peptides, there is an increase of intracellular Ca+2, which stimulates the secretion of Cl− (Yu and Beyenbach, 2002). A basolateral Na+/Ca2+ exchanger could provide a mechanism for lowering intracellular [Ca2+] after the influx of Ca2+ or maintaining low intracellular [Ca2+] before activation by kinins.
These previously unstudied transporters in A. aegypti represent potential mechanisms for disrupting diuresis, and exploiting as targets for novel insecticides development (Raphemot et al., 2013).
It is still unclear how the MT transport water from the hemolymph to the tubule lumen, but paracellular transport through septate junctions and transcellular transport through aquaporins are the most likely routes. RNAi studies by our group suggest that aquaporins play a major role (Drake et al., 2010, 2012; Benoit et al., 2014). In the present study, we found reads aligned against all six currently known aquaporins in the A. aegypti genome, as well as the B and C isoforms of aquaporin 5 (Eglp2) in both samples (Table 7). The pattern of expression observed for each aquaporin in the MT was very similar to the results of Drake et al. with aquaporins 1(Drip), 2 (Prip), 4 (Eglp1), and 5 (Eglp2) being the most highly expressed, while 3 (Bib) and 6 (Aqp12L) are expressed at very low levels (Drake et al., 2010). We found that aquaporins are generally upregulated in larvae with the exception of AQP6. However, our statistical analysis did not show significant changes except for AQP 1 and 3 (P < 0.05). These results may indicate that AQP expression overall is higher in MT of larvae vs. adults and that transcellular water transport is more important in MT of larvae than in those of adults.
Waste Excretion and Detoxification
In addition to maintaining fluid balance, the Malpighian tubules are a major site of detoxification and excretion of waste products, from pesticides to excess nitrogen and iron, through mechanisms such as cytochrome P450's and glutathione S-transferases (Folwell et al., 2006; Yang et al., 2007; Dow, 2009; Chahine and O'Donnell, 2011). We investigated which of these detoxification mechanisms are present by examining the expression of transcripts and comparing adult and larval expression.
Cytochrome P450 and Glutathione
We found reads aligning to 181 identified cytochrome P450 transcripts and 115 of them were differentially expressed with statistical significance (P < 0.05); of these we found 35 of 115 significant transcripts to be expressed at low levels (<100 reads) (Table S2) in both larvae and adults. Of the transcripts expressed at a low level, we found 17 (all exhibiting at least a 2 fold induction) to be more highly expressed in adults compared to larvae and 18 (17 exhibiting at least a 2 fold induction) to be more highly expressed in the larvae compared to adults. The most highly expressed cytochrome P450 transcript in adults was AAEL004054 and in larvae AAEL017539. We compared the changes in expression to the expression profiles of cytochrome P450's found to be involved in xenobiotic metabolism and pyrethroid resistance (Poupardin et al., 2010; Bariami et al., 2012). We found that many cytochrome P450's that are up-regulated in resistant strains were upregulated in larvae. However, there were many exceptions and it does not appear that the complement of larval up-regulation carry over to adult expression.
There were reads aligning to 34 transcripts which have been categorized as glutathione-interacting enzymes and are likely to be involved in detoxification pathways (Table S3). Eighteen of these were expressed above a low level (both larvae and adult >1,000). 14 were found to have statistically significant differential expression, with 10 more abundant in larvae, and four more abundant in adults. The most highly expressed detoxification-related transcript in MT of larvae and adults was AAEL001071 (GSTD5). We found the epsilon class of GSTs, which are capable of detoxification of DDT, to be much more highly expressed in adults (2–1,500 fold higher) (Lumjuan et al., 2007, 2011).
Our results indicate that mosquito larvae have a much more active detoxification system compared to adults, which makes physiological sense, because larvae are more likely to encounter persistent exposure to harmful wastes and xenobiotics in closed aquatic larval habitats.
Three transcripts relating to uric acid production, aldehyde oxidase, xanthine dehydrogenase and uricase in particular, were found to be significantly more abundant in MT of larvae vs. adults, by 13.9, 7.6, and 7.5 fold, respectively (Table S4). Arginase, which produces urea, was found to be evenly expressed in larvae and adults (Isoe and Scaraffia, 2013). However, expression of transcripts encoding arginase in MT of both larvae and adults is relatively low (>1,000 reads) while other enzymes in the urea cycle are much more highly expressed, and significantly more abundant in larvae (Figure 3). Results indicate that adults may be producing more urea as a waste product, while larvae are producing uric acid; representing possible differential nitrogen excretion between adult females and larvae (Kuzhivelil and Mohamed, 1998; von Dungern and Briegel, 2001; Scaraffia et al., 2008; Isoe and Scaraffia, 2013). This could be due to an aquatic vs. terrestrial lifestyle or larvae consuming cat food with protein while the adults consume sugar water. However, we did not detect significant differences in the amount of uric acid between MT of adults (non-blood fed) and larvae (data not shown).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/243161/fphys-08-00283-HTML/image_m/fphys-08-00283-g003.jpg)
Figure 3. Conversion of nitrogen in waste products. Log scale expression of transcripts for enzymes involved in the synthesis of nitrogen waste products in A. aegypti MTs. In general, in the classical urea cycle, transcripts are lower in larvae compared to adults, with the exception of argininosuccinate lyase. In the alternative urea pathway, transcripts are higher in larvae compared to adults with the exception of allantoicase. Shown is the log scale transcript expression of each transporter we found to be expressed, larvae in red, adults in blue.
We found 8 of 13 ferritin transcripts, including a ferritin precursor, to be highly expressed within the MT (>1,000 in both adult and larval MT, Table S5). Overall, ferritins were more abundant in the larval MT (1.4x). However, we would expect this to be reverse after the adults blood feed and this seems to be the case according to data in aeGEPUCI (Dissanayake et al., 2010). That is, in whole mosquitoes, several ferritins expressed at a low level increase rapidly post-blood meal and some of the highly expressed transcripts increase.
From these results, it appears that MT of adults potentially utilize stored iron to create heme peroxidases for immunity, dealing with oxidative stress or for other uses prior to blood meal. Using some stored iron in this way could assist in dealing with the influx of iron following a blood meal.
Conclusions
This study found significant changes between the MT transcriptomes of 4th instar larvae and adult female A. aegypti 3 days post-eclosion and identified several genes of potential interest for further functional analysis. Complete results of DEseq and Cuffdiff analysis are shown in Tables S6–S8.
Alterations in the expression of classical diuretic genes, such as ion and water transporters, reflect differing challenges in undertaking diuresis in aquatic and terrestrial environments, while changes in the abundance of glutathione S-transferases and cytochrome P450s likely reflect changes in the environment and diet. Changes in expression of enzymes of the nitrogen cycle and iron metabolism may reflect the dramatic change in diets between larval and adult mosquitoes.
We also have identified candidates for sodium and chloride channels predicted by the models of KCl and NaCl excretion of Hine et al. (2014).
Author Contributions
Performed the experiments: HD, FS, and YL; Analyzed the data: YL, PP, FS, and CE; Wrote the paper: IH and PP; Edited the manuscript: YL, PP, CE, FS, and IH.
Funding
This project was supported by the NIH grants SC1AI109055 and P20GM103451, and pilot funds from the National Center for Genome Resources (NCGR). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. David Price for his support of the early stages of this project.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fphys.2017.00283/full#supplementary-material
Table S1. Overall sequencing results. General statistics pertaining to the sequencing run. Number of reads and % aligning is the number of reads from the library which aligned to the reference transcripts using Bowtie2 and Tophat2.
Table S2. Transcripts related to Cytochrome P450. Expression of transcripts related to Cytochrome P450 in adult and larval MTs.
Table S3. Transcripts related to glutathione. Expression of transcripts related to glutathione in adult and larval MTs.
Table S4. Transcripts related to uric acid production. Expression of transcripts related to uric acid production in adult and larval MTs.
Table S5. Transcripts related to iron sequestration. Expression of transcripts related to iron sequestration and metabolism in adult and larval MTs.
Table S6. Completed DEseq result. The completed DEseq result table with annotation for each gene.
Table S7. Completed Cuffdiff Gene result. The complete Cuffdiff result table of each gene.
Table S8. Completed Cuffdiff isofrom result. The complete Cuffdiff result table for different isoforms.
References
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Anders, S., Pyl, P. T., and Huber, W. (2015). HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics 2, 166–169. doi: 10.1093/bioinformatics/btu638
Bariami, V., Jones, C. M., Poupardin, R., Vontas, J., and Ranson, H. (2012). Gene amplification, ABC transporters and cytochrome P450s: unraveling the molecular basis of pyrethroid resistance in the dengue vector, Aedes aegypti. PLoS Negl. Trop. Dis. 6:e1692. doi: 10.1371/journal.pntd.0001692
Benelli, G., and Mehlhorn, H. (2016). Declining malaria, rising of dengue and Zika virus: insights for mosquito vector control. Parasitol. Res. 115, 1747–1754. doi: 10.1007/s00436-016-4971-z
Benoit, J. B., Hansen, I. A., Szuter, E. M., Drake, L. L., Burnett, D. L., and Attardo, G. M. (2014). Emerging roles of aquaporins in relation to the physiology of blood-feeding arthropods. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 184, 811–825. doi: 10.1007/s00360-014-0836-x
Beyenbach, K. W. (2003). Transport mechanisms of diuresis in Malpighian tubules of insects. J. Exp. Biol. 206, 3845–3856. doi: 10.1242/jeb.00639
Beyenbach, K. W., and Piermarini, P. M. (2011). Transcellular and paracellular pathways of transepithelial fluid secretion in Malpighian (renal) tubules of the yellow fever mosquito Aedes aegypti. Acta Physiol. 202, 387–407. doi: 10.1111/j.1748-1716.2010.02195.x
Beyenbach, K. W., Skaer, H., and Dow, J. A. (2010). The developmental, molecular, and transport biology of Malpighian tubules. Annu. Rev. Entomol. 55, 351–374. doi: 10.1146/annurev-ento-112408-085512
Bonizzoni, M., Dunn, W. A., Campbell, C. L., Olson, K. E., Dimon, M. T., Marinotti, O., et al. (2011). RNA-seq analyses of blood-induced changes in gene expression in the mosquito vector species, Aedes aegypti. BMC Genomics 12:82. doi: 10.1186/1471-2164-12-82
Bradley, T. J. (1987). Physiology of osmoregulation in mosquitoes. Annu. Rev. Entomol. 32, 439–462. doi: 10.1146/annurev.en.32.010187.002255
Cabrero, P., Terhzaz, S., Romero, M. F., Davies, S. A., Blumenthal, E. M., and Dow, J. A. (2014). Chloride channels in stellate cells are essential for uniquely high secretion rates in neuropeptide-stimulated Drosophila diuresis. Proc. Natl. Acad. Sci. U.S.A. 111, 14301–14306. doi: 10.1073/pnas.1412706111
CDC (2009). Dengue Fever Fact Sheet. Available online at: https://www.cdc.gov/dengue/index.html
Chahine, S., and O'Donnell, M. J. (2011). Interactions between detoxification mechanisms and excretion in Malpighian tubules of Drosophila melanogaster. J. Exp. Biol. 214, 462–468. doi: 10.1242/jeb.048884
Chomczynski, P. (1993). A reagent for the single-step simultaneous isolation of RNA, DNA and proteins from cell and tissue samples. Biotechniques 15, 532–534, 536–537.
Dissanayake, S. N., Ribeiro, J. M., Wang, M. H., Dunn, W. A., Yan, G., James, A. A., et al. (2010). aeGEPUCI: a database of gene expression in the dengue vector mosquito, Aedes aegypti. BMC Res. Notes 3:248. doi: 10.1186/1756-0500-3-248
Dow, J. A. T. (2009). Insights into the Malpighian tubule from functional genomics. J. Exp. Biol. 212, 435–445. doi: 10.1242/jeb.024224
Drake, L. L., Boudko, D. Y., Marinotti, O., Carpenter, V. K., Dawe, A. L., and Hansen, I. A. (2010). The Aquaporin gene family of the yellow fever mosquito, Aedes aegypti. PLoS ONE 5:e15578. doi: 10.1371/journal.pone.0015578
Drake, L. L., Price, D. P., Aguirre, S. E., and Hansen, I. A. (2012). RNAi-mediated gene knockdown and in vivo diuresis assay in adult female Aedes aegypti mosquitoes. J. Vis. Exp. e3479. doi: 10.3791/3479
Drake, L. L., Rodriguez, S. D., and Hansen, I. A. (2015). Functional characterization of aquaporins and aquaglyceroporins of the yellow fever mosquito, Aedes aegypti. Sci. Rep. 5:7795. doi: 10.1038/srep07795
Esquivel, C. J., Cassone, B. J., and Piermarini, P. M. (2014). Transcriptomic evidence for a dramatic functional transition of the malpighian tubules after a blood meal in the Asian tiger mosquito Aedes albopictus. PLoS Negl. Trop. Dis. 8:e2929. doi: 10.1371/journal.pntd.0002929
Esquivel, C. J., Cassone, B. J., and Piermarini, P. M. (2016). A de novo transcriptome of the Malpighian tubules in non-blood-fed and blood-fed Asian tiger mosquitoes Aedes albopictus: insights into diuresis, detoxification, and blood meal processing. PeerJ 4:e1784. doi: 10.7717/peerj.1784
Folwell, J. L., Barton, C. H., and Shepherd, D. (2006). Immunolocalisation of the D. melanogaster Nramp homologue Malvolio to gut and Malpighian tubules provides evidence that Malvolio and Nramp2 are orthologous. J. Exp. Biol. 209, 1988–1995. doi: 10.1242/jeb.02193
Gade, G. (2004). Regulation of intermediary metabolism and water balance of insects by neuropeptides. Annu. Rev. Entomol. 49, 93–113. doi: 10.1146/annurev.ento.49.061802.123354
Galton, F. (1877). Typical laws of heredity. Nature 15, 492–495, 512–514, 532–533. doi: 10.1038/015492a0
Hine, R. M., Rouhier, M. F., Park, S. T., Qi, Z., Piermarini, P. M., and Beyenbach, K. W. (2014). The excretion of NaCl and KCl loads in mosquitoes. 1. Control data. Am. J. Physiol. Regul. Integr. Comp. Physiol. 307, R837–R849. doi: 10.1152/ajpregu.00105.2014
Huang, D. W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi: 10.1038/nprot.2008.211
Ingham, V. A., Jones, C. M., Pignatelli, P., Balabanidou, V., Vontas, J., Wagstaff, S. C., et al. (2014). Dissecting the organ specificity of insecticide resistance candidate genes in Anopheles gambiae: known and novel candidate genes. BMC Genomics 15:1018. doi: 10.1186/1471-2164-15-1018
Isoe, J., and Scaraffia, P. Y. (2013). Urea synthesis and excretion in Aedes aegypti mosquitoes are regulated by a unique cross-talk mechanism. PLoS ONE 8:e65393. doi: 10.1371/journal.pone.0065393
Kuzhivelil, B. T., and Mohamed, U. V. (1998). Allantoin and allantoic acid titre in the faeces and tissues of the developing larva of the moth, Orthaga exvinacea Hampson. Insect Biochem. Mol. Biol. 28, 979–986. doi: 10.1016/S0965-1748(98)00086-1
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Lawson, D., Arensburger, P., Atkinson, P., Besansky, N. J., Bruggner, R. V., Butler, R., et al. (2009). VectorBase: a data resource for invertebrate vector genomics. Nucleic Acids Res. 37, D583–D587. doi: 10.1093/nar/gkn857
Lazear, H. M., Stringer, E. M., and de Silva, A. M. (2016). The emerging Zika virus epidemic in the americas: research priorities. JAMA 315, 1945–1946. doi: 10.1001/jama.2016.2899
Leinonen, R., Sugawara, H., and Shumway, M. (2011). The sequence read archive. Nucleic Acids Res. 39, D19–D21. doi: 10.1093/nar/gkq1019
Liu, K., Tsujimoto, H., Cha, S. J., Agre, P., and Rasgon, J. L. (2011). Aquaporin water channel AgAQP1 in the malaria vector mosquito Anopheles gambiae during blood feeding and humidity adaptation. Proc. Natl. Acad. Sci. U.S.A. 108, 6062–6066. doi: 10.1073/pnas.1102629108
Lumjuan, N., Rajatileka, S., Changsom, D., Wicheer, J., Leelapat, P., Prapanthadara, L. A., et al. (2011). The role of the Aedes aegypti Epsilon glutathione transferases in conferring resistance to DDT and pyrethroid insecticides. Insect Biochem. Mol. Biol. 41, 203–209. doi: 10.1016/j.ibmb.2010.12.005
Lumjuan, N., Stevenson, B. J., Prapanthadara, L. A., Somboon, P., Brophy, P. M., Loftus, B. J., et al. (2007). The Aedes aegypti glutathione transferase family. Insect Biochem. Mol. Biol. 37, 1026–1035. doi: 10.1016/j.ibmb.2007.05.018
Marusalin, J., Matier, B. J., Rheault, M. R., and Donini, A. (2012). Aquaporin homologs and water transport in the anal papillae of the larval mosquito, Aedes aegypti. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 182, 1047–1056. doi: 10.1007/s00360-012-0679-2
Nene, V., Wortman, J. R., Lawson, D., Haas, B., Kodira, C., Tu, Z. J., et al. (2007). Genome sequence of Aedes aegypti, a major arbovirus vector. Science 316, 1718–1723. doi: 10.1126/science.1138878
Niaid, N. (2006). Database for Annotation, Visualization and Integrated Discovery (DAVID). Available online at: https://david.ncifcrf.gov/
O'Donnell, M. J., Maddrell, S. H. P., and Gardiner, B. O. C. (1983). Transport of uric acid by the Malpighian tubules of Rhodnius prolixus and other insects. J. Exp. Biol. 103, 169–184.
Overend, G., Cabrero, P., Halberg, K. A., Ranford-Cartwright, L. C., Woods, D. J., Davies, S. A., et al. (2015). A comprehensive transcriptomic view of renal function in the malaria vector, Anopheles gambiae. Insect Biochem. Mol. Biol. 67, 47–58. doi: 10.1016/j.ibmb.2015.05.007
Patrick, M. L., Aimanova, K., Sanders, H. R., and Gill, S. S. (2006). P-type Na+/K+-ATPase and V-type H+-ATPase expression patterns in the osmoregulatory organs of larval and adult mosquito Aedes aegypti. J. Exp. Biol. 209, 4638–4651. doi: 10.1242/jeb.02551
Patrick, M. L., Gonzalez, R. J., and Bradley, T. J. (2001). Sodium and chloride regulation in freshwater and osmoconforming larvae of Culex mosquitoes. J. Exp. Biol. 204, 3345–3354. Available online at: http://jeb.biologists.org/content/jexbio/204/19/3345.full.pdf
Pearson, K. (1895). Notes on regression and inheritance in the case of two parents. Proc. R. Soc. Lond. 58, 240–242. doi: 10.1098/rspl.1895.0041
Piermarini, P. (2016). Renal excretory processes in mosquitoes. Adv. Insect Phys. 51, 393–433. doi: 10.1016/bs.aiip.2016.04.003
Piermarini, P. M., Dunemann, S. M., Rouhier, M. F., Calkins, T. L., Raphemot, R., Denton, J. S., et al. (2015). Localization and role of inward rectifier K channels in Malpighian tubules of the yellow fever mosquito Aedes aegypti. Insect Biochem. Mol. Biol. 67, 59–73. doi: 10.1016/j.ibmb.2015.06.006
Piermarini, P. M., Hine, R. M., Schepel, M., Miyauchi, J., and Beyenbach, K. W. (2011). Role of an apical K,Cl cotransporter in urine formation by renal tubules of the yellow fever mosquito (Aedes aegypti). Am. J. Physiol. Regul. Integr. Comp. Physiol. 301, R1318–R1337. doi: 10.1152/ajpregu.00223.2011
Pietrantonio, P. V., Jagge, C., Keeley, L. L., and Ross, L. S. (2000). Cloning of an aquaporin-like cDNA and in situ hybridization in adults of the mosquito Aedes aegypti (Diptera: Culicidae). Insect Mol. Biol. 9, 407–418. doi: 10.1046/j.1365-2583.2000.00201.x
Poupardin, R., Riaz, M. A., Vontas, J., David, J. P., and Reynaud, S. (2010). Transcription profiling of eleven cytochrome P450s potentially involved in xenobiotic metabolism in the mosquito Aedes aegypti. Insect Mol. Biol. 19, 185–193. doi: 10.1111/j.1365-2583.2009.00967.x
Price, D. P., Nagarajan, V., Churbanov, A., Houde, P., Milligan, B., Drake, L. L., et al. (2011). The fat body transcriptomes of the yellow fever mosquito Aedes aegypti, pre- and post- blood meal. PLoS ONE 6:e22573. doi: 10.1371/journal.pone.0022573
Raphemot, R., Rouhier, M. F., Hopkins, C. R., Gogliotti, R. D., Lovell, K. M., Hine, R. M., et al. (2013). Eliciting renal failure in mosquitoes with a small-molecule inhibitor of inward-rectifying potassium channels. PLoS ONE 8:e64905. doi: 10.1371/journal.pone.0064905
Rodriguez-Morales, A. J. (2015). Zika: the new arbovirus threat for Latin America. J. Infect. Dev. Ctries 9, 684–685. doi: 10.3855/jidc.7230
Scanlon, V. C., and Tina, S. (2007). Essentials of Anatomy and Physiology. Version eBook 5th Edn: Philladelphia, PA: FA Davis Co.
Scaraffia, P. Y., Tan, G., Isoe, J., Wysocki, V. H., Wells, M. A., and Miesfeld, R. L. (2008). Discovery of an alternate metabolic pathway for urea synthesis in adult Aedes aegypti mosquitoes. Proc. Natl. Acad. Sci. U.S.A. 105, 518–523. doi: 10.1073/pnas.0708098105
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Tsujimoto, H., Hanley, K. A., Sundararajan, A., Devitt, N. P., Schilkey, F. D., and Hansen, I. A. (2017). Dengue virus serotype 2 infection alters midgut and carcass gene expression in the Asian tiger mosquito, Aedes albopictus. PLoS ONE 12:e0171345. doi: 10.1371/journal.pone.0171345
Tsujimoto, H., Liu, K., Linser, P. J., Agre, P., and Rasgon, J. L. (2013). Organ-specific splice variants of aquaporin water channel AgAQP1 in the malaria vector Anopheles gambiae. PLoS ONE 8:e75888. doi: 10.1371/journal.pone.0075888
von Dungern, P., and Briegel, H. (2001). Protein catabolism in mosquitoes: ureotely and uricotely in larval and imaginal Aedes aegypti. J. Insect Phys. 47, 131–141. doi: 10.1016/S0022-1910(00)00096-2
Wang, J., Kean, L., Yang, J., Allan, A. K., Davies, S. A., Herzyk, P., et al. (2004). Function-informed transcriptome analysis of Drosophila renal tubule. Genome Biol. 5:R69. doi: 10.1186/gb-2004-5-9-r69
WHO (2011). Yellow Fever Fact Sheet. Available online at: http://www.who.int/mediacentre/factsheets/fs100/en/
Williams, J. C., Hagedorn, H. H., and Beyenbach, K. W. (1983). Dynamic changes in flow rate and composition of urine during the post blood meal diuresis in Aedes aegypti. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 153, 257–266. doi: 10.1007/BF00689629
Xiang, M. A., Linser, P. J., Price, D. A., and Harvey, W. R. (2012). Localization of two Na+- or K+-H+ antiporters, AgNHA1 and AgNHA2, in Anopheles gambiae larval Malpighian tubules and the functional expression of AgNHA2 in yeast. J. Insect Phys. 58, 570–579. doi: 10.1016/j.jinsphys.2011.12.009
Yang, J., McCart, C., Woods, D. J., Terhzaz, S., Greenwood, K. G., ffrench-Constant, R. H., et al. (2007). A Drosophila systems approach to xenobiotic metabolism. Physiol. Genomics 30, 223–231. doi: 10.1152/physiolgenomics.00018.2007
Yang, Z., Statler, B.-M., Calkins, T. L., Alfaro, E., Esquivel, C. J., Rouhier, M. F., et al. (2017). Dynamic expression of genes encoding subunits of inward rectifier potassium (Kir) channels in the yellow fever mosquito Aedes aegypti. Comp. Biochem. Phys. B Biochem. Mol. Biol. 204, 35–44. doi: 10.1016/j.cbpb.2016.11.003
Keywords: mosquito, Aedes aegypti, Malpighian tubules, RNAseq, diuresis, detoxification
Citation: Li Y, Piermarini PM, Esquivel CJ, Drumm HE, Schilkey FD and Hansen IA (2017) RNA-Seq Comparison of Larval and Adult Malpighian Tubules of the Yellow Fever Mosquito Aedes aegypti Reveals Life Stage-Specific Changes in Renal Function. Front. Physiol. 8:283. doi: 10.3389/fphys.2017.00283
Received: 16 November 2016; Accepted: 19 April 2017;
Published: 09 May 2017.
Edited by:
Jan Adrianus Veenstra, University of Bordeaux 1, FranceReviewed by:
Nicolas Durand, University of Orléans, FranceZainulabeuddin Syed, University of Notre Dame, USA
Copyright © 2017 Li, Piermarini, Esquivel, Drumm, Schilkey and Hansen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Immo A. Hansen, immoh@nmsu.edu