Exome sequencing improves the molecular diagnostics of paediatric unexplained neurodevelopmental disorders.

BACKGROUND: Neurodevelopmental disorders (NDDs) and/or associated multiple congenital abnormalities (MCAs) represent a genetically heterogeneous group of conditions with an adverse prognosis for the quality of intellectual and social abilities and common daily functioning. The rapid development of exome sequencing (ES) techniques, together with trio-based analysis, nowadays leads to up to 50% diagnostic yield. Therefore, it is considered as the state-of-the-art approach in these diagnoses. RESULTS: In our study, we present the results of ES in a cohort of 85 families with 90 children with severe NDDs and MCAs. The interconnection of the in-house bioinformatic pipeline and a unique algorithm for variant prioritization resulted in a diagnostic yield of up to 48.9% (44/90), including rare and novel causative variants (41/90) and intragenic copy-number variations (CNVs) (3/90). Of the total number of 47 causative variants, 53.2% (25/47) were novel, highlighting the clinical benefit of ES for unexplained NDDs. Moreover, trio-based ES was verified as a reliable tool for the detection of rare CNVs, ranging from intragenic exon deletions (GRIN2A, ZC4H2 genes) to a 6-Mb duplication. The functional analysis using PANTHER Gene Ontology confirmed the involvement of genes with causative variants in a wide spectrum of developmental processes and molecular pathways, which form essential structural and functional components of the central nervous system. CONCLUSION: Taken together, we present one of the first ES studies of this scale from the central European region. Based on the high diagnostic yield for paediatric NDDs in this study, 48.9%, we confirm trio-based ES as an effective and reliable first-tier diagnostic test in the genetic evaluation of children with NDDs.

Broad H3K4me3 domains: Maintaining cellular identity and their implication in super-enhancer hijacking.

The human and mouse genomes are complex from a genomic standpoint. Each cell has the same genomic sequence, yet a wide array of cell types exists due to the presence of a plethora of regulatory elements in the non-coding genome. Recent advances in epigenomic profiling have uncovered non-coding gene proximal promoters and distal enhancers of transcription genome-wide. Extension of promoter-associated H3K4me3 histone mark across the gene body, known as a broad H3K4me3 domain (H3K4me3-BD), is a signature of constitutive expression of cell-type-specific regulation and of tumour suppressor genes in healthy cells. Recently, it has been discovered that the presence of H3K4me3-BDs over oncogenes is a cancer-specific feature associated with their dysregulated gene expression and tumourigenesis. Moreover, it has been shown that the hijacking of clusters of enhancers, known as super-enhancers (SE), by proto-oncogenes results in the presence of H3K4me3-BDs over the gene body. Therefore, H3K4me3-BDs and SE crosstalk in healthy and cancer cells therefore represents an important mechanism to identify future treatments for patients with SE driven cancers.

A unique case of Bloom syndrome with a combination of genetic hits: A lesson from trio­based exome sequencing: A case report.

Pathogenic variants affecting the BLM gene are responsible for the manifestation of extremely rare cancer­predisposing Bloom syndrome. The present study reports on a case of an infant with a congenital hypotrophy, short stature and abnormal facial appearance. Initially she was examined using a routine molecular diagnostic algorithm, including the cytogenetic analysis of her karyotype, microarray analysis and methylation­specific MLPA, however, she remained undiagnosed on a molecular level. Therefore, she and her parents were enrolled in the project of trio­based exome sequencing (ES) using Human Core Exome kit. She was revealed as a carrier of an extremely rare combination of causative sequence variants altering the BLM gene (NM_000057.4), c.1642C>T and c.2207_2212delinsTAGATTC in the compound heterozygosity, resulting in a diagnosis of Bloom syndrome. Simultaneously, a mosaic loss of heterozygosity of chromosome 11p was detected and then confirmed as a borderline imprinting center 1 hypermethylation on chromosome 11p15. The diagnosis of Bloom syndrome and mosaic copy­number neutral loss of heterozygosity of chromosome 11p increases a lifetime risk to develop any types of malignancy. This case demonstrates the trio­based ES as a complex approach for the molecular diagnostics of rare pediatric diseases.

