X-linked DNAme variation in placental cells is driven by cell type and sex.
To identify the primary variables associated with X-linked DNAme variation in the cell-type specific data from GSE159526, we first performed principal component analysis (PCA) on DNAme at 14,766 X-chromosome probes in all 94 samples. We tested for association between DNAme variation described by the principal components (PC scores) and sample variables (via linear models of the form PC ~ variables of interest) (Figure 1A, 1B). Cell type of the individual DNAme sample was strongly associated with the top 4 PCs (all nominal p-values < 0.001), which together explained > 70% of the X-chromosome DNAme variation in the dataset (PC1: 30%, PC2: 23%, PC3: 15%, PC4: 6%). Sex of the sample (XX or XY) was associated with PC1 and PC3 (p-value < 0.001). Although PC1 primarily separated samples by sex (Figure 1C), there was further subdivision by cell type within each sex along PC1, with XX stromal and endothelial cells clustering more closely to XY cells, while XX cytotrophoblasts and villi fell furthest away from XY samples along PC1. PC2 predominantly separated Hofbauer cells from all other cell types. Inferred ancestry probabilities and the technical chip variable (Sentrix ID) had only weak or non-significant associations with any of the top 10 PCs (Fig 1A, 1B).
To further visualize the relationship between sex and cell type, we performed hierarchical clustering on the same 14,766 X-chromosome probes (Figure 1D). This revealed three major clusters based on cell type: (i) cytotrophoblasts and villi, (ii) Hofbauer cells, and (iii) endothelial and stromal cells. Trophoblasts are the predominant cell type in whole chorionic villus tissue, which explains the clustering of cytotrophoblasts with villi in this analysis. Samples further subdivided by sex within each of the three main cell type clusters. Together with the PCA results, these hierarchical clustering results showed that (i) both sex and cell type are major drivers of X-chromosome DNAme in this cohort, and (ii) cell-specific patterns of DNAme are similar between cells with shared developmental origins (i.e. endothelial and stromal cells).
The distribution of X-linked DNAme of placental cells shows three distinct patterns.
To better understand what differentiates cell-type specific patterns of X-chromosome DNAme, we compared the distributions of X-chromosome DNAme across all cell types, separately in each sex (Figure 2A, 2B, Supplementary table 2). In XX cells, three general distributions of DNAme were observed. Hofbauer cells displayed a trimodal distribution, as is typical for somatic cells: a low methylated peak (0.2 ≤ β), a high methylated peak (β < 0.8), and a distinct peak of intermediate DNAme (0.2 < β ≤ 0.8). The intermediate DNAme peak reflects allele-specific DNAme associated with XCI, where promoters are fully methylated on the Xi and fully unmethylated on the Xa [43]. By contrast, cytotrophoblast and whole chorionic villi showed relatively few highly-methylated sites and lacked distinct intermediate methylation peaks. Both endothelial and stromal cells showed a distinct peak of low DNAme and a smaller peak of high DNAme, but had no clear intermediate DNAme peak, similar to the distribution observed amongst XY samples across all cell types.
DNAme of X-linked promoters show few sex differences in endothelial and stromal cells.
In somatic tissues, most X-chromosome gene promoters are roughly 50% methylated (β » 0.5) in XX cells, and unmethylated in XY cells [44]. Although whole chorionic villi show lower DNAme at X-linked promoters relative to somatic tissues (XX), genes subject to XCI in placenta do tend to show relatively higher levels of X-linked promoter DNAme in XX versus XY samples [18]. To compare the level of X promoter DNAme amongst the different cell types we calculated the sex difference in DNAme (|∆β| = XX - XY) at the 1,393 CpGs in X-linked promoter-associated CGI, based on CGI regions defined by Bizet et al. (2022) (Figure 2C). Sex differences in DNAme (i.e. DNAme |∆β| > 0.1) were observed at most X-linked CGIs in Hofbauer cells (76%) and at many CGI loci in cytotrophoblasts (45%). However, in endothelial and stromal cells, very few X-linked CGIs (3% and 7% respectively) had a sex difference in DNAme of |∆β| > 0.1, consistent with low/absent DNAme in both sexes at X-linked CGI promoters in endothelial and stromal cells. The lack of X-linked promoter DNAme in endothelial and stromal cells is also illustrated by a high correlation of X-linked DNAme between these cell types in both XX and XY samples (Supplementary Figure 2A, 2B), and our PCA showing close clustering of these cell types (Figure 1C).
In somatic cells, unlike promoters, gene bodies and intergenic regions on the X-chromosome show lower DNAme in XX compared to XY tissues [45]. We wanted to determine if the same patterns can be found in placental cells. To further evaluate DNAme by X-chromosome genomic region, we identified sex-differentially methylated CpGs (sex-DMCs) in each cell type using linear regression models with thresholds of |∆β| > |0.1| and FDR < 0.05) (Supplementary Table 3). As expected, most (45 - 75%) of the X promoter DMCs had higher XX relative to XY DNAme in Hofbauer cells, cytotrophoblasts and chorionic villi, but not in endothelial or stromal cells where few loci (4% or 7%, respectively) had sex-differential DNAme (Figure 2D). However, all cell types showed higher DNAme in XY relative to XX cells at DMCS located in X-linked enhancers, intergenic regions, and gene bodies, with cytotrophoblast showing the greatest number of sex-differentially methylated sites.
DNAme profile of placental endothelial cells differs from umbilical cord endothelial cells
To further characterize the unique X-linked DNAme patterns we observed in placental endothelial and stromal cells, we compared them to endothelial cells derived from other gestational tissues to determine whether we were detecting a specific endothelial X-linked DNAme signature. The endothelial cells evaluated in this work were FACS-isolated and were derived from placental microvessels within the chorionic villi. Blood flows between these placental microvessels and larger arteries and veins within the placental chorionic plate; these larger vessels then connect to the fetal vasculature via the umbilical cord. While, the placental villous endothelial cells derive from extraembryonic mesoderm, endothelial cells within the fetal compartment (fetal vessels) like the umbilical cord vessels, derive from embryonic precursors [46]. We thus sought to determine if low X-linked promoter DNAme was a shared property of all endothelial cells including those derived from the umbilical cord and large vessels of the placenta. To evaluate this, we compared DNAme patterns of our placental microvascular endothelial cells (henceforth called “pME”), to public Illumina HumanMethylation 450K data derived from (i) cultured human placental arterial and venous endothelial cells (pAE/pVE) obtained from the chorionic plate (GSE106099) (npAE = 12, XX = 7, XY = 5 / npVE = 17, XX = 9, XY = 8), and EPIC data derived from (ii) cultured human umbilical venous endothelial cells (uVE) (GSE144804) (nuVE = 36, XX = 22, XY = 14) (Supplementary table 4).
PCA and hierarchical clustering on the 8,012 X-chromosome CpGs common to all endothelial cell datasets (total sample n = 65) demonstrated DNAme differences by both sex and endothelial cell sampling location (Figure 3A-3C). PC1 separated samples by sex and cell type, with the greatest separation between XX and XY cells observed in uVEs and least separation in pAEs (Figure 3B). Hierarchical clustering also showed separation by both sex and cell type, although XX pVE clustered with XY samples (Figure 3C). The X-chromosome DNAme distributions also showed similar trends (Figure 3D), with XX uVEs showing a distinct intermediate methylated peak (57% at 0.2 < ∆β ≤ 0.8) characteristic of XX somatic cells, while in pAEs and pVEs this peak was largely absent and a large portion of X-linked CpGs had low DNAme (41%, and 57% respectively at ∆β ≤ 0.2), similar to what we observed in pMEs. In contrast, X-chromosome DNAme in XY samples was similar across all endothelial cell types.
Most promoter-CGI CpGs showed a difference in DNAme by sex (∆β > |0.1|) in uVEs (76%), but not pVEs (9%), while pAEs (44%) showed an intermediate result (Figure 3E). We next tested the cell-type correlation in DNAme at all X-chromosome CpGs (Supplementary Figure 3A, 3B). The strongest correlation was observed between XX pVE and pME (r = 0.78), and uVE and pAE (r = 0.75). The weakest correlation was observed between pME and uVE samples. In other words, placental microvessels (pME) studied in our first set of analyses were most similar to placental venous endothelial cells (pVE), while the umbilical cord vein (uVE) was most similar to placental artery (pAE). These results overall suggest considerable heterogeneity in the developmental origin of endothelial cells amongst these tissues, as inter-cell correlations in X DNAme between cells of similar origin is typically over 0.90 [47].
Y-chromosome DNAme varies by cell type.
Like the X, the Y-chromosome is also under-studied in epigenome-wide DNAme analyses. As the Y contains few genes, most of which function in the testes, we did not anticipate many DNAme differences by cell type in placenta. Nonetheless, PCA and hierarchical clustering of XY samples based on Y-chromosome DNAme (nCpGs = 293) showed three distinct clusters by cell type, parallelling those observed for the X-chromosome (Figure 4A-4C). In linear models comparing the average Y-chromosome DNAme of each cell type to the average of all other cell types, we identified many cell-influenced Y-chromosome DMCs (78, 84, 105, and 116 DMCs in endothelial, stromal, Hofbuaer cells and cytotrophoblasts, respectively) (Supplementary Figure 4). These cell-influenced Y DMCs overlapped 13 genes (see Supplementary Table 5). Five of these 13 genes (DDX3Y, EIF1AY, RPS4Y1, USP9Y, and ZFY), were previously reported to be expressed in XY term placentas [5]. These genes all possess X-linked homologs, which are also expressed in XX term placentas [5]. To confirm that Y-chromosome DNAme attributed to these X-Y homologs in our dataset was not arising from DNAme array probes cross-hybridization to their X-linked pairs, we performed a Command-line nucleotide BLAST (blastn) for short sequences on the 50-nucleotide probe sequences (probeSeqA/B) of significant Y-linked cell-DMCs of all cell types to exclude any possibility of cross-hybridization on the X chromosome due to its sequence between X and Y chromosomes. None of the Y DMC probes filtered by the selected criteria matched, suggesting that our results mostly reflect true Y-chromosome DNA methylation patterns.
Sex-influenced autosomal DNA methylation differs by cell type
We previously reported 145 CpGs on the autosomes that show sex-influenced DNAme in whole chorionic villi [28]. To determine if these sex-influenced autosomal CpGs were consistent across placental cell types we performed PCA using these 145 sex-influenced CpGs in the sorted placental cells. In scatterplots of PC1 versus PC2, samples predominantly separated by cell type and not sex; only the cytotrophoblasts/chorionic villi separated by sex (Figure 5A). These results likely reflect that the 145 CpGs were originally identified as having sex differences in data from bulk villi, and may not show a sex difference in other cell types.
Although our samples sizes were small, we wanted to evaluate the existence of sex-influenced autosomal DNAme in the individual placental cell types. By comparing DNAme in each cell type by sex using linear models, we identified multiple significant sex-influenced DMCs in endothelial cells (nCpGs = 35 DMCs at FDR < 0.05 and Δβ > |0.1|) and whole chorionic villi (nCpGs = 7 DMCs at FDR < 0.05 and Δβ > |0.1|) (Figure 5B, Supplementary figure 5), but none in the other cell types. Of the 35 endothelial sex-associated DMCs, the majority (89 %) had higher DNAme in XY as compared to XX cells, and were located in promoter or enhancer regions. Among the endothelial-associated sex-DMCs were CpGs in genes including LDB3, INHBB, NSD1, RAB7A, and ZNF300; of these genes LDB3 and ZNF300 had 2 or more DMCs each (nCpGs = 2/4, respectively).
As we were likely underpowered to detect sex differences in all cell types, to evaluate whether the sex differences in DNAme identified in endothelial cells were truly cell-specific, we performed PCA on the endothelial sex-associated DMCs in all placental cell types. The scatterplot of PC1 versus PC2 (Figure 5C) showed separation by cell type on PC1, with separation by sex observed in all cell types along PC2. Not surprisingly, the sex difference in endothelial cells was greatest, however, these results suggest that some of the identified DMCs show sex differences across multiple cell types, and that we are likely underpowered and missing significance in detecting sex-associated DMCs in other sorted placental cell types.
Finally, considering sex differences in DNAme at specific genes, the ZNF300 gene was previously reported to be DM by sex in Inkster et al., (2021), and was shown to be associated with placental morphology and development [48]. In endothelial cells, we identified 4 sex-DMCs at ZNF300 (DNAme XY > XX in all cell types except Hofbauer cells), in the same promoter region reported to be sex-DMCs in Inkster et al., (2019). Among the other sex-DMCs we identified in endothelial cells, several were associated with NSD1 and LDB3: 2 DMCs were observed in NSD1 (Nuclear Receptor Binding SET Domain Protein 1) and LDB3 (LIM domain binding 3), which plays a role as histone lysine methyltransferase and generate proteins maintaining the stability of the muscle structure, with higher DNAme in XX than XY cells (both endothelial and cytotrophoblast). Figures exemplifying the sex-differential DNAme patterns at these genes are shown in Supplementary Figure 6.