Evaluation of enzyme activity predictions for variants of unknown significance in Arylsulfatase A.

Continued advances in variant effect prediction are necessary to demonstrate the ability of machine learning methods to accurately determine the clinical impact of variants of unknown significance (VUS). Towards this goal, the ARSA Critical Assessment of Genome Interpretation (CAGI) challenge was designed to characterize progress by utilizing 219 experimentally assayed missense VUS in the Arylsulfatase A (ARSA) gene to assess the performance of community-submitted predictions of variant functional effects. The challenge involved 15 teams, and evaluated additional predictions from established and recently released models. Notably, a model developed by participants of a genetics and coding bootcamp, trained with standard machine-learning tools in Python, demonstrated superior performance among submissions. Furthermore, the study observed that state-of-the-art deep learning methods provided small but statistically significant improvement in predictive performance compared to less elaborate techniques. These findings underscore the utility of variant effect prediction, and the potential for models trained with modest resources to accurately classify VUS in genetic and clinical research.

3D genome mapping identifies subgroup-specific chromosome conformations and tumor-dependency genes in ependymoma.

Ependymoma is a tumor of the brain or spinal cord. The two most common and aggressive molecular groups of ependymoma are the supratentorial ZFTA-fusion associated and the posterior fossa ependymoma group A. In both groups, tumors occur mainly in young children and frequently recur after treatment. Although molecular mechanisms underlying these diseases have recently been uncovered, they remain difficult to target and innovative therapeutic approaches are urgently needed. Here, we use genome-wide chromosome conformation capture (Hi-C), complemented with CTCF and H3K27ac ChIP-seq, as well as gene expression and DNA methylation analysis in primary and relapsed ependymoma tumors, to identify chromosomal conformations and regulatory mechanisms associated with aberrant gene expression. In particular, we observe the formation of new topologically associating domains ('neo-TADs') caused by structural variants, group-specific 3D chromatin loops, and the replacement of CTCF insulators by DNA hyper-methylation. Through inhibition experiments, we validate that genes implicated by these 3D genome conformations are essential for the survival of patient-derived ependymoma models in a group-specific manner. Thus, this study extends our ability to reveal tumor-dependency genes by 3D genome conformations even in tumors that lack targetable genetic alterations.

A Genetics-First Approach Revealed Monogenic Disorders in Patients With ARM and VACTERL Anomalies.

Background: The VATER/VACTERL association (VACTERL) is defined as the non-random occurrence of the following congenital anomalies: Vertebral, Anal, Cardiac, Tracheal-Esophageal, Renal, and Limb anomalies. As no unequivocal candidate gene has been identified yet, patients are diagnosed phenotypically. The aims of this study were to identify patients with monogenic disorders using a genetics-first approach, and to study whether variants in candidate genes are involved in the etiology of VACTERL or the individual features of VACTERL: Anorectal malformation (ARM) or esophageal atresia with or without trachea-esophageal fistula (EA/TEF). Methods: Using molecular inversion probes, a candidate gene panel of 56 genes was sequenced in three patient groups: VACTERL (n = 211), ARM (n = 204), and EA/TEF (n = 95). Loss-of-function (LoF) and additional likely pathogenic missense variants, were prioritized and validated using Sanger sequencing. Validated variants were tested for segregation and patients were clinically re-evaluated. Results: In 7 out of the 510 patients (1.4%), pathogenic or likely pathogenic variants were identified in SALL1, SALL4, and MID1, genes that are associated with Townes-Brocks, Duane-radial-ray, and Opitz-G/BBB syndrome. These syndromes always include ARM or EA/TEF, in combination with at least two other VACTERL features. We did not identify LoF variants in the remaining candidate genes. Conclusions: None of the other candidate genes were identified as novel unequivocal disease genes for VACTERL. However, a genetics-first approach allowed refinement of the clinical diagnosis in seven patients, in whom an alternative molecular-based diagnosis was found with important implications for the counseling of the families.

Hi-C Identifies Complex Genomic Rearrangements and TAD-Shuffling in Developmental Diseases.

Genome-wide analysis methods, such as array comparative genomic hybridization (CGH) and whole-genome sequencing (WGS), have greatly advanced the identification of structural variants (SVs) in the human genome. However, even with standard high-throughput sequencing techniques, complex rearrangements with multiple breakpoints are often difficult to resolve, and predicting their effects on gene expression and phenotype remains a challenge. Here, we address these problems by using high-throughput chromosome conformation capture (Hi-C) generated from cultured cells of nine individuals with developmental disorders (DDs). Three individuals had previously been identified as harboring duplications at the SOX9 locus and six had been identified with translocations. Hi-C resolved the positions of the duplications and was instructive in interpreting their distinct pathogenic effects, including the formation of new topologically associating domains (neo-TADs). Hi-C was very sensitive in detecting translocations, and it revealed previously unrecognized complex rearrangements at the breakpoints. In several cases, we observed the formation of fused-TADs promoting ectopic enhancer-promoter interactions that were likely to be involved in the disease pathology. In summary, we show that Hi-C is a sensible method for the detection of complex SVs in a clinical setting. The results help interpret the possible pathogenic effects of the SVs in individuals with DDs.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations.

