Here we have used a systems biology approach to study innate and adaptive responses to vaccination against influenza in humans during three consecutive influenza seasons. We studied healthy adults vaccinated with trivalent inactivated influenza vaccine (TIV) or live attenuated influenza vaccine (LAIV). TIV induced higher antibody titers and more plasmablasts than LAIV did. In subjects vaccinated with TIV, early molecular signatures correlated with and could be used to accurately predict later antibody titers in two independent trials. Notably, expression of the kinase CaMKIV at day 3 was inversely correlated with later antibody titers. Vaccination of CaMKIV-deficient mice with TIV induced enhanced antigen-specific antibody titers, which demonstrated an unappreciated role for CaMKIV in the regulation of antibody responses. Thus, systems approaches can be used to predict immunogenicity and provide new mechanistic insights about vaccines.

miR-137 plays critical roles in the nervous system and tumor development; an increase in its expression is required for neuronal differentiation while its reduction is implicated in gliomagenesis. To evaluate the potential of miR-137 in glioblastoma therapy, we conducted genome-wide target mapping in glioblastoma cells by measuring the level of association between PABP and mRNAs in cells transfected with miR-137 mimics vs. controls via RIPSeq. Impact on mRNA levels was also measured by RNASeq. By combining the results of both experimental approaches, 1468 genes were found to be negatively impacted by miR-137--among them, 595 (40%) contain miR-137 predicted sites. The most relevant targets include oncogenic proteins and key players in neurogenesis like c-KIT, YBX1, AKT2, CDC42, CDK6 and TGFβ2. Interestingly, we observed that several identified miR-137 targets are also predicted to be regulated by miR-124, miR-128 and miR-7, which are equally implicated in neuronal differentiation and gliomagenesis. We suggest that the concomitant increase of these four miRNAs in neuronal stem cells or their repression in tumor cells could produce a robust regulatory effect with major consequences to neuronal differentiation and tumorigenesis.

Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC.

Composting operations are a rich source for prospection of biomass degradation enzymes. We have analyzed the microbiomes of two composting samples collected in a facility inside the São Paulo Zoo Park, in Brazil. All organic waste produced in the park is processed in this facility, at a rate of four tons/day. Total DNA was extracted and sequenced with Roche/454 technology, generating about 3 million reads per sample. To our knowledge this work is the first report of a composting whole-microbial community using high-throughput sequencing and analysis. The phylogenetic profiles of the two microbiomes analyzed are quite different, with a clear dominance of members of the Lactobacillus genus in one of them. We found a general agreement of the distribution of functional categories in the Zoo compost metagenomes compared with seven selected public metagenomes of biomass deconstruction environments, indicating the potential for different bacterial communities to provide alternative mechanisms for the same functional purposes. Our results indicate that biomass degradation in this composting process, including deconstruction of recalcitrant lignocellulose, is fully performed by bacterial enzymes, most likely by members of the Clostridiales and Actinomycetales orders.

Systems approaches have been used to describe molecular signatures driving immunity to influenza vaccination in humans. Whether such signatures are similar across multiple seasons and in diverse populations is unknown. We applied systems approaches to study immune responses in young, elderly, and diabetic subjects vaccinated with the seasonal influenza vaccine across five consecutive seasons. Signatures of innate immunity and plasmablasts correlated with and predicted influenza antibody titers at 1 month after vaccination with >80% accuracy across multiple seasons but were not associated with the longevity of the response. Baseline signatures of lymphocyte and monocyte inflammation were positively and negatively correlated, respectively, with antibody responses at 1 month. Finally, integrative analysis of microRNAs and transcriptomic profiling revealed potential regulators of vaccine immunity. These results identify shared vaccine-induced signatures across multiple seasons and in diverse populations and might help guide the development of next-generation vaccines that provide persistent immunity against influenza.

