Single-cell assessment of primary and stem cell-derived human trophoblast organoids as placenta-modeling platforms

Human trophoblast stem cells (hTSCs) and related trophoblast organoids are state-of-the-art culture systems that facilitate the study of trophoblast development and human placentation. Using single-cell transcriptomics, we evaluate how organoids derived from freshly isolated first-trimester trophoblasts or from established hTSC cell lines reproduce developmental cell trajectories and transcriptional regulatory processes defined in vivo. Although organoids from primary trophoblasts and hTSCs overall model trophoblast differentiation with accuracy, specific features related to trophoblast composition, trophoblast differentiation, and transcriptional drivers of trophoblast development show levels of misalignment. This is best illustrated by the identification of an expanded progenitor state in stem cell-derived organoids that is nearly absent in vivo and transcriptionally shares both villous cytotrophoblast and extravillous trophoblast characteristics. Together, this work provides a comprehensive resource that identifies strengths and limitations of current trophoblast organoid platforms.

Sex-dependent placental methylation quantitative trait loci provide insight into the prenatal origins of childhood onset traits and conditions

Molecular quantitative trait loci (QTLs) allow us to understand the biology captured in genome-wide association studies (GWASs). The placenta regulates fetal development and shows sex differences in DNA methylation. We therefore hypothesized that placental methylation QTL (mQTL) explain variation in genetic risk for childhood onset traits, and that effects differ by sex. We analyzed 411 term placentas from two studies and found 49,252 methylation (CpG) sites with mQTL and 2,489 CpG sites with sex-dependent mQTL. All mQTL were enriched in regions that typically affect gene expression in prenatal tissues. All mQTL were also enriched in GWAS results for growth- and immune-related traits, but male- and female-specific mQTL were more enriched than cross-sex mQTL. mQTL colocalized with trait loci at 777 CpG sites, with 216 (28%) specific to males or females. Overall, mQTL specific to male and female placenta capture otherwise overlooked variation in childhood traits.

The application of epiphenotyping approaches to DNA methylation array studies of the human placenta

BACKGROUND

Genome-wide DNA methylation (DNAme) profiling of the placenta with Illumina Infinium Methylation bead arrays is often used to explore the connections between in utero exposures, placental pathology, and fetal development. However, many technical and biological factors can lead to signals of DNAme variation between samples and between cohorts, and understanding and accounting for these factors is essential to ensure meaningful and replicable data analysis. Recently, "epiphenotyping" approaches have been developed whereby DNAme data can be used to impute information about phenotypic variables such as gestational age, sex, cell composition, and ancestry. These epiphenotypes offer avenues to compare phenotypic data across cohorts, and to understand how phenotypic variables relate to DNAme variability. However, the relationships between placental epiphenotyping variables and other technical and biological variables, and their application to downstream epigenome analyses, have not been well studied.

RESULTS

Using DNAme data from 204 placentas across three cohorts, we applied the PlaNET R package to estimate epiphenotypes gestational age, ancestry, and cell composition in these samples. PlaNET ancestry estimates were highly correlated with independent polymorphic ancestry-informative markers, and epigenetic gestational age, on average, was estimated within 4 days of reported gestational age, underscoring the accuracy of these tools. Cell composition estimates varied both within and between cohorts, as well as over very long placental processing times. Interestingly, the ratio of cytotrophoblast to syncytiotrophoblast proportion decreased with increasing gestational age, and differed slightly by both maternal ethnicity (lower in white vs. non-white) and genetic ancestry (lower in higher probability European ancestry). The cohort of origin and cytotrophoblast proportion were the largest drivers of DNAme variation in this dataset, based on their associations with the first principal component.

CONCLUSIONS

This work confirms that cohort, array (technical) batch, cell type proportion, self-reported ethnicity, genetic ancestry, and biological sex are important variables to consider in any analyses of Illumina DNAme data. We further demonstrate the specific utility of epiphenotyping tools developed for use with placental DNAme data, and show that these variables (i) provide an independent check of clinically obtained data and (ii) provide a robust approach to compare variables across different datasets. Finally, we present a general framework for the processing and analysis of placental DNAme data, integrating the epiphenotype variables discussed here.

The application of epiphenotyping approaches to DNA methylation array studies of the human placenta

Background : Genome-wide DNA methylation (DNAme) profiling of the placenta with Illumina Infinium Methylation bead arrays is often used to explore the connections between in utero exposures, placental pathology, and fetal development. However, many technical and biological factors can lead to signals of DNAme variation between samples and between cohorts, and understanding and accounting for these factors is essential to ensure meaningful and replicable data analysis. Recently, "epiphenotyping" approaches have been developed whereby DNAme data can be used to impute information about phenotypic variables such as gestational age, sex, cell composition, and ancestry. These epiphenotypes offer avenues to compare phenotypic data across cohorts, and to understand how phenotypic variables relate to DNAme variability. However, the relationships between placental epiphenotyping variables and other technical and biological variables, and their application to downstream epigenome analyses, have not been well studied. Results : Using DNAme data from 204 placentas across three cohorts, we applied the PlaNET R package to estimate epiphenotypes gestational age, ancestry, and cell composition in these samples. PlaNET ancestry estimates were highly correlated with independent polymorphic ancestry informative markers, and epigenetic gestational age, on average, was estimated within 4 days of reported gestational age, underscoring the accuracy of these tools. Cell composition estimates varied both within and between cohorts, but reassuringly were robust to placental processing time. Interestingly, the ratio of cytotrophoblast to syncytiotrophoblast proportion decreased with increasing gestational age, and differed slightly by both maternal ethnicity (lower in white vs. non-white) and genetic ancestry (lower in higher probability European ancestry). The cohort of origin and cytotrophoblast proportion were the largest drivers of DNAme variation in this dataset, based on their associations with the first principal component. Conclusions : This work confirms that cohort, array (technical) batch, cell type proportion, self-reported ethnicity, genetic ancestry, and biological sex are important variables to consider in any analyses of Illumina DNAme data. Further, we demonstrate that estimating epiphenotype variables from the DNAme data itself, when possible, provides both an independent check of clinically-obtained data and can provide a robust approach to compare variables across different datasets.

Who's afraid of the X? Incorporating the X and Y chromosomes into the analysis of DNA methylation array data

BACKGROUND

Many human disease phenotypes manifest differently by sex, making the development of methods for incorporating X and Y-chromosome data into analyses vital. Unfortunately, X and Y chromosome data are frequently excluded from large-scale analyses of the human genome and epigenome due to analytical complexity associated with sex chromosome dosage differences between XX and XY individuals, and the impact of X-chromosome inactivation (XCI) on the epigenome. As such, little attention has been given to considering the methods by which sex chromosome data may be included in analyses of DNA methylation (DNAme) array data.

RESULTS

With Illumina Infinium HumanMethylation450 DNAme array data from 634 placental samples, we investigated the effects of probe filtering, normalization, and batch correction on DNAme data from the X and Y chromosomes. Processing steps were evaluated in both mixed-sex and sex-stratified subsets of the analysis cohort to identify whether including both sexes impacted processing results. We found that identification of probes that have a high detection p-value, or that are non-variable, should be performed in sex-stratified data subsets to avoid over- and under-estimation of the quantity of probes eligible for removal, respectively. All normalization techniques investigated returned X and Y DNAme data that were highly correlated with the raw data from the same samples. We found no difference in batch correction results after application to mixed-sex or sex-stratified cohorts. Additionally, we identify two analytical methods suitable for XY chromosome data, the choice between which should be guided by the research question of interest, and we performed a proof-of-concept analysis studying differential DNAme on the X and Y chromosome in the context of placental acute chorioamnionitis. Finally, we provide an annotation of probe types that may be desirable to filter in X and Y chromosome analyses, including probes in repetitive elements, the X-transposed region, and cancer-testis gene promoters.

CONCLUSION

