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.

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 CRL2pVHL 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 CRL2pVHL 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 CRL2pVHL 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 CRL2pVHL 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.

Examining the Vanishing Twin Hypothesis of Neural Tube Defects: Application of an Epigenetic Predictor for Monozygotic Twinning.

Strong associations between neural tube defects (NTDs) and monozygotic (MZ) twinning have long been noted, and it has been suggested that NTD cases who do not present as MZ twins may be the survivors of MZ twinning events. We have recently shown that MZ twins carry a strong, distinctive DNA methylation signature and have developed an algorithm based on genomewide DNA methylation array data that distinguishes MZ twins from dizygotic twins and other relatives at well above chance level. We have applied this algorithm to published methylation data from five fetal tissues (placental chorionic villi, kidney, spinal cord, brain and muscle) collected from spina bifida cases (n = 22), anencephalic cases (n = 15) and controls (n = 19). We see no difference in signature between cases and controls, providing no support for a common etiological role of MZ twinning in NTDs. The strong associations therefore continue to await elucidation.

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.

Mistakes Are Common; Should We Worry about Them?

Mutations arising early in human development are surprisingly common, but most often are confined to the placenta. These mutations provide clues to the normal developmental processes leading to a healthy placenta, despite these features being shared in common with cancer.

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.

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.

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.

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.

Nucleated red blood cells impact DNA methylation and expression analyses of cord blood hematopoietic cells.

BACKGROUND: Genome-wide DNA methylation (DNAm) studies have proven extremely useful to understand human hematopoiesis. Due to their active DNA content, nucleated red blood cells (nRBCs) contribute to epigenetic and transcriptomic studies derived from whole cord blood. Genomic studies of cord blood hematopoietic cells isolated by fluorescence-activated cell sorting (FACS) may be significantly altered by heterotopic interactions with nRBCs during conventional cell sorting. RESULTS: We report that cord blood T cells, and to a lesser extent monocytes and B cells, physically engage with nRBCs during FACS. These heterotopic interactions resulted in significant cross-contamination of genome-wide epigenetic and transcriptomic data. Formal exclusion of erythroid lineage-specific markers yielded DNAm profiles (measured by the Illumina 450K array) of cord blood CD4 and CD8 T lymphocytes, B lymphocytes, natural killer (NK) cells, granulocytes, monocytes, and nRBCs that were more consistent with expected hematopoietic lineage relationships. Additionally, we identified eight highly differentially methylated CpG sites in nRBCs (false detection rate <5 %, |Δß| >0.50) that can be used to detect nRBC contamination of purified hematopoietic cells or to assess the impact of nRBCs on whole cord blood DNAm profiles. Several of these erythroid markers are located in or near genes involved in erythropoiesis (ZFPM1, HDAC4) or immune function (MAP3K14, IFIT1B), reinforcing a possible immune regulatory role for nRBCs in early life. CONCLUSIONS: Heterotopic interactions between erythroid cells and white blood cells can result in contaminated cell populations if not properly excluded during cell sorting. Cord blood nRBCs have a distinct DNAm profile that can significantly skew epigenetic studies. Our findings have major implications for the design and interpretation of genome-wide epigenetic and transcriptomic studies using human cord blood.

Genome-wide DNA methylation identifies trophoblast invasion-related genes: Claudin-4 and Fucosyltransferase IV control mobility via altering matrix metalloproteinase activity.

Previously we showed that extravillous cytotrophoblast (EVT) outgrowth and migration on a collagen gel explant model were affected by exposure to decidual natural killer cells (dNK). This study investigates the molecular causes behind this phenomenon. Genome wide DNA methylation of exposed and unexposed EVT was assessed using the Illumina Infinium HumanMethylation450 BeadChip array (450 K array). We identified 444 differentially methylated CpG loci in dNK-treated EVT compared with medium control (P < 0.05). The genes associated with these loci had critical biological roles in cellular development, cellular growth and proliferation, cell signaling, cellular assembly and organization by Ingenuity Pathway Analysis (IPA). Furthermore, 23 mobility-related genes were identified by IPA from dNK-treated EVT. Among these genes, CLDN4 (encoding claudin-4) and FUT4 (encoding fucosyltransferase IV) were chosen for follow-up studies because of their biological relevance from research on tumor cells. The results showed that the mRNA and protein expressions of both CLDN4 and FUT4 in dNK-treated EVT were significantly reduced compared with control (P < 0.01 for both CLDN4 and FUT4 mRNA expression; P < 0.001 for CLDN4 and P < 0.01 for FUT4 protein expression), and were inversely correlated with DNA methylation. Knocking down CLDN4 and FUT4 by small interfering RNA reduced trophoblast invasion, possibly through the altered matrix metalloproteinase (MMP)-2 and/or MMP-9 expression and activity. Taken together, dNK alter EVT mobility at least partially in association with an alteration of DNA methylation profile. Hypermethylation of CLDN4 and FUT4 reduces protein expression. CLDN4 and FUT4 are representative genes that participate in modulating trophoblast mobility.

Quantification of cell-free DNA in normal and complicated pregnancies: overcoming biological and technical issues.