Chronic viral infections are characterized by a state of CD8+ T-cell dysfunction that is associated with expression of the programmed cell death 1 (PD-1) inhibitory receptor. A better understanding of the mechanisms that regulate CD8+ T-cell responses during chronic infection is required to improve immunotherapies that restore function in exhausted CD8+ T cells. Here we identify a population of virus-specific CD8+ T cells that proliferate after blockade of the PD-1 inhibitory pathway in mice chronically infected with lymphocytic choriomeningitis virus (LCMV). These LCMV-specific CD8+ T cells expressed the PD-1 inhibitory receptor, but also expressed several costimulatory molecules such as ICOS and CD28. This CD8+ T-cell subset was characterized by a unique gene signature that was related to that of CD4+ T follicular helper (TFH) cells, CD8+ T cell memory precursors and haematopoietic stem cell progenitors, but that was distinct from that of CD4+ TH1 cells and CD8+ terminal effectors. This CD8+ T-cell population was found only in lymphoid tissues and resided predominantly in the T-cell zones along with naive CD8+ T cells. These PD-1+CD8+ T cells resembled stem cells during chronic LCMV infection, undergoing self-renewal and also differentiating into the terminally exhausted CD8+ T cells that were present in both lymphoid and non-lymphoid tissues. The proliferative burst after PD-1 blockade came almost exclusively from this CD8+ T-cell subset. Notably, the transcription factor TCF1 had a cell-intrinsic and essential role in the generation of this CD8+ T-cell subset. These findings provide a better understanding of T-cell exhaustion and have implications in the optimization of PD-1-directed immunotherapy in chronic infections and cancer.

Schistosoma mansoni is the major causative agent of schistosomiasis. The parasite takes advantage of host signals to complete its development in the human body. Tumor necrosis factor-alpha (TNF-alpha) is a human cytokine involved in skin inflammatory responses, and although its effect on the adult parasite's metabolism and egg-laying process has been previously described, a comprehensive assessment of the TNF-alpha pathway and its downstream molecular effects is lacking.

Transcription of large numbers of non-coding RNAs originating from intronic regions of human genes has been recently reported, but mechanisms governing their biosynthesis and biological functions are largely unknown. In this work, we evaluated the existence of a common mechanism of transcription regulation shared by protein-coding mRNAs and intronic RNAs by measuring the effect of androgen on the transcriptional profile of a prostate cancer cell line.

Congenital Zika syndrome (CZS) causes early brain development impairment by affecting neural progenitor cells (NPCs). Here, we analyze NPCs from three pairs of dizygotic twins discordant for CZS. We compare by RNA-Seq the NPCs derived from CZS-affected and CZS-unaffected twins. Prior to Zika virus (ZIKV) infection the NPCs from CZS babies show a significantly different gene expression signature of mTOR and Wnt pathway regulators, key to a neurodevelopmental program. Following ZIKV in vitro infection, cells from affected individuals have significantly higher ZIKV replication and reduced cell growth. Whole-exome analysis in 18 affected CZS babies as compared to 5 unaffected twins and 609 controls excludes a monogenic model to explain resistance or increased susceptibility to CZS development. Overall, our results indicate that CZS is not a stochastic event and depends on NPC intrinsic susceptibility, possibly related to oligogenic and/or epigenetic mechanisms.

The mammalian genome is transcribed into large numbers of long noncoding RNAs (lncRNAs), but the definition of functional lncRNA groups has proven difficult, partly due to their low sequence conservation and lack of identified shared properties. Here we consider promoter conservation and positional conservation as indicators of functional commonality.

The analysis of modular gene co-expression networks is a well-established method commonly used for discovering the systems-level functionality of genes. In addition, these studies provide a basis for the discovery of clinically relevant molecular pathways underlying different diseases and conditions.

Understanding the mechanisms of vaccine-elicited protection contributes to the development of new vaccines. The emerging field of systems vaccinology provides detailed information on host responses to vaccination and has been successfully applied to study the molecular mechanisms of several vaccines. Long noncoding RNAs (lncRNAs) are crucially involved in multiple biological processes, but their role in vaccine-induced immunity has not been explored. We performed an analysis of over 2,000 blood transcriptome samples from 17 vaccine cohorts to identify lncRNAs potentially involved with antibody responses to influenza and yellow fever vaccines. We have created an online database where all results from this analysis can be accessed easily. We found that lncRNAs participate in distinct immunological pathways related to vaccine-elicited responses. Among them, we showed that the expression of lncRNA FAM30A was high in B cells and correlates with the expression of immunoglobulin genes located in its genomic vicinity. We also identified altered expression of these lncRNAs in RNA-sequencing (RNA-seq) data from a cohort of children following immunization with intranasal live attenuated influenza vaccine, suggesting a common role across several diverse vaccines. Taken together, these findings provide evidence that lncRNAs have a significant impact on immune responses induced by vaccination.

