RNA-seq validation: software for selection of reference and variable candidate genes for RT-qPCR.

BACKGROUND: Real-time quantitative PCR (RT-qPCR) is one of the most widely used gene expression analyses for validating RNA-seq data. This technique requires reference genes that are stable and highly expressed, at least across the different biological conditions present in the transcriptome. Reference and variable candidate gene selection is often neglected, leading to misinterpretation of the results. RESULTS: We developed a software named "Gene Selector for Validation" (GSV), which identifies the best reference and variable candidate genes for validation within a quantitative transcriptome. This tool also filters the candidate genes concerning the RT-qPCR assay detection limit. GSV was compared with other software using synthetic datasets and performed better, removing stable low-expression genes from the reference candidate list and creating the variable-expression validation list. GSV software was used on a real case, an Aedes aegypti transcriptome. The top GSV reference candidate genes were selected for RT-qPCR analysis, confirming that eiF1A and eiF3j were the most stable genes tested. The tool also confirmed that traditional mosquito reference genes were less stable in the analyzed samples, highlighting the possibility of inappropriate choices. A meta-transcriptome dataset with more than ninety thousand genes was also processed successfully. CONCLUSION: The GSV tool is a time and cost-effective tool that can be used to select reference and validation candidate genes from the biological conditions present in transcriptomic data.

Structure and expression of Rhodnius prolixus GH18 chitinases and chitinase-like proteins: Characterization of the physiological role of RpCht7, a gene from subgroup VIII, in vector fitness and reproduction.

Chitinases are enzymes responsible for the hydrolysis of glycosidic linkages within chitin chains. In insects, chitinases are typically members of the multigenic glycoside hydrolase family 18 (GH18). They participate in the relocation of chitin during development and molt, and in digestion in detritivores and predatory insects, and they control the peritrophic membrane thickness. Chitin metabolism is a promising target for developing vector control strategies, and knowledge of the roles of chitinases may reveal new targets and illuminate unique aspects of their physiology and interaction with microorganisms. Rhodnius prolixus is an important vector of Chagas disease, which is caused by the parasite Trypanosoma cruzi. In this study, we performed annotation and structural characterization of nine chitinase and chitinase-like protein genes in the R. prolixus genome. The roles of their corresponding transcripts were studied in more depth; their physiological roles were studied through RNAi silencing. Phylogenetic analysis of coding sequences showed that these genes belong to different subfamilies of GH18 chitinases already described in other insects. The expression patterns of these genes in different tissues and developmental stages were initially characterized using RT-PCR. RNAi screening showed silencing of the gene family members with very different efficiencies. Based on the knockdown results and the general lack of information about subgroup VIII of GH18, the RpCht7 gene was chosen for phenotype analysis. RpCht7 knockdown doubled the mortality in starving fifth-instar nymphs compared to dsGFP-injected controls. However, it did not alter blood intake, diuresis, digestion, molting rate, molting defects, sexual ratio, percentage of hatching, or average hatching time. Nevertheless, female oviposition was reduced by 53% in RpCht7-silenced insects, and differences in oviposition occurred within 14-20 days after a saturating blood meal. These results suggest that RpCht7 may be involved in the reproductive physiology and vector fitness of R. prolixus.

Evolution of Toll, Spatzle and MyD88 in insects: the problem of the Diptera bias.

BACKGROUND: Arthropoda, the most numerous and diverse metazoan phylum, has species in many habitats where they encounter various microorganisms and, as a result, mechanisms for pathogen recognition and elimination have evolved. The Toll pathway, involved in the innate immune system, was first described as part of the developmental pathway for dorsal-ventral differentiation in Drosophila. Its later discovery in vertebrates suggested that this system was extremely conserved. However, there is variation in presence/absence, copy number and sequence divergence in various genes along the pathway. As most studies have only focused on Diptera, for a comprehensive and accurate homology-based approach it is important to understand gene function in a number of different species and, in a group as diverse as insects, the use of species belonging to different taxonomic groups is essential. RESULTS: We evaluated the diversity of Toll pathway gene families in 39 Arthropod genomes, encompassing 13 different Insect Orders. Through computational methods, we shed some light into the evolution and functional annotation of protein families involved in the Toll pathway innate immune response. Our data indicates that: 1) intracellular proteins of the Toll pathway show mostly species-specific expansions; 2) the different Toll subfamilies seem to have distinct evolutionary backgrounds; 3) patterns of gene expansion observed in the Toll phylogenetic tree indicate that homology based methods of functional inference might not be accurate for some subfamilies; 4) Spatzle subfamilies are highly divergent and also pose a problem for homology based inference; 5) Spatzle subfamilies should not be analyzed together in the same phylogenetic framework; 6) network analyses seem to be a good first step in inferring functional groups in these cases. We specifically show that understanding Drosophila's Toll functions might not indicate the same function in other species. CONCLUSIONS: Our results show the importance of using species representing the different orders to better understand insect gene content, origin and evolution. More specifically, in intracellular Toll pathway gene families the presence of orthologues has important implications for homology based functional inference. Also, the different evolutionary backgrounds of Toll gene subfamilies should be taken into consideration when functional studies are performed, especially for TOLL9, TOLL, TOLL2_7, and the new TOLL10 clade. The presence of Diptera specific clades or the ones lacking Diptera species show the importance of overcoming the Diptera bias when performing functional characterization of Toll pathways.

