Comprehensive transcriptome mining identified the gene expression signature and differentially regulated pathways of the late-onset preeclampsia

Construction of a pathway-level model for preeclampsia based on gene expression data

Preeclampsia (PE) is a heterogeneous disease that seriously affects the health of mothers and fetuses. Lack of detection assays, its diagnosis and intervention are often delayed when the clinical symptoms are atypical. Using personalized pathway-based analysis and machine learning algorithms, we built a PE diagnosis model consisting of nine core pathways using multiple cohorts from the Gene Expression Omnibus database. The model showed an area under the receiver operating characteristic (AUROC) curve of 0.959 with the data from the placental tissue samples in the development cohort. In the two validation cohorts, the AUROCs were 0.898 and 0.876, respectively. The model also performed well with the maternal plasma data in another validation cohort (AUROC: 0.815). Moreover, we identified tyrosine-protein kinase Lck (LCK) as the hub gene in this model and found that LCK and pLCK proteins were downregulated in placentas from PE patients. The pathway-level model for PE can provide a novel direction to develop molecular diagnostic assay and investigate potential mechanisms of PE in future studies.

Pre-eclampsia

Pre-eclampsia is a life-threatening disease of pregnancy unique to humans and a leading cause of maternal and neonatal morbidity and mortality. Women who survive pre-eclampsia have reduced life expectancy, with increased risks of stroke, cardiovascular disease and diabetes, while babies from a pre-eclamptic pregnancy have increased risks of preterm birth, perinatal death and neurodevelopmental disability and cardiovascular and metabolic disease later in life. Pre-eclampsia is a complex multisystem disease, diagnosed by sudden-onset hypertension (>20 weeks of gestation) and at least one other associated complication, including proteinuria, maternal organ dysfunction or uteroplacental dysfunction. Pre-eclampsia is found only when a placenta is or was recently present and is classified as preterm (delivery <37 weeks of gestation), term (delivery ≥37 weeks of gestation) and postpartum pre-eclampsia. The maternal syndrome of pre-eclampsia is driven by a dysfunctional placenta, which releases factors into maternal blood causing systemic inflammation and widespread maternal endothelial dysfunction. Available treatments target maternal hypertension and seizures, but the only 'cure' for pre-eclampsia is delivery of the dysfunctional placenta and baby, often prematurely. Despite decades of research, the aetiology of pre-eclampsia, particularly of term and postpartum pre-eclampsia, remains poorly defined. Significant advances have been made in the prediction and prevention of preterm pre-eclampsia, which is predicted in early pregnancy through combined screening and is prevented with daily low-dose aspirin, starting before 16 weeks of gestation. By contrast, the prediction of term and postpartum pre-eclampsia is limited and there are no preventive treatments. Future research must investigate the pathogenesis of pre-eclampsia, in particular of term and postpartum pre-eclampsia, and evaluate new prognostic tests and treatments in adequately powered clinical trials.

Bioinformatics analysis combined with clinical sample screening reveals that leptin may be a biomarker of preeclampsia

Introduction: Preeclampsia (PE) is a gestational hypertensive disease with unclear pathogenesis. This study aimed to identify the genes that play an important role in determining the pathogenesis of PE using bioinformatics analysis and fundamental researches. Materials and methods: Datasets from the Gene Expression Omnibus (GEO) database were used to screen for differentially expressed genes (DEGs). The NCBI, SangerBox, and other databases were used to analyze the functions of the DEGs. Targetscan7, miRWalk, ENCORI, DIANA TOOLS, CircBank databases, and the Cytoscape tool were used to construct the lncRNA/circRNA-miRNA- LEP network. SRAMP, RPISeq, RBPsuite, and catRPAID were used to analyze the RNA modifications of LEP. Immune cell infiltration was analyzed using the dataset GSE75010. Placental tissues from normal pregnant women and PE patients were collected, screened for gene expression using reverse transcription quantitative polymerase chain reaction (RT-qPCR) and western blotting. The results were further verified in HTR-8/SVneo cell line hypoxia model and PE mouse model. Results: Our analyses revealed that LEP was significantly upregulated in eight datasets. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses indicated that LEP was involved in the JAK/STAT signaling pathway, angiogenesis, and placental development. Immune cell infiltration analysis showed that M1 and M2 macrophages differed between normal pregnancies and those in PE patients. A competing endogenous RNA (ceRNA) network was constructed, and proteins interacting with LEP were identified. RNA modification sites of LEP were also identified. Finally, the overexpression of LEP in PE was confirmed in clinical samples, HTR-8/SVneo cell line and PE mouse model. Conclusion: Our results indicate that LEP overexpression is associated with PE and may be a potential diagnostic marker and therapeutic target.

Genome-Wide Copy Number Variant and High-Throughput Transcriptomics Analyses of Placental Tissues Underscore Persisting Child Susceptibility in At-Risk Pregnancies Cleared in Standard Genetic Testing