While there may be no single "best" approach for analyzing DNAme array data from the X and Y chromosome, analysts must consider key factors during processing and analysis of sex chromosome data to accommodate the underlying biology of these chromosomes, and the technical limitations of DNA methylation arrays.

Profiling placental DNA methylation associated with maternal SSRI treatment during pregnancy

Selective serotonin reuptake inhibitors (SSRIs) for treatment of prenatal maternal depression have been associated with neonatal neurobehavioral disturbances, though the molecular mechanisms remain poorly understood.  In utero exposure to SSRIs may affect DNA methylation (DNAme) in the human placenta, an epigenetic mark that is established during development and is associated with gene expression. Chorionic villus samples from 64 human placentas were profiled with the Illumina MethylationEPIC BeadChip; clinical assessments of maternal mood and SSRI treatment records were collected at multiple time points during pregnancy. Case distribution was 20 SSRI-exposed cases and 44 SSRI non-exposed cases. Maternal depression was defined using a mean maternal Hamilton Depression score > 8 to indicate symptomatic depressed mood ("maternally-depressed"), and we further classified cases into SSRI-exposed, maternally-depressed (n = 14); SSRI-exposed, not maternally-depressed (n = 6); SSRI non-exposed, maternally-depressed (n = 20); and SSRI non-exposed, not maternally-depressed (n = 24). For replication, Illumina 450K DNAme profiles were obtained from 34 additional cases from an independent cohort (n = 17 SSRI-exposed, n = 17 SSRI non-exposed). No CpGs were differentially methylated at FDR < 0.05 comparing SSRI-exposed to non-exposed placentas, in a model adjusted for mean maternal Hamilton Depression score, or in a model restricted to maternally-depressed cases with and without SSRI exposure. However, at a relaxed threshold of FDR < 0.25, five CpGs were differentially methylated (|Δβ| > 0.03) by SSRI exposure status. Four were covered by the replication cohort measured by the 450K array, but none replicated. No CpGs were differentially methylated (FDR < 0.25) comparing maternally depressed to not depressed cases. In sex-stratified analyses for SSRI-exposed versus non-exposed cases (females n = 31; males n = 33), three additional CpGs in females, but none in males, were differentially methylated at the relaxed FDR < 0.25 cut-off. We did not observe large-scale alterations of DNAme in placentas exposed to maternal SSRI treatment, as compared to placentas with no SSRI exposure. We also found no evidence for altered DNAme in maternal depression-exposed versus depression non-exposed placentas. This novel work in a prospectively-recruited cohort with clinician-ascertained SSRI exposure and mood assessments would benefit from future replication.

Genetic variation in placental insufficiency: What have we learned over time?

Genetic variation shapes placental development and function, which has long been known to impact fetal growth and pregnancy outcomes such as miscarriage or maternal pre-eclampsia. Early epidemiology studies provided evidence of a strong heritable component to these conditions with both maternal and fetal-placental genetic factors contributing. Subsequently, cytogenetic studies of the placenta and the advent of prenatal diagnosis to detect chromosomal abnormalities provided direct evidence of the importance of spontaneously arising genetic variation in the placenta, such as trisomy and uniparental disomy, drawing inferences that remain relevant to this day. Candidate gene approaches highlighted the role of genetic variation in genes influencing immune interactions at the maternal-fetal interface and angiogenic factors. More recently, the emergence of molecular techniques and in particular high-throughput technologies such as Single-Nucleotide Polymorphism (SNP) arrays, has facilitated the discovery of copy number variation and study of SNP associations with conditions related to placental insufficiency. This review integrates past and more recent knowledge to provide important insights into the role of placental function on fetal and perinatal health, as well as into the mechanisms leading to genetic variation during development.

Genetic and Epigenetic Mechanisms Deregulate the CRL2pVHL Complex in Hepatocellular Carcinoma

Dysregulation of ubiquitin-proteasome pathway genes through copy number alteration, promoter hypomethylation, and miRNA deregulation is involved in cancer development and progression. Further characterizing alterations in these genes may uncover novel drug targets across a range of diseases in which druggable alterations are uncommon, including hepatocellular carcinoma (HCC). We analyzed 377 HCC and 59 adjacent non-malignant liver tissue samples, focusing on alterations to component genes of the widely studied CRL2^pVHL E3 ubiquitin ligase complex. mRNA upregulation of the component genes was common, and was correlated with DNA hypomethylation and copy number increase, but many tumours displayed overexpression that was not explained by either mechanism. Interestingly, we found 66 miRNAs, including 39 previously unannotated miRNAs, that were downregulated in HCC and predicted to target one or more CRL2^pVHL components. Several miRNAs, including hsa-miR-101-3p and hsa-miR-139-5p, were negatively correlated with multiple component genes, suggesting that miRNA deregulation may contribute to CRL2^pVHL overexpression. Combining miRNA and mRNA expression, DNA copy number, and methylation status into one multidimensional survival analysis, we found a significant association between greater numbers of alterations and poorer overall survival for multiple component genes. While the intricacies of CRL2^pVHL complex gene regulation require additional research, it is evident that multiple causes for the deregulation of these genes must be considered in HCC, including non-traditional mechanisms.

Small Non-Coding RNAs in the Human Placenta: Regulatory Roles and Clinical Utility

The placenta is a vital organ formed during pregnancy, and being the interface between the mother and fetus, it is paramount that placental functioning is strictly controlled. Gene expression in the placenta is finely tuned-with aberrant expression causing placental pathologies and inducing stress on both mother and fetus. Gene regulation is brought upon by several mechanisms, and small non-coding RNAs (sncRNAs) have recently been appreciated for their contribution in gene repression. Their dysregulation has been implicated in a range of somatic and inherited disorders, highlighting their importance in maintaining healthy organ function. Their specific roles within the placenta, however, are not well understood, and require further exploration. To this end, we summarize the mechanisms of microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), small nuclear RNAs (snRNAs), small nucleolar RNAs (snoRNAs), and transfer RNAs (tRNAs), their known contributions to human placental health and disease, the relevance of sncRNAs as promising biomarkers throughout pregnancy, and the current challenges faced by placental sncRNA studies.

Are sex differences in cognitive impairment reflected in epigenetic age acceleration metrics?

Sex differences are well-established in Alzheimer's disease (AD) frequency and pathogenesis, but are not mechanistically understood. Accelerated epigenetic age has been associated with both cognitive aging and AD pathophysiology, but has not been studied by sex in AD or related cognitive impairment. Using the ADNI cohort, we found that none of sex, cognitive impairment diagnosis, nor load of APOEε4 alleles (strongest genetic AD risk factor) were associated with epigenetic age acceleration (DNAmAge, Intrinsic DNAmAge, PhenoAge, or GrimAge), although females exhibit more accelerated epigenetic aging using the Skin & Blood clock in the transition from normal cognition to cognitive impairment than males.

Human placental piwi-interacting RNA transcriptome is characterized by expression from the DLK1-DIO3 imprinted region

The placenta is vital to embryonic development and requires a finely-tuned pattern of gene expression, achieved in part by its unique epigenetic landscape. Piwi-interacting RNAs (piRNAs) are a class of small-non-coding RNA with established roles as epigenetic regulators of gene expression, largely via methylation of targeted DNA sequences. The expression of piRNAs have mainly been described in germ cells, but a fraction have been shown to retain expression in adult somatic tissues. To aid in understanding the contribution of these regulators in the placenta, we provide the first description of the piRNA transcriptome in human placentas. We find 297 piRNAs to be preferentially expressed in the human placenta, a subset of which are expressed at higher levels relative to testes samples. We also observed a large proportion of placental piRNAs to be expressed from a single locus, as distinct from canonical cluster locations associated with transposable element silencing. Finally, we find that 15 of the highest-expressed placental piRNAs maps to the DLK1-DIO3 locus, suggesting a link to placental biology. Our findings suggest that piRNAs could contribute to the molecular networks defining placental function in humans, and a biological impact of piRNA expression beyond germ cells.