SerpinB2 (plasminogen activator inhibitor type 2) has been called the "undecided serpin" with no clear consensus on its physiological role, although it is well described as an inhibitor of urokinase plasminogen activator (uPA). In macrophages, pro-inflammatory stimuli usually induce SerpinB2; however, expression is constitutive in Gata6+ large peritoneal macrophages (LPM). Interrogation of expression data from human macrophages treated with a range of stimuli using a new bioinformatics tool, CEMiTool, suggested that SerpinB2 is most tightly co- and counter-regulated with genes associated with cell movement. Using LPM from SerpinB2-/- and SerpinB2R380A (active site mutant) mice, we show that migration on Matrigel was faster than for their wild-type controls. Confocal microscopy illustrated that SerpinB2 and F-actin staining overlapped in focal adhesions and lamellipodia. Genes associated with migration and extracellular matrix interactions were also identified by RNA-Seq analysis of migrating RPM from wild-type and SerpinB2R380A mice. Subsequent gene set enrichment analyses (GSEA) suggested SerpinB2 counter-regulates many Gata6-regulated genes associated with migration. These data argue that the role of SerpinB2 in macrophages is inhibition of uPA-mediated plasmin generation during cell migration. GSEA also suggested that SerpinB2 expression (likely via ensuing modulation of uPA-receptor/integrin signaling) promotes the adoption of a resolution phase signature.

Next Generation Sequencing (NGS) strategies, like RNA-Seq, have revealed the transcription of a wide variety of long non-coding RNAs (lncRNAs) in the genomes of several organisms. In the present work we assessed the lncRNAs complement of Schistosoma mansoni, the blood fluke that causes schistosomiasis, ranked among the most prevalent parasitic diseases worldwide. We focused on the long intergenic/intervening ncRNAs (lincRNAs), hidden within the large amount of information obtained through RNA-Seq in S. mansoni (88 libraries). Our computational pipeline identified 7029 canonically-spliced putative lincRNA genes on 2596 genomic loci (at an average 2.7 isoforms per lincRNA locus), as well as 402 spliced lncRNAs that are antisense to protein-coding (PC) genes. Hundreds of lincRNAs showed traits for being functional, such as the presence of epigenetic marks at their transcription start sites, evolutionary conservation among other schistosome species and differential expression across five different life-cycle stages of the parasite. Real-time qPCR has confirmed the differential life-cycle stage expression of a set of selected lincRNAs. We have built PC gene and lincRNA co-expression networks, unraveling key biological processes where lincRNAs might be involved during parasite development. This is the first report of a large-scale identification and structural annotation of lncRNAs in the S. mansoni genome.

Schistosoma mansoni is a flatworm that causes schistosomiasis, a neglected tropical disease that affects more than 200 million people worldwide. There is only one drug indicated for treatment, praziquantel, which may lead to parasite resistance emergence. The ribonucleoside analogue 5-azacytidine (5-AzaC) is an epigenetic drug that inhibits S. mansoni oviposition and ovarian development through interference with parasite transcription, translation and stem cell activities. Therefore, studying the downstream pathways affected by 5-AzaC in S. mansoni may contribute to the discovery of new drug targets. Long non-coding RNAs (lncRNAs) are transcripts longer than 200 nucleotides with low or no protein coding potential that have been involved in reproduction, stem cell maintenance and drug resistance. We have recently published a catalog of lncRNAs expressed in S. mansoni life-cycle stages, tissues and single cells. However, it remains largely unknown if lncRNAs are responsive to epigenetic drugs in parasites. Here, we show by RNA-Seq re-analyses that hundreds of lncRNAs are differentially expressed after in vitro 5-AzaC treatment of S. mansoni females, including intergenic, antisense and sense lncRNAs. Many of these lncRNAs belong to co-expression network modules related to male metabolism and are also differentially expressed in unpaired compared with paired females and ovaries. Half of these lncRNAs possess histone marks at their genomic loci, indicating regulation by histone modification. Among a selected set of 8 lncRNAs, half of them were validated by RT-qPCR as differentially expressed in females, and some of them also in males. Interestingly, these lncRNAs are also expressed in other life-cycle stages. This study demonstrates that many lncRNAs potentially involved with S. mansoni reproductive biology are modulated by 5-AzaC and sheds light on the relevance of exploring lncRNAs in response to drug treatments in parasites.

