Germline pathogenic variants in HNRNPU are associated with alterations in blood methylome.

HNRNPU encodes a multifunctional RNA-binding protein that plays critical roles in regulating pre-mRNA splicing, mRNA stability, and translation. Aberrant expression and dysregulation of HNRNPU have been implicated in various human diseases, including cancers and neurological disorders. We applied a next generation sequencing based assay (EPIC-NGS) to investigate genome-wide methylation profiling for >2 M CpGs for 7 individuals with a neurodevelopmental disorder associated with HNRNPU germline pathogenic loss-of-function variants. Compared to healthy individuals, 227 HNRNPU-associated differentially methylated positions were detected. Both hyper- and hypomethylation alterations were identified but the former predominated. The identification of a methylation episignature for HNRNPU-associated neurodevelopmental disorder (NDD) implicates HNPRNPU-related chromatin alterations in the aetiopathogenesis of this disorder and suggests that episignature profiling should have clinical utility as a predictor for the pathogenicity of HNRNPU variants of uncertain significance. The detection of a methylation episignaure for HNRNPU-associated NDD is consistent with a recent report of a methylation episignature for HNRNPK-associated NDD.

Epigenotype-genotype-phenotype correlations in SETD1A and SETD2 chromatin disorders.

Germline pathogenic variants in two genes encoding the lysine-specific histone methyltransferase genes SETD1A and SETD2 are associated with neurodevelopmental disorders (NDDs) characterized by developmental delay and congenital anomalies. The SETD1A and SETD2 gene products play a critical role in chromatin-mediated regulation of gene expression. Specific methylation episignatures have been detected for a range of chromatin gene-related NDDs and have impacted clinical practice by improving the interpretation of variant pathogenicity. To investigate if SETD1A and/or SETD2-related NDDs are associated with a detectable episignature, we undertook targeted genome-wide methylation profiling of > 2 M CpGs using a next-generation sequencing-based assay. A comparison of methylation profiles in patients with SETD1A variants (n = 6) did not reveal evidence of a strong methylation episignature. A review of the clinical and genetic features of the SETD2 patient group revealed that, as reported previously, there were phenotypic differences between patients with truncating mutations (n = 4, Luscan-Lumish syndrome; MIM:616831) and those with missense codon 1740 variants [p.Arg1740Trp (n = 4) and p.Arg1740Gln (n = 2)]. Both SETD2 subgroups demonstrated a methylation episignature, which was characterized by hypomethylation and hypermethylation events, respectively. Within the codon 1740 subgroup, both the methylation changes and clinical phenotype were more severe in those with p.Arg1740Trp variants. We also noted that two of 10 cases with a SETD2-NDD had developed a neoplasm. These findings reveal novel epigenotype-genotype-phenotype correlations in SETD2-NDDs and predict a gain-of-function mechanism for SETD2 codon 1740 pathogenic variants.

Molecular characterisation of 36 multilocus imprinting disturbance (MLID) patients: a comprehensive approach.

BACKGROUND: Imprinting disorders (ImpDis) comprise diseases which are caused by aberrant regulation of monoallelically and parent-of-origin-dependent expressed genes. A characteristic molecular change in ImpDis patients is aberrant methylation signatures at disease-specific loci, without an obvious DNA change at the specific differentially methylated region (DMR). However, there is a growing number of reports on multilocus imprinting disturbances (MLIDs), i.e. aberrant methylation at different DMRs in the same patient. These MLIDs account for a significant number of patients with specific ImpDis, and several reports indicate a central role of pathogenic maternal effect variants in their aetiology by affecting the maturation of the oocyte and the early embryo. Though several studies on the prevalence and the molecular causes of MLID have been conducted, homogeneous datasets comprising both genomic and methylation data are still lacking. RESULTS: Based on a cohort of 36 MLID patients, we here present both methylation data obtained from next-generation sequencing (NGS, ImprintSeq) approaches and whole-exome sequencing (WES). The compilation of methylation data did not reveal a disease-specific MLID episignature, and a predisposition for the phenotypic modification was not obvious as well. In fact, this lack of epigenotype-phenotype correlation might be related to the mosaic distribution of imprinting defects and their functional relevance in specific tissues. CONCLUSIONS: Due to the higher sensitivity of NGS-based approaches, we suggest that ImprintSeq might be offered at reference centres in case of ImpDis patients with unusual phenotypes but MLID negative by conventional tests. By WES, additional MLID causes than the already known maternal effect variants could not be identified, neither in the patients nor in the maternal exomes. In cases with negative WES results, it is currently unclear to what extent either environmental factors or undetected genetic variants contribute to MLID.