Sex Differences Are Here to Stay: Relevance to Prenatal Care

Sex differences exist in the incidence and presentation of many pregnancy complications, including but not limited to pregnancy loss, spontaneous preterm birth, and fetal growth restriction. Sex differences arise very early in development due to differential gene expression from the X and Y chromosomes, and later may also be influenced by the action of gonadal steroid hormones. Though offspring sex is not considered in most prenatal diagnostic or therapeutic strategies currently in use, it may be beneficial to consider sex differences and the associated mechanisms underlying pregnancy complications. This review will cover (i) the prevalence and presentation of sex differences that occur in perinatal complications, particularly with a focus on the placenta; (ii) possible mechanisms underlying the development of sex differences in placental function and pregnancy phenotypes; and (iii) knowledge gaps that should be addressed in the development of diagnostic or risk prediction tools for such complications, with an emphasis on those for which it would be important to consider sex.

Reactivation of Multiple Fetal miRNAs in Lung Adenocarcinoma

MicroRNAs (miRNAs) play vital roles in the regulation of normal developmental pathways. However, cancer cells can co-opt these miRNAs, and the pathways that they regulate, to drive pro-tumourigenic phenotypes. Characterization of the miRNA transcriptomes of fetal organs is essential for identifying these oncofetal miRNAs, but it has been limited by fetal sample availability. As oncofetal miRNAs are absent from healthy adult lungs, they represent ideal targets for developing diagnostic and therapeutic strategies. We conducted small RNA sequencing of a rare collection of 25 human fetal lung (FL) samples and compared them to two independent cohorts (n = 140, n = 427), each comprised of adult non-neoplastic lung (ANL) and lung adenocarcinoma (LUAD) samples. We identified 13 oncofetal miRNAs that were expressed in FL and LUAD but not in ANL. These oncofetal miRNAs are potential biomarkers for LUAD detection (AUC = 0.963). Five of these miRNAs are derived from the imprinted C14MC miRNA cluster at the 14q32 locus, which has been associated with cancer development and abnormal fetal and placental development. Additionally, we observed the pulmonary expression of 44 previously unannotated miRNAs. The sequencing of these fetal lung samples also provides a baseline resource against which aberrant samples can be compared.

A cross-cohort analysis of autosomal DNA methylation sex differences in the term placenta

BACKGROUND

Human placental DNA methylation (DNAme) data is a valuable resource for studying sex differences during gestation, as DNAme profiles after delivery reflect the cumulative effects of gene expression patterns and exposures across gestation. Here, we present an analysis of sex differences in autosomal DNAme in the uncomplicated term placenta (n = 343) using the Illumina 450K array.

RESULTS

At a false discovery rate < 0.05 and a mean sex difference in DNAme beta value of > 0.10, we identified 162 autosomal CpG sites that were differentially methylated by sex and replicated in an independent cohort of samples (n = 293). Several of these differentially methylated CpG sites were part of larger correlated regions of sex differential DNAme. Although global DNAme levels did not differ by sex, the majority of significantly differentially methylated CpGs were more highly methylated in male placentae, the opposite of what is seen in differential methylation analyses of somatic tissues. Patterns of autosomal DNAme at these 162 CpGs were significantly associated with maternal age (in males) and newborn birthweight standard deviation (in females).

CONCLUSIONS

Our results provide a comprehensive analysis of sex differences in autosomal DNAme in the term human placenta. We report a list of high-confidence autosomal sex-associated differentially methylated CpGs and identify several key features of these loci that suggest their relevance to sex differences observed in normative and complicated pregnancies.

Confined placental mosaicism involving multiple de novo copy number variants associated with fetal growth restriction: A case report

The presence of multiple large (>1 Mb) copy number variants (CNVs) in non‐malignant tissue is rare in human genetics. We present a liveborn male with a birth weight below the first percentile associated with placental mosaicism involving eight 2.4–3.9 Mb de novo duplications. We found that the duplications likely co‐localized to the same cells, were mosaic in the placenta, and impacted maternal and paternal chromosomes. In addition, 27.4 Mb and 240 genes were duplicated in affected cells, including candidate placental genes KISS1 and REN. We ruled out involvement of homologous recombination‐based mechanisms or an altered epigenome in generating the CNVs. This case highlights the diversity of genetic abnormalities in the human placenta and the gaps in our knowledge of how such errors arise.

Genomic imbalances in the placenta are associated with poor fetal growth

Background

Fetal growth restriction (FGR) is associated with increased risks for complications before, during, and after birth, in addition to risk of disease through to adulthood. Although placental insufficiency, failure to supply the fetus with adequate nutrients, underlies most cases of FGR, its causes are diverse and not fully understood. One of the few diagnosable causes of placental insufficiency in ongoing pregnancies is the presence of large chromosomal imbalances such as trisomy confined to the placenta; however, the impact of smaller copy number variants (CNVs) has not yet been adequately addressed. In this study, we confirm the importance of placental aneuploidy, and assess the potential contribution of CNVs to fetal growth.

Methods

We used molecular-cytogenetic approaches to identify aneuploidy in placentas from 101 infants born small-for-gestational age (SGA), typically used as a surrogate for FGR, and from 173 non-SGA controls from uncomplicated pregnancies. We confirmed aneuploidies and assessed mosaicism by microsatellite genotyping. We then profiled CNVs using high-resolution microarrays in a subset of 53 SGA and 61 control euploid placentas, and compared the load, impact, gene enrichment and clinical relevance of CNVs between groups. Candidate CNVs were confirmed using quantitative PCR.

Results

Aneuploidy was over tenfold more frequent in SGA-associated placentas compared to controls (11.9% vs. 1.1%; p = 0.0002, OR = 11.4, 95% CI 2.5–107.4), was confined to the placenta, and typically involved autosomes, whereas only sex chromosome abnormalities were observed in controls. We found no significant difference in CNV load or number of placental-expressed or imprinted genes in CNVs between SGA and controls, however, a rare and likely clinically-relevant germline CNV was identified in 5.7% of SGA cases. These CNVs involved candidate genes INHBB, HSD11B2, CTCF, and CSMD3.

Conclusions

We conclude that placental genomic imbalances at the cytogenetic and submicroscopic level may underlie up to ~ 18% of SGA cases in our population. This work contributes to the understanding of the underlying causes of placental insufficiency and FGR, which is important for counselling and prediction of long term outcomes for affected cases.

Cell-specific characterization of the placental methylome

BACKGROUND

DNA methylation (DNAm) profiling has emerged as a powerful tool for characterizing the placental methylome. However, previous studies have focused primarily on whole placental tissue, which is a mixture of epigenetically distinct cell populations. Here, we present the first methylome-wide analysis of first trimester (n = 9) and term (n = 19) human placental samples of four cell populations: trophoblasts, Hofbauer cells, endothelial cells, and stromal cells, using the Illumina EPIC methylation array, which quantifies DNAm at > 850,000 CpGs.

RESULTS

The most distinct DNAm profiles were those of placental trophoblasts, which are central to many pregnancy-essential functions, and Hofbauer cells, which are a rare fetal-derived macrophage population. Cell-specific DNAm occurs at functionally-relevant genes, including genes associated with placental development and preeclampsia. Known placental-specific methylation marks, such as those associated with genomic imprinting, repetitive element hypomethylation, and placental partially methylated domains, were found to be more pronounced in trophoblasts and often absent in Hofbauer cells. Lastly, we characterize the cell composition and cell-specific DNAm dynamics across gestation.

CONCLUSIONS

