One in five pregnant women suffer from gestational complications, prevalently driven by placental malfunction. Using RNASeq, we analyzed differential placental gene expression in cases of normal gestation, late-onset preeclampsia (LO-PE), gestational diabetes (GD) and pregnancies ending with the birth of small-for-gestational-age (SGA) or large-for-gestational-age (LGA) newborns (n = 8/group). In all groups, the highest expression was detected for small noncoding RNAs and genes specifically implicated in placental function and hormonal regulation. The transcriptome of LO-PE placentas was clearly distinct, showing statistically significant (after FDR) expressional disturbances for hundreds of genes. Taqman RT-qPCR validation of 45 genes in an extended sample (n = 24/group) provided concordant results. A limited number of transcription factors including LRF, SP1 and AP2 were identified as possible drivers of these changes. Notable differences were detected in differential expression signatures of LO-PE subtypes defined by the presence or absence of intrauterine growth restriction (IUGR). LO-PE with IUGR showed higher correlation with SGA and LO-PE without IUGR with LGA placentas. Whereas changes in placental transcriptome in SGA, LGA and GD cases were less prominent, the overall profiles of expressional disturbances overlapped among pregnancy complications providing support to shared placental responses. The dataset represent a rich catalogue for potential biomarkers and therapeutic targets.

CDH13 encodes T-cadherin, a receptor for high molecular weight (HMW) adiponectin and low-density lipoprotein, promoting proliferation and migration of endothelial cells. Genome-wide association studies have mapped multiple variants in CDH13 associated with cardiometabolic traits (CMT) with variable effects across studies. We hypothesized that this heterogeneity might reflect interplay with DNA methylation within the region. Resequencing and EpiTYPER™ assay were applied for the HYPertension in ESTonia/Coronary Artery Disease in Czech (HYPEST/CADCZ; n = 358) samples to identify CDH13 promoter SNPs acting as methylation Quantitative Trait Loci (meQTLs) and to investigate their associations with CMT. In silico data were extracted from genome-wide DNA methylation and genotype datasets of the population-based sample Estonian Genome Center of the University of Tartu (EGCUT; n = 165). HYPEST–CADCZ meta-analysis identified a rare variant rs113460564 as highly significant meQTL for a 134-bp distant CpG site (P = 5.90 × 10−6; β = 3.19 %). Four common SNPs (rs12443878, rs12444338, rs62040565, rs8060301) exhibited effect on methylation level of up to 3 neighboring CpG sites in both datasets. The strongest association was detected in EGCUT between rs8060301 and cg09415485 (false discovery rate corrected P value = 1.89 × 10−30). Simultaneously, rs8060301 showed association with diastolic blood pressure, serum high-density lipoprotein and HMW adiponectin (P < 0.005). Novel strong associations were identified between rare CDH13 promoter meQTLs (minor allele frequency <5 %) and HMW adiponectin: rs2239857 (P = 5.50 × 10−5, β = −1,841.9 ng/mL) and rs77068073 (P = 2.67 × 10−4, β = −2,484.4 ng/mL). Our study shows conclusively that CDH13 promoter harbors meQTLs associated with CMTs. It paves the way to deeper understanding of the interplay between DNA variation and methylation in susceptibility to common diseases.Electronic supplementary materialThe online version of this article (doi:10.1007/s00439-014-1521-6) contains supplementary material, which is available to authorized users.