Accurate Measurement of DNA Methylation: Challenges and Bias Correction.

DNA methylation is a key epigenetic modification involved in gene regulation whose contribution to disease susceptibility is still not fully understood. As the cost of genome sequencing technologies continues to drop, it will soon become commonplace to perform genome-wide quantification of DNA methylation at a single base-pair resolution. However, the demand for its accurate quantification might vary across studies. When the scope of the analysis is to detect regions of the genome with different methylation patterns between two or more conditions, e.g., case vs control; treatments vs placebo, accuracy is not crucial. This is the case in epigenome-wide association studies used as genome-wide screening of methylation changes to detect new candidate genes and regions associated with a specific disease or condition. If the aim of the analysis is to use DNA methylation measurements as a biomarker for diseases diagnosis and treatment (Laird, Nat Rev Cancer 3:253-266, 2003; Bock, Epigenomics 1:99-110, 2009), it is instead recommended to produce accurate methylation measurements. Furthermore, if the objective is the detection of DNA methylation in subclonal tumor cell populations or in circulating tumor DNA or in any case of mosaicism, the importance of accuracy becomes critical. The aim of this chapter is to describe the factors that could affect the precise measurement of methylation levels and a recent Bayesian statistical method called MethylCal and its extension that have been proposed to minimize this problem.

Comparison of methylation episignatures in KMT2B- and KMT2D-related human disorders.

Aim & methods: To investigate peripheral blood methylation episignatures in KMT2B-related dystonia (DYT-KMT2B), the authors undertook genome-wide methylation profiling of ∼2 M CpGs using a next-generation sequencing-based assay and compared the findings with those in controls and patients with KMT2D-related Kabuki syndrome type 1 (KS1). Results: A total of 1812 significantly differentially methylated CpG positions (false discovery rate < 0.05) were detected in DYT-KMT2B samples compared with controls. Multi-dimensional scaling analysis showed that the 10 DYT-KMT2B samples clustered together and separately from 29 controls and 10 with pathogenic variants in KMT2D. The authors found that most differentially methylated CpG positions were specific to one disorder and that all (DYT-KMT2B) and most (Kabuki syndrome type 1) methylation alterations in CpG islands were gain of methylation events. Conclusion: Using sensitive methylation profiling methodology, the authors replicated recent reports of a methylation episignature for DYT-KMT2B. These findings will facilitate the development of episignature-based assays to improve diagnostic accuracy.

The authors compared the DNA methylation patterns in blood from individuals with two rare neurodevelopmental disorders (childhood-onset dystonia [DYT-KMT2B] and Kabuki syndrome type 1) and healthy control samples. These two disorders are associated with pathogenic variants in KMT2B and KMT2D, which encode proteins with related functions but cause distinct inherited disorders. Comparison of the methylation patterns in the two disorders showed that most DNA regions with altered methylation patterns differed between the two disorders and controls. These findings suggest that analyzing DNA methylation patterns could improve diagnostic testing for these disorders and might provide insights into how the clinical features of these disorders are caused.

ImprintSeq, a novel tool to interrogate DNA methylation at human imprinted regions and diagnose multilocus imprinting disturbance.

