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

Comparing the Catalytic and Structural Characteristics of a 'Short' Unspecific Peroxygenase (UPO) Expressed in Pichia†pastoris and Escherichia†coli.

Unspecific peroxygenases (UPOs) have emerged as valuable tools for the oxygenation of non-activated carbon atoms, as they exhibit high turnovers, good stability and depend only on hydrogen peroxide as the external oxidant for activity. However, the isolation of UPOs from their natural fungal sources remains a barrier to wider application. We have cloned the gene encoding an 'artificial' peroxygenase (artUPO), close in sequence to the 'short' UPO from Marasmius rotula (MroUPO), and expressed it in both the yeast Pichia pastoris and E. coli to compare the catalytic and structural characteristics of the enzymes produced in each system. Catalytic efficiency for the UPO substrate 5-nitro-1,3-benzodioxole (NBD) was largely the same for both enzymes, and the structures also revealed few differences apart from the expected glycosylation of the yeast enzyme. However, the glycosylated enzyme displayed greater stability, as determined by nano differential scanning fluorimetry (nano-DSF) measurements. Interestingly, while artUPO hydroxylated ethylbenzene derivatives to give the (R)-alcohols, also given by a variant of the 'long' UPO from Agrocybe aegerita (AaeUPO), it gave the opposite (S)-series of sulfoxide products from a range of sulfide substrates, broadening the scope for application of the enzymes. The structures of artUPO reveal substantial differences to that of AaeUPO, and provide a platform for investigating the distinctive activity of this and related'short' UPOs.

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

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.

Profiling the small non-coding RNA transcriptome of the human placenta.

Proper functioning of the human placenta is critical for maternal and fetal health. While microRNAs (miRNAs) are known to impact placental gene expression, the effects of other small non-coding RNAs (sncRNAs) on the placental transcriptome are not well-established, and are emerging topics in the study of environmental influence on fetal development and reproductive health. Here, we assembled a cohort of 30 placental chorionic villi samples of varying gestational ages (M ± SD = 23.7 ± 11.3 weeks) to delineate the human placental sncRNA transcriptome through small RNA sequence analysis. We observed expression of 1544 sncRNAs, which include 48 miRNAs previously unannotated in humans. Additionally, 18,003 miRNA variants (isomiRs) were identified from the 654 observed miRNA species. This characterization of the term and pre-term placental sncRNA transcriptomes provides data fundamental to future investigations of their regulatory functions in the human placenta, and the baseline expression pattern needed for identifying changes in response to environmental factors, or under disease conditions.

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.

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.

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.

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.

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.