Patients infected with the Dengue virus (DENV) often present with a massive generation of DENV-specific antibody-secreting cells (ASCs) in the blood. In some cases, these ASCs represent more than 50% of the circulating B cells, a higher magnitude than those induced by other infections, vaccinations, and plasma cell lymphomas. However, it remains unclear how the DENV infection elicits this colossal response. To address this issue, we utilised an in vitro strategy to induce human PBMCs of healthy individuals incubated with DENV particles (DENV4 TVP/360) to differentiate into ASCs. As controls, PBMCs were incubated with a mitogen cocktail or supernatants of uninfected C6/36 cells (mock). The ASC phenotype and function were increasingly detected in the DENV and mitogen-cultured PBMCs as compared to mock-treated cells. In contrast to the in vivo condition, secreted IgG derived from the PBMC-DENV culture was not DENV-specific. Lower ASC numbers were observed when inactivated viral particles or purified B cells were added to the cultures. The physical contact was essential between B cells and the remaining PBMCs for the DENV-mediated ASC response. Considering the evidence for the activation of the tryptophan metabolism detected in the serum of Dengue patients, we assessed its relevance in the DENV-mediated ASC differentiation. For this, tryptophan and its respective metabolites were quantified in the supernatants of cell cultures through mass spectrophotometry. Tryptophan depletion and kynurenine accumulation were found in the supernatants of PBMC-DENV cultures, which presented enhanced detection of indoleamine 2,3-dioxygenase 1 and 2 transcripts as compared to controls. In PBMC-DENV cultures, tryptophan and kynurenine levels strongly correlated to the respective ASC numbers, while the kynurenine levels were directly proportional to the secreted IgG titers. Contrastingly, PBMCs incubated with Zika or attenuated Yellow Fever viruses showed no correlation between their kynurenine concentrations and ASC numbers. Therefore, our data revealed the existence of distinct pathways for the DENV-mediated ASC differentiation and suggest the involvement of the tryptophan metabolism in this cellular process triggered by flavivirus infections.

Cell signaling events triggered by androgen hormone in prostate cells is dependent on activation of the androgen receptor (AR) transcription factor. Androgen hormone binding to AR promotes its displacement from the cytoplasm to the nucleus and AR binding to DNA motifs, thus inducing activatory and inhibitory transcriptional programs through a complex regulatory mechanism not yet fully understood. In this work, we performed RNA-seq deep-sequencing of LNCaP prostate cancer cells and found over 7000 expressed long intergenic non-coding RNAs (lincRNAs), of which ∼4000 are novel lincRNAs, and 258 lincRNAs have their expression activated by androgen. Immunoprecipitation of AR, followed by large-scale sequencing of co-immunoprecipitated RNAs (RIP-Seq) has identified in the LNCaP cell line a total of 619 lincRNAs that were significantly enriched (FDR < 10%, DESeq2) in the anti-Androgen Receptor (antiAR) fraction in relation to the control fraction (non-specific IgG), and we named them Androgen-Receptor-Associated lincRNAs (ARA-lincRNAs). A genome-wide analysis showed that protein-coding gene neighbors to ARA-lincRNAs had a significantly higher androgen-induced change in expression than protein-coding genes neighboring lincRNAs not associated to AR. To find relevant epigenetic signatures enriched at the ARA-lincRNAs' transcription start sites (TSSs) we used a machine learning approach and identified that the ARA-lincRNA genomic loci in LNCaP cells are significantly enriched with epigenetic marks that are characteristic of in cis enhancer RNA regulators, and that the H3K27ac mark of active enhancers is conspicuously enriched at the TSS of ARA-lincRNAs adjacent to androgen-activated protein-coding genes. In addition, LNCaP topologically associating domains (TADs) that comprise chromatin regions with ARA-lincRNAs exhibit transcription factor contents, epigenetic marks and gene transcriptional activities that are significantly different from TADs not containing ARA-lincRNAs. This work highlights the possible involvement of hundreds of lincRNAs working in synergy with the AR on the genome-wide androgen-induced gene regulatory program in prostate cells.