PURPOSE: Disruptions of genomic imprinting are associated with congenital imprinting disorders (CIDs) and other disease states, including cancer. CIDs are most often associated with altered methylation at imprinted differentially methylated regions (iDMRs). In some cases, multiple iDMRs are affected causing multilocus imprinting disturbances (MLIDs). The availability of accurate, quantitative, and scalable high-throughput methods to interrogate multiple iDMRs simultaneously would enhance clinical diagnostics and research. METHODS: We report the development of a custom targeted methylation sequencing panel that covered most relevant 63 iDMRs for CIDs and the detection of MLIDs. We tested it in 70 healthy controls and 147 individuals with CIDs. We distinguished loss and gain of methylation per differentially methylated region and classified high and moderate methylation alterations. RESULTS: Across a range of CIDs with a variety of molecular mechanisms, ImprintSeq performed at 98.4% sensitivity, 99.9% specificity, and 99.9% accuracy (when compared with previous diagnostic testing). ImprintSeq was highly sensitive for detecting MLIDs and enabled diagnostic criteria for MLID to be proposed. In a child with extreme MLID profile a probable genetic cause was identified. CONCLUSION: ImprintSeq provides a novel assay for clinical diagnostic and research studies of CIDs, MLIDs, and the role of disordered imprinting in human disease states.

Alteration of Genomic Imprinting after Assisted Reproductive Technologies and Long-Term Health.

Assisted reproductive technologies (ART) are the treatment of choice for some infertile couples and even though these procedures are generally considered safe, children conceived by ART have shown higher reported risks of some perinatal and postnatal complications such as low birth weight, preterm birth, and childhood cancer. In addition, the frequency of some congenital imprinting disorders, like Beckwith-Wiedemann Syndrome and Silver-Russell Syndrome, is higher than expected in the general population after ART. Experimental evidence from animal studies suggests that ART can induce stress in the embryo and influence gene expression and DNA methylation. Human epigenome studies have generally revealed an enrichment of alterations in imprinted regions in children conceived by ART, but no global methylation alterations. ART procedures occur simultaneously with the establishment and maintenance of imprinting during embryonic development, so this may underlie the apparent sensitivity of imprinted regions to ART. The impact in adulthood of imprinting alterations that occurred during early embryonic development is still unclear, but some experimental evidence in mice showed higher risk to obesity and cardiovascular disease after the restriction of some imprinted genes in early embryonic development. This supports the hypothesis that imprinting alterations in early development might induce epigenetic programming of metabolism and affect long-term health. Given the growing use of ART, it is important to determine the impact of ART in genomic imprinting and long-term health.

A cross-disease meta-GWAS identifies four new susceptibility loci shared between systemic sclerosis and Crohn's disease.

Genome-wide association studies (GWASs) have identified a number of genetic risk loci associated with systemic sclerosis (SSc) and Crohn's disease (CD), some of which confer susceptibility to both diseases. In order to identify new risk loci shared between these two immune-mediated disorders, we performed a cross-disease meta-analysis including GWAS data from 5,734 SSc patients, 4,588 CD patients and 14,568 controls of European origin. We identified 4 new loci shared between SSc and CD, IL12RB2, IRF1/SLC22A5, STAT3 and an intergenic locus at 6p21.31. Pleiotropic variants within these loci showed opposite allelic effects in the two analysed diseases and all of them showed a significant effect on gene expression. In addition, an enrichment in the IL-12 family and type I interferon signaling pathways was observed among the set of SSc-CD common genetic risk loci. In conclusion, through the first cross-disease meta-analysis of SSc and CD, we identified genetic variants with pleiotropic effects on two clinically distinct immune-mediated disorders. The fact that all these pleiotropic SNPs have opposite allelic effects in SSc and CD reveals the complexity of the molecular mechanisms by which polymorphisms affect diseases.

SDHC epi-mutation testing in gastrointestinal stromal tumours and related tumours in clinical practice.