Placenta is a unique organ that serves its own function, and contributes to maternal gestational adaptation and fetal development. Coordination of its transcriptome to satisfy all the maternal-fetal needs across gestation is not fully understood. MicroRNAs are powerful transcriptome modulators capable to adjust rapidly the expression level and dynamics of large gene sets. This MiR-Seq based study presents a multi-omics investigation of the human placental miRNome and its synergy with the transcriptome. The analysis included 52 placentas representing three trimesters of normal pregnancy, and term cases of late-onset preeclampsia (LO-PE), gestational diabetes and affected fetal growth. Gestational-age dependent differential expression (FDR &lt; 0.05) was detected for 319 of 417 tested miRNAs (76.5%). A shared list of target genes of dynamic miRNAs suggested their coordinated action. The most abundant miR-143-3p revealed as a marker for pregnancy progression. The data suggested critical, but distinct roles of placenta-specific imprinted C19MC and C14MC miRNA clusters. Paternally encoded primate-specific C19MC was highly transcribed during first trimester, potentially fine-tuning the early placental transcriptome in dosage-sensitive manner. Maternally encoded eutherian C14MC showed high expression until term, underlining its key contribution across gestation. A major shift in placental miRNome (16% miRNAs) was observed in LO-PE, but not in other term pregnancy complications. Notably, 13/38 upregulated miRNAs were transcribed from C19MC and only one from C14MC, whereas 11/28 downregulated miRNAs represented C14MC and none C19MC. miR-210-3p, miR-512-5p, miR-32-5p, miR-19a-3p, miR-590-3p, miR-379-5p were differentially expressed in LO-PE and cases of small-for-gestational-age newborns, supporting a shared etiology. Expression correlation analysis with the RNA-Seq data (16,567 genes) of the same samples clustered PE-linked miRNAs into five groups. Large notable clusters of miRNA–gene pairs showing directly and inversely correlated expression dynamics suggested potential functional relationships in both scenarios. The first genome-wide study of placental miR-eQTLs identified 66 placental SNVs associated with the expression of neighboring miRNAs, including PE-linked miRNAs miR-30a-5p, miR-210-3p, miR-490-3p and miR-518-5p. This study provided a rich catalog of miRNAs for further in-depth investigations of their individual and joint effect on placental transcriptome. Several highlighted miRNAs may serve as potential biomarkers for pregnancy monitoring and targets to prevent or treat gestational disorders.

The variant rs4769613 T/C within the enhancer element near FLT1 , an acknowledged gene in preeclampsia, was previously identified as a risk factor for preeclampsia in the genome-wide association study (GWAS) targeting placental genotypes. We aimed to test the robustness of this association in 2 Estonian cohorts. Both placental sample sets HAPPY PREGNANCY (Development of novel non-invasive biomarkers for fertility and healthy pregnancy; preeclampsia, n=44 versus nonpreeclampsia, n=1724) and REPROMETA (REPROgrammed fetal and/or maternal METAbolism; 52/277) exhibited suggestive association between rs4769613[C] variant and preeclampsia (logistic regression adjusted for gestational age and fetal sex, nominal P <0.05). Meta-analysis across 2 samples (96/2001) replicated the genome-wide association study outcome (Bonferroni corrected P =4×10 -3 ; odds ratio, 1.75 [95% CI, 1.23–2.49]). No association was detected with gestational diabetes mellitus, preterm birth, and newborn parameters. Also, neither maternal nor paternal rs4769613 genotypes predisposed to preeclampsia. The exact role of placental rs4769613 genotype in the preeclampsia pathogenesis is to be clarified as no effect was detected on maternal baseline serum sFlt-1 (soluble fms-related receptor tyrosine kinase 1) levels. However, when placental FLT1 gene expression and maternal serum sFlt-1 measurements were stratified by placental rs4769613 genotypes, significantly higher transcript and biomarker levels were detected in preeclampsia versus nonpreeclampsia cases in the CC- and CT- (Student t test, P ≤0.02), but not in the TT-genotype subgroup. We suggest that rs4769613 represents a conditional expression Quantitative Trait Locus, whereby only the enhancer with the C-allele reacts to promote the FLT1 expression in unfavorable placental conditions. The study highlighted that the placental FLT1 rs4769613 C-allele is a preeclampsia-specific risk factor. It may contribute to early identification of high-risk women, for example, when genotyped in the cffDNA available in maternal blood plasma.

MicroRNAs are among the key modulators of placental transcriptome dynamics. This study aimed at comparative profiling of urinary (sampling at 228-230 gestational days, g.days), serum (217-230 g.days), and placental (279-286 g.days) microRNAs in three healthy pregnant women using miR-Nome sequencing. Placenta showed significant enrichment of microRNAs compared with serum and urine (1174, 341, and 193, respectively; P < 10 À5 ). All sample types shared 153 microRNAs, representing candidate biomarkers for placental health. Urine samples contained eight of 56 transcripts from the placenta-specific chromosome 19 microRNA cluster C19MC and one of 91 transcripts (miR-432-5p) from the chromosome 14 cluster C14MC. These data suggest an active filtering at the maternal-fetal interface in passing through only selected microRNAs. Urine represents a valid source to monitor the signature of placenta-expressed microRNAs that are differentially expressed in pregnancy complications.