The characterization of cell-free DNA (cfDNA) originating from placental trophoblast in maternal plasma provides a powerful tool for non-invasive diagnosis of fetal and obstetrical complications. Due to its placental origin, the specific epigenetic features of this DNA (commonly known as cell-free fetal DNA) can be utilized in creating universal 'fetal' markers in maternal plasma, thus overcoming the limitations of gender- or rhesus-specific ones. The goal of this study was to compare the performance of relevant approaches and assays evaluating the amount of cfDNA in maternal plasma throughout gestation (7.2-39.5 weeks). Two fetal- or placental-specific duplex assays (RPP30/SRY and RASSF1A/ß-Actin) were applied using two technologies, real-time quantitative PCR (qPCR) and droplet digital PCR (ddPCR). Both methods revealed similar performance parameters within the studied dynamic range. Data obtained using qPCR and ddPCR for these assays were positively correlated (total cfDNA (RPP30): R = 0.57, p = 0.001/placental cfDNA (SRY): R = 0.85, p<0.0001; placental cfDNA (RASSF1A): R = 0.75, p<0.0001). There was a significant correlation in SRY and RASSF1A results measured with qPCR (R = 0.68, p = 0.013) and ddPCR (R = 0.56, p = 0.039). Different approaches also gave comparable results with regard to the correlation of the placental cfDNA concentration with gestational age and pathological outcome. We conclude that ddPCR is a practical approach, adaptable to existing qPCR assays and well suited for analysis of cell-free DNA in plasma. However, it may need further optimization to surpass the performance of qPCR.

Global analysis of DNA methylation changes during progression of oral cancer.

OBJECTIVES: Earlier studies involving a priori gene selection have identified promoter regions deregulated by DNA methylation changes in oral squamous cell cancers (OSCCs) and precancers. Interrogation of global DNA methylation patterns for such specimens has not been reported, though such analyses are needed to uncover novel molecular factors driving disease. MATERIALS AND METHODS: We evaluated global DNA methylation patterns for 30 biopsies obtained from 10 patients undergoing surgical removal of an OSCC or carcinoma in situ (CIS). From a disease field in each patient, we collected (i) dysplastic, (ii) CIS or OSCC, and (iii) adjacent normal biopsies. DNA isolated from each biopsy was profiled for methylation status using the Illumina HumanMethylation27K platform. RESULTS: Our data demonstrate that aberrant methylation of promoter CpG islands exists across oral precancer and OSCC genomes. Non-hierarchical clustering of all methylation data revealed distinct methylation patterns between the normal and the CIS/OSCC tissues (with results for dysplastic biopsies split between groups). Multiple genes exhibiting recurrent aberrant DNA methylation were found for both dysplastic and CIS/OSCC groups, and included enrichment for genes found in the WNT and MAPK signaling pathways. CONCLUSION: In identifying aberrant DNA methylation at the earliest stages of oral precancer and finding recurring epigenetic disruption of specific genes/pathways across our analyzed cohort, we see evidence that CpG methylation changes may play a role in oral cancer progression and that global DNA methylation analyses may have significant utility in wider studies that seek to derive biomarkers or potentially druggable targets to improve oral cancer outcomes.

Widespread DNA hypomethylation at gene enhancer regions in placentas associated with early-onset pre-eclampsia.

Pre-eclampsia is a serious complication of pregnancy that can affect both maternal and fetal outcomes. Early-onset pre-eclampsia (EOPET) is a severe form of pre-eclampsia that is associated with altered physiological characteristics and gene expression in the placenta. DNA methylation is a relatively stable epigenetic modification to DNA that can reflect gene expression, and can provide insight into the mechanisms underlying such expression changes. This case-control study focused on DNA methylation and gene expression of whole chorionic villi samples from 20 EOPET placentas and 20 gestational age-matched controls from pre-term births. DNA methylation was also assessed in placentas affected by late-onset pre-eclampsia (LOPET) and normotensive intrauterine growth restriction (nIUGR). The Illumina HumanMethylation450 BeadChip was used to assess DNA methylation at >480 000 cytosine-guanine dinucleotide (CpG) sites. The Illumina HT-12v4 Expression BeadChip was used to assess gene expression of >45 000 transcripts in a subset of cases and controls. DNA methylation analysis by pyrosequencing was used to follow-up the initial findings in four genes with a larger cohort of cases and controls, including nIUGR and LOPET placentas. Bioinformatic analysis was used to identify overrepresentation of gene ontology categories and transcription factor binding motifs. We identified 38 840 CpG sites with significant (false discovery rate <0.01) DNA methylation alterations in EOPET, of which 282 had >12.5% methylation difference compared with the controls. Significant sites were enriched at the enhancers and low CpG density regions of the associated genes and the majority (74.5%) of these sites were hypomethylated in EOPET. EOPET, but not associated clinical features, such as intrauterine growth restriction (IUGR), presented a distinct DNA methylation profile. CpG sites from four genes relevant to pre-eclampsia (INHBA, BHLHE40, SLC2A1 and ADAM12) showed different extent of changes in LOPET and nIUGR. Genome-wide expression in a subset of samples showed that some of the gene expression changes were negatively correlated with DNA methylation changes, particularly for genes that are responsible for angiogenesis (such as EPAS1 and FLT1). Results could be confounded by altered cell populations in abnormal placentas. Larger sample sizes are needed to fully address the possibility of sub-profiles of methylation within the EOPET cohort. Based on DNA methylation profiling, we conclude that there are widespread DNA methylation alterations in EOPET that may be associated with changes in placental function. This property may provide a useful tool for early screening of such placentas. This study identifies DNA methylation changes at many loci previously reported to have altered gene expression in EOPET placentas, as well as in novel biologically relevant genes we confirmed to be differentially expressed. These results may be useful for DNA- methylation-based non-invasive prenatal diagnosis of at-risk pregnancies.