The enzyme succinate dehydrogenase (SDH) functions in the citric acid cycle and loss of function predisposes to the development of phaeochromocytoma/paraganglioma (PPGL), wild type gastrointestinal stromal tumour (wtGIST) and renal cell carcinoma. SDH-deficient tumours are most commonly associated with a germline SDH subunit gene (SDHA/B/C/D) mutation but can also be associated with epigenetic silencing of the SDHC gene. However, clinical diagnostic testing for an SDHC epimutation is not widely available. The objective of this study was to investigate the indications for and the optimum diagnostic pathways for the detection of SDHC epimutations in clinical practice. SDHC promoter methylation analysis of 32 paraffin embedded tumours (including 15 GIST and 17 PPGL) was performed using a pyrosequencing technique and correlated with SDHC gene expression. SDHC promoter methylation was identified in 6 (18.7%) tumours. All 6 SDHC epimutation cases presented with SDH deficient wtGIST and 3/6 cases had multiple primary tumours. No case of constitutional SDHC promoter hypermethylation was detected. Whole genome sequencing of germline DNA from three wtGIST cases with an SDHC epimutation, did not reveal any causative sequence anomalies. Herein, we recommend a diagnostic workflow for the detection of an SDHC epimutation in a service setting.

MethylCal: Bayesian calibration of methylation levels.

Bisulfite amplicon sequencing has become the primary choice for single-base methylation quantification of multiple targets in parallel. The main limitation of this technology is a preferential amplification of an allele and strand in the PCR due to methylation state. This effect, known as 'PCR bias', causes inaccurate estimation of the methylation levels and calibration methods based on standard controls have been proposed to correct for it. Here, we present a Bayesian calibration tool, MethylCal, which can analyse jointly all CpGs within a CpG island (CGI) or a Differentially Methylated Region (DMR), avoiding 'one-at-a-time' CpG calibration. This enables more precise modeling of the methylation levels observed in the standard controls. It also provides accurate predictions of the methylation levels not considered in the controlled experiment, a feature that is paramount in the derivation of the corrected methylation degree. We tested the proposed method on eight independent assays (two CpG islands and six imprinting DMRs) and demonstrated its benefits, including the ability to detect outliers. We also evaluated MethylCal's calibration in two practical cases, a clinical diagnostic test on 18 patients potentially affected by Beckwith-Wiedemann syndrome, and 17 individuals with celiac disease. The calibration of the methylation levels obtained by MethylCal allows a clearer identification of patients undergoing loss or gain of methylation in borderline cases and could influence further clinical or treatment decisions.

Erratum: PrimerPooler: automated primer pooling to prepare library for targeted sequencing.

[This corrects the article DOI: 10.1093/biomethods/bpx006.][This corrects the article DOI: 10.1093/biomethods/bpx006.].

Angioimmunoblastic T cell lymphoma: novel molecular insights by mutation profiling.

Angioimmunoblastic T cell lymphoma (AITL) originates from follicular helper T-cells and is characterised by a polymorphic infiltrate with the neoplastic T-cells forming small clusters around the follicle and high endothelial venules. Despite the recent advances in its phenotypic characterisation, the genetics and molecular mechanisms underlying AITL are not fully understood. In the present study, we performed whole exome sequencing in 9 cases of AITL from Taiwan (n = 6) and U.K. (n = 3). We confirmed frequent mutations in TET2 (9/9), DNMT3A (3/9), IDH2 (3/9), RHOA (3/9) and PLCG1 (2/9) as recently reported by others. More importantly, we identified mutations in TNFRSF21 (1/9), CCND3 (1/9) and SAMSN1 (1/9), which are not yet seen or strongly implicated in the pathogenesis of AITL. Among the pathogenic mutations identified in AITL, mutations in DNA methylation regulators TET2 and DNMT3A occur early in hematopoietic stem cells as shown by previous studies, and these genetic events enhance the self-renewal of hematopoietic stem cells, but are unlikely to have any major impact on T-cell differentiation. Mutations in RHOA, PLCG1 and TNFRSF21 (DR6), which encode proteins critical for T-cell biology, most likely promote T-cell differentiation and malignant transformation, consequently generating the malignant phenotype. Our findings extend the molecular insights into the multistage development of AITL.

PrimerPooler: automated primer pooling to prepare library for targeted sequencing.