Molecular characterization and clinical outcome of B-cell precursor acute lymphoblastic leukemia with IG-MYC rearrangement.

Rarely, immunophenotypically immature B-cell precursor acute lymphoblastic leukemia (BCP-ALL) carries an immunoglobulin- MYC rearrangement (IG-MYC-r). This can result in diagnostic confusion with Burkitt lymphoma/leukemia and use of individualized treatment schedules of unproven efficacy. Here we compare the molecular characteristics of these conditions and investigate historic clinical outcome data. We identified 90 cases registered in a national BCP-ALL clinical trial/registry. When present, diagnostic material underwent cytogenetic, exome, methylome and transcriptome analyses. The outcomes analyzed were 3-year event-free survival and overall survival. IG-MYC-r was identified in diverse cytogenetic backgrounds, co-existing with either established BCP-ALL-specific abnormalities (high hyperdiploidy, n=3; KMT2A-rearrangement, n=6; iAMP21, n=1; BCR-ABL1, n=1); BCL2/BCL6-rearrangements (n=15); or, most commonly, as the only defining feature (n=64). Within this final group, precursor-like V(D)J breakpoints predominated (8/9) and KRAS mutations were common (5/11). DNA methylation identified a cluster of V(D)J-rearranged cases, clearly distinct from Burkitt leukemia/lymphoma. Children with IG-MYC-r within that subgroup had a 3-year event-free survival of 47% and overall survival of 60%, representing a high-risk BCP-ALL. To develop effective management strategies this group of patients must be allowed access to contemporary, minimal residual disease-adapted, prospective clinical trial protocols.

High-resolution simulations of chromatin folding at genomic rearrangements in malignant B cells provide mechanistic insights into proto-oncogene deregulation.

Genomic rearrangements are known to result in proto-oncogene deregulation in many cancers, but the link to 3D genome structure remains poorly understood. Here, we used the highly predictive heteromorphic polymer (HiP-HoP) model to predict chromatin conformations at the proto-oncogene CCND1 in healthy and malignant B cells. After confirming that the model gives good predictions of Hi-C data for the nonmalignant human B cell-derived cell line GM12878, we generated predictions for two cancer cell lines, U266 and Z-138. These possess genome rearrangements involving CCND1 and the immunoglobulin heavy locus (IGH), which we mapped using targeted genome sequencing. Our simulations showed that a rearrangement in U266 cells where a single IGH super-enhancer is inserted next to CCND1 leaves the local topologically associated domain (TAD) structure intact. We also observed extensive changes in enhancer-promoter interactions within the TAD, suggesting that it is the downstream chromatin remodeling which gives rise to the oncogene activation, rather than the presence of the inserted super-enhancer DNA sequence per se. Simulations of the IGH-CCND1 reciprocal translocation in Z-138 cells revealed that an oncogenic fusion TAD is created, encompassing CCND1 and the IGH super-enhancers. We predicted how the structure and expression of CCND1 changes in these different cell lines, validating this using qPCR and fluorescence in situ hybridization microscopy. Our work demonstrates the power of polymer simulations to predict differences in chromatin interactions and gene expression for different translocation breakpoints.

Structural variants shape the genomic landscape and clinical outcome of multiple myeloma.

Deciphering genomic architecture is key to identifying novel disease drivers and understanding the mechanisms underlying myeloma initiation and progression. In this work, using the CoMMpass dataset, we show that structural variants (SV) occur in a nonrandom fashion throughout the genome with an increased frequency in the t(4;14), RB1, or TP53 mutated cases and reduced frequency in t(11;14) cases. By mapping sites of chromosomal rearrangements to topologically associated domains and identifying significantly upregulated genes by RNAseq we identify both predicted and novel putative driver genes. These data highlight the heterogeneity of transcriptional dysregulation occurring as a consequence of both the canonical and novel structural variants. Further, it shows that the complex rearrangements chromoplexy, chromothripsis and templated insertions are common in MM with each variant having its own distinct frequency and impact on clinical outcome. Chromothripsis is associated with a significant independent negative impact on clinical outcome in newly diagnosed cases consistent with its use alongside other clinical and genetic risk factors to identify prognosis.

