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.

In-depth transcriptomic analysis of Anopheles gambiae hemocytes uncovers novel genes and the oenocytoid developmental lineage.

BACKGROUND: Hemocytes are immune cells that patrol the mosquito hemocoel and mediate critical cellular defense responses against pathogens. However, despite their importance, a comprehensive transcriptome of these cells was lacking because they constitute a very small fraction of the total cells in the insect, limiting the study of hemocyte differentiation and immune function. RESULTS: In this study, an in-depth hemocyte transcriptome was built by extensive bulk RNA sequencing and assembly of hemocyte RNAs from adult A. gambiae female mosquitoes, based on approximately 2.4 billion short Illumina and about 9.4 million long PacBio high-quality reads that mapped to the A. gambiae PEST genome (P4.14 version). A total of 34,939 transcripts were annotated including 4,020 transcripts from novel genes and 20,008 novel isoforms that result from extensive differential splicing of transcripts from previously annotated genes. Most hemocyte transcripts identified (89.8%) are protein-coding while 10.2% are non-coding RNAs. The number of transcripts identified in the novel hemocyte transcriptome is twice the number in the current annotation of the A. gambiae genome (P4.14 version). Furthermore, we were able to refine the analysis of a previously published single-cell transcriptome (scRNAseq) data set by using the novel hemocyte transcriptome as a reference to re-define the hemocyte clusters and determine the path of hemocyte differentiation. Unsupervised pseudo-temporal ordering using the Tools for Single Cell Analysis software uncovered a novel putative prohemocyte precursor cell type that gives rise to prohemocytes. Pseudo-temporal ordering with the Monocle 3 software, which analyses changes in gene expression during dynamic biological processes, determined that oenocytoids derive from prohemocytes, a cell population that also gives rise to the granulocyte lineage. CONCLUSION: A high number of mRNA splice variants are expressed in hemocytes, and they may account for the plasticity required to mount efficient responses to many different pathogens. This study highlights the importance of a comprehensive set of reference transcripts to perform robust single-cell transcriptomic data analysis of cells present in low abundance. The detailed annotation of the hemocyte transcriptome will uncover new facets of hemocyte development and function in adult dipterans and is a valuable community resource for future studies on mosquito cellular immunity.

Redox imbalance induces remodeling of glucose metabolism in Rhipicephalus microplus embryonic cell line.

Carbohydrate metabolism not only functions in supplying cellular energy but also has an important role in maintaining physiological homeostasis and in preventing oxidative damage caused by reactive oxygen species. Previously, we showed that arthropod embryonic cell lines have high tolerance to H2O2 exposure. Here, we describe that Rhipicephalus microplus tick embryonic cell line (BME26) employs an adaptive glucose metabolism mechanism that confers tolerance to hydrogen peroxide at concentrations too high for other organisms. This adaptive mechanism sustained by glucose metabolism remodeling promotes cell survival and redox balance in BME26 cell line after millimolar H2O2 exposure. The present work shows that this tick cell line could tolerate high H2O2 concentrations by initiating a carbohydrate-related adaptive response. We demonstrate that gluconeogenesis was induced as a compensation strategy that involved, among other molecules, the metabolic enzymes NADP-ICDH, G6PDH, and PEPCK. We also found that this phenomenon was coupled to glycogen accumulation and glucose uptake, supporting the pentose phosphate pathway to sustain NADPH production and leading to cell survival and proliferation. Our findings suggest that the described response is not atypical, being also observed in cancer cells, which highlights the importance of this model to all proliferative cells. We propose that these results will be useful in generating basic biological information to support the development of new strategies for disease treatment and parasite control.

Aedes fluviatilis cell lines as new tools to study metabolic and immune interactions in mosquito-Wolbachia symbiosis.

In the present work, we established two novel embryonic cell lines from the mosquito Aedes fluviatilis containing or not the naturally occurring symbiont bacteria Wolbachia, which were called wAflu1 and Aflu2, respectively. We also obtained wAflu1 without Wolbachia after tetracycline treatment, named wAflu1.tet. Morphofunctional characterization was performed to help elucidate the symbiont-host interaction in the context of energy metabolism regulation and molecular mechanisms of the immune responses involved. The presence of Wolbachia pipientis improves energy performance in A. fluviatilis cells; it affects the regulation of key energy sources such as lipids, proteins, and carbohydrates, making the distribution of actin more peripheral and with extensions that come into contact with neighboring cells. Additionally, innate immunity mechanisms were activated, showing that the wAflu1 and wAflu1.tet cells are responsive after the stimulus using Gram negative bacteria. Therefore, this work confirms the natural, mutually co-regulating symbiotic relationship between W. pipientis and A. fluviatilis, modulating the host metabolism and immune pathway activation. The results presented here add important resources to the current knowledge of Wolbachia-arthropod interactions.

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.