Comprehensive Quantitative Proteome Analysis of Aedes aegypti Identifies Proteins and Pathways Involved in Wolbachia pipientis and Zika Virus Interference Phenomenon.

Zika virus (ZIKV) is a global public health emergency due to its association with microcephaly, Guillain-Barré syndrome, neuropathy, and myelitis in children and adults. A total of 87 countries have had evidence of autochthonous mosquito-borne transmission of ZIKV, distributed across four continents, and no antivirus therapy or vaccines are available. Therefore, several strategies have been developed to target the main mosquito vector, Aedes aegypti, to reduce the burden of different arboviruses. Among such strategies, the use of the maternally-inherited endosymbiont Wolbachia pipientis has been applied successfully to reduce virus susceptibility and decrease transmission. However, the mechanisms by which Wolbachia orchestrate resistance to ZIKV infection remain to be elucidated. In this study, we apply isobaric labeling quantitative mass spectrometry (MS)-based proteomics to quantify proteins and identify pathways altered during ZIKV infection; Wolbachia infection; co-infection with Wolbachia/ZIKV in the A. aegypti heads and salivary glands. We show that Wolbachia regulates proteins involved in reactive oxygen species production, regulates humoral immune response, and antioxidant production. The reduction of ZIKV polyprotein in the presence of Wolbachia in mosquitoes was determined by MS and corroborates the idea that Wolbachia helps to block ZIKV infections in A. aegypti. The present study offers a rich resource of data that may help to elucidate mechanisms by which Wolbachia orchestrate resistance to ZIKV infection in A. aegypti, and represents a step further on the development of new targeted methods to detect and quantify ZIKV and Wolbachia directly in complex tissues.

A physiologic overview of the organ-specific transcriptome of the cattle tick Rhipicephalus microplus.

To further obtain insights into the Rhipicephalus microplus transcriptome, we used RNA-seq to carry out a study of expression in (i) embryos; (ii) ovaries from partially and fully engorged females; (iii) salivary glands from partially engorged females; (iv) fat body from partially and fully engorged females; and (v) digestive cells from partially, and (vi) fully engorged females. We obtained > 500 million Illumina reads which were assembled de novo, producing > 190,000 contigs, identifying 18,857 coding sequences (CDS). Reads from each library were mapped back into the assembled transcriptome giving a view of gene expression in different tissues. Transcriptomic expression and pathway analysis showed that several genes related in blood digestion and host-parasite interaction were overexpressed in digestive cells compared with other tissues. Furthermore, essential genes for the cell development and embryogenesis were overexpressed in ovaries. Taken altogether, these data offer novel insights into the physiology of production and role of saliva, blood digestion, energy metabolism, and development with submission of 10,932 novel tissue/cell specific CDS to the NCBI database for this important tick species.

Polyphenol-Rich Diets Exacerbate AMPK-Mediated Autophagy, Decreasing Proliferation of Mosquito Midgut Microbiota, and Extending Vector Lifespan.

BACKGROUND: Mosquitoes feed on plant-derived fluids such as nectar and sap and are exposed to bioactive molecules found in this dietary source. However, the role of such molecules on mosquito vectorial capacity is unknown. Weather has been recognized as a major determinant of the spread of dengue, and plants under abiotic stress increase their production of polyphenols. RESULTS: Here, we show that including polyphenols in mosquito meals promoted the activation of AMP-dependent protein kinase (AMPK). AMPK positively regulated midgut autophagy leading to a decrease in bacterial proliferation and an increase in vector lifespan. Suppression of AMPK activity resulted in a 6-fold increase in midgut microbiota. Similarly, inhibition of polyphenol-induced autophagy induced an 8-fold increase in bacterial proliferation. Mosquitoes maintained on the polyphenol diet were readily infected by dengue virus. CONCLUSION: The present findings uncover a new direct route by which exacerbation of autophagy through activation of the AMPK pathway leads to a more efficient control of mosquito midgut microbiota and increases the average mosquito lifespan. Our results suggest for the first time that the polyphenol content and availability of the surrounding vegetation may increase the population of mosquitoes prone to infection with arboviruses.

