Genome wide DNA methylation analysis of alveolar capillary dysplasia lung tissue reveals aberrant methylation of genes involved in development including the FOXF1 locus
Clinical Epigenetics volume 13, Article number: 148 (2021)
Alveolar capillary dysplasia with or without misalignment of the pulmonary veins (ACD/MPV) is a lethal congenital lung disorder associated with a variety of heterozygous genomic alterations in the FOXF1 gene or its 60 kb enhancer. Cases without a genomic alteration in the FOXF1 locus have been described as well. The mechanisms responsible for FOXF1 haploinsufficiency and the cause of ACD/MPV in patients without a genomic FOXF1 variant are poorly understood, complicating the search for potential therapeutic targets for ACD/MPV. To investigate the contribution of aberrant DNA methylation, genome wide methylation patterns of ACD/MPV lung tissues were compared with methylation patterns of control lung tissues using the recently developed technique Methylated DNA sequencing (MeD-seq).
Eight ACD/MPV lung tissue samples and three control samples were sequenced and their mutual comparison resulted in identification of 319 differentially methylated regions (DMRs) genome wide, involving 115 protein coding genes. The potentially upregulated genes were significantly enriched in developmental signalling pathways, whereas potentially downregulated genes were mainly enriched in O-linked glycosylation. In patients with a large maternal deletion encompassing the 60 kb FOXF1 enhancer, DNA methylation patterns in this FOXF1 enhancer were not significantly different compared to controls. However, two hypermethylated regions were detected in the 60 kb FOXF1 enhancer of patients harbouring a FOXF1 point mutation. Lastly, a large hypermethylated region overlapping the first FOXF1 exon was found in one of the ACD/MPV patients without a known pathogenic FOXF1 variation.
This is the first study providing genome wide methylation data on lung tissue of ACD/MPV patients. DNA methylation analyses in the FOXF1 locus excludes maternal imprinting of the 60 kb FOXF1 enhancer. Hypermethylation at the 60 kb FOXF1 enhancer might contribute to FOXF1 haploinsufficiency caused by heterozygous mutations in the FOXF1 coding region. Interestingly, DNA methylation analyses of patients without a genomic FOXF1 variant suggest that abnormal hypermethylation of exon 1 might play a role in some ACD/MPV in patients.
Alveolar capillary dysplasia with or without misalignment of the pulmonary veins (ACD/MPV) is a lethal congenital lung disorder associated with heterozygous genomic alterations in the Forkhead Box F1 (FOXF1) gene locus. These genomic alterations vary from point mutations to large copy number variations (CNVs) but all seem to cause haploinsufficiency leading to ACD/MPV [1, 2]. The molecular mechanism responsible for the haploinsufficiency are poorly understood which complicates the search for potential therapeutic targets for ACD/MPV. Furthermore, in approximately 30% of ACD/MPV patients no genomic variants in the FOXF1 locus could be found, challenging the diagnostic process [3, 4].
The majority of ACD/MPV patients die within the first weeks of life due to insufficient gas exchange and therapy-resistant pulmonary hypertension. Histologically, ACD/MPV is characterized by reduced numbers of pulmonary capillaries, increased medial wall thickening in pulmonary arterioles and, in most cases, displaced pulmonary veins (reviewed in Slot et al. ). Genomic alterations associated with ACD/MPV are either small indel mutations in the FOXF1 gene or large deletions that include the 60 kb enhancer region [4,5,6]. This enhancer is located 250 kb upstream of the FOXF1 transcription start site and lies within the same topologically associated domain (TAD) as FOXF1 . Furthermore, this FOXF1 enhancer region contains GLI2 binding sites that upregulate FOXF1 expression upon Sonic hedgehog activation [4, 5, 7].
Heterozygous FOXF1 variants seem to have a dominant negative effect that cause haploinsufficiency of FOXF1, associated with ACD/MPV, and it has been suggested that FOXF1 is subjected to parental imprinting [1, 8,9,10]. One of the well-established mechanisms of parental imprinting is allele specific DNA methylation of gene promotors . However, several studies have shown that the FOXF1 promoter is not methylated in both normal and ACD/MPV lung tissue, excluding the possibility of allele specific silencing of FOXF1 through promotor methylation [10, 12, 13]. Furthermore, FOXF1 seems to be similarly transcribed from both alleles, which opposes imprinting of the FOXF1 transcription region itself [13, 14]. However, this does not completely exclude the 60 kb enhancer region to be subject to parental imprinting on the paternal allele (Additional file 3: Figure S 1), as suggested previously  based on the finding that all but one of the large deletions detected in ACD/MPV patients involved the 60 kb enhancer on the maternal chromosome [3, 10]. Furthermore, although very unlikely, parental imprinting of the FOXF1 enhancer on the maternal allele has never been excluded either. Up to now, studies investigating allele specific methylation in the 60 kb enhancer are limited and mainly use DNA isolated from blood or skin tissue which might not be representative for methylation patterns in lung tissue [15, 16].
Although most patients have a similar clinical course, atypical cases have been described that deviate in time of disease presentation and progression or reliance on life supportive care. So far, the phenotypical differences could not be correlated with the presence, type or location of the FOXF1 variant. Considering the complex genotype–phenotype correlation and the absence of genomic FOXF1 alterations in some of the patients, abnormal DNA methylation might contribute to the pathogenesis of ACD/MPV.
In 2018, we developed Methylated DNA Sequencing (MeD-seq) which enables genome wide DNA-methylation profiling at single-nucleotide resolution without the need for deep sequencing . Using this technique, this study aimed at detection of genome wide methylation patterns in lung tissue of eight ACD/MPV patients, with specific focus on the FOXF1 locus, to identify potential differences that might play a role in the pathogenesis of ACD/MPV.
Differential methylation of developmental genes in ACD/MPV lung tissue
Formalin-fixed and paraffin-embedded (FFPE) lung tissue samples of eight ACD/MPV patients and three age-matched controls were available to study lung specific methylation (Table 1). All ACD/MPV patients have been previously described and tested for FOXF1 alterations extensively . Three patients with a point mutation in the first FOXF1 exon (ACD-mut), three patients with a large deletion overlapping the FOXF1 enhancer on the maternal allele (ACD-del) and two patients without a known pathogenic variant in the FOXF1 locus (ACD-none) (Additional file 4: FigureS 2) were selected. Control patients were age-matched individuals with a medical condition unrelated to lung development or underlying pulmonary pathology.
Genome wide DNA methylation patterns of ACD/MPV lung tissues were compared with methylation patterns of control lung tissues. Differentially methylated regions (DMRs) located on the X- and Y-chromosome were removed because samples were not gender-matched, as well as DMRs with a fold change below 2. This resulted in 319 DMRs of which 184 were hypermethylated and 135 were hypomethylated in ACD/MPV lung samples (Fig. 1a). Approximately half of the DMRs (43%) were intergenic, the other 57% of DMRs fully or partly overlapped putative gene promoters (TSS) or gene bodies (Fig. 1b). Depending on the location, methylation acts differently on gene regulation. In general, promotor methylation is associated with gene silencing, whereas gene body methylation is mostly associated with gene activation [18,19,20,21,22,23,24,25]. Based on these assumptions, we labelled genes as potentially up- or downregulated in ACD/MPV lung tissue compared to control lung tissue. Exclusion of pseudogenes, long-non coding RNAs and duplicates resulted in 79 potentially upregulated and 36 potentially downregulated protein encoding genes (Fig. 1b, Additional file 1: table S 1). To study if these potentially up- and downregulated genes are specifically involved in certain biological pathways, gene enrichment analyses were performed using Metascape . Gene enrichment analysis of the potentially upregulated genes revealed 16 statistically significant enriched gene clusters of which the top cluster (named ‘embryonic morphogenesis’ (GO:0,048,598)) (Fig. 1c) included different signalling pathways involved in developmental processes such as ‘embryonic organ morphogenesis’ (GO:0,048,562) and ‘blood vessel morphogenesis’ (GO: 0,048,514) (Additional file 2: table S 2). The genes included in this top cluster were: COL4A2, GATA2, HOXA3, HOXB3, HOXB6, HOXD3, NOTCH1, TGFB1, SOCS3, KDM6B, TENM4, CCDC40, TIE1, NXN, DHRS3, RXRA. Gene enrichment analysis of potentially downregulated genes revealed three significantly overrepresented gene clusters, with the most significant pathway being ‘O-linked glycosylation’ which included MUC5AC, ADAMTS2 and GALNT15.
The FOXF1 enhancer is not maternally imprinted
Next, we focused on methylation patterns in the FOXF1 locus (Fig. 2a). Comparison of all ACD/MPV samples with controls revealed no significant DMRs in the FOXF1 gene, the 60 kb FOXF1 enhancer nor the 250 kb between these two regions (Table 2). In order to identify potential correlations between methylation patterns and patient-specific FOXF1 variants, we separated the ACD/MPV samples with known genomic FOXF1 variants into two groups. The first group (ACD-del) contained the three samples with a large deletion involving the FOXF1 enhancer on the maternal chromosome (Fig. 2a), the second group (ACD-mut) contained the three ACD/MPV patients with a point mutation in the first exon of FOXF1 (ACD-mut). Both groups were compared with control samples and analysed for significant DMRs in the FOXF1 locus. In addition, we separately compared DNA methylation patterns of the two ACD/MPV patients without a pathogenic FOXF1 variant (ACD-none1 and ACD-none2) with control samples.
We analysed DNA methylation patterns in ACD-del samples. Although methylated regions were present in both ACD-del samples and controls, no significant differences were detected between the groups (Table 2, Fig. 2b). Since one of the ACD-del patients (ACD-del3) harboured a deletion starting in the middle of the 60 kb enhancer, ACD-del samples ACD-del1 and ACD-del2 were also separately compared with control samples but again revealed no significant DMRs. These results indicate that the large maternal deletions in ACD-del samples have no effect on methylation patterns in the 60 kb FOXF1 enhancer. Szafranski and colleagues previously suggested that the 60 kb FOXF1 enhancer is subjected to allele specific imprinting through DNA methylation at the paternal allele . However, they did not study allele specific methylation and therefore could not exclude allele specific imprinting through DNA methylation at the maternal allele. If the 60 kb FOXF1 enhancer normally contains regions that are specifically methylated on the maternal allele, these would be lost in the ACD-del samples studied by us, resulting in significant DMRs between ACD-del and control samples (Additional file 5: Figure S 3A, left panel). Therefore, the absence of significant DMRs between ACD-del and control samples confirms that the 60 kb FOXF1 enhancer is devoid of allele specific methylation on the maternal allele. Consequently, this means that if this FOXF1 enhancer is subjected to allele specific imprinting through DNA methylation, this concerns the paternal allele (Additional file 5: Figure S3A, right panel).
Using bisulphite sequencing, Schulze and colleagues recently identified a paternally imprinted region in the 60 kb FOXF1 enhancer 10 kb upstream of LINC01082 (Hg38 chr16: 86,186,428–86,186,443) . Since our ACD-del and control samples carried an intact paternal allele, we investigated whether this region was methylated in our samples. We did not detect abundant CpG methylation in this region, even though MeD-seq would be able to detect methylation at half of the CpG sites in this region if methylated, questioning the validity of paternal imprinting in this region (Additional file 5: Figure S 3B).
Two hypermethylated regions in the 60 kb FOXF1 enhancer of ACD/MPV patients with FOXF1 point mutations
Next, we compared methylation patterns between ACD-mut and control samples and found seven DMRs distributed over the FOXF1 enhancer and the 250 kb region between the enhancer and FOXF1 (Table 2; Fig. 3a). The two DMRs located inside the 60 kb FOXF1 enhancer were both hypermethylated in ACD-mut samples compared to controls, but also compared to ACD-del and ACD-none samples (Fig. 3a, b). Interestingly, both DMRs were located within the 15 kb of the enhancer region which is thought to be the critical region that needs to be deleted in order to cause the typical ACD/MPV phenotype (Fig. 3c). Moreover, one of these DMRs (DMR2) encompasses the 250 bp that have been shown to physically interact with the FOXF1 promoter . Since ACD-mut samples do not contain single nucleotide polymorphisms (SNPs) in the identified DMRs, it is unclear whether only one or both parental alleles of ACD-mut samples are hypermethylated.
Hypermethylation of the first FOXF1 exon in lung tissue of an ACD/MPV patient without a genomic FOXF1 variant
In our previous study, we described four patients with typical severe ACD/MPV phenotypes in which we could not detect genomic FOXF1 alterations classified as pathogenic or likely pathogenic . To investigate whether aberrant DNA methylation could contribute to the phenotypes of these patients, we compared methylation patterns of the two patients from which we had lung tissue available (ACD-none) with methylation patterns of control lung samples (Table 2). As we only had two samples of this group, normalized read counts in the FOXF1 locus were manually reviewed without bioinformatical analysis. In contrast to the ACD-mut samples, we did not observe DMRs located in the FOXF1 enhancer region. However, we found a highly methylated region in sample ACD-none2 that covered exon 1 of FOXF1 (Fig. 4a). This hypermethylated region was not observed in any of the other lung samples tested in this study. Finally, since we had bowel and lymph node tissue of patient ACD-none1 available from autopsy, we compared the methylation pattern at the FOXF1 locus in bowel, lymph node and lung tissues to explore the presence of tissue specific methylation patterns in the FOXF1 locus. The DNA methylation patterns in bowel tissue showed no evident differences compared to the methylation pattern in lung tissue. However, lymph node tissue of ACD-none1 contained a highly methylated region covering exon 1 that was not present in lung tissue of ACD-none1 (Fig. 4b). As it was remarkable that this region was similarly methylated in the lung sample of ACD-none2 and not in any of the other lung tissues tested, the origin of each tissue block was re-confirmed by the department of pathology at the Erasmus MC.
The aim of the present study was to identify genome wide methylation patterns in lung tissue of ACD/MPV patients, with specific focus on the FOXF1 locus to identify potential differences that might play a role in the pathogenesis of ACD/MPV. Comparison of ACD/MPV lung tissues with control lung tissues resulted in 319 DMRs genome wide, possibly affecting gene regulation of 115 protein coding genes. Interestingly, pathway analysis showed that the potentially upregulated genes are mainly involved in developmental processes. This conforms the developmental delay that is indicated by histological aspects of ACD/MPV lungs, such as reduced numbers of alveolar capillaries and increased numbers of type 2 pneumocytes . For example, the most significantly enriched gene cluster of potentially upregulated genes included TIE1 and NOTCH1, which are implicated in vascular development [27,28,29,30,31,32]. Furthermore, this cluster included COL4A2 and TGFB1, which are involved in alveolar and airway branching morphogenesis respectively [33, 34]. SOCS3, an inhibitor of the JAK/STAT signalling pathway [35, 36] was among the potentially upregulated genes as well, in line with the previously demonstrated reduction of STAT3 protein in lung tissue of ACD/MPV patients . The most significantly enriched gene cluster of potentially downregulated genes included MUC5AC, ADAMTS2 and GALNT15 which are all involved in O-linked glycosylation, a post-translational process regulating protein function. Of these, MUC5AC is secreted by goblet cells of the lower respiratory tract where it contributes to the mucocilliary clearance as part of the innate immune system [38,39,40]. However, the relevance of downregulation of O-linked glycosylation in relation to ACD/MPV is unclear. In summary, the results of our genome-wide methylation analysis implicate up- and downregulation of a range of genes in ACD/MPV. Whether these genes are indeed dysregulated and what their specific role is in ACD/MPV, requires further investigation. Additionally, similar to the patients included in the current study, genetic testing of ACD/MPV patients has been mainly directed at the FOXF1 locus, potentially obscuring detection of alterations in other genes. Therefore, our findings could serve as a basis to expand genetic testing in ACD/MPV patients to study whether these genes harbour genomic alterations that contribute to ACD/MPV.
Through comparison of all included ACD/MPV lung samples together with control samples we did not find significant DMRs in the FOXF1 locus. Similarly, no significant DMRs were found in the FOXF1 locus when we compared the group of ACD-del samples with control samples. The three ACD-del samples included in this study harbour a large deletion (partly) encompassing the FOXF1 enhancer on the maternal allele. And thus, the absence of DMRs in comparison indicates that the methylation pattern of the unaffected paternal allele is similar to the combined methylation patterns of both alleles in control samples. In line with previous studies, these results indicate that if the FOXF1 enhancer is imprinted through DNA methylation, this involves the paternal allele (Additional file 4: Figure S 2A, right panel) [4, 6, 9, 16]. Although MeD-seq revealed multiple regions within the 60 kb enhancer that are methylated on the paternal allele of ACD-del samples, additional methylation studies and analyses are needed to define if these regions are paternally imprinted. For instance, a detailed SNP analysis on methylation data from multiple control samples could indicate whether certain regions at the 60 kb enhancer favour methylation at one allele. Using parental SNP data, one could identify on which parental allele the specific methylation occurs. Furthermore, bisulphite sequencing could be applied to identify paternally methylated regions using ACD-del samples. Since this technique also sequences unmethylated regions, it is possible to perform a quantitative comparison between the methylation in ACD-del samples harbouring a maternal deletion, and control samples. If specific regions in the 60 kb FOXF1 enhancer are paternally imprinted, then 100% of the reads of the ACD-del samples would be methylated, compared to 50% of the reads in control samples. Due to DNA degradation in the FFPE tissue blocks that have been stored for many years, we were unable to perform bisulphite sequencing with the current collected samples .
When we compared the three ACD-mut samples with control samples, two significant hypermethylated regions were detected in the FOXF1 enhancer, potentially interfering with normal gene regulation. Recently, it was shown that rare genomic variants in transcription factor binding sites can influence DNA methylation patterns, presumably leading to altered expression of nearby genes . Although the mechanisms are unclear yet, it has been proposed that mutations could cause altered transcription factor binding, leading to formation of protein complexes with the DNA methylation machinery and thereby, alter DNA methylation. To investigate whether FOXF1 point mutations influence methylation at the FOXF1 enhancer, it would be interesting to study the DNA methylation status in primary human lung cells before and after introduction of a site-specific mutation in the FOXF1 gene. Endothelial colony forming cells, for instance, have a high FOXF1 expression, are easily isolated from lung tissue, and may therefore, be a suitable cell type to investigate their methylation profile in the 60 kb enhancer before and after introduction of a FOXF1 mutation. However, it would be ideal to study multiple cell types so that cell type specific differences in methylation patterns can be recognized. One of the mechanisms by which FOXF1 enhancer methylation could interfere with FOXF1 regulation is by changing the interaction between the enhancer and the FOXF1 promoter region. Therefore, it would be an interesting next step to use the ACD-mut samples to perform assay for transposase-accessible chromatin with sequencing (ATAC-Seq)  or universal NicE-seq (UniNicE-seq) . This could clarify whether enhancer methylation alters the chromatin state in the TAD harbouring the FOXF1 enhancer and gene, leading to an altered interaction.
In one of the ACD-none samples (ACD-none2) we detected a highly methylated region completely overlapping exon 1 of FOXF1. According to the Eukaryotic Promoter Database (EPD), the FOXF1 core-promoter sequence ends approximately 40 bp downstream of the transcription start site (TSS) and 30 bp upstream of the start of exon 1 which means that the hypermethylated region does not overlap with the FOXF1 promoter . Currently, very little is known about the exact consequences of DNA methylation of exon 1 methylation on gene regulation. Therefore, it should be further investigated whether this hypermethylated region is associated with gene activation similar to other methylated regions overlapping gene bodies, or if this methylated region is close enough to the promoter to cause gene silencing, for instance by disrupting polymerase binding. Although only one lymph node sample was tested, the observation that exon 1 is hypermethylated in lymph tissue but not in the gastro-intestinal tissue or any of the lung tissues samples other than ACD-none2, could indicate that exon 1 methylation is indeed associated with reduced FOXF1 expression since lymph node tissue is also the only tested tissue that normally does not express FOXF1 . However, Szafranski and colleagues also found methylation of exon 1 in a normal lung sample . Therefore, expression analyses in relation to exon 1 methylation are necessary before any conclusions about the relevance of the hypermethylated region in ACD-none2 can be drawn. In our previous study we showed that patient ACD-none2 carried a duplication in the 3’ UTR of FOXF1 that was also present in the healthy father of the patient . Although this variant has been classified as likely benign according to the ACMG classification , we cannot exclude that this variant contributed to the ACD/MPV phenotype, for instance through altering DNA methylation as suggested by our findings in ACD-mut samples. As multiple ACD/MPV patients have been described without a known genomic FOXF1 variant, it would be interesting to investigate their methylation status to get a better indication of the contribution of abnormal FOXF1 methylation in the pathogenesis of ACD/MPV .
In this study we investigated methylation patterns that led to many new hypotheses about the role of DNA methylation in the pathogenesis of ACD/MPV. To fully comprehend these results, additional expression analysis would be required. However, a limitation of this study was the quality of FFPE tissue blocks which was poor and highly variable, complicating reliable expression analyses. In our experience, these type of blocks frequently have partially degraded RNA, hampering RNA in situ hybridisations or qPCR analyses. Additional studies with FFPE blocks that have been processed similarly and stored for less than 5 years, or with fresh or frozen lung tissue, are needed to confirm our hypotheses.
Because DNA methylation is tissue specific, we used lung tissues to study methylation patterns that could contribute to the lung abnormalities observed in ACD/MPV patients [15, 21]. However, lung tissue contains a heterogeneous cell population and thus, the DMRs identified are not cell type specific. In human and mice lung tissue, FOXF1 is mainly expressed in mesenchymal, smooth muscle and endothelial cells [47, 49,50,51]. To study whether aberrant methylation of the FOXF1 locus is specific for one of these cell types, it would be interesting to study DNA methylation patterns in sorted cells from fresh ACD/MPV lung tissues.
This study provides genome wide methylation data on lung tissue of ACD/MPV patients. A detailed investigation of the FOXF1 locus revealed hypermethylation of the FOXF1 enhancer in ACD/MPV patients harbouring a point mutation in the FOXF1 gene and abnormal hypermethylation of exon 1 in a patient without a genomic FOXF1 variant. These abnormal methylation patterns potentially change the regulation of FOXF1 contributing to the pathogenesis of ACD/MPV. This is the first study presenting FOXF1 specific methylation data in relation to different genomic FOXF1 variants associated with ACD/MPV and serves as important starting point for further research.
ACD/MPV and control FFPE lung tissue samples were collected from the Erasmus MC Sophia Children’s Hospital in Rotterdam, the VU University Medical Centre in Amsterdam and the Hospital for Sick Children in Toronto. All specimen were obtained as part of routine autopsy following informed consent of the parents/legal representatives or as part of a diagnostic procedure to diagnose ACD/MPV based on clinical suspicion. Samples were anonymized before they were subjected to MeD-seq. The research proposal was reviewed and approved by the Daily Board of the Medical Ethics Committee (METC) Erasmus University Medical Centre Rotterdam, The Netherlands.
DNA was isolated from FFPE lung tissue using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. Eight to ten 10 µm sections were pre-treated with xylene to remove the paraffin and with sodium thiocyanate to permeabilize the tissue. Incubation with (ATL) lysis buffer and proteinase K was extended to 36 h to allow complete lysis. DNA concentrations were measured with the Quant-iT Picogreen assay kit (ThermoFisher Scientific) according the manufacturer’s instructions.
MeD-seq sample preparation
DNA samples were prepared for MeD-seq as described previously . In brief, DNA samples were digested with LpnPI (New England Biolabs, Ipswhich, MA, USA) and resulted in fragments of 32 bp with the methylated cytosine in the centre. Fragments were either purified on 10% TBE gel before preparation or purified by Pippin system gel after preparation. The 32 bp DNA fragments were prepared for sequencing using a ThruPlex DNA-seq 96D kit (Takara Bio Inc, Kusatasu, Japan) according to manufacturer’s protocol. To include dual indexed barcodes, stem–loop adaptors were blunt-end ligated to repaired input DNA and amplified (4 + 10 cycles) using a high-fidelity DNA polymerase. Multiplexed samples were sequenced on Illumina HiSeq2500 systems for single reads of 50 bp according to the manufacturer’s instructions. Dual indexed samples were demultiplexed using bcl2fastq software (Illumina).
Data processing and analysis
MeD-seq data were processed and analysed with Python 2.7.5 using specifically created scripts as described previously . In short, before mapping of the reads to the Hg38 genome using bowtie 2.1.0., the raw FASTQ files were subjected to Illumina adaptor trimming and filtered for the presence of LpnPI restriction sites 13-17 bp from the 3’ or 5’ end. For visualization of the mapped reads, BAM files were generated using SAMtools. LpnPI site scores were used to produce read count scores for the transcription start sites (TSS) [1 kb before and 1 kb after], gene bodies (1 kb after the TSS until the transcription end site) and CpG islands. Gene and CpG island annotations were downloaded from UCSC (hg38). To detect DMRs between two data sets, genome wide read counts were compared using the Chi-Squared test. Significance was set at p < 0.05 and was called with a Bonferroni correction or FDR using the Benjamini–Hochberg procedure. In addition, a genome-wide sliding window was used to detect sequentially differentially methylated LpnPI sites. Statistical significance was called between LpnPI sites in predetermined groups using the Chi-square test with a Bonferroni correction. Neighbouring significantly called LpnPI sites were binned and reported. Overlap of genome wide detected DMRs was reported for TSS, CpGisland or gene body regions using the annotations of UCSC (Hg38). Gene enrichment analyses were carried out using Metascape . Metascape identifies significant enriched terms such as Gene Ontology (GO) terms and canonical pathways, and clusters them into a tree, based on Kappa-statistical similarities among their gene memberships .
Alveolar capillary dysplasia with or without misalignment of the pulmonary veins
Methylation DNA sequencing
Forkhead Box F1
Copy number variation
Formalin-fixed and paraffin-embedded
Differentially methylated region
Transcription start site
Stankiewicz P, Sen P, Bhatt SS, Storer M, Xia Z, Bejjani BA, et al. Genomic and genic deletions of the FOX gene cluster on 16q24.1 and inactivating mutations of FOXF1 cause alveolar capillary dysplasia and other malformations. Am J Hum Genet. 2009;84(6):780–91.
Slot E, Edel G, Cutz E, van Heijst A, Post M, Schnater M, et al. Alveolar capillary dysplasia with misalignment of the pulmonary veins: clinical, histological, and genetic aspects. Pulm Circ. 2018;8(3):2045894018795143.
Slot E, von der Thusen JH, van Heijst A, van Marion R, Magielsen F, Dubbink HJ, et al. Fast detection of FOXF1 variants in patients with alveolar capillary dysplasia with misalignment of pulmonary veins using targeted sequencing. Pediatr Res. 2021;89(3):518–525.
Szafranski P, Gambin T, Dharmadhikari AV, Akdemir KC, Jhangiani SN, Schuette J, et al. Pathogenetics of alveolar capillary dysplasia with misalignment of pulmonary veins. Hum Genet. 2016;135(5):569–86.
Szafranski P, Dharmadhikari AV, Wambach JA, Towe CT, White FV, Grady RM, et al. Two deletions overlapping a distant FOXF1 enhancer unravel the role of lncRNA LINC01081 in etiology of alveolar capillary dysplasia with misalignment of pulmonary veins. Am J Med Genet A. 2014;164A(8):2013–9.
Szafranski P, Herrera C, Proe LA, Coffman B, Kearney DL, Popek E, et al. Narrowing the FOXF1 distant enhancer region on 16q24.1 critical for ACDMPV. Clin Epigenetics. 2016;8:112.
Mahlapuu M, Enerback S, Carlsson P. Haploinsufficiency of the forkhead gene Foxf1, a target for sonic hedgehog signaling, causes lung and foregut malformations. Development. 2001;128(12):2397–406.
Dharmadhikari AV, Szafranski P, Kalinichenko VV, Stankiewicz P. Genomic and epigenetic complexity of the FOXF1 locus in 16q24.1: implications for development and disease. Curr Genomics. 2015;16(2):107–16.
Sen P, Gerychova R, Janku P, Jezova M, Valaskova I, Navarro C, et al. A familial case of alveolar capillary dysplasia with misalignment of pulmonary veins supports paternal imprinting of FOXF1 in human. Eur J Hum Genet. 2013;21(4):474–7.
Szafranski P, Dharmadhikari AV, Brosens E, Gurha P, Kolodziejska KE, Zhishuo O, et al. Small noncoding differentially methylated copy-number variants, including lncRNA genes, cause a lethal lung developmental disorder. Genome Res. 2013;23(1):23–33.
Deaton AM, Bird A. CpG islands and the regulation of transcription. Genes Dev. 2011;25(10):1010–22.
Lo PK, Lee JS, Liang X, Han L, Mori T, Fackler MJ, et al. Epigenetic inactivation of the potential tumor suppressor gene FOXF1 in breast cancer. Cancer Res. 2010;70(14):6047–58.
Alsina Casanova M, Monteagudo-Sanchez A, RodiguezGuerineau L, Court F, Gazquez Serrano I, Martorell L, et al. Maternal mutations of FOXF1 cause alveolar capillary dysplasia despite not being imprinted. Hum Mutat. 2017;38(6):615–20.
Szafranski P, Gambin T, Karolak JA, Popek E, Stankiewicz P. Lung-specific distant enhancer cis regulates expression of FOXF1 and lncRNA FENDRR. Hum Mutat. 2021;42(6):694–8.
Schultz MD, He Y, Whitaker JW, Hariharan M, Mukamel EA, Leung D, et al. Human body epigenome maps reveal noncanonical DNA methylation variation. Nature. 2015;523(7559):212–6.
Schulze KV, Szafranski P, Lesmana H, Hopkin RJ, Hamvas A, Wambach JA, et al. Novel parent-of-origin-specific differentially methylated loci on chromosome 16. Clin Epigenet. 2019;11(1):60.
Boers R, Boers J, de Hoon B, Kockx C, Ozgur Z, Molijn A, et al. Genome-wide DNA methylation profiling using the methylation-dependent restriction enzyme LpnPI. Genome Res. 2018;28(1):88–99.
Arechederra M, Daian F, Yim A, Bazai SK, Richelme S, Dono R, et al. Hypermethylation of gene body CpG islands predicts high dosage of functional oncogenes in liver cancer. Nat Commun. 2018;9(1):3164.
Bird A. DNA methylation patterns and epigenetic memory. Genes Dev. 2002;16(1):6–21.
Fan L, Zhang A, Deng P. LMX1B mRNA expression and its gene body CpG methylation are valuable prognostic biomarkers for laryngeal squamous cell carcinoma. Biomed Pharmacother. 2019;117:109174.
Greenberg MVC, Bourc’his D. The diverse roles of DNA methylation in mammalian development and disease. Nat Rev Mol Cell Biol. 2019;20(10):590–607.
Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.
Jones PA, Baylin SB. The epigenomics of cancer. Cell. 2007;128(4):683–92.
Ponnaluri VK, Ehrlich KC, Zhang G, Lacey M, Johnston D, Pradhan S, et al. Association of 5-hydroxymethylation and 5-methylation of DNA cytosine with tissue-specific gene expression. Epigenetics. 2017;12(2):123–38.
Wang YW, Ma X, Zhang YA, Wang MJ, Yatabe Y, Lam S, et al. ITPKA gene body methylation regulates gene expression and serves as an early diagnostic marker in lung and other cancers. J Thorac Oncol. 2016;11(9):1469–81.
Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523.
Partanen J, Puri MC, Schwartz L, Fischer KD, Bernstein A, Rossant J. Cell autonomous functions of the receptor tyrosine kinase TIE in a late phase of angiogenic capillary growth and endothelial cell survival during murine development. Development. 1996;122(10):3013–21.
Iljin K, Petrova TV, Veikkola T, Kumar V, Poutanen M, Alitalo K. A fluorescent Tie1 reporter allows monitoring of vascular development and endothelial cell isolation from transgenic mouse embryos. FASEB J. 2002;16(13):1764–74.
Augustin HG, Koh GY, Thurston G, Alitalo K. Control of vascular morphogenesis and homeostasis through the angiopoietin-Tie system. Nat Rev Mol Cell Biol. 2009;10(3):165–77.
Benedito R, Roca C, Sorensen I, Adams S, Gossler A, Fruttiger M, et al. The notch ligands Dll4 and Jagged1 have opposing effects on angiogenesis. Cell. 2009;137(6):1124–35.
Hasan SS, Tsaryk R, Lange M, Wisniewski L, Moore JC, Lawson ND, et al. Endothelial Notch signalling limits angiogenesis via control of artery formation. Nat Cell Biol. 2017;19(8):928–40.
Benedito R, Rocha SF, Woeste M, Zamykal M, Radtke F, Casanovas O, et al. Notch-dependent VEGFR3 upregulation allows angiogenesis without VEGF-VEGFR2 signalling. Nature. 2012;484(7392):110–4.
Loscertales M, Nicolaou F, Jeanne M, Longoni M, Gould DB, Sun Y, et al. Type IV collagen drives alveolar epithelial-endothelial association and the morphogenetic movements of septation. BMC Biol. 2016;14:59.
Noe N, Shim A, Millette K, Luo Y, Azhar M, Shi W, et al. Mesenchyme-specific deletion of Tgf-beta1 in the embryonic lung disrupts branching morphogenesis and induces lung hypoplasia. Lab Invest. 2019;99(9):1363–75.
Babon JJ, Kershaw NJ, Murphy JM, Varghese LN, Laktyushin A, Young SN, et al. Suppression of cytokine signaling by SOCS3: characterization of the mode of inhibition and the basis of its specificity. Immunity. 2012;36(2):239–50.
Carow B, Rottenberg ME. SOCS3, a Major Regulator of Infection and Inflammation. Front Immunol. 2014;5:58.
Pradhan A, Dunn A, Ustiyan V, Bolte C, Wang G, Whitsett JA, et al. The S52F FOXF1 Mutation Inhibits STAT3 Signaling and Causes Alveolar Capillary Dysplasia. Am J Respir Crit Care Med. 2019;200(8):1045–1056.
Holmen JM, Karlsson NG, Abdullah LH, Randell SH, Sheehan JK, Hansson GC, et al. Mucins and their O-Glycans from human bronchial epithelial cell cultures. Am J Physiol Lung Cell Mol Physiol. 2004;287(4):L824–34.
Rose MC, Voynow JA. Respiratory tract mucin genes and mucin glycoproteins in health and disease. Physiol Rev. 2006;86(1):245–78.
Ridley C, Thornton DJ. Mucins: the frontline defence of the lung. Biochem Soc Trans. 2018;46(5):1099–106.
Ludgate JL, Wright J, Stockwell PA, Morison IM, Eccles MR, Chatterjee A. A streamlined method for analysing genome-wide DNA methylation patterns from low amounts of FFPE DNA. BMC Med Genomics. 2017;10(1):54.
Martin-Trujillo A, Patel N, Richter F, Jadhav B, Garg P, Morton SU, et al. Rare genetic variation at transcription factor binding sites modulates local DNA methylation profiles. PLoS Genet. 2020;16(11):e1009189.
Sturtzel C, Lipnik K, Hofer-Warbinek R, Testori J, Ebner B, Seigner J, et al. FOXF1 Mediates Endothelial Progenitor Functions and Regulates Vascular Sprouting. Front Bioeng Biotechnol. 2018;6:76.
Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position. Nat Methods. 2013;10(12):1213–8.
Chin HG, Sun Z, Vishnu US, Hao P, Cejas P, Spracklin G, et al. Universal NicE-seq for high-resolution accessible chromatin profiling for formaldehyde-fixed and FFPE tissues. Clin Epigenetics. 2020;12(1):143.
Dreos R, Ambrosini G, Perier RC, Bucher P. The Eukaryotic Promoter Database: expansion of EPDnew and new promoter analysis tools. Nucleic Acids Res. 2015;43(Database issue):D92–6.
Uhlen M, Fagerberg L, Hallstrom BM, Lindskog C, Oksvold P, Mardinoglu A, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347(6220):1260419.
Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17(5):405–24.
Kalinichenko VV, Lim L, Stolz DB, Shin B, Rausa FM, Clark J, et al. Defects in pulmonary vasculature and perinatal lung hemorrhage in mice heterozygous null for the Forkhead Box f1 transcription factor. Dev Biol. 2001;235(2):489–506.
Mahlapuu M, Pelto-Huikko M, Aitola M, Enerback S, Carlsson P. FREAC-1 contains a cell-type-specific transcriptional activation domain and is expressed in epithelial-mesenchymal interfaces. Dev Biol. 1998;202(2):183–95.
Steiner LA, Getman M, Schiralli Lester GM, Iqbal MA, Katzman P, Szafranski P, et al. Disruption of normal patterns of FOXF1 expression in a lethal disorder of lung development. J Med Genet. 2020;57(5):296–300.
We thank the VU University Medical Center Amsterdam and the Hospital for Sick Children Toronto for their contribution to our sample collection. Further, we thank the Friends from Sophia Foundation for their support and, we thank the Royal Netherlands Academy of Arts and Sciences for the Ter Meulen Grant.
Funding was provided by Stichting Vrienden van het Sophia (Grant No. 2017-010/WO) and Koninklijke Nederlandse Akademie van Wetenschappen (Grant No. 1453/WetF/TMB17/329).
The authors declare no conflict of interest or financial interests except for RB, JB, WVI and JG, who report being shareholder in Methylomics B.V., a commercial company that applies MeD-seq to develop methylation markers for cancer staging.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. Table S1. Lists of potentially up- and downregulated genes based on methylation status and TSS/GB overlap TSS: promotor, GB: Gene body. ↑: hypermethylated in ACD/MPV samples, ↓: hypomethylated in ACD/MPV samples.
. Table S 2. Results of gene ontology enrichment analysis using Metascape (26).
. Figure S 1. Simplistic illustration of the two models for parental imprinting of the 60kb FOXF1 enhancer Based on the suggested models by Szafranski et al. (4). Continuous arrow: full function of the FOXF1 enhancer. Dotted arrow: reduced function of the FOXF1 enhancer.
. Figure S 2. Haematoxylin and eosin staining of lung tissues of patients ACD-none1 and ACD-none2 demonstrating the main characteristics of ACD/MPV (2). Arrow: thickened alveolar wall. V: misaligned pulmonary vein. B: bronchiole. Scale bar: 200µm.
. Figure S 3. A: Models of maternal and paternal imprinting of the 60kb FOXF1 enhancers through DNA methylation. The absence of DMRs in ACD-del samples makes maternal imprinting of the FOXF1 enhancer unlikely and supports a paternal imprinting model. B: Methylation status of CpG sites of previously suggested paternally methylated region (rectangular frame) (16) in ACD-del and control samples. Read counts of each individual samples is demonstrated with IGV viewer (red: ACD-del, grey: control).
About this article
Cite this article
Slot, E., Boers, R., Boers, J. et al. Genome wide DNA methylation analysis of alveolar capillary dysplasia lung tissue reveals aberrant methylation of genes involved in development including the FOXF1 locus. Clin Epigenet 13, 148 (2021). https://doi.org/10.1186/s13148-021-01134-1
- DNA methylation
- Lung development