Targeted next-generation sequencing based on PCR amplification involves pooling of hundreds to thousands of primers, for preamplification and subsequent parallel single/multiplex PCR. It is often necessary to allocate the set of primers into subpools, a common issue being potential cross-hybridization. For smaller numbers of primers, pool division can be done manually using trial and error to minimize potential hybridization, but this becomes inefficient and time consuming with increasing numbers of primer pairs. We developed PrimerPooler that automates swapping of primer pairs between any user-defined number of subpools to obtain combinations with low-potential interactions. PrimerPooler performs inter-/intra-primer hybridization analysis to identify the adverse interactions, as well as simultaneous mapping of all primers onto a genome sequence in a single run without requiring a prior index of the genome. This allows detection of overlapping primer pairs and allocation of these primer pairs into separate subpools where tiling approaches are used. Using PrimerPooler, 1153 primer pairs were assigned to three preamplification pools (388, 389 and 376 primer pairs each), then 144 subpools of six- to nine-plex PCR for Fluidigm Access Array PCR, followed by Illumina MiSeq sequencing. With optimized experimental protocols, an average of 3269 reads was achieved for the targeted regions, with 95% of targets covered by at least 50 reads, the minimal depth of reads for confident variant calling. PrimerPooler provides a fast and highly efficient stratification of primer pairs for targeted enrichment, thus ensuring representative amplification of the targeted sequences. PrimerPooler is also able to analyse degenerate primers, and is thus also useful for microbiological identification and related target sequencing.

LDLR and PCSK9 Are Associated with the Presence of Antiphospholipid Antibodies and the Development of Thrombosis in aPLA Carriers.

INTRODUCTION: The identification of the genetic risk factors that could discriminate non- thrombotic from thrombotic antiphospholipid antibodies (aPLA) carriers will improve prognosis of these patients. Several human studies have shown the presence of aPLAs associated with atherosclerotic plaque, which is a known risk factor for thrombosis. Hence, in order to determine the implication of atherosclerosis in the risk of developing thrombosis in aPLA positive patients, we performed a genetic association study with 3 candidate genes, APOH, LDLR and PCSK9. MATERIAL & METHODS: For genetic association study we analyzed 190 aPLA carriers -100 with non-thrombotic events and 90 with thrombotic events- and 557 healthy controls. Analyses were performed by χ2 test and were corrected by false discovery rate. To evaluate the functional implication of the newly established susceptibility loci, we performed expression analyses in 86 aPLA carrier individuals (43 with thrombotic manifestations and 43 without it) and in 45 healthy controls. RESULTS: Our results revealed significant associations after correction in SNPs located in LDLR gene with aPLA carriers and thrombotic aPLA carriers, when compared with healthy controls. The most significant association in LDLR gene was found between SNP rs129083082 and aPLA carriers in recessive model (adjusted P-value = 2.55 x 10-3; OR = 2.18; 95%CI = 1.49-3.21). Furthermore, our work detected significant allelic association after correction between thrombotic aPLA carriers and healthy controls in SNP rs562556 located in PCSK9 gene (adjusted P-value = 1.03 x 10-2; OR = 1.60; 95%CI = 1.24-2.06). Expression level study showed significantly decreased expression level of LDLR gene in aPLA carriers (P-value <0.0001; 95%CI 0.16-2.10; SE 0.38-1.27) in comparison to the control group. DISCUSSION: Our work has identified LDLR gene as a new susceptibility gene associated with the development of thrombosis in aPLA carriers, describing for the first time the deregulation of LDLR expression in individuals with aPLAs. Besides, thrombotic aPLA carriers also showed significant association with PCSK9 gene, a regulator of LDLR plasma levels. These results highlight the importance of atherosclerotic processes in the development of thrombosis in patients with aPLA.

Confirmation of CCR6 as a risk factor for anti-topoisomerase I antibodies in systemic sclerosis.

