Genome-wide multilocus imprinting disturbance analysis in Temple syndrome and Kagami-Ogata syndrome.

PURPOSE: Recent studies have identified multilocus imprinting disturbances (MLIDs) in a subset of patients with imprinting diseases (IDs) caused by epimutations. We examined MLIDs in patients with Temple syndrome (TS14) and Kagami-Ogata syndrome (KOS14).
METHODS: We studied four TS14 patients (patients 1-4) and five KOS14 patients (patients 5-9) with epimutations. We performed HumanMethylation450 BeadChip (HM450k) analysis for 43 differentially methylated regions (DMRs) (753 CpG sites) and pyrosequencing for 12 DMRs (62 CpG sites) using leukocyte genomic DNA (Leu-gDNA) of patients 1-9, and performed HM450k analysis for 43 DMRs (a slightly different set of 753 CpG sites) using buccal cell gDNA (Buc-gDNA) of patients 1, 3, and 4. We also performed mutation analysis for six causative and candidate genes for MLIDs and quantitative expression analysis using immortalized lymphocytes in MLID-positive patients.
RESULTS: Methylation analysis showed hypermethylated ZDBF2-DMR and ZNF597/NAA60-DMR, hypomethylated ZNF597-DMR in both Leu-gDNA and Buc-gDNA, and hypomethylated PPIEL-DMR in Buc-gDNA of patient 1, and hypermethylated GNAS-A/B-DMR in Leu-gDNA of patient 3. No mutations were detected in the six genes for MLIDs. Expression patterns of ZDBF2, ZNF597, and GNAS-A/B were consistent with the identified MLIDs.
CONCLUSION: This study indicates the presence of MLIDs in TS14 patients but not in KOS14 patients.Genet Med 19 4, 476-482.

Introduction

The human chromosome 14q32.2 imprinted region harbors paternally and maternally expressed genes, together with the germ line–derived primary DLK1-MEG3 intergenic differentially methylated region (IG-DMR) and the postfertilization-derived secondary MEG3-DMR that function as the imprinting control centers in the placenta and the body, respectively.[1] Consistent with this, maternal uniparental disomy 14 (UPD(14)mat) and loss-of-methylation (LOM)-type epimutations and microdeletions affecting the paternally derived imprinted region lead to Temple syndrome (TS14, OMIM 616222), which is associated with pre- and postnatal growth deficiency, hypotonia, small hands, and precocious puberty.[2] Similarly, paternal uniparental disomy 14 (UPD(14)pat) and gain-of-methylation (GOM)-type epimutations and microdeletions involving the maternally inherited imprinted region result in Kagami-Ogata syndrome (KOS14, OMIM 608149), which is characterized by unique facial features; a small, bell-shaped thorax with a coat-hanger appearance of the ribs; abdominal-wall defects; placentomegaly; and polyhydramnios.[1]

Recent studies have identified multilocus imprinting disturbances (MLIDs) in a subset of patients with imprinting diseases (IDs) caused by epimutations.[3] Indeed, although the frequency of MLIDs is variable among studies, MLIDs have been detected in approximately 50% of patients with transient neonatal diabetes mellitus caused by LOM of the PLAGL1-DMR, approximately 12% of patients with Beckwith-Wiedemann syndrome caused by LOM of the Kv-DMR or GOM of the H19-DMR, 8–10% of patients with Silver-Russell syndrome (SRS) caused by LOM of the H19-DMR, and 8–10% of patients with pseudohypoparathyroidism type 1b caused by LOM of the GNAS-A/B-DMR, as well as in a single patient with Angelman syndrome caused by LOM of the SNRPN-DMR.[3] In addition, although ID-related DMRs have usually been examined to detect MLIDs by conventional methylation analyses including pyrosequencing, genome-wide MLID analysis has become possible by array-based methods.[4] The underlying factors for the MLIDs remain to be clarified in most patients, but the causative genes for MLIDs (ZFP57, NLRP2, and NLRP5), as well as candidate genes for MLIDs (NLRP7, KHDC3L, and TRIM28), have been identified in a few patients with MLIDs.[4,5]

MLIDs have been poorly studied in patients with TS14 and KOS14, primarily because of their rarity. To our knowledge, TS14 has been identified in 58 patients, including 10 patients with epimutations of the IG-DMR and the MEG3-DMR, and KOS14 has been identified in 57 patients, including 7 patients with epimutations of the two DMRs.[1,2,6,7,8,9] Thus, we examined MLIDs in TS14 and KOS14 caused by epimutations.

Materials and Methods

Ethical approval

This study was approved by the Institutional Review Board Committee at the National Center for Child Health and Development and was performed after obtaining written informed consent to publish the clinical and molecular information.

Patients

This study consisted of four TS14 patients with LOM-type epimutations (patients 1–4) and five KOS14 patients with GOM-type epimutations (patients 5–9) (). LOMs/GOMs were confirmed by the conventional methylation analysis for IG-DMR and MEG3-DMR using leukocyte genomic DNA (Leu-gDNA). UPD(14)mat/pat was ruled out by microsatellite analysis using parental Leu-gDNA, and microdeletions were excluded by fluorescence in situ hybridization analysis of the two DMRs and a custom-built array comparative genomic hybridization analysis for the 14q32.2 imprinted region.[1,10]