Our results provide a comprehensive analysis of DNAm in human placental cell types from first trimester and term pregnancies. This data will serve as a useful DNAm reference for future placental studies, and we provide access to this data via download from GEO (GSE159526), through interactive exploration from the web browser ( https://robinsonlab.shinyapps.io/Placental_Methylome_Browser/ ), and through the R package planet, which allows estimation of cell composition directly from placental DNAm data.

Low oxygen enhances trophoblast column growth by potentiating differentiation of the extravillous lineage and promoting LOX activity

Early placental development and the establishment of the invasive trophoblast lineage take place within a low oxygen environment. However, conflicting and inconsistent findings have obscured the role of oxygen in regulating invasive trophoblast differentiation. In this study, the effect of hypoxic, normoxic and atmospheric oxygen on invasive extravillous pathway progression was examined using a human placental explant model. Here, we show that exposure to low oxygen enhances extravillous column outgrowth and promotes the expression of genes that align with extravillous trophoblast (EVT) lineage commitment. By contrast, supra-physiological atmospheric levels of oxygen promote trophoblast proliferation while simultaneously stalling EVT progression. Low oxygen-induced EVT differentiation coincided with elevated transcriptomic levels of lysyl oxidase (LOX) in trophoblast anchoring columns, in which functional experiments established a role for LOX activity in promoting EVT column outgrowth. The findings of this work support a role for low oxygen in potentiating the differentiation of trophoblasts along the extravillous pathway. In addition, these findings generate insight into new molecular processes controlled by oxygen during early placental development.

Renpenning syndrome in a female

Renpenning syndrome (OMIM: 309500) is a rare X-linked disorder that causes intellectual disability, microcephaly, short stature, a variety of eye anomalies, and characteristic craniofacial features. This condition results from pathogenic variation of PQBP1, a polyglutamine-binding protein involved in transcription and pre-mRNA splicing. Renpenning syndrome has only been reported in affected males. Carrier females do not usually have clinical features, and in reported families with Renpenning syndrome, most female carriers exhibit favorable skewing of X-chromosome inactivation. We describe a female with syndromic features typical of Renpenning syndrome. She was identified by exome sequencing to have a de novo heterozygous c.459_462delAGAG mutation in PQBP1 (Xp11.23), affecting the AG hexamer in exon 4, which is the most common causative mutation in this syndrome. Streaky hypopigmentation of the skin was observed, supporting a hypothesized presence of an actively expressed, PQBP1 mutation-bearing X-chromosome in some cells. X-inactivation studies on peripheral blood cells demonstrated complete skewing in both the proband and her mother with preferential inactivation of the maternal X chromosome in the child. We demonstrated expression of the PQBP1 mutant transcript in leukocytes of the affected girl. Therefore, it is highly likely that the PQBP1 mutation arose from the paternal X chromosome.

The significance of the placental genome and methylome in fetal and maternal health

The placenta is a crucial organ for supporting a healthy pregnancy, and defective development or function of the placenta is implicated in a number of complications of pregnancy that affect both maternal and fetal health, including maternal preeclampsia, fetal growth restriction, and spontaneous preterm birth. In this review, we highlight the role of the placental genome in mediating fetal and maternal health by discussing the impact of a variety of genetic alterations, from large whole-chromosome aneuploidies to single-nucleotide variants, on placental development and function. We also discuss the placental methylome in relation to its potential applications for refining diagnosis, predicting pathology, and identifying genetic variants with potential functional significance. We conclude that understanding the influence of the placental genome on common placental-mediated pathologies is critical to improving perinatal health outcomes.

Considerations when processing and interpreting genomics data of the placenta

The application of genomic approaches to placental research has opened exciting new avenues to help us understand basic biological properties of the placenta, improve prenatal screening/diagnosis, and measure effects of in utero exposures on child health outcomes. In the last decade, such large-scale genomic data (including epigenomics and transcriptomics) have become more easily accessible to researchers from many disciplines due to the increasing ease of obtaining such data and the rapidly evolving computational tools available for analysis. While the potential of large-scale studies has been widely promoted, less attention has been given to some of the challenges associated with processing and interpreting such data. We hereby share some of our experiences in assessing data quality, reproducibility, and interpretation in the context of genome-wide studies of the placenta, with the aim to improve future studies. There is rarely a single "best" approach, as that can depend on the study question and sample cohort. However, being consistent, thoroughly assessing potential confounders in the data, and communicating key variables in the methods section of the manuscript are critically important to help researchers to collaborate and build on each other's work.

Accurate ethnicity prediction from placental DNA methylation data

Background

The influence of genetics on variation in DNA methylation (DNAme) is well documented. Yet confounding from population stratification is often unaccounted for in DNAme association studies. Existing approaches to address confounding by population stratification using DNAme data may not generalize to populations or tissues outside those in which they were developed. To aid future placental DNAme studies in assessing population stratification, we developed an ethnicity classifier, PlaNET (Placental DNAme Elastic Net Ethnicity Tool), using five cohorts with Infinium Human Methylation 450k BeadChip array (HM450k) data from placental samples that is also compatible with the newer EPIC platform.

Results

Data from 509 placental samples were used to develop PlaNET and show that it accurately predicts (accuracy = 0.938, kappa = 0.823) major classes of self-reported ethnicity/race (African: n = 58, Asian: n = 53, Caucasian: n = 389), and produces ethnicity probabilities that are highly correlated with genetic ancestry inferred from genome-wide SNP arrays (> 2.5 million SNP) and ancestry informative markers (n = 50 SNPs). PlaNET’s ethnicity classification relies on 1860 HM450K microarray sites, and over half of these were linked to nearby genetic polymorphisms (n = 955). Our placental-optimized method outperforms existing approaches in assessing population stratification in placental samples from individuals of Asian, African, and Caucasian ethnicities.

Conclusion

PlaNET provides an improved approach to address population stratification in placental DNAme association studies. The method can be applied to predict ethnicity as a discrete or continuous variable and will be especially useful when self-reported ethnicity information is missing and genotyping markers are unavailable.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0296-3) contains supplementary material, which is available to authorized users.

Placental epigenetic clocks: estimating gestational age using placental DNA methylation levels

The human pan-tissue epigenetic clock is widely used for estimating age across the entire lifespan, but it does not lend itself well to estimating gestational age (GA) based on placental DNAm methylation (DNAm) data. We replicate previous findings demonstrating a strong correlation between GA and genome-wide DNAm changes. Using substantially more DNAm arrays (n=1,102 in the training set) than a previous study, we present three new placental epigenetic clocks: 1) a robust placental clock (RPC) which is unaffected by common pregnancy complications (e.g., gestational diabetes, preeclampsia), and 2) a control placental clock (CPC) constructed using placental samples from pregnancies without known placental pathology, and 3) a refined RPC for uncomplicated term pregnancies. These placental clocks are highly accurate estimators of GA based on placental tissue; e.g., predicted GA based on RPC is highly correlated with actual GA (r>0.95 in test data, median error less than one week). We show that epigenetic clocks derived from cord blood or other tissues do not accurately estimate GA in placental samples. While fundamentally different from Horvath's pan-tissue epigenetic clock, placental clocks closely track fetal age during development and may have interesting applications.

Altered levels of placental miR-338-3p and miR-518b are associated with acute chorioamnionitis and IL6 genotype

Placental-derived miRNAs are attractive candidates as biomarkers of placental health, but their associations with specific pathologies, such as acute chorioamnionitis (aCA), are not well explored. Samples of chorionic villi from 57 placentas (33 aCA and 24 non-aCA) were analyzed. Expression was quantified for six candidate miRNAs (miR-146a, miR-210, miR-223, miR-338-3p, miR-411, and miR-518b), using quantitative real-time PCR. miR-518b and miR-338-3p were differentially expressed between aCA cases and non-aCA cases (Bonferroni-corrected p < 0.05). Further, we observed that placental miR-518b expression was associated with an IL6 SNP (rs1800796), a polymorphism we previously reported as a risk-conferring variant for aCA.