OBJECTIVES: The current knowledge of the influence of systemic sclerosis (SSc) risk loci in the clinical sub-phenotypes is still limited. The main limitation lies in the low frequency of some sub-phenotypes which could be solved by replication studies in independent cohorts and meta-analysis between studies. In this regard, CCR6 gene variants have been recently associated with anti-topoisomerase I positive (ATA+) production in SSc patients in a candidate gene study. This gene has been proposed to have a critical role in IL-17-driven autoimmunity in human diseases. METHODS: In order to confirm the association between CCR6 and ATA+ SSc patients, we performed an independent replication study in populations of European ancestry. We studied two CCR6 genetic variants (rs968334 and rs3093024) in a total of 901 ATA+ SSc cases, 3,258 ATA- SSc cases and 7,865 healthy controls and compared allelic frequencies for those SNPs in ATA+ SSc with healthy controls and also with ATA- SSc patients. RESULTS: The comparison performed between ATA+ SSc patients and healthy controls showed significant association with SNP rs968334 (p=4.88x10(-2), OR=1.11). When we compared ATA+ SSc cases with ATA- SSc, both SNPs, rs3093024 and rs968334, showed significant associations (p=2.89x10(-2), OR=1.13; p=1.69x10(-2), OR=1.15). Finally, in order to increase even more sample size and statistical power, we meta-analysed our study with the previous reported and found a significant association between SNP rs3093024 and ATA+ SSc patients (p=1.00x10(-4), OR=1.16) comparing with healthy controls. CONCLUSIONS: Our work confirms the association of CCR6 gene and ATA+ SSc patients.

Immunochip analysis identifies multiple susceptibility loci for systemic sclerosis.

In this study, 1,833 systemic sclerosis (SSc) cases and 3,466 controls were genotyped with the Immunochip array. Classical alleles, amino acid residues, and SNPs across the human leukocyte antigen (HLA) region were imputed and tested. These analyses resulted in a model composed of six polymorphic amino acid positions and seven SNPs that explained the observed significant associations in the region. In addition, a replication step comprising 4,017 SSc cases and 5,935 controls was carried out for several selected non-HLA variants, reaching a total of 5,850 cases and 9,401 controls of European ancestry. Following this strategy, we identified and validated three SSc risk loci, including DNASE1L3 at 3p14, the SCHIP1-IL12A locus at 3q25, and ATG5 at 6q21, as well as a suggested association of the TREH-DDX6 locus at 11q23. The associations of several previously reported SSc risk loci were validated and further refined, and the observed peak of association in PXK was related to DNASE1L3. Our study has increased the number of known genetic associations with SSc, provided further insight into the pleiotropic effects of shared autoimmune risk factors, and highlighted the power of dense mapping for detecting previously overlooked susceptibility loci.

Thrombotic antiphospholipid syndrome shows strong haplotypic association with SH2B3-ATXN2 locus.

BACKGROUND: Thrombotic antiphospholipid syndrome is defined as a complex form of thrombophilia that is developed by a fraction of antiphospholipid antibody (aPLA) carriers. Little is known about the genetic risk factors involved in thrombosis development among aPLA carriers. METHODS: To identify new loci conferring susceptibility to thrombotic antiphospholipid syndrome, a two-stage genotyping strategy was performed. In stage one, 19,000 CNV loci were genotyped in 14 thrombotic aPLA+ patients and 14 healthy controls by array-CGH. In stage two, significant CNV loci were fine-mapped in a larger cohort (85 thrombotic aPLA+, 100 non-thrombotic aPLA+ and 569 healthy controls). RESULTS: Array-CGH and fine-mapping analysis led to the identification of 12q24.12 locus as a new susceptibility locus for thrombotic APS. Within this region, a TAC risk haplotype comprising one SNP in SH2B3 gene (rs3184504) and two SNPs in ATXN2 gene (rs10774625 and rs653178) exhibited the strongest association with thrombotic antiphospholipid syndrome (p-value = 5,9 × 10(-4) OR 95% CI 1.84 (1.32-2.55)). CONCLUSION: The presence of a TAC risk haplotype in ATXN2-SH2B3 locus may contribute to increased thrombotic risk in aPLA carriers.