Dynamics of broad H3K4me3 domains uncover an epigenetic switch between cell identity and cancer-related genes.

Broad domains of H3K4 methylation have been associated with consistent expression of tissue-specific, cell identity, and tumor suppressor genes. Here, we identified broad domain-associated genes in healthy human thymic T cell populations and a collection of T cell acute lymphoblastic leukemia (T-ALL) primary samples and cell lines. We found that broad domains are highly dynamic throughout T cell differentiation, and their varying breadth allows the distinction between normal and neoplastic cells. Although broad domains preferentially associate with cell identity and tumor suppressor genes in normal thymocytes, they flag key oncogenes in T-ALL samples. Moreover, the expression of broad domain-associated genes, both coding and noncoding, is frequently deregulated in T-ALL. Using two distinct leukemic models, we showed that the ectopic expression of T-ALL oncogenic transcription factor preferentially impacts the expression of broad domain-associated genes in preleukemic cells. Finally, an H3K4me3 demethylase inhibitor differentially targets T-ALL cell lines depending on the extent and number of broad domains. Our results show that the regulation of broad H3K4me3 domains is associated with leukemogenesis, and suggest that the presence of these structures might be used for epigenetic prioritization of cancer-relevant genes, including long noncoding RNAs.

Novel de novo pathogenic variant in the GNAI1 gene as a cause of severe disorders of intellectual development.

Pathogenic sequence variant in the GNAI1 gene were recently introduced as a cause of novel syndrome with a manifestation of variable developmental delay and autistic features. In our study, we report a case of monozygotic twins with severe intellectual disability and motor delay and developmental dysphasia. Both probands and their parents were examined using multi-step molecular diagnostic algorithm including whole-exome sequencing (WES), resulting in the identification of a novel, de novo pathogenic sequence variant in the GNAI1 gene, NM_002069.6:c.815 A>G, p.(Asp272Gly) in probands. Using WES we also verified the microarray findings of a familial 8q24.23q24.3 duplication and heterozygous 5q13.2 deletion, not associated with clinical symptoms in probands. Our results confirmed the role of the GNAI1 gene in the pathogenesis of syndromic neurodevelopmental disorders. They support trio- or quatro-based WES as a suitable molecular diagnostics method for the simultaneous detection of clinically relevant sequence variants and CNVs in individuals with neurodevelopmental disorders and rare diseases.

Epigenomic translocation of H3K4me3 broad domains over oncogenes following hijacking of super-enhancers.

Chromosomal translocations are important drivers of haematological malignancies whereby proto-oncogenes are activated by juxtaposition with enhancers, often called enhancer hijacking We analyzed the epigenomic consequences of rearrangements between the super-enhancers of the immunoglobulin heavy locus (IGH) and proto-oncogene CCND1 that are common in B cell malignancies. By integrating BLUEPRINT epigenomic data with DNA breakpoint detection, we characterized the normal chromatin landscape of the human IGH locus and its dynamics after pathological genomic rearrangement. We detected an H3K4me3 broad domain (BD) within the IGH locus of healthy B cells that was absent in samples with IGH-CCND1 translocations. The appearance of H3K4me3-BD over CCND1 in the latter was associated with overexpression and extensive chromatin accessibility of its gene body. We observed similar cancer-specific H3K4me3-BDs associated with hijacking of super-enhancers of other common oncogenes in B cell (MAF, MYC, and FGFR3/NSD2) and T cell malignancies (LMO2, TLX3, and TAL1). Our analysis suggests that H3K4me3-BDs can be created by super-enhancers and supports the new concept of epigenomic translocation, in which the relocation of H3K4me3-BDs from cell identity genes to oncogenes accompanies the translocation of super-enhancers.

Case Report: Contiguous Xq22.3 Deletion Associated with ATS-ID Syndrome: From Genotype to Further Delineation of the Phenotype.