Protein 2DE reference map of the anterior midgut of the blood-sucking bug Rhodnius prolixus.

Rhodnius prolixus is an important vector of Trypanosoma cruzi, the causative agent of Chagas' disease, an illness that affects 20% of Latin America population. The obligatory course of the parasite in the vector digestive tract has made it an important target for investigation in order to control the parasite transmission and thus interrupt its biological cycle in the insect vector. Therefore, an insight into the vector midgut physiology is valuable for insect control as well as to provide potential novel targets for drugs and vaccines development and thus disease treatment. In this study, the first 2DE map of R. prolixus anterior midgut is described. Proteins were separated by 2DE and analyzed by nano-LC MS/MS. The results yielded 489 proteins from 475 spots. These proteins were classified into 28 functional groups and their physiological roles in the insect midgut are discussed. All MS data have been deposited in the ProteomeXchange with identifiers PXD001488 and PXD001489 (http://proteomecentral.proteomexchange.org/dataset/PXD001488, http://proteomecentral.proteomexchange.org/dataset/PXD001489).

Molecular analysis of Aedes aegypti classical protein tyrosine phosphatases uncovers an ortholog of mammalian PTP-1B implicated in the control of egg production in mosquitoes.

BACKGROUND: Protein Tyrosine Phosphatases (PTPs) are enzymes that catalyze phosphotyrosine dephosphorylation and modulate cell differentiation, growth and metabolism. In mammals, PTPs play a key role in the modulation of canonical pathways involved in metabolism and immunity. PTP1B is the prototype member of classical PTPs and a major target for treating human diseases, such as cancer, obesity and diabetes. These signaling enzymes are, hence, targets of a wide array of inhibitors. Anautogenous mosquitoes rely on blood meals to lay eggs and are vectors of the most prevalent human diseases. Identifying the mosquito ortholog of PTP1B and determining its involvement in egg production is, therefore, important in the search for a novel and crucial target for vector control. METHODOLOGY/PRINCIPAL FINDINGS: We conducted an analysis to identify the ortholog of mammalian PTP1B in the Aedes aegypti genome. We identified eight genes coding for classical PTPs. In silico structural and functional analyses of proteins coded by such genes revealed that four of these code for catalytically active enzymes. Among the four genes coding for active PTPs, AAEL001919 exhibits the greatest degree of homology with the mammalian PTP1B. Next, we evaluated the role of this enzyme in egg formation. Blood feeding largely affects AAEL001919 expression, especially in the fat body and ovaries. These tissues are critically involved in the synthesis and storage of vitellogenin, the major yolk protein. Including the classical PTP inhibitor sodium orthovanadate or the PTP substrate DiFMUP in the blood meal decreased vitellogenin synthesis and egg production. Similarly, silencing AAEL001919 using RNA interference (RNAi) assays resulted in 30% suppression of egg production. CONCLUSIONS/SIGNIFICANCE: The data reported herein implicate, for the first time, a gene that codes for a classical PTP in mosquito egg formation. These findings raise the possibility that this class of enzymes may be used as novel targets to block egg formation in mosquitoes.

Glycoinositolphospholipids from Trypanosomatids subvert nitric oxide production in Rhodnius prolixus salivary glands.

BACKGROUND: Rhodnius prolixus is a blood-sucking bug vector of Trypanosoma cruzi and T. rangeli. T. cruzi is transmitted by vector feces deposited close to the wound produced by insect mouthparts, whereas T. rangeli invades salivary glands and is inoculated into the host skin. Bug saliva contains a set of nitric oxide-binding proteins, called nitrophorins, which deliver NO to host vessels and ensure vasodilation and blood feeding. NO is generated by nitric oxide synthases (NOS) present in the epithelium of bug salivary glands. Thus, T. rangeli is in close contact with NO while in the salivary glands. METHODOLOGY/PRINCIPAL FINDINGS: Here we show by immunohistochemical, biochemical and molecular techniques that inositolphosphate-containing glycolipids from trypanosomatids downregulate NO synthesis in the salivary glands of R. prolixus. Injecting insects with T. rangeli-derived glycoinositolphospholipids (Tr GIPL) or T. cruzi-derived glycoinositolphospholipids (Tc GIPL) specifically decreased NO production. Salivary gland treatment with Tc GIPL blocks NO production without greatly affecting NOS mRNA levels. NOS protein is virtually absent from either Tr GIPL- or Tc GIPL-treated salivary glands. Evaluation of NO synthesis by using a fluorescent NO probe showed that T. rangeli-infected or Tc GIPL-treated glands do not show extensive labeling. The same effect is readily obtained by treatment of salivary glands with the classical protein tyrosine phosphatase (PTP) inhibitor, sodium orthovanadate (SO). This suggests that parasite GIPLs induce the inhibition of a salivary gland PTP. GIPLs specifically suppressed NO production and did not affect other anti-hemostatic properties of saliva, such as the anti-clotting and anti-platelet activities. CONCLUSIONS/SIGNIFICANCE: Taken together, these data suggest that trypanosomatids have overcome NO generation using their surface GIPLs. Therefore, these molecules ensure parasite survival and may ultimately enhance parasite transmission.