Several studies have shown that children from pregnancies with estimated first-trimester risk based on fetal nuchal translucency thickness and abnormal maternal serum pregnancy protein and hormone levels maintain a higher likelihood of adverse outcomes, even if initial testing for known genetic conditions is negative. We used the Finnish InTraUterine cohort (ITU), which is a comprehensively characterized perinatal cohort consisting of 943 mothers and their babies followed throughout pregnancy and 18 months postnatally, including mothers shortlisted for prenatal genetic testing but cleared for major aneuploidies (cases: n = 544, 57.7%) and control pregnancies (n = 399, 42.3%). Using genome-wide genotyping and RNA sequencing of first-trimester and term placental tissue, combined with medical information from registry data and maternal self-report data, we investigated potential negative medical outcomes and genetic susceptibility to disease and their correlates in placenta gene expression. Case mothers did not present with higher levels of depression, perceived stress, or anxiety during pregnancy. Case children were significantly diagnosed more often with congenital malformations of the circulatory system (4.12 (95% CI [1.22−13.93]) higher hazard) and presented with significantly more copy number duplications as compared to controls (burden analysis, based on all copy number variants (CNVs) with at most 10% frequency, 823 called duplications in 297 cases versus 626 called duplications in 277 controls, p = 0.01). Fifteen genes showed differential gene expression (FDR < 0.1) in association with congenital malformations in first-trimester but not term placenta. These were significantly enriched for genes associated with placental dysfunction. In spite of normal routine follow-up prenatal testing results in early pregnancy, case children presented with an increased likelihood of negative outcomes, which should prompt vigilance in follow-up during pregnancy and after birth.

An integrated approach to understand fluid shear stress-driven and reactive oxygen species-mediated metastasis of colon adenocarcinoma through mRNA-miRNA-lncRNA-circRNA networks

Development of colon adenocarcinoma (COAD) metastasis involves several mediators including fluid shear stress (FSS), intracellular ROS levels, and non-coding RNAs. In our present study, we identified and investigated the role of regulatory non-coding RNA molecules specifically involved in COAD metastasis and their association with FSS and ROS. Interactions between the mRNAs associated with FSS and ROS, the corresponding microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) in COAD metastasis were used to generate the mRNA-miRNA-lncRNA-circRNA network. Experimental validation of the identified RNA hubs using quantitative real-time PCR demonstrated a direct effect of the FSS on their expression levels in cancer cells. FSS resulted in the downregulation of HMGA1 and RAN, as well as the upregulation of HSP90AA1, PMAIP1 and BIRC5. Application of shear stress also led to downregulation of hsa-miR-26b-5p and hsa-miR-34a-5p levels in HCT116 cells. Further, functional enrichment and survival analysis of the significant miRNAs, as well as the OncoPrint and the survival analyses of the selected mRNAs were performed. Subsequently, their functional role was also corroborated with existing literature. Ten significant miRNA hubs were identified, out of which hsa-miR-17-5p and hsa-miR-20a-5p were found to interact with lncRNA (CCAT2) while hsa-miR-335 was found to interact with four circRNAs. Fifteen significant miRNAs were identified in 10 different modules suggesting their importance in FSS and ROS-mediated COAD metastasis. Finally, 10 miRNAs and 3 mRNAs associated with FSS and/or ROS were identified as significant overall survival markers; 33 mRNAs were also identified as metastasis-free survival markers whereas 15 mRNAs showed > 10% gene alterations in TCGA-COAD data and may serve as promising therapeutic biomarkers in the COAD metastasis.

Stress Reactivity, Susceptibility to Hypertension, and Differential Expression of Genes in Hypertensive Compared to Normotensive Patients

Although half of hypertensive patients have hypertensive parents, known hypertension-related human loci identified by genome-wide analysis explain only 3% of hypertension heredity. Therefore, mainstream transcriptome profiling of hypertensive subjects addresses differentially expressed genes (DEGs) specific to gender, age, and comorbidities in accordance with predictive preventive personalized participatory medicine treating patients according to their symptoms, individual lifestyle, and genetic background. Within this mainstream paradigm, here, we determined whether, among the known hypertension-related DEGs that we could find, there is any genome-wide hypertension theranostic molecular marker applicable to everyone, everywhere, anytime. Therefore, we sequenced the hippocampal transcriptome of tame and aggressive rats, corresponding to low and high stress reactivity, an increase of which raises hypertensive risk; we identified stress-reactivity-related rat DEGs and compared them with their known homologous hypertension-related animal DEGs. This yielded significant correlations between stress reactivity-related and hypertension-related fold changes (log2 values) of these DEG homologs. We found principal components, PC1 and PC2, corresponding to a half-difference and half-sum of these log2 values. Using the DEGs of hypertensive versus normotensive patients (as the control), we verified the correlations and principal components. This analysis highlighted downregulation of β-protocadherins and hemoglobin as whole-genome hypertension theranostic molecular markers associated with a wide vascular inner diameter and low blood viscosity, respectively.