INTRODUCTION: Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model. METHODS: We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing. RESULTS: When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing. CONCLUSION: Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Serial genomic inversions induce tissue-specific architectural stripes, gene misexpression and congenital malformations.

Balanced chromosomal rearrangements such as inversions and translocations can cause congenital disease or cancer by inappropriately rewiring promoter-enhancer contacts1,2. To study the potentially pathogenic consequences of balanced chromosomal rearrangements, we generated a series of genomic inversions by placing an active limb enhancer cluster from the Epha4 regulatory domain at different positions within a neighbouring gene-dense region and investigated their effects on gene regulation in vivo in mice. Expression studies and high-throughput chromosome conformation capture from embryonic limb buds showed that the enhancer cluster activated several genes downstream that are located within asymmetric regions of contact, the so-called architectural stripes3. The ectopic activation of genes led to a limb phenotype that could be rescued by deleting the CCCTC-binding factor (CTCF) anchor of the stripe. Architectural stripes appear to be driven by enhancer activity, because they do not form in mouse embryonic stem cells. Furthermore, we show that architectural stripes are a frequent feature of developmental three-dimensional genome architecture often associated with active enhancers. Therefore, balanced chromosomal rearrangements can induce ectopic gene expression and the formation of asymmetric chromatin contact patterns that are dependent on CTCF anchors and enhancer activity.

Somatic variants in autosomal dominant genes are a rare cause of sporadic Alzheimer's disease.

INTRODUCTION: A minority of patients with sporadic early-onset Alzheimer's disease (AD) exhibit de novo germ line mutations in the autosomal dominant genes such as APP, PSEN1, or PSEN2. We hypothesized that negatively screened patients may harbor somatic variants in these genes. METHODS: We applied an ultrasensitive approach based on single-molecule molecular inversion probes followed by deep next generation sequencing of 11 genes to 100 brain and 355 blood samples from 445 sporadic patients with AD (>80% exhibited an early onset, <66 years). RESULTS: We identified and confirmed nine somatic variants (allele fractions: 0.2%-10.8%): two APP, five SORL1, one NCSTN, and one MARK4 variants by independent amplicon-based deep sequencing. DISCUSSION: Two of the SORL1 variant might have contributed to the disease, the two APP variants were interpreted as likely benign and the other variants remained of unknown significance. Somatic variants in the autosomal dominant AD genes may not be a common cause of sporadic AD, including early onset cases.

Expanding the clinical spectrum of recessive truncating mutations of KLHL7 to a Bohring-Opitz-like phenotype.

BACKGROUND: Bohring-Opitz syndrome (BOS) is a rare genetic disorder characterised by a recognisable craniofacial appearance and a typical 'BOS' posture. BOS is caused by sporadic mutations ofASXL1. However, several typical patients with BOS have no molecular diagnosis, suggesting clinical and genetic heterogeneity. OBJECTIVES: To expand the phenotypical spectrum of autosomal recessive variants of KLHL7, reported as causing Crisponi syndrome/cold-induced sweating syndrome type 1 (CS/CISS1)-like syndrome. METHODS: We performed whole-exome sequencing in two families with a suspected recessive mode of inheritance. We used the Matchmaker Exchange initiative to identify additional patients. RESULTS: Here, we report six patients with microcephaly, facial dysmorphism, including exophthalmos, nevus flammeus of the glabella and joint contractures with a suspected BOS posture in five out of six patients. We identified autosomal recessive truncating mutations in the KLHL7 gene. KLHL7 encodes a BTB-kelch protein implicated in the cell cycle and in protein degradation by the ubiquitin-proteasome pathway. Recently, biallelic mutations in the KLHL7 gene were reported in four families and associated with CS/CISS1, characterised by clinical features overlapping with our patients. CONCLUSION: We have expanded the clinical spectrum of KLHL7 autosomal recessive variants by describing a syndrome with features overlapping CS/CISS1 and BOS.

The spectrum of DNMT3A variants in Tatton-Brown-Rahman syndrome overlaps with that in hematologic malignancies.