Association of a placental Interleukin-6 genetic variant (rs1800796) with DNA methylation, gene expression and risk of acute chorioamnionitis

BACKGROUND

Acute chorioamnionitis (aCA), inflammation of the placenta and fetal membranes, is a frequently reported lesion in preterm deliveries. Genetic variants in innate immune system genes such as Interleukin-6 (IL6) may contribute to the placenta's inflammatory response, thus predisposing some pregnancies to aCA. These genetic variants may modulate molecular processes such as DNA methylation and gene expression, and in turn might affect susceptibility to aCA. Currently, there is remarkably little research on the role of fetal (placental) genetic variation in aCA. We aimed to explore the associations between genetic variants in candidate immune-system genes and susceptibility towards inflammatory responses in the placenta, which is linked to a strong inflammatory response in the newborn.

METHODS

DNA samples from 269 placentas (72 aCA cases, 197 non-aCA cases) were collected for this study. Samples were genotyped at 55 ancestry informative markers (AIMs) and 16 additional single nucleotide polymorphisms (SNPs) in 12 candidate innate immune system genes using the Sequenom iPLEX Gold Assay. Publicly available datasets were used to obtain DNA methylation (GSE100197, GSE74738, GSE115508, GSE44667, GSE98224) and gene expression data (GSE44711, GSE98224).

RESULTS

Differences in IL6 placental allele frequencies were associated with aCA (rs1800796, p = 0.04) with the CC genotype specifically implicated (OR = 3.1; p = 0.02). In a subset of the placental samples (n = 67; chorionic villi), we showed that the IL6 SNP (rs1800796) was associated with differential DNA methylation in five IL6-related CpG sites (cg01770232, cg02335517, cg07998387, cg13104385, and cg0526589), where individuals with a CC genotype showed higher DNA methylation levels than individuals carrying the GG genotype. Using two publicly available datasets, we observed that the DNA methylation levels at cg01770232 negatively correlated with IL6 gene expression in the placenta (r = - 0.67, p < 0.004; r = - 0.56, p < 2.937e-05).

CONCLUSIONS

We demonstrated that the minor C allele at the IL6 SNP (rs1800796), which is largely limited to East Asian populations, is associated with the presence of aCA. This SNP was associated with increased DNA methylation at a nearby MEPC2 binding site, which was also associated with decreased expression of IL6 in the placenta. Decreased expression of IL6 may increase vulnerability to microbial infection. Additional studies are required to confirm this association in Asian populations with larger sample sizes.

Mining DNA methylation alterations towards a classification of placental pathologies

Placental health is a key component to a successful pregnancy. Placental insufficiency (PI), inadequate nutrient delivery to the fetus, is associated with preeclampsia (PE), a maternal hypertensive disorder, and intrauterine growth restriction (IUGR), pathologically poor fetal growth. PI is more common in early-onset PE (EOPE) than late-onset PE (LOPE). However, the relationship between these disorders remains unclear. While DNA methylation (DNAm) alterations have been identified in PE and IUGR, these entities can overlap and few studies have analysed them separately. This study aims to utilize DNAm profiling to better understand the underlying placental variation associated with PE and IUGR. Placental samples from a discovery (43 controls, 22 EOPE, 18 LOPE, 11 IUGR) and validation cohort (15 controls, 22 EOPE, 11 LOPE) were evaluated using the Illumina HumanMethylation450 array. To account for gestational age (GA) effects, EOPE samples were compared with pre-term births of varying etiologies (GA <37 weeks). LOPE and IUGR were compared with term controls (GA >37 weeks). While 1703 sites were differentially methylated (DM) (FDR < 0.05, Δβ > 0.1) in EOPE, few changes were associated with LOPE (N = 5), or IUGR (N = 0). Of the 1703 EOPE sites, 599 validated in the second cohort. Using these 599 sites, both cohorts clustered into three distinct groups. Interestingly, LOPE samples diagnosed between 34 and 36 weeks with co-occurring IUGR clustered with the EOPE. DNAm profiling may provide an independent tool to refine clinical/pathological diagnoses into subgroups with more uniform pathology. Despite large changes observed in EOPE, there were challenges in reproducing genome-wide DNAm hits that are discussed.

CHIPS-Child: Testing the developmental programming hypothesis in the offspring of the CHIPS trial

OBJECTIVES

As a follow-up to the CHIPS trial (Control of Hypertension In Pregnancy Study) of 'less tight' (versus 'tight') control of maternal blood pressure in pregnancy, CHIPS-Child investigated potential developmental programming of maternal blood pressure control in pregnancy, by examining measures of postnatal growth rate and hypothalamic-pituitary adrenal (HPA) axis activation.

METHODS

CHIPS follow-up was extended to 12 ± 2 months corrected post-gestational age for anthropometry (weight, length, head/waist circumference). For eligible children with consent for a study visit, we collected biological samples (hair/buccal samples) to evaluate HPA axis function (hair cortisol levels) and epigenetic change (DNA methylation analysis of buccal cells). The primary outcome was 'change in z-score for weight' between birth and 12 ± 2 mos. Secondary outcomes were hair cortisol and genome-wide DNA methylation status.

RESULTS

Of 683 eligible babies, 183 (26.8%) were lost to follow-up, 83 (12.2%) declined, 3 (0.4%) agreed only to ongoing contact, and 414 (60.6%) consented. 372/414 (89.9%) had weight measured at 12mos. In 'less tight' (vs. 'tight') control, the primary outcome was similar [-0.26 (-0.53, +0.01); p = 0.14, p_adjusted = 0.06]; median (95% confidence interval) hair cortisol (N = 35 samples) was lower [-496 (-892, -100) ng/g; p = 0.02], and buccal swab DNA methylation (N = 16 samples) was similar. No differences in growth rate could be demonstrated up to 5 years.

CONCLUSIONS

Results demonstrate no compelling evidence for developmental programming of growth or the HPA axis. Clinicians should look to the clinical findings of CHIPS to guide practice. Researchers should seek to replicate these findings and extend outcomes to paediatric blood pressure and neurodevelopment.

DNA methylation profiling of acute chorioamnionitis-associated placentas and fetal membranes: insights into epigenetic variation in spontaneous preterm births

BACKGROUND

Placental inflammation, often presenting as acute chorioamnionitis (aCA), is commonly associated with preterm birth. Preterm birth can have both immediate and long-term adverse effects on the health of the baby. Developing biomarkers of inflammation in the placenta can help to understand its effects and potentially lead to new approaches for rapid prenatal diagnosis of aCA. We aimed to characterize epigenetic variation associated with aCA in placenta (chorionic villi) and fetal membranes (chorion and amnion) to better understand how aCA may impact processes that lead to preterm birth. This study lays the groundwork for development of novel biomarkers for aCA.

METHODS

Samples from 44 preterm placentas (chorionic villi) as well as matched chorion and amnion for 16 of these cases were collected for this study. These samples were profiled using the Illumina Infinium HumanMethylation850 BeadChip to measure DNA methylation (DNAm) at 866,895 CpGs across the genome. An additional 78 placental samples were utilized to independently validate the array findings by pyrosequencing.

RESULTS

Using a false discovery rate of < 0.15 and average group difference in DNAm of > 0.05, 66 differentially methylated (DM) CpG sites were identified between aCA cases and non-aCA cases in chorionic villi. For the majority of these 66 DM CpGs, the DNAm profile of the aCA cases as compared to the non-aCA cases trended in the direction of the blood cell DNAm. Interestingly, neutrophil-specific DNAm signatures, but not those associated with other immune cell types, were capable of separating aCA cases from the non-aCA cases.

CONCLUSIONS

Our results suggest that aCA-associated placentas showed altered DNAm signatures that were not observed in the absence of aCA. This DNAm profile is consistent with the activation of the innate immune response in the placenta and/or reflect increase in neutrophils as a response to inflammation.

