Differential gene expression elicited by ZIKV infection in trophoblasts from congenital Zika syndrome discordant twins.

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.

Weighted gene co-expression analyses point to long non-coding RNA hub genes at different Schistosoma mansoni life-cycle stages

Long non-coding RNAs (lncRNAs) (>200 nt) are expressed at levels lower than those of the protein-coding mRNAs, and in all eukaryotic model species where they have been characterized, they are transcribed from thousands of different genomic loci. In humans, some four dozen lncRNAs have been studied in detail, and they have been shown to play important roles in transcriptional regulation, acting in conjunction with transcription factors and epigenetic marks to modulate the tissue-type specific programs of transcriptional gene activation and repression. In Schistosoma mansoni, around 10,000 lncRNAs have been identified in previous works. However, the limited number of RNA-sequencing (RNA-seq) libraries that had been previously assessed, together with the use of old and incomplete versions of the S. mansoni genome and protein-coding transcriptome annotations, have hampered the identification of all lncRNAs expressed in the parasite. Here we have used 633 publicly available S. mansoni RNA-seq libraries from whole worms at different stages (n = 121), from isolated tissues (n = 24), from cell-populations (n = 81), and from single-cells (n = 407). We have assembled a set of 16,583 lncRNA transcripts originated from 10,024 genes, of which 11,022 are novel S. mansoni lncRNA transcripts, whereas the remaining 5,561 transcripts comprise 120 lncRNAs that are identical to and 5,441 lncRNAs that have gene overlap with S. mansoni lncRNAs already reported in previous works. Most importantly, our more stringent assembly and filtering pipeline has identified and removed a set of 4,293 lncRNA transcripts from previous publications that were in fact derived from partially processed mRNAs with intron retention. We have used weighted gene co-expression network analyses and identified 15 different gene co-expression modules. Each parasite life-cycle stage has at least one highly correlated gene co-expression module, and each module is comprised of hundreds to thousands lncRNAs and mRNAs having correlated co-expression patterns at different stages. Inspection of the top most highly connected genes within the modules’ networks has shown that different lncRNAs are hub genes at different life-cycle stages, being among the most promising candidate lncRNAs to be further explored for functional characterization.