De novo, germline variants in DNMT3A cause Tatton-Brown-Rahman syndrome (TBRS). This condition is characterized by overgrowth, distinctive facial appearance, and intellectual disability. Somatic DNMT3A variants frequently occur in hematologic malignances, particularly acute myeloid leukemia. The Arg882 residue is the most common site of somatic DNMT3A variants, and has also been altered in patients with TBRS. Here we present three additional patients with this disorder attributed to DNMT3A germline variants that disrupt the Arg882 codon, suggesting that this codon may be a germline mutation hotspot in this disorder. Furthermore, based on the investigation of previously reported variants in patients with TBRS, we found overlap in the spectrum of DNMT3A variants observed in this disorder and somatic variants in hematological malignancies.

Ultra-sensitive Sequencing Identifies High Prevalence of Clonal Hematopoiesis-Associated Mutations throughout Adult Life.

Clonal hematopoiesis results from somatic mutations in hematopoietic stem cells, which give an advantage to mutant cells, driving their clonal expansion and potentially leading to leukemia. The acquisition of clonal hematopoiesis-driver mutations (CHDMs) occurs with normal aging and these mutations have been detected in more than 10% of individuals ≥65 years. We aimed to examine the prevalence and characteristics of CHDMs throughout adult life. We developed a targeted re-sequencing assay combining high-throughput with ultra-high sensitivity based on single-molecule molecular inversion probes (smMIPs). Using smMIPs, we screened more than 100 loci for CHDMs in more than 2,000 blood DNA samples from population controls between 20 and 69 years of age. Loci screened included 40 regions known to drive clonal hematopoiesis when mutated and 64 novel candidate loci. We identified 224 somatic mutations throughout our cohort, of which 216 were coding mutations in known driver genes (DNMT3A, JAK2, GNAS, TET2, and ASXL1), including 196 point mutations and 20 indels. Our assay's improved sensitivity allowed us to detect mutations with variant allele frequencies as low as 0.001. CHDMs were identified in more than 20% of individuals 60 to 69 years of age and in 3% of individuals 20 to 29 years of age, approximately double the previously reported prevalence despite screening a limited set of loci. Our findings support the occurrence of clonal hematopoiesis-associated mutations as a widespread mechanism linked with aging, suggesting that mosaicism as a result of clonal evolution of cells harboring somatic mutations is a universal mechanism occurring at all ages in healthy humans.

Overlapping SETBP1 gain-of-function mutations in Schinzel-Giedion syndrome and hematologic malignancies.

Schinzel-Giedion syndrome (SGS) is a rare developmental disorder characterized by multiple malformations, severe neurological alterations and increased risk of malignancy. SGS is caused by de novo germline mutations clustering to a 12bp hotspot in exon 4 of SETBP1. Mutations in this hotspot disrupt a degron, a signal for the regulation of protein degradation, and lead to the accumulation of SETBP1 protein. Overlapping SETBP1 hotspot mutations have been observed recurrently as somatic events in leukemia. We collected clinical information of 47 SGS patients (including 26 novel cases) with germline SETBP1 mutations and of four individuals with a milder phenotype caused by de novo germline mutations adjacent to the SETBP1 hotspot. Different mutations within and around the SETBP1 hotspot have varying effects on SETBP1 stability and protein levels in vitro and in in silico modeling. Substitutions in SETBP1 residue I871 result in a weak increase in protein levels and mutations affecting this residue are significantly more frequent in SGS than in leukemia. On the other hand, substitutions in residue D868 lead to the largest increase in protein levels. Individuals with germline mutations affecting D868 have enhanced cell proliferation in vitro and higher incidence of cancer compared to patients with other germline SETBP1 mutations. Our findings substantiate that, despite their overlap, somatic SETBP1 mutations driving malignancy are more disruptive to the degron than germline SETBP1 mutations causing SGS. Additionally, this suggests that the functional threshold for the development of cancer driven by the disruption of the SETBP1 degron is higher than for the alteration in prenatal development in SGS. Drawing on previous studies of somatic SETBP1 mutations in leukemia, our results reveal a genotype-phenotype correlation in germline SETBP1 mutations spanning a molecular, cellular and clinical phenotype.

New insights into the generation and role of de novo mutations in health and disease.

Aside from inheriting half of the genome of each of our parents, we are born with a small number of novel mutations that occurred during gametogenesis and postzygotically. Recent genome and exome sequencing studies of parent-offspring trios have provided the first insights into the number and distribution of these de novo mutations in health and disease, pointing to risk factors that increase their number in the offspring. De novo mutations have been shown to be a major cause of severe early-onset genetic disorders such as intellectual disability, autism spectrum disorder, and other developmental diseases. In fact, the occurrence of novel mutations in each generation explains why these reproductively lethal disorders continue to occur in our population. Recent studies have also shown that de novo mutations are predominantly of paternal origin and that their number increases with advanced paternal age. Here, we review the recent literature on de novo mutations, covering their detection, biological characterization, and medical impact.