Alport syndrome with intellectual disability (ATS-ID, AMME complex; OMIM #300194) is an X-linked contiguous gene deletion syndrome associated with an Xq22.3 locus mainly characterized by hematuria, renal failure, hearing loss/deafness, neurodevelopmental disorder (NDD), midface retrusion, and elliptocytosis. It is thought that ATS-ID is caused by the loss of function of COL4A5 (ATS) and FACL4 (ACSL4) genes through the interstitial (micro)deletion of chromosomal band Xq22.3. We report detailed phenotypic description and results from genome-wide screening of a Czech family with diagnosis ATS-ID (proband, maternal uncle, and two female carriers). Female carriers showed mild clinical features of microscopic hematuria only, while affected males displayed several novel clinical features associated with ATS-ID. Utilization of whole-exome sequencing discovered the presence of approximately 3 Mb of deletion in the Xq23 area, which affected 19 genes from TSC22D3 to CHRDL1. We compared the clinical phenotype with previously reported three ATS-ID families worldwide and correlated their clinical manifestations with the incidence of genes in both telomeric and centromeric regions of the deleted chromosomal area. In addition to previously described phenotypes associated with aberrations in AMMECR1 and FACL4, we identified two genes, members of tripartite motif family MID2 and subunit of the proteasome PA700/19S complex (PSMD10), respectively, as prime candidate genes responsible for additional clinical features observed in our patients with ATS-ID. Overall, our findings further improve the knowledge about the clinical impact of Xq23 deletions and bring novel information about phenotype/genotype association of this chromosomal aberration.

GATA2 deficiency phenotype associated with tandem duplication of GATA2 and overexpression of GATA2-AS1.

A 3-year-old girl of nonconsanguineous healthy parents presented with cervical and mediastinal lymphadenopathy due to Mycobacterium fortuitum infection. Routine blood analysis showed normal hemoglobin, neutrophils, and platelets but profound mononuclear cell deficiency (monocytes < 0.1 × 109/L; B cells 78/µL; NK cells 48/µL). A 548 902-bp region containing GATA2 was sequenced by targeted capture and deep sequencing. This revealed a de novo 187-kb duplication of the entire GATA2 locus, containing a maternally inherited copy number variation deletion of 25 kb (GRCh37: esv2725896 and nsv513733). Many GATA2-associated phenotypes have been attributed to amino acid substitution, frameshift/deletion, loss of intronic enhancer function, or aberrant splicing. Gene deletion has been described, but other structural variation has not been reported in the germline configuration. In this case, duplication of the GATA2 locus was paradoxically associated with skewed diminished expression of GATA2 messenger RNA and loss of GATA2 protein. Chimeric RNA fusion transcripts were not detected. A possible mechanism involves increased transcription of the anti-sense long noncoding RNA GATA2-AS1 (RP11-472.220), which was increased several fold. This case further highlights that evaluation of the allele count is essential in any case of suspected GATA2-related syndrome.

Whole-genome sequencing reveals progressive versus stable myeloma precursor conditions as two distinct entities.

Multiple myeloma (MM) is consistently preceded by precursor conditions recognized clinically as monoclonal gammopathy of undetermined significance (MGUS) or smoldering myeloma (SMM). We interrogate the whole genome sequence (WGS) profile of 18 MGUS and compare them with those from 14 SMMs and 80 MMs. We show that cases with a non-progressing, clinically stable myeloma precursor condition (n = 15) are characterized by later initiation in the patient's life and by the absence of myeloma defining genomic events including: chromothripsis, templated insertions, mutations in driver genes, aneuploidy, and canonical APOBEC mutational activity. This data provides evidence that WGS can be used to recognize two biologically and clinically distinct myeloma precursor entities that are either progressive or stable.

Germline TET2 loss of function causes childhood immunodeficiency and lymphoma.

Molecular dissection of inborn errors of immunity can help to elucidate the nonredundant functions of individual genes. We studied 3 children with an immune dysregulation syndrome of susceptibility to infection, lymphadenopathy, hepatosplenomegaly, developmental delay, autoimmunity, and lymphoma of B-cell (n = 2) or T-cell (n = 1) origin. All 3 showed early autologous T-cell reconstitution following allogeneic hematopoietic stem cell transplantation. By whole-exome sequencing, we identified rare homozygous germline missense or nonsense variants in a known epigenetic regulator of gene expression: ten-eleven translocation methylcytosine dioxygenase 2 (TET2). Mutated TET2 protein was absent or enzymatically defective for 5-hydroxymethylating activity, resulting in whole-blood DNA hypermethylation. Circulating T cells showed an abnormal immunophenotype including expanded double-negative, but depleted follicular helper, T-cell compartments and impaired Fas-dependent apoptosis in 2 of 3 patients. Moreover, TET2-deficient B cells showed defective class-switch recombination. The hematopoietic potential of patient-derived induced pluripotent stem cells was skewed toward the myeloid lineage. These are the first reported cases of autosomal-recessive germline TET2 deficiency in humans, causing clinically significant immunodeficiency and an autoimmune lymphoproliferative syndrome with marked predisposition to lymphoma. This disease phenotype demonstrates the broad role of TET2 within the human immune system.

Microhomology-mediated end joining drives complex rearrangements and overexpression of MYC and PVT1 in multiple myeloma.

MYC is a widely acting transcription factor and its deregulation is a crucial event in many human cancers. MYC is important biologically and clinically in multiple myeloma, but the mechanisms underlying its dysregulation are poorly understood. We show that MYC rearrangements are present in 36.0% of newly diagnosed myeloma patients, as detected in the largest set of next generation sequencing data to date (n=1,267). Rearrangements were complex and associated with increased expression of MYC and PVT1, but not other genes at 8q24. The highest effect on gene expression was detected in cases where the MYC locus is juxtaposed next to super-enhancers associated with genes such as IGH, IGK, IGL, TXNDC5/BMP6, FAM46C and FOXO3 We identified three hotspots of recombination at 8q24, one of which is enriched for IGH-MYC translocations. Breakpoint analysis indicates primary myeloma rearrangements involving the IGH locus occur through non-homologous end joining, whereas secondary MYC rearrangements occur through microhomology-mediated end joining. This mechanism is different to lymphomas, where non-homologous end joining generates MYC rearrangements. Rearrangements resulted in overexpression of key genes and chromatin immunoprecipitation-sequencing identified that HK2, a member of the glucose metabolism pathway, is directly over-expressed through binding of MYC at its promoter.

Severe type I interferonopathy and unrestrained interferon signaling due to a homozygous germline mutation in STAT2.

Excessive type I interferon (IFNα/ß) activity is implicated in a spectrum of human disease, yet its direct role remains to be conclusively proven. We investigated two siblings with severe early-onset autoinflammatory disease and an elevated IFN signature. Whole-exome sequencing revealed a shared homozygous missense Arg148Trp variant in STAT2, a transcription factor that functions exclusively downstream of innate IFNs. Cells bearing STAT2R148W in homozygosity (but not heterozygosity) were hypersensitive to IFNα/ß, which manifest as prolonged Janus kinase-signal transducers and activators of transcription (STAT) signaling and transcriptional activation. We show that this gain of IFN activity results from the failure of mutant STAT2R148W to interact with ubiquitin-specific protease 18, a key STAT2-dependent negative regulator of IFNα/ß signaling. These observations reveal an essential in vivo function of STAT2 in the regulation of human IFNα/ß signaling, providing concrete evidence of the serious pathological consequences of unrestrained IFNα/ß activity and supporting efforts to target this pathway therapeutically in IFN-associated disease.

The spectrum of somatic mutations in monoclonal gammopathy of undetermined significance indicates a less complex genomic landscape than that in multiple myeloma.

Monoclonal gammopathy of undetermined significance is a pre-malignant precursor of multiple myeloma with a 1% risk of progression per year. Although targeted analyses have shown the presence of specific genetic abnormalities such as IGH translocations, RB1 deletion, 1q gain, hyperdiploidy or RAS gene mutations, little is known about the molecular mechanism of malignant transformation. We performed whole exome sequencing together with comparative genomic hybridization plus single nucleotide polymorphism array analysis in 33 flow-cytometry-separated abnormal plasma cell samples from patients with monoclonal gammopathy of undetermined significance to describe somatic gene mutations and chromosome changes at the genome-wide level. Non-synonymous mutations and copy-number alterations were present in 97.0% and in 60.6% of cases, respectively. Importantly, the number of somatic mutations was significantly lower in monoclonal gammopathy of undetermined significance than in myeloma (P<10-4) and we identified six genes that were significantly mutated in myeloma (KRAS, NRAS, DIS3, HIST1H1E, EGR1 and LTB) within the monoclonal gammopathy of undetermined significance dataset. We also found a positive correlation with increasing chromosome changes and somatic gene mutations. IGH translocations, comprising t(4;14), t(11;14), t(14;16) and t(14;20), were present in 27.3% of cases and in a similar frequency to myeloma, consistent with the primary lesion hypothesis. MYC translocations and TP53 deletions or mutations were not detected in samples from patients with monoclonal gammopathy of undetermined significance, indicating that they may be drivers of progression to myeloma. Data from this study show that monoclonal gammopathy of undetermined significance is genetically similar to myeloma, however overall genetic abnormalities are present at significantly lower levels in monoclonal gammopathy of undetermined significant than in myeloma.

A first Czech analysis of 1887 cases with monoclonal gammopathy of undetermined significance.

INTRODUCTION: Monoclonal gammopathy of undetermined significance (MGUS) is a premalignant condition with a risk of malignant conversion. PATIENTS AND METHODS: With the aim to estimate the cumulative risk MGUS progression to hematologic malignancies, we analyzed a nationwide population-based cohort of 1887 MGUS patients from the Czech Registry of Monoclonal Gammopathies (RMG) between 2007 and 2013. RESULTS: During the follow-up period (median 4 years; range 0.6-34.8), progression to hematologic malignancies was observed in 8.6% (162 of 1887) of patients. Factors associated with progression were as follows: M-protein concentration ≥1.5 g/dL, pathological sFLC (<0.26 or >1.65) ratio, bone marrow plasma cells (BMPCs) in cytology >5%, immunoparesis, age ≥69 years, and the level of serum hemoglobin at baseline <12.0 g/dL. Combining these factors, we propose a new risk model (CMG model). The risk of progression at 10 years was 1.6%, 16.9%, 22.9%, 39.4%, and 52.3%, respectively, if 0 (reference group), one, two, three, or four to five risk factors are present (P<.001) with HR 63 times higher compared to the reference MGUS group. CONCLUSION: The new CMG model was established with an advantage for better identification of MGUS patients at low risk.

Genomewide profiling of copy-number alteration in monoclonal gammopathy of undetermined significance.

Monoclonal gammopathy of undetermined significance (MGUS) is a benign condition with an approximate 1% annual risk of symptomatic plasma cell disorder development, mostly to multiple myeloma (MM). We performed genomewide screening of copy-number alterations (CNAs) in 90 MGUS and 33 MM patients using high-density DNA microarrays. We identified CNAs in a smaller proportion of MGUS (65.6%) than in MM (100.0%, P = 1.31 × 10-5 ) and showed median number of CNAs is lower in MGUS (3, range 0-22) than in MM (13, range 4-38, P = 1.82 × 10-10 ). In the MGUS cohort, the most frequent losses were located at 1p (5.6%), 6q (6.7%), 13q (30.0%), 14q (14.4%), 16q (8.9%), 21q (5.6%), and gains at 1q (23.3%), 2p (6.7%), 6p (13.3%), and Xq (7.8%). Hyperdiploidy was detected in 38.9% of MGUS cases, and the most frequent whole chromosome gains were 3 (25.6%), 5 (23.3%), 9 (37.8%), 15 (23.3%), and 19 (32.2%). We also identified CNAs such as 1p, 6q, 8p, 12p, 13q, 16q losses, 1q gain and hypodiploidy, which are potentially associated with an adverse prognosis in MGUS. In summary, we showed that MGUS is similar to MM in that it is a genetically heterogeneous disorder, but overall cytogenetic instability is lower than in MM, which confirms that genetic abnormalities play important role in monoclonal gammopathies.