Adjusting for Batch Effects in DNA Methylation Microarray Data, a Lesson Learned

It is well-known, but frequently overlooked, that low- and high-throughput molecular data may contain batch effects, i.e., systematic technical variation. Confounding of experimental batches with the variable(s) of interest is especially concerning, as a batch effect may then be interpreted as a biologically significant finding. An integral step toward reducing false discovery in molecular data analysis includes inspection for batch effects and accounting for this signal if present. In a 30-sample pilot Illumina Infinium HumanMethylation450 (450k array) experiment, we identified two sources of batch effects: row and chip. Here, we demonstrate two approaches taken to process the 450k data in which an R function, ComBat, was applied to adjust for the non-biological signal. In the "initial analysis," the application of ComBat to an unbalanced study design resulted in 9,612 and 19,214 significant (FDR < 0.05) DNA methylation differences, despite none present prior to correction. Suspicious of this dramatic change, a "revised processing" included changes to our analysis as well as a greater number of samples, and successfully reduced batch effects without introducing false signal. Our work supports conclusions made by an article previously published in this journal: though the ultimate antidote to batch effects is thoughtful study design, every DNA methylation microarray analysis should inspect, assess and, if necessary, account for batch effects. The analysis experience presented here can serve as a reminder to the broader community to establish research questions a priori, ensure that they match with study design and encourage communication between technicians and analysts.

No evidence for association of MTHFR 677C>T and 1298A>C variants with placental DNA methylation

Background

5,10-Methylenetetrahydrofolate reductase (MTHFR) is a key enzyme in one-carbon metabolism that ensures the availability of methyl groups for methylation reactions. Two single-nucleotide polymorphisms (SNPs) in the MTHFR gene, 677C>T and 1298A>C, result in a thermolabile enzyme with reduced function. These variants, in both the maternal and/or fetal genes, have been associated with pregnancy complications including miscarriage, neural tube defects (NTDs), and preeclampsia (PE), perhaps due to altered capacity for DNA methylation (DNAm). In this study, we assessed the association between MTHFR 677TT and 1298CC genotypes and risk of NTDs, PE, or normotensive intrauterine growth restriction (nIUGR). Additionally, we assessed whether these high-risk genotypes are associated with altered DNAm in the placenta.

Results

In 303 placentas screened for this study, we observed no significant association between the occurrence of NTDs (N = 55), PE (early-onset: N = 28, late-onset: N = 20), or nIUGR (N = 21) and placental (fetal) MTHFR 677TT or 1298CC genotypes compared to healthy pregnancies (N = 179), though a trend of increased 677TT genotype in PE/IUGR together was observed (OR 2.53, p = 0.048). DNAm was profiled in 10 high-risk 677 (677TT + 1298AA), 10 high-risk 1298 (677CC + 1298CC), and 10 reference (677CC + 1298AA) genotype placentas. Linear modeling identified no significantly differentially methylated sites between high-risk 677 or 1298 and reference placentas at a false discovery rate < 0.05 and Δβ ≥ 0.05 using the Illumina Infinium HumanMethylation450 BeadChip. Using a differentially methylated region analysis or separating cytosine-guanine dinucleotides (CpGs) by CpG density to reduce multiple comparisons also did not identify differential methylation. Additionally, there was no consistent evidence for altered methylation of repetitive DNA between high-risk and reference placentas.

Conclusions

We conclude that large-scale, genome-wide disruption in DNAm does not occur in placentas with the high-risk MTHFR 677TT or 1298CC genotypes. Furthermore, there was no evidence for an association of the 1298CC genotype and only a tendency to higher 677TT in pregnancy complications of PE/IUGR. This may be due to small sample sizes or folate repletion in our Canadian population attenuating effects of the high-risk MTHFR variants. However, given our results and the conflicting results in the literature, investigations into alternative mechanisms that may explain the link between MTHFR variants and pregnancy complications, or in populations at risk of folate deficiencies, are warranted.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0468-1) contains supplementary material, which is available to authorized users.

Epigenetic regulation of placental gene expression in transcriptional subtypes of preeclampsia

Background

Preeclampsia (PE) is a heterogeneous, hypertensive disorder of pregnancy, with no robust biomarkers or effective treatments. We hypothesized that this heterogeneity is due to the existence of multiple subtypes of PE and, in support of this hypothesis, we recently identified five clusters of placentas within a large gene expression microarray dataset (N = 330), of which four (clusters 1, 2, 3, and 5) contained a substantial number of PE samples. However, while transcriptional analysis of placentas can subtype patients, we propose that the addition of epigenetic information could discern gene regulatory mechanisms behind the distinct PE pathologies, as well as identify clinically useful potential biomarkers.

Results

We subjected 48 of our samples from transcriptional clusters 1, 2, 3, and 5 to Infinium HumanMethylation450 arrays. Samples belonging to transcriptional clusters 1–3 still showed visible relationships to each other by methylation, but cluster 5, with known chromosomal abnormalities, no longer formed a cohesive group. Within transcriptional clusters 2 and 3, controlling for fetal sex and gestational age in the identification of differentially methylated sites, compared to the healthier cluster 1, dramatically reduced the number of significant sites, but increased the percentage that demonstrated a strong linear correlation with gene expression (from 5% and 2% to 9% and 8%, respectively). Locations exhibiting a positive relationship between methylation and gene expression were most frequently found in CpG open sea enhancer regions within the gene body, while those with a significant negative correlation were often annotated to the promoter in a CpG shore region. Integrated transcriptome and epigenome analysis revealed modifications in TGF-beta signaling, cell adhesion, oxidative phosphorylation, and metabolism pathways in cluster 2 placentas, and aberrations in antigen presentation, allograft rejection, and cytokine-cytokine receptor interaction in cluster 3 samples.

Conclusions

Overall, we have established DNA methylation alterations underlying a portion of the transcriptional development of “canonical” PE in cluster 2 and “immunological” PE in cluster 3. However, a significant number of the observed methylation changes were not associated with corresponding changes in gene expression, and vice versa, indicating that alternate methods of gene regulation will need to be explored to fully comprehend these PE subtypes.

Electronic supplementary material

The online version of this article (10.1186/s13148-018-0463-6) contains supplementary material, which is available to authorized users.

A de novo mosaic mutation in SPAST with two novel alternative alleles and chromosomal copy number variant in a boy with spastic paraplegia and autism spectrum disorder

Here we report a 12 year old male with an extreme presentation of spastic paraplegia along with autism and dysmorphisms. Whole exome sequencing identified a predicted pathogenic pair of missense variants in SPAST at the same chromosomal location, each with a different alternative allele, while a chromosome microarray identified a 1.73 Mb paternally inherited copy gain of 1q21.1q21.2 resulting in a blended phenotype of both Spastic paraplegia 4 and 1q21.1 microduplication syndrome. We believe that the extreme phenotype observed is likely caused by the presence of cells which contain only mutant SPAST, but that the viability of the patient is possible due mosaicism of mutant alleles observed in different proportions across tissues.

Child mortality, hypothalamic-pituitary-adrenal axis activity and cellular aging in mothers

Psychological challenges, including traumatic events, have been hypothesized to increase the age-related pace of biological aging. Here we test the hypothesis that psychological challenges can affect the pace of telomere attrition, a marker of cellular aging, using data from an ongoing longitudinal-cohort study of Kaqchikel Mayan women living in a population with a high frequency of child mortality, a traumatic life event. Specifically, we evaluate the associations between child mortality, maternal telomere length and the mothers’ hypothalamic-pituitary-adrenal axis (HPAA), or stress axis, activity. Child mortality data were collected in 2000 and 2013. HPAA activity was assessed by quantifying cortisol levels in first morning urinary specimens collected every other day for seven weeks in 2013. Telomere length (TL) was quantified using qPCR in 55 women from buccal specimens collected in 2013. Results: Shorter TL with increasing age was only observed in women who experienced child mortality (p = 0.015). Women with higher average basal cortisol (p = 0.007) and greater within-individual variation (standard deviation) in basal cortisol (p = 0.053) presented shorter TL. Non-parametric bootstrapping to estimate mediation effects suggests that HPAA activity mediates the effect of child mortality on TL. Our results are, thus, consistent with the hypothesis that traumatic events can influence cellular aging and that HPAA activity may play a mediatory role. Future large-scale longitudinal studies are necessary to confirm our results and further explore the role of the HPAA in cellular aging, as well as to advance our understanding of the underlying mechanisms involved.

Cord blood hematopoietic cells from preterm infants display altered DNA methylation patterns

Background

Premature infants are highly vulnerable to infection. This is partly attributable to the preterm immune system, which differs from that of the term neonate in cell composition and function. Multiple studies have found differential DNA methylation (DNAm) between preterm and term infants’ cord blood; however, interpretation of these studies is limited by the confounding factor of blood cell composition. This study evaluates the epigenetic impact of preterm birth in isolated hematopoietic cell populations, reducing the concern of cell composition differences.

Methods

Genome-wide DNAm was measured using the Illumina 450K array in T cells, monocytes, granulocytes, and nucleated red blood cells (nRBCs) isolated from cord blood of 5 term and 5 preterm (<31 weeks gestational age) newborns. DNAm of hematopoietic cells was compared globally across the 450K array and through site-specific linear modeling.

Results

Nucleated red blood cells (nRBCs) showed the most extensive changes in DNAm, with 9258 differentially methylated (DM) sites (FDR < 5%, |Δβ| > 0.10) discovered between preterm and term infants compared to the <1000 prematurity-DM sites identified in white blood cell populations. The direction of DNAm change with gestational age at these prematurity-DM sites followed known patterns of hematopoietic differentiation, suggesting that term hematopoietic cell populations are more epigenetically mature than their preterm counterparts. Consistent shifts in DNAm between preterm and term cells were observed at 25 CpG sites, with many of these sites located in genes involved in growth and proliferation, hematopoietic lineage commitment, and the cytoskeleton. DNAm in preterm and term hematopoietic cells conformed to previously identified DNAm signatures of fetal liver and bone marrow, respectively.

Conclusions

This study presents the first genome-wide mapping of epigenetic differences in hematopoietic cells across the late gestational period. DNAm differences in hematopoietic cells between term and <31 weeks were consistent with the hematopoietic origin of these cells during ontogeny, reflecting an important role of DNAm in their regulation. Due to the limited sample size and the high coincidence of prematurity and multiple births, the relationship between cause of preterm birth and DNAm could not be evaluated. These findings highlight gene regulatory mechanisms at both cell-specific and systemic levels that may be involved in fetal immune system maturation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0339-1) contains supplementary material, which is available to authorized users.

An empirically driven data reduction method on the human 450K methylation array to remove tissue specific non-variable CpGs

Background

Population based epigenetic association studies of disease and exposures are becoming more common with the availability of economical genome-wide technologies for interrogation of the methylome, such as the Illumina 450K Human Methylation Array (450K). Often, the expected small number of differentially methylated cytosine-guanine pairs (CpGs) in studies of the human methylome presents a statistical challenge, as the large number of CpGs measured on the 450K necessitates careful multiple test correction. While the 450K is a highly useful tool for population epigenetic studies, many of the CpGs tested are not variable and thus of limited information content in the context of the study and tissue. CpGs with observed lack of variability in the tissue under study could be removed to reduce the data dimensionality, limit the severity of multiple test correction and allow for improved detection of differential DNA methylation.

Methods

Here, we performed a meta-analysis of 450K data from three commonly studied human tissues, namely blood (605 samples), buccal epithelial cells (121 samples) and placenta (157 samples). We developed lists of CpGs that are non-variable in each tissue.

Results

These lists are surprisingly large (blood 114,204 CpGs, buccal epithelial cells 120,009 CpGs and placenta 101,367 CpGs) and thus will be valuable filters for epigenetic association studies, considerably reducing the dimensionality of the 450K and subsequently the multiple testing correction severity.

Conclusions

We propose this empirically derived method for data reduction to allow for more power in detecting differential DNA methylation associated with exposures in studies on the human methylome.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-017-0320-z) contains supplementary material, which is available to authorized users.

Placental telomere length decline with gestational age differs by sex and TERT, DNMT1, and DNMT3A DNA methylation

INTRODUCTION

Telomere length (TL) has been suggested to be influenced by inherited genetic and epigenetic variation, hormonal effects, oxidative stress and age. However, the dynamics of TL during in utero development have not been well explored. This study investigates the relationship between placental TL and sex, gestational age (GA), and DNA methylation (DNAm). Placental TL is further evaluated in pregnancies complicated by preeclampsia (PE) and intrauterine growth restriction (IUGR), conditions hypothesized to lead to decreased placental TL due to increased oxidative stress.

METHODS

Average TL in 21 early-onset PE (EOPE), 18 late-onset PE (LOPE), 9 IUGR, 59 viable and 33 non-viable control placentas were measured by qPCR. Of these, 13 control, 20 EOPE, 17 LOPE, and 8 IUGR samples were also run on the Illumina 450K array. ANOVA was used to compare TL between controls and EOPE, LOPE, and IUGR. Linear regression correcting for GA and sex, assessed the association between TL and DNAm in biologically-relevant genes (TERC, TERT, DNMT1, DNMT3a, DNMT3b), and array-wide.

RESULTS

Male sex and increasing GA were associated with shorter placental TL. Correcting for these factors, no significant difference in TL was observed between EOPE, LOPE, and IUGR placentas compared to controls. Targeted analysis revealed TL was associated with DNAm at TERT, DNMT1, and DNMT3a. An array-wide approach found no additional sites associated with TL.

CONCLUSION

Variability in placental TL is associated with alterations in DNAm at TERT, DNMT1, and DNMT3a. Placental TL is not strongly influenced by EOPE, LOPE, or IUGR.

Characterizing the hypomethylated DNA methylation profile of nucleated red blood cells from cord blood

AIM

To provide insight into fetal nucleated red blood cell (nRBC) development using genome-wide DNA methylation (DNAm) profiling.

MATERIALS & METHODS

The DNAm profile (Illumina 450K array) of cord blood (n = 7) derived nRBCs was compared with B cells, CD4 and CD8 T cells, natural killer cells, granulocytes, monocytes and placenta (n = 5).

RESULTS

nRBCs and placenta had similarly low array-wide DNAm compared with white blood cells, but their patterns of hypomethylation differed at biologically relevant subsets of the array. High interindividual variability in nRBC DNAm was driven by a negative association between DNAm and nRBC count.

CONCLUSION

nRBC hypomethylation is likely an epigenetic signature of erythropoiesis rather than of early development. Variability in nRBC DNAm may stem from differences in the cell population's maturity or hematopoietic source.

Developmental transcription factor NFIB is a putative target of oncofetal miRNAs and is associated with tumour aggressiveness in lung adenocarcinoma

Genes involved in fetal lung development are thought to play crucial roles in the malignant transformation of adult lung cells. Consequently, the study of lung tumour biology in the context of lung development has the potential to reveal key developmentally relevant genes that play critical roles in lung cancer initiation/progression. Here, we describe for the first time a comprehensive characterization of miRNA expression in human fetal lung tissue, with subsequent identification of 37 miRNAs in non-small cell lung cancer (NSCLC) that recapitulate their fetal expression patterns. Nuclear factor I/B (NFIB), a transcription factor essential for lung development, was identified as a potential frequent target for these 'oncofetal' miRNAs. Concordantly, analysis of NFIB expression in multiple NSCLC independent cohorts revealed its recurrent underexpression (in ∼40-70% of tumours). Interrogation of NFIB copy number, methylation, and mutation status revealed that DNA level disruption of this gene is rare, and further supports the notion that oncofetal miRNAs are likely the primary mechanism responsible for NFIB underexpression in NSCLC. Reflecting its functional role in regulating lung differentiation, low expression of NFIB was significantly associated with biologically more aggressive subtypes and, ultimately, poorer survival in lung adenocarcinoma patients. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Profiling placental and fetal DNA methylation in human neural tube defects

Background

The incidence of neural tube defects (NTDs) declined by about 40 % in Canada with the introduction of a national folic acid (FA) fortification program. Despite the fact that few Canadians currently exhibit folate deficiency, NTDs are still the second most common congenital abnormality. FA fortification may have aided in reducing the incidence of NTDs by overcoming abnormal one carbon metabolism cycling, the process which provides one carbon units for methylation of DNA. We considered that NTDs persisting in a folate-replete population may also occur in the context of FA-independent compromised one carbon metabolism, and that this might manifest as abnormal DNA methylation (DNAm). Second trimester human placental chorionic villi, kidney, spinal cord, brain, and muscle were collected from 19 control, 22 spina bifida, and 15 anencephalic fetuses in British Columbia, Canada. DNA was extracted, assessed for methylenetetrahydrofolate reductase (MTHFR) genotype and for genome-wide DNAm using repetitive elements, in addition to the Illumina Infinium HumanMethylation450 (450k) array.

Results

No difference in repetitive element DNAm was noted between NTD status groups. Using a false discovery rate <0.05 and average group difference in DNAm ≥0.05, differentially methylated array sites were identified only in (1) the comparison of anencephaly to controls in chorionic villi (n = 4 sites) and (2) the comparison of spina bifida to controls in kidney (n = 3342 sites).

Conclusions

We suggest that the distinctive DNAm of spina bifida kidneys may be consequent to the neural tube defect or reflective of a common etiology for abnormal neural tube and renal development. Though there were some small shifts in DNAm in the other tested tissues, our data do not support the long-standing hypothesis of generalized altered genome-wide DNAm in NTDs. This finding may be related to the fact that most Canadians are not folate deficient, but it importantly opens the field to the investigation of other epigenetic and non-epigenetic mechanisms in the etiology of NTDs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0054-8) contains supplementary material, which is available to authorized users.

Pervasive polymorphic imprinted methylation in the human placenta

The maternal and paternal copies of the genome are both required for mammalian development, and this is primarily due to imprinted genes, those that are monoallelically expressed based on parent-of-origin. Typically, this pattern of expression is regulated by differentially methylated regions (DMRs) that are established in the germline and maintained after fertilization. There are a large number of germline DMRs that have not yet been associated with imprinting, and their function in development is unknown. In this study, we developed a genome-wide approach to identify novel imprinted DMRs in the human placenta and investigated the dynamics of these imprinted DMRs during development in somatic and extraembryonic tissues. DNA methylation was evaluated using the Illumina HumanMethylation450 array in 134 human tissue samples, publicly available reduced representation bisulfite sequencing in the human embryo and germ cells, and targeted bisulfite sequencing in term placentas. Forty-three known and 101 novel imprinted DMRs were identified in the human placenta by comparing methylation between diandric and digynic triploid conceptions in addition to female and male gametes. Seventy-two novel DMRs showed a pattern consistent with placental-specific imprinting, and this monoallelic methylation was entirely maternal in origin. Strikingly, these DMRs exhibited polymorphic imprinted methylation between placental samples. These data suggest that imprinting in human development is far more extensive and dynamic than previously reported and that the placenta preferentially maintains maternal germline-derived DNA methylation.

Placental DNA methylation at term reflects maternal serum levels of INHA and FN1, but not PAPPA, early in pregnancy

Background

Early detection of pregnancies at risk of complications, such as intrauterine growth restriction (IUGR) and preeclampsia (PE), is critical for improved monitoring and preventative treatment to optimize health outcomes. We predict that levels of placental-derived proteins circulating in maternal blood reflect placental gene expression, which is associated with placental DNA methylation (DNAm) profiles. As such, placental DNAm profiling may be useful to distinguish pregnancies at risk of developing complications and correlation between DNAm and protein levels in maternal blood may give further evidence for a protein’s use as a biomarker. However, few studies investigate all clinical parameters that may influence DNAm and/or protein expression, which can significantly affect the relationship between these measures.

Results

Candidate genes were chosen based on i) reported alterations of protein levels in maternal blood and ii) observed changes in placental DNAm (∆β > 0.05 and False Discovery Rate (FDR) <0.05) in pregnancies complicated by PE/IUGR. Fibronectin (FN1) enhancer DNAm and placental gene expression were inversely correlated (r = −0.88 p < 0.01). The same trend was observed between promoter DNAm and gene expression for INHBA and PAPPA, though not significant. INHBA and FN1 DNAm was associated with gestational–age corrected birth weight, while INHA levels were associated with fetal: placental weight ratio and FN1 level was associated with maternal body mass index (BMI).

DNAm at the INHBA promoter in the term placenta was negatively correlated with second trimester maternal serum levels (r = −0.50 p = 0.01) and DNAm at the FN1 enhancer was negatively associated with third trimester maternal serum levels (r = −0.38, p = 0.009). However, a similar correlation was not found for PAPPA.

Conclusions

These results show that establishing a correlation between altered DNAm in the term placenta and altered maternal serum levels of the corresponding protein, is affected by a number of factors. Nonetheless, the correlation between placental DNAm of INHBA/FN1 and maternal serum INHA/FN1 levels indicate that DNAm may be a useful tool to identify novel biomarkers for adverse pregnancy outcomes in some cases.

Electronic supplementary material

The online version of this article (doi:10.1186/s12881-015-0257-z) contains supplementary material, which is available to authorized users.

Noninvasive nucleic acid-based approaches to monitor placental health and predict pregnancy-related complications

During pregnancy, the placenta releases a variety of nucleic acids (including deoxyribonucleic acid, messenger ribonucleic acid, or microribonucleic acids) either as a result of cell turnover or as an active messaging system between the placenta and cells in the maternal body. The profile of released nucleic acids changes with the gestational age and has been associated with maternal and fetal parameters. It also can directly reflect pathological changes in the placenta. Nucleic acids may therefore provide a rich source of novel biomarkers for the prediction of pregnancy complications. However, their utility in the clinical setting depends, first, on overcoming some technical considerations in their quantification, and, second, on developing a better understanding of the factors that influence their function and abundance.

Defects in fatty acid amide hydrolase 2 in a male with neurologic and psychiatric symptoms

Background

Fatty acid amide hydrolase 2 (FAAH2) is a hydrolase that mediates the degradation of endocannabinoids in man. Alterations in the endocannabinoid system are associated with a wide variety of neurologic and psychiatric conditions, but the phenotype and biochemical characterization of patients with genetic defects of FAAH2 activity have not previously been described. We report a male with autistic features with an onset before the age of 2 years who subsequently developed additional features including anxiety, pseudoseizures, ataxia, supranuclear gaze palsy, and isolated learning disabilities but was otherwise cognitively intact as an adult.

Methods and results

Whole exome sequencing identified a rare missense mutation in FAAH2, hg19: g.57475100G > T (c.1372G > T) resulting in an amino acid change (p.Ala458Ser), which was Sanger confirmed as maternally inherited and absent in his healthy brother. Alterations in lipid metabolism with abnormalities of the whole blood acyl carnitine profile were found. Biochemical and molecular modeling studies confirmed that the p.Ala458Ser mutation results in partial inactivation of FAAH2. Studies in patient derived fibroblasts confirmed a defect in FAAH2 activity resulting in altered levels of endocannabinoid metabolites.

Conclusions

We propose that genetic alterations in FAAH2 activity contribute to neurologic and psychiatric disorders in humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s13023-015-0248-3) contains supplementary material, which is available to authorized users.