Detailed phenotypes are shown in Supplementary Tables S1 and S2 online. In brief, patient 1 exhibited Prader-Willi syndrome (PWS)-like phenotype in the absence of GOM at the SNRPN-DMR, and patients 2–4 manifested SRS-compatible phenotype in the absence of LOM at the H19-DMR and UPD(7)mat. Prenatal growth failure was absent from patient 1 and present in patients 2–4, whereas postnatal growth failure was present in patients 1–4. KOS14 patients were identified by the typical phenotype.

Molecular studies

We conducted methylation, mutation, and expression analyses. Detailed methods are described in the Supplementary Methods online. In brief, methylation analysis was performed by array-based HumanMethylation450 BeadChip (HM450k) (Illumina, San Diego, CA) analysis for 753 CpG sites on 43 DMRs using Leu-gDNA of patients 1–9 and for a slightly different set of 753 CpG sites on 43 DMRs using epithelial buccal cell gDNA (Buc-gDNA) of patients 1, 3, and 4 (Supplementary Table S3 online), and by pyrosequencing analysis for 50 CpG sites on nine ID-related DMRs and for 12 CpG sites on three DMRs that were assessed to be abnormally methylated by the HM450k analysis () using Leu-gDNA of patients 1–9 (for primer information, see Supplementary Table S4 online). Mutation analysis was performed for the causative and candidate genes for MLIDs (ZFP57, NLRP2, NLRP7, KHDC3L, NLRP5, and TRIM28)[4,5] by whole-exome sequencing using Leu-gDNA of MLID-positive patients and their mothers. Expression analysis was performed by quantitative real-time polymerase chain reaction for imprinted genes regulated by the abnormally methylated DMRs using mRNA samples extracted from immortalized lymphocytes derived from MLID-positive patients. We utilized the nomenclature for DMRs and imprinted genes proposed by the European Network of Human Congenital Imprinting Disorders (http://www.imprinting-disorders.eu/).

Results

Methylation analysis

MLID was identified in patients 1 and 3 with TS14 (, Supplementary Table S3 online, and ) (actual methylation levels of each patient and the average ± 3 SD values of control subjects obtained by HM450k analysis are shown in Supplementary Tables S5 and S6 online). In patient 1, GOM at ZDBF2-DMR, LOM at ZNF597-DMR, and GOM at ZNF597/NAA60-DMR were identified by HM450k analysis for Leu-gDNA and Buc-gDNA, with more obvious aberrant methylation patterns in Buc-gDNA than in Leu-gDNA, and by pyrosequencing analysis for Leu-gDNA. In addition, LOM at PPIEL-DMR was found only in Buc-gDNA by HM450k analysis. In patient 3, GOM at GNAS-A/B-DMR was detected only in Leu-gDNA by HM450k analysis, but not by pyrosequencing analysis.

No other MLID was identified. Although pyrosequencing analysis showed slight GOMs at the #17 CpG site within PEG10-DMR in patient 2 and the #27 CpG site within MEST-DMR in patient 3 (), other CpG sites within these DMRs were normal. Similarly, although HM450k analysis revealed several abnormally methylated CpG sites (, Supplementary Table S3 online), they accounted for <20% of CpGs within corresponding DMRs.

Of the 11 CpG sites examined by HM450k and pyrosequencing analyses, discordant methylation data were obtained for the #17 CpG site in patient 2 and for the #44 and #45 CpG sites in patient 3. Furthermore, HM450k analysis indicated (i) more severe LOMs at MEG3-DMR of Leu-gDNA in patient 1 than in patients 2–4; (ii) more severe LOMs at MEG3-DMR in Buc-gDNA than in Leu-gDNA of patients 3 and 4; (iii) invariable GOMs and LOMs at MEG8-DMR in TS14 and KOS14 patients, respectively; (iv) normal methylation patterns at H19-DMR and SNRPN-DMR in Buc-gDNA of patients 1, 3, and 4, as well as in Leu-gDNA of patients 1–9; and (v) hypermethylated GPR1-AS-DMR in Leu-gDNA of patients 1–9 and control subjects (Supplementary Table S5 online).

Mutation analysis

No pathogenic mutation was identified in the six genes of MLID-positive patients 1 and 3 and their mothers.

Expression analysis

Expression analysis revealed significantly increased ZDBF2 expression and decreased ZNF597 expressions in patient 1 and significantly decreased GNAS-A/B expression in patient 3 ().

Discussion

We performed genome-wide MLID analysis for four TS14 and five KOS14 patients with epimutations. MLIDs were found in Leu-gDNA and Buc-gDNA of patient 1 and in Leu-gDNA of patient 3 with TS14. In particular, the MLIDs identified in Leu-gDNA (the GOMs at ZDBF2-DMR and ZNF597/NAA60-DMR in patient 1 and the GOM at GNAS-A/B-DMR in patient 3) would be consistent with the expression data using immortalized lymphocytes. Indeed, under normal conditions, ZDBF2 is expressed from a paternally derived chromosome with methylated ZDBF2-DMR,[11] ZNF597 is expressed from a maternally inherited chromosome with unmethylated ZNF597/NAA60-DMR,[12] and GNAS-A/B is expressed from a paternally derived chromosome with unmethylated GNAS-A/B-DMR.[13] By contrast, because MLID was absent from KOS14 patients, the GOM at IG-DMR and MEG3-DMR may occur as an isolated epimutation.

Several findings are notable with regard to the identified MLIDs. First, the degrees of GOMs at ZBDF2-DMR and ZNF597/NAA60-DMR in patient 1 differed between Leu-gDNA and Buc-gDNA, and the LOM at PPIEL-DMR in patient 1 and GOM at GNAS-A/B-DMR in patient 3 were found only in Buc-gDNA and Leu-gDNA, respectively. The extents of LOMs at MEG3-DMR also differed between Buc-gDNA and Leu-gDNA of patients 3 and 4. These findings imply the occurrence of somatic mosaicism caused by defective maintenance of methylation patterns at these DMRs during development.[14] This notion assumes that the associated GOM-type dominant MLIDs in patients 1 and 3 with TS14 have taken place independently of the LOM-type epimutations at MEG3-DMR. Consistent with this, no pathologic mutation was identified in the six genes for MLIDs. However, the possibility of a mutation in a hitherto-unknown gene for MLIDs has not been excluded formally, and further studies are required to elucidate the underlying factors for the development of MLIDs. Second, GOM at GNAS-A/B-DMR in Leu-gDNA of patient 3 was identified by HM450k analysis but not by pyrosequencing analysis, with discordant results for the same #44 and #45 CpG sites (). This would be due primarily to the differences in the control subjects and the normal ranges (average ± 3 SD range for HM450k analysis as used in the previous study,[4] and the minimum–maximum range for pyrosequencing analysis as used in our previous study[10]) between the two analyses. Age-, sex-, and race-specific normal ranges would need to be established to permit data comparison among different studies and different methods. Finally, MLIDs in patient 1 were detected in ID-unrelated nonclassic DMRs. Recent array-based studies have also detected MLIDs in such DMRs.[5] Thus, application of array-based analysis in epimutation-positive patients, including those who have been studied by conventional methylation analyses only, would facilitate the identification of MLIDs affecting such DMRs.

Several additional matters are also worth pointing out. First, MEG8-DMR exhibited methylation patterns opposite to those of IG-DMR and MEG3-DMR, as previously reported.[15,16] This exemplifies the complexity of the methylation patterns of DMRs within a given imprinted region. Second, GPR1-AS-DMR was hypermethylated in Leu-gDNA of all patients and control subjects. Thus, although GPR1-AS-DMR remains a maternally methylated DMR and appears to control the methylation pattern of the paternally methylated secondary ZDBF2-DMR in the placenta,[11,15,17] it is unlikely that such an interaction underlies in the development of GOM at ZDBF2-DMR in Leu-gDNA and Buc-gDNA of patient 1. Third, LOM at ZNF597-DMR co-existed with GOM at ZNF597/NAA60-DMR in Leu-gDNA of patient 1. This is consistent with the maternally methylated ZNF597-DMR functioning as an upstream regulator for the methylation pattern of the secondary paternally methylated ZNF597/NAA60-DMR.[12,15]

Clinical features of patients 1 and 3 appear to be within the phenotypic spectrum of TS14, despite the presence of MLIDs. Indeed, TS14 patients have been identified primarily by PWS- or SRS-like phenotypes.[2] However, it is known that UPD(16)mat with GOM of ZNF597-DMR and excessive ZNF597 expression is associated with pre- and postnatal growth failure,[18] as is UPD(20)mat with GOM at GNAS-A/B-DMR and no GNAS-A/B expression.[19] Thus, it may be possible that the LOM of ZNF597-DMR is relevant to the normal fetal growth in patient 1, and that the GOM at GNAS-A/B -DMR is involved in the growth failure in patient 3. In addition, because ZDBF2 is expressed in the fetal mouse,[11] the GOM of ZDBF2-DMR also may have played a certain role in the growth pattern of patient 1. Furthermore, although it remains unknown why TS14 patients show PWS- or SRS-like features, more severe LOMs at MEG3-DMR in patient 1 than in patients 2–4 might be involved in the phenotypic difference. In addition, although the methylation patterns of the PWS- and SRS-related DMRs were normal in the examined tissues of patients 1–4, there might be abnormally methylated DMR(s) hidden in the critical tissues.[20]

In summary, we identified MLIDs in two of four patients with TS14 and none of five patients with KOS14 with epimutations. Further studies will permit a better assessment of the frequency of MLIDs and the phenotypic consequences of MLIDs in patients with IDs.

Disclosure

The authors declare no conflict of interest.

Table 1

Methylation indices (%) for CpG dinucleotides determined by pyrosequencing analysis for leukocyte genomic DNA samples