The dengue vector Aedes aegypti contains a functional high mobility group box 1 (HMGB1) protein with a unique regulatory C-terminus.

The mosquito Aedes aegypti can spread the dengue, chikungunya and yellow fever viruses. Thus, the search for key molecules involved in the mosquito survival represents today a promising vector control strategy. High Mobility Group Box (HMGB) proteins are essential nuclear factors that maintain the high-order structure of chromatin, keeping eukaryotic cells viable. Outside the nucleus, secreted HMGB proteins could alert the innate immune system to foreign antigens and trigger the initiation of host defenses. In this work, we cloned and functionally characterized the HMGB1 protein from Aedes aegypti (AaHMGB1). The AaHMGB1 protein typically consists of two HMG-box DNA binding domains and an acidic C-terminus. Interestingly, AaHMGB1 contains a unique alanine/glutamine-rich (AQ-rich) C-terminal region that seems to be exclusive of dipteran HMGB proteins. AaHMGB1 is localized to the cell nucleus, mainly associated with heterochromatin. Circular dichroism analyses of AaHMGB1 or the C-terminal truncated proteins revealed α-helical structures. We showed that AaHMGB1 can effectively bind and change the topology of DNA, and that the AQ-rich and the C-terminal acidic regions can modulate its ability to promote DNA supercoiling, as well as its preference to bind supercoiled DNA. AaHMGB1 is phosphorylated by PKA and PKC, but not by CK2. Importantly, phosphorylation of AaHMGB1 by PKA or PKC completely abolishes its DNA bending activity. Thus, our study shows that a functional HMGB1 protein occurs in Aedes aegypt and we provide the first description of a HMGB1 protein containing an AQ-rich regulatory C-terminus.

Functional properties of Schistosoma mansoni single-stranded DNA-binding protein SmPUR-alpha. Description of the interaction between SmPUR-alpha and SMYB1.

PUR-alpha is a highly conserved protein in eukaryotes belonging to the family of single-stranded DNA-binding proteins. Because PUR-alpha is a multifunctional protein that participates in several regulatory events at the level of gene transcription, it became relevant to investigate the structural features of Schistosoma mansoni PUR-alpha (SmPUR-alpha) that could be correlated to its mode of action. Using deletion constructs on a dot blot assay we mapped the domains of GST-SmPUR-alpha fusion protein involved in the interactions with DNA and RNA. Individually, the N-terminal amino acid residues 1-26 and the C-terminal residues 196-276 of GST-SmPUR-alpha which did not contain nucleic acid-binding domains, did not bind ssDNA or RNA. In contrast, domains encompassing the N-terminal and Class I and C-terminal plus Class I exhibited the highest binding affinity. Seemingly, the latter (GST-SmPUR-alpha 174-276) played a major role in nucleic acid interaction as judged by affinity alone. Other combinations of the deletion constructs displayed either intermediary or no binding affinity to the DNA or RNA probes. Gel shift competition assay showed that GST-SmPUR-alpha bound to ssDNA with higher affinity than to RNA. Because SmPUR-alpha contains two putative phosphorylation sites the protein was tested as a substrate to casein kinase II. GST-SmPUR-alpha could be phosphorylated in vitro by casein kinase II at both, the N- and C-terminus of the protein. The multifunctional nature of SmPUR-alpha was demonstrated by experiments measuring the physical interaction between SmPUR-alpha and the transcription factor SMYB1. This was determined in vivo (yeast two hybrid) and in vitro (GST-pull down). Furthermore, we showed that SmPUR-alpha and SMYB1 acted synergistically to bind preferentially to pyrimidine-rich sequences.