Post-zygotic Point Mutations Are an Underrecognized Source of De Novo Genomic Variation.

De novo mutations are recognized both as an important source of genetic variation and as a prominent cause of sporadic disease in humans. Mutations identified as de novo are generally assumed to have occurred during gametogenesis and, consequently, to be present as germline events in an individual. Because Sanger sequencing does not provide the sensitivity to reliably distinguish somatic from germline mutations, the proportion of de novo mutations that occur somatically rather than in the germline remains largely unknown. To determine the contribution of post-zygotic events to de novo mutations, we analyzed a set of 107 de novo mutations in 50 parent-offspring trios. Using four different sequencing techniques, we found that 7 (6.5%) of these presumed germline de novo mutations were in fact present as mosaic mutations in the blood of the offspring and were therefore likely to have occurred post-zygotically. Furthermore, genome-wide analysis of "de novo" variants in the proband led to the identification of 4/4,081 variants that were also detectable in the blood of one of the parents, implying parental mosaicism as the origin of these variants. Thus, our results show that an important fraction of de novo mutations presumed to be germline in fact occurred either post-zygotically in the offspring or were inherited as a consequence of low-level mosaicism in one of the parents.

Thyroid hormone resistance syndrome due to mutations in the thyroid hormone receptor α gene (THRA).

BACKGROUND: Resistance to thyroid hormone is characterised by a lack of response of peripheral tissues to the active form of thyroid hormone (triiodothyronine, T3). In about 85% of cases, a mutation in THRB, the gene coding for thyroid receptor ß (TRß), is the cause of this disorder. Recently, individual reports described the first patients with thyroid hormone receptor α gene (THRA) defects. METHODS: We used longitudinal clinical assessments over a period of 18 years at one hospital setting combined with biochemical and molecular studies to characterise a novel thyroid hormone resistance syndrome in a cohort of six patients from five families. FINDINGS: Using whole exome sequencing and subsequent Sanger sequencing, we identified truncating and missense mutations in the THRA gene in five of six individuals and describe a distinct and consistent phenotype of mild hypothyroidism (growth retardation, relatively high birth length and weight, mild-to-moderate mental retardation, mild skeletal dysplasia and constipation), specific facial features (round, somewhat coarse and flat face) and macrocephaly. Laboratory investigations revealed anaemia and slightly elevated cholesterol, while the thyroid profile showed low free thyroxine (fT4) levels coupled with high free T3 (fT3), leading to an altered T4 : T3 ratio, along with normal thyroid-stimulating hormone levels. We observed a genotype-phenotype correlation, with milder outcomes for missense mutations and more severe phenotypical effects for truncating mutations. INTERPRETATION: THRA mutations may be more common than expected. In patients with clinical symptoms of mild hypothyreosis without confirmation in endocrine studies, a molecular study of THRA defects is strongly recommended.

Neu-Laxova syndrome is a heterogeneous metabolic disorder caused by defects in enzymes of the L-serine biosynthesis pathway.

Neu-Laxova syndrome (NLS) is a rare autosomal-recessive disorder characterized by a recognizable pattern of severe malformations leading to prenatal or early postnatal lethality. Homozygous mutations in PHGDH, a gene involved in the first and limiting step in L-serine biosynthesis, were recently identified as the cause of the disease in three families. By studying a cohort of 12 unrelated families affected by NLS, we provide evidence that NLS is genetically heterogeneous and can be caused by mutations in all three genes encoding enzymes of the L-serine biosynthesis pathway. Consistent with recently reported findings, we could identify PHGDH missense mutations in three unrelated families of our cohort. Furthermore, we mapped an overlapping homozygous chromosome 9 region containing PSAT1 in four consanguineous families. This gene encodes phosphoserine aminotransferase, the enzyme for the second step in L-serine biosynthesis. We identified six families with three different missense and frameshift PSAT1 mutations fully segregating with the disease. In another family, we discovered a homozygous frameshift mutation in PSPH, the gene encoding phosphoserine phosphatase, which catalyzes the last step of L-serine biosynthesis. Interestingly, all three identified genes have been previously implicated in serine-deficiency disorders, characterized by variable neurological manifestations. Our findings expand our understanding of NLS as a disorder of the L-serine biosynthesis pathway and suggest that NLS represents the severe end of serine-deficiency disorders, demonstrating that certain complex syndromes characterized by early lethality could indeed be the extreme end of the phenotypic spectrum of already known disorders.