Acute respiratory infection (ARI) is the most frequent cause for hospitalization in infants and young children. Using multiplexed nCounter technology to digitally quantify 600 human mRNAs in parallel with 14 virus- and 5 bacterium-specific RNAs, we characterized viral and bacterial presence in nasopharyngeal aspirates (NPA) of 58 children with ARI and determined the corresponding in situ immune profiles. NPA contained different groups of organisms and these were classified into bacterial (n = 27), viral (n = 5), codetection [containing both viral and bacterial transcripts (n = 21), or indeterminate intermediate where microbial load is below threshold (n = 5)]. We then identified differentially expressed immune transcripts (DEITs) comparing NPAs from symptomatic children vs. healthy controls, and comparing children presenting NPAs with detectable microbial load vs. indeterminate. We observed a strong innate immune response in NPAs, due to the presence of evolutionarily conserved type I Interferon (IFN)-stimulated genes (ISG), which was correlated with total bacterial and/or viral load. In comparison with indeterminate NPAs, adaptive immunity transcripts discriminated among viral, bacterial, and codetected microbial profiles. In viral NPAs, B cell transcripts were significantly enriched among DEITs, while only type III IFN was correlated with viral load. In bacterial NPAs, myeloid cells and coinhibitory transcripts were enriched and significantly correlated with bacterial load. In conclusion, digital nCounter transcriptomics provide a microbial and immunological in situ "snapshot" of the nasopharyngeal interface in children with ARI. This enabled discrimination among viral, bacterial, codetection, and indeterminate transcripts in the samples using non-invasive sampling.

Zika virus (ZIKV) causes congenital Zika syndrome (CZS), which is characterized by fetal demise, microcephaly and other abnormalities. ZIKV in the pregnant woman circulation must cross the placental barrier that includes fetal endothelial cells and trophoblasts, in order to reach the fetus. CZS occurs in ~1-40% of cases of pregnant women infected by ZIKV, suggesting that mothers' infection by ZIKV during pregnancy is not deterministic for CZS phenotype in the fetus. Therefore, other susceptibility factors might be involved, including the host genetic background. We have previously shown that in three pairs of dizygotic twins discordant for CZS, neural progenitor cells (NPCs) from the CZS-affected twins presented differential in vitro ZIKV susceptibility compared with NPCs from the non-affected. Here, we analyzed human-induced-pluripotent-stem-cell-derived (hiPSC-derived) trophoblasts from these twins and compared by RNA-Seq the trophoblasts from CZS-affected and non-affected twins. Following in vitro exposure to a Brazilian ZIKV strain (ZIKVBR), trophoblasts from CZS-affected twins were significantly more susceptible to ZIKVBR infection when compared with trophoblasts from the non-affected. Transcriptome profiling revealed no differences in gene expression levels of ZIKV candidate attachment factors, IFN receptors and IFN in the trophoblasts, either before or after ZIKVBR infection. Most importantly, ZIKVBR infection caused, only in the trophoblasts from CZS-affected twins, the downregulation of genes related to extracellular matrix organization and to leukocyte activation, which are important for trophoblast adhesion and immune response activation. In addition, only trophoblasts from non-affected twins secreted significantly increased amounts of chemokines RANTES/CCL5 and IP10 after infection with ZIKVBR. Overall, our results showed that trophoblasts from non-affected twins have the ability to more efficiently activate genes that are known to play important roles in cell adhesion and in triggering the immune response to ZIKV infection in the placenta, and this may contribute to predict protection from ZIKV dissemination into fetuses' tissues.

Trypanosoma cruzi-the causative agent of Chagas disease-like other kinetoplastids, relies mostly on post-transcriptional mechanisms for regulation of gene expression. However, trypanosomatids undergo drastic changes in nuclear architecture and chromatin structure along their complex life cycle which, combined with a remarkable set of reversible histone post-translational modifications, indicate that chromatin is also a target for control of gene expression and differentiation signals in these organisms. Chromatin-modifying enzymes have a direct impact on gene expression programs and DNA metabolism. In this work, we have investigated the function of T. cruzi histone deacetylase 4 (TcHDAC4). We show that, although TcHDAC4 is not essential for viability, metacyclic trypomastigote TcHDAC4 null mutants show a thin cell body and a round and less condensed nucleus located very close to the kinetoplast. Sixty-four acetylation sites were quantitatively evaluated, which revealed H2AT85ac, H4K10ac and H4K78ac as potential target sites of TcHDAC4. Gene expression analyses identified three chromosomes with overrepresented regions of differentially expressed genes in the TcHDAC4 knockout mutant compared with the wild type, showing clusters of either up or downregulated genes. The adjacent chromosomal location of some of these genes indicates that TcHDAC4 participates in gene expression regulation during T. cruzi differentiation.