DNA methylation profiling of placental villi from karyotypically normal miscarriage and recurrent miscarriage.

Miscarriage occurs in 15% of clinical pregnancies. Although chromosomal errors are observed in >50%, causes of karyotypically normal losses are poorly understood. DNA methylation undergoes reprogramming during development and must be appropriately set to maintain a healthy pregnancy. We hypothesize that aberrant DNA methylation may cause karyotypically normal miscarriage, particularly among women experiencing recurrent miscarriage (RM). DNA methylation in first-trimester chorionic villi was assessed in chromosomally normal miscarriages from women with RM (N = 33) or isolated miscarriage (M; N = 21) and elective terminations (TA; N = 16). Differentially methylated candidate loci were identified using the Illumina Infinium HumanMethylation27 BeadChip array. Follow-up bisulfite pyrosequencing at promoter regions showed an increase in methylation in M compared with TA at cytochrome P450, subfamily 1A, polypeptide 2 (CYP1A2; P = 0.002) and RM compared with TA at AXL receptor tyrosine kinase (P = 0.02), and a decrease in RM and M compared with TA at defensin ß 1 (DEFB1; P = 0.008). Gene ontology analysis showed an enrichment of imprinted genes (P = 9.53 × 10(-10)) and genes previously associated with RM (P = 9.51 × 10(-6)). An increase of outliers at seven imprinted loci was observed in RM (3.9%) compared with M (0%) and TA (0.9%) (P = 0.02), with increased average methylation at H19/IGF2 ICR1 in M samples (P < 0.0001). Altered DNA methylation in the placenta at specific loci, and global dysregulation in specific cases, may contribute to or be a consequence of poor placental function in karyotypically normal miscarriage.

The human placenta methylome.

Tissue-specific DNA methylation is found at promoters, enhancers, and CpG islands but also over larger genomic regions. In most human tissues, the vast majority of the genome is highly methylated (>70%). Recently, sequencing of bisulfite-treated DNA (MethylC-seq) has revealed large partially methylated domains (PMDs) in some human cell lines. PMDs cover up to 40% of the genome and are associated with gene repression and inactive chromatin marks. However, to date, only cultured cells and cancers have shown evidence for PMDs. Here, we performed MethylC-seq in full-term human placenta and demonstrate it is the first known normal tissue showing clear evidence of PMDs. We found that PMDs cover 37% of the placental genome, are stable throughout gestation and between individuals, and can be observed with lower sensitivity in Illumina 450K Infinium data. RNA-seq analysis confirmed that genes in PMDs are repressed in placenta. Using a hidden Markov model to map placental PMDs genome-wide and compare them to PMDs in other cell lines, we found that genes within placental PMDs have tissue-specific functions. For regulatory regions, methylation levels in promoter CpG islands are actually higher for genes within placental PMDs, despite the lower overall methylation of surrounding regions. Similar to PMDs, polycomb-regulated regions are hypomethylated but smaller and distinct from PMDs, with some being hypermethylated in placenta compared with other tissues. These results suggest that PMDs are a developmentally dynamic feature of the methylome that are relevant for understanding both normal development and cancer and may be of use as epigenetic biomarkers.

Beckwith-Wiedemann syndrome in sibs discordant for IC2 methylation.

Genetically heterogeneous imprinting disorders include Beckwith-Wiedemann syndrome (BWS) and multiple maternal hypomethylation syndrome (MMHS). Using DNA sequencing, quantitative PCR, SNuPE, pyrosequencing, and hybridization to the Illumina GoldenGate Methylation Cancer Panel 1 array, we characterized the genomic DNA of two brothers with BWS who were discordant for loss of methylation at several differentially methylated regions (DMR), including imprinting center 2 (IC2) on chromosome band 11p15.5, which is often hypomethylated in BWS. In keeping with MMHS, the elder child had hypomethylation of SGCE and PLAGL1 as well as of IC2, whereas the younger brother demonstrated no loss of methylation at these DMRs. Although this discordance is consistent with the observation that 15-20% of individuals with BWS do not have detectable genetic or epigenetic alterations of 11p15.5, this is the first report of familial recurrence of BWS with discordance for chromosomal 11p15.5 alterations. We hypothesize that this apparent discordance arises either from mosaicism precluding identification of IC2 hypomethylation in blood or buccal mucosa DNA of the younger child, or from hypomethylation at a site not interrogated by our molecular studies.