Loss-of-function variants affecting the STAGA complex component SUPT7L cause a developmental disorder with generalized lipodystrophy.

Generalized lipodystrophy is a feature of various hereditary disorders, often leading to a progeroid appearance. In the present study we identified a missense and a frameshift variant in a compound heterozygous state in SUPT7L in a boy with intrauterine growth retardation, generalized lipodystrophy, and additional progeroid features. SUPT7L encodes a component of the transcriptional coactivator complex STAGA. By transcriptome sequencing, we showed the predicted missense variant to cause aberrant splicing, leading to exon truncation and thereby to a complete absence of SUPT7L in dermal fibroblasts. In addition, we found altered expression of genes encoding DNA repair pathway components. This pathway was further investigated and an increased rate of DNA damage was detected in proband-derived fibroblasts and genome-edited HeLa cells. Finally, we performed transient overexpression of wildtype SUPT7L in both cellular systems, which normalizes the number of DNA damage events. Our findings suggest SUPT7L as a novel disease gene and underline the link between genome instability and progeroid phenotypes.

AXIN1 bi-allelic variants disrupting the C-terminal DIX domain cause craniometadiaphyseal osteosclerosis with hip dysplasia.

Sclerosing skeletal dysplasias result from an imbalance between bone formation and resorption. We identified three homozygous, C-terminally truncating AXIN1 variants in seven individuals from four families affected by macrocephaly, cranial hyperostosis, and vertebral endplate sclerosis. Other frequent findings included hip dysplasia, heart malformations, variable developmental delay, and hematological anomalies. In line with AXIN1 being a central component of the ß-catenin destruction complex, analyses of primary and genome-edited cells harboring the truncating variants revealed enhanced basal canonical Wnt pathway activity. All three AXIN1-truncating variants resulted in reduced protein levels and impaired AXIN1 polymerization mediated by its C-terminal DIX domain but partially retained Wnt-inhibitory function upon overexpression. Addition of a tankyrase inhibitor attenuated Wnt overactivity in the AXIN1-mutant model systems. Our data suggest that AXIN1 coordinates the action of osteoblasts and osteoclasts and that tankyrase inhibitors can attenuate the effects of AXIN1 hypomorphic variants.

Aorto-aortic bypass in an infant with middle aortic syndrome and Marfan syndrome: a 15-year follow-up.

We present a 15-year follow-up after aorto-aortic bypass surgery in a 7-month-old infant with middle aortic syndrome and confirmed Marfan syndrome. In anticipation of her growth, the length of the graft was adjusted to the anticipated length of the narrowed aorta in her adolescence. In addition, her height was controlled by oestrogen, and her growth was stopped at 178 cm. To date, the patient is free from aortic reoperation and lower limb malperfusion.

Biallelic variants in ADAMTS15 cause a novel form of distal arthrogryposis.

PURPOSE: We aimed to identify the underlying genetic cause for a novel form of distal arthrogryposis. METHODS: Rare variant family-based genomics, exome sequencing, and disease-specific panel sequencing were used to detect ADAMTS15 variants in affected individuals. Adamts15 expression was analyzed at the single-cell level during murine embryogenesis. Expression patterns were characterized using in situ hybridization and RNAscope. RESULTS: We identified homozygous rare variant alleles of ADAMTS15 in 5 affected individuals from 4 unrelated consanguineous families presenting with congenital flexion contractures of the interphalangeal joints and hypoplastic or absent palmar creases. Radiographic investigations showed physiological interphalangeal joint morphology. Additional features included knee, Achilles tendon, and toe contractures, spinal stiffness, scoliosis, and orthodontic abnormalities. Analysis of mouse whole-embryo single-cell sequencing data revealed a tightly regulated Adamts15 expression in the limb mesenchyme between embryonic stages E11.5 and E15.0. A perimuscular and peritendinous expression was evident in in situ hybridization in the developing mouse limb. In accordance, RNAscope analysis detected a significant coexpression with Osr1, but not with markers for skeletal muscle or joint formation. CONCLUSION: In aggregate, our findings provide evidence that rare biallelic recessive trait variants in ADAMTS15 cause a novel autosomal recessive connective tissue disorder, resulting in a distal arthrogryposis syndrome.

Autosomal recessive cutis laxa type IIIA: Report of a patient with severe phenotype and review of the literature.

Autosomal recessive cutis laxa type IIIA is a very rare genetic condition, caused by pathogenic variants in ALDH18A1, encoding delta-1-pyrroline-5-carboxylate synthase (P5CS). This enzyme catalyzes the reduction of glutamic acid to delta1-pyrroline-5-carboxylate, playing a key role in the de novo biosynthesis of proline, ornithine, and arginine. Autosomal recessive cutis laxa type IIIA is characterized by abundant and wrinkled skin, skeletal anomalies, cataract or corneal clouding and neuro-developmental disorders of variable degree. We report on a patient with autosomal recessive cutis laxa type IIIA, due to a homozygous missense c.1273C > T; p. (Arg425Cys) pathogenic variant in ALDH18A1. The patient presented a severe phenotype with serious urological involvement, peculiar cerebro-vascular abnormalities and neurodevelopmental compromise. This description contributes to better characterize the phenotypic spectrum associated with ALDH18A1 pathogenic variants, confirming the systemic involvement as a typical feature of autosomal recessive cutis laxa type IIIA.

ClearCNV: CNV calling from NGS panel data in the presence of ambiguity and noise.

MOTIVATION: While the identification of small variants in panel sequencing data can be considered a solved problem, the identification of larger, multi-exon copy number variants (CNVs) still poses a considerable challenge. Thus, CNV calling has not been established in all laboratories performing panel sequencing. At the same time, such laboratories have accumulated large datasets and thus have the need to identify CNVs on their data to close the diagnostic gap. RESULTS: In this article, we present our method clearCNV that addresses this need in two ways. First, it helps laboratories to properly assign datasets to enrichment kits. Based on homogeneous subsets of data, clearCNV identifies CNVs affecting the targeted regions. Using real-world datasets and validation, we show that our method is highly competitive with previous methods and preferable in terms of specificity. AVAILABILITY AND IMPLEMENTATION: The software is available for free under a permissible license at https://github.com/bihealth/clear-cnv. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Alternative splicing of BUD13 determines the severity of a developmental disorder with lipodystrophy and progeroid features.

PURPOSE: In this study we aimed to identify the molecular genetic cause of a progressive multisystem disease with prominent lipodystrophy. METHODS: In total, 5 affected individuals were investigated using exome sequencing. Dermal fibroblasts were characterized using RNA sequencing, proteomics, immunoblotting, immunostaining, and electron microscopy. Subcellular localization and rescue studies were performed. RESULTS: We identified a lipodystrophy phenotype with a typical facial appearance, corneal clouding, achalasia, progressive hearing loss, and variable severity. Although 3 individuals showed stunted growth, intellectual disability, and died within the first decade of life (A1, A2, and A3), 2 are adults with normal intellectual development (A4 and A5). All individuals harbored an identical homozygous nonsense variant affecting the retention and splicing complex component BUD13. The nucleotide substitution caused alternative splicing of BUD13 leading to a stable truncated protein whose expression positively correlated with disease expression and life expectancy. In dermal fibroblasts, we found elevated intron retention, a global reduction of spliceosomal proteins, and nuclei with multiple invaginations, which were more pronounced in A1, A2, and A3. Overexpression of both BUD13 isoforms normalized the nuclear morphology. CONCLUSION: Our results define a hitherto unknown syndrome and show that the alternative splice product converts a loss-of-function into a hypomorphic allele, thereby probably determining the severity of the disease and the survival of affected individuals.

Different ascending aortic phenotypes with similar mutations in 2 patients with Loeys-Dietz syndrome type 2.

Our goal was to present 2 infants with confirmed Loeys-Dietz syndrome. The missense mutations in exon 7 of the TGFBR2 gene are only 5 codons apart (c.1597T>C and c.1582C>G). Phenotypically, the aneurysms of the ascending aorta were restricted to different segments of the aorta: the suprajunctional segment in 1 patient and the aortic root in another. These cases highlight the complexity of signaling pathways and gene expression in the pathogenesis of aortic aneurysms.

AutozygosityMapper: Identification of disease-mutations in consanguineous families.

With the shift from SNP arrays to high-throughput sequencing, most researchers studying diseases in consanguineous families do not rely on linkage analysis any longer, but simply search for deleterious variants which are homozygous in all patients. AutozygosityMapper allows the fast and convenient identification of disease mutations in patients from consanguineous pedigrees by focussing on homozygous segments shared by all patients. Users can upload multi-sample VCF files, including WGS data, without any pre-processing. Genome-wide runs of homozygosity and the underlying genotypes are presented in graphical interfaces. AutozygosityMapper extends the functions of its predecessor, HomozygosityMapper, to the search for autozygous regions, in which all patients share the same homozygous genotype. We provide export of VCF files containing only the variants found in homozygous regions, this usually reduces the number of variants by two orders of magnitude. These regions can also directly be analysed with our disease mutation identification tool MutationDistiller. The application comes with simple and intuitive graphical interfaces for data upload, analysis, and results. We kept the structure of HomozygosityMapper so that previous users will find it easy to switch. With AutozygosityMapper, we provide a fast web-based way to identify disease mutations in consanguineous families. AutozygosityMapper is freely available at https://www.genecascade.org/AutozygosityMapper/.

Xq27.1 palindrome mediated interchromosomal insertion likely causes familial congenital bilateral laryngeal abductor paralysis (Plott syndrome).

Bilateral laryngeal abductor paralysis is a rare entity and the second most common cause of stridor in newborns. So far, no conclusive genetic or chromosomal aberration has been reported for X-linked isolated bilateral vocal cord paralysis, also referred to as Plott syndrome. Via whole genome sequencing (WGS), we identified a complex interchromosomal insertion in a large family with seven affected males. The 404 kb inserted fragment originates from chromosome 10q21.3, contains no genes and is inserted inversionally into the intergenic chromosomal region Xq27.1, 82 kb centromeric to the nearest gene SOX3. The patterns found at the breakpoint junctions resemble typical characteristics that arise in replication-based mechanisms with long-distance template switching. Non protein-coding insertions into the same genomic region have been described to result in different phenotypes, indicating that the phenotypic outcome likely depends on the introduction of regulatory elements. In conclusion, our data adds Plott syndrome as another entity, likely caused by the insertion of non-coding DNA into the intergenic chromosomal region Xq27.1. In this regard, we demonstrate the importance of WGS as a powerful diagnostic test in unsolved genetic diseases, as this genomic rearrangement has not been detected by current first-line diagnostic tests, i.e., exome sequencing and chromosomal microarray analysis.

Biallelic truncating variants in ATP9A cause a novel neurodevelopmental disorder involving postnatal microcephaly and failure to thrive.

BACKGROUND: Genes implicated in the Golgi and endosomal trafficking machinery are crucial for brain development, and mutations in them are particularly associated with postnatal microcephaly (POM). METHODS: Exome sequencing was performed in three affected individuals from two unrelated consanguineous families presenting with delayed neurodevelopment, intellectual disability of variable degree, POM and failure to thrive. Patient-derived fibroblasts were tested for functional effects of the variants. RESULTS: We detected homozygous truncating variants in ATP9A. While the variant in family A is predicted to result in an early premature termination codon, the variant in family B affects a canonical splice site. Both variants lead to a substantial reduction of ATP9A mRNA expression. It has been shown previously that ATP9A localises to early and recycling endosomes, whereas its depletion leads to altered gene expression of components from this compartment. Consistent with previous findings, we also observed overexpression of ARPC3 and SNX3, genes strongly interacting with ATP9A. CONCLUSION: In aggregate, our findings show that pathogenic variants in ATP9A cause a novel autosomal recessive neurodevelopmental disorder with POM. While the physiological function of endogenous ATP9A is still largely elusive, our results underline a crucial role of this gene in endosomal transport in brain tissue.

Complete lung agenesis caused by complex genomic rearrangements with neo-TAD formation at the SHH locus.

During human organogenesis, lung development is a timely and tightly regulated developmental process under the control of a large number of signaling molecules. Understanding how genetic variants can disturb normal lung development causing different lung malformations is a major goal for dissecting molecular mechanisms during embryogenesis. Here, through exome sequencing (ES), array CGH, genome sequencing (GS) and Hi-C, we aimed at elucidating the molecular basis of bilateral isolated lung agenesis in three fetuses born to a non-consanguineous family. We detected a complex genomic rearrangement containing duplicated, triplicated and deleted fragments involving the SHH locus in fetuses presenting complete agenesis of both lungs and near-complete agenesis of the trachea, diagnosed by ultrasound screening and confirmed at autopsy following termination. The rearrangement did not include SHH itself, but several regulatory elements for lung development, such as MACS1, a major SHH lung enhancer, and the neighboring genes MNX1 and NOM1. The rearrangement incorporated parts of two topologically associating domains (TADs) including their boundaries. Hi-C of cells from one of the affected fetuses showed the formation of two novel TADs each containing SHH enhancers and the MNX1 and NOM1 genes. Hi-C together with GS indicate that the new 3D conformation is likely causative for this condition by an inappropriate activation of MNX1 included in the neo-TADs by MACS1 enhancer, further highlighting the importance of the 3D chromatin conformation in human disease.

A CRISPR-Cas9-engineered mouse model for GPI-anchor deficiency mirrors human phenotypes and exhibits hippocampal synaptic dysfunctions.

Pathogenic germline mutations in PIGV lead to glycosylphosphatidylinositol biosynthesis deficiency (GPIBD). Individuals with pathogenic biallelic mutations in genes of the glycosylphosphatidylinositol (GPI)-anchor pathway exhibit cognitive impairments, motor delay, and often epilepsy. Thus far, the pathophysiology underlying the disease remains unclear, and suitable rodent models that mirror all symptoms observed in human patients have not been available. Therefore, we used CRISPR-Cas9 to introduce the most prevalent hypomorphic missense mutation in European patients, Pigv:c.1022C > A (p.A341E), at a site that is conserved in mice. Mirroring the human pathology, mutant Pigv341E mice exhibited deficits in motor coordination, cognitive impairments, and alterations in sociability and sleep patterns, as well as increased seizure susceptibility. Furthermore, immunohistochemistry revealed reduced synaptophysin immunoreactivity in Pigv341E mice, and electrophysiology recordings showed decreased hippocampal synaptic transmission that could underlie impaired memory formation. In single-cell RNA sequencing, Pigv341E-hippocampal cells exhibited changes in gene expression, most prominently in a subtype of microglia and subicular neurons. A significant reduction in Abl1 transcript levels in several cell clusters suggested a link to the signaling pathway of GPI-anchored ephrins. We also observed elevated levels of Hdc transcripts, which might affect histamine metabolism with consequences for circadian rhythm. This mouse model will not only open the doors to further investigation into the pathophysiology of GPIBD, but will also deepen our understanding of the role of GPI-anchor-related pathways in brain development.

Expanding the clinical and molecular spectrum of ATP6V1A related metabolic cutis laxa.

Several inborn errors of metabolism show cutis laxa as a highly recognizable feature. One group of these metabolic cutis laxa conditions is autosomal recessive cutis laxa type 2 caused by defects in v-ATPase components or the mitochondrial proline cycle. Besides cutis laxa, muscular hypotonia and cardiac abnormalities are hallmarks of autosomal recessive cutis laxa type 2D (ARCL2D) due to pathogenic variants in ATP6V1A encoding subunit A of the v-ATPase. Here, we report on three affected individuals from two families with ARCL2D in whom we performed whole exome and Sanger sequencing. We performed functional studies in fibroblasts from one individual, summarized all known probands' clinical, molecular, and biochemical features and compared them, also to other metabolic forms of cutis laxa. We identified novel missense and the first nonsense variant strongly affecting ATP6V1A expression. All six ARCL2D affected individuals show equally severe cutis laxa and dysmorphism at birth. While for one no information was available, two died in infancy and three are now adolescents with mild or absent intellectual disability. Muscular weakness, ptosis, contractures, and elevated muscle enzymes indicated a persistent myopathy. In cellular studies, a fragmented Golgi compartment, a delayed Brefeldin A-induced retrograde transport and glycosylation abnormalities were present in fibroblasts from two individuals. This is the second and confirmatory report on pathogenic variants in ATP6V1A as the cause of this extremely rare condition and the first to describe a nonsense allele. Our data highlight the tremendous clinical variability of ATP6V1A related phenotypes even within the same family.

An intronic splice site alteration in combination with a large deletion affecting VPS13B (COH1) causes Cohen syndrome.

Cohen syndrome (CS) is a rare, autosomal recessive disorder characterized by intellectual disability, postnatal microcephaly, facial abnormalities, abnormal truncal fat distribution, myopia, and pigmentary retinopathy. It is often considered an underdiagnosed condition, especially in children with developmental delay and intellectual disability. Here we report on four individuals from a large Jordanian family clinically diagnosed with CS. Using Trio Exome Sequencing (Trio-WES) and MLPA analyses we identified a maternally inherited novel intronic nucleotide substitution c.3446-23T>G leading to the activation of a cryptic splice site and a paternally inherited multi-exon deletion in VPS13B (previously termed COH1) in the index patient. Expression analysis showed a strong decrease of VPS13B mRNA levels and direct sequencing of cDNA confirmed splicing at a cryptic upstream splice acceptor site, resulting in the inclusion of 22 intronic bases. This extension results in a frameshift and a premature stop of translation (p.Gly1149Valfs*9). Segregation analysis revealed that three affected maternal cousins were homozygous for the intronic splice site variant. Our data show causality of both alterations and strongly suggest the expansion of the diagnostic strategy to search for intronic splice variants in molecularly unconfirmed patients affected by CS.

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.

VarFish: comprehensive DNA variant analysis for diagnostics and research.

VarFish is a user-friendly web application for the quality control, filtering, prioritization, analysis, and user-based annotation of DNA variant data with a focus on rare disease genetics. It is capable of processing variant call files with single or multiple samples. The variants are automatically annotated with population frequencies, molecular impact, and presence in databases such as ClinVar. Further, it provides support for pathogenicity scores including CADD, MutationTaster, and phenotypic similarity scores. Users can filter variants based on these annotations and presumed inheritance pattern and sort the results by these scores. Variants passing the filter are listed with their annotations and many useful link-outs to genome browsers, other gene/variant data portals, and external tools for variant assessment. VarFish allows users to create their own annotations including support for variant assessment following ACMG-AMP guidelines. In close collaboration with medical practitioners, VarFish was designed for variant analysis and prioritization in diagnostic and research settings as described in the software's extensive manual. The user interface has been optimized for supporting these protocols. Users can install VarFish on their own in-house servers where it provides additional lab notebook features for collaborative analysis and allows re-analysis of cases, e.g. after update of genotype or phenotype databases.

Biallelic variants in KYNU cause a multisystemic syndrome with hand hyperphalangism.

Catel-Manzke syndrome is characterized by the combination of Pierre Robin sequence and radial deviation, shortening as well as clinodactyly of the index fingers, due to an accessory ossification center. Mutations in TGDS have been identified as one cause of Catel-Manzke syndrome, but cannot be found as causative in every patient with the clinical diagnosis. We performed a chromosome microarray and/or exome sequencing in three patients with hand hyperphalangism, heart defect, short stature, and mild to severe developmental delay, all of whom were initially given a clinical diagnosis of Catel-Manzke syndrome. In one patient, we detected a large deletion of exons 1-8 and the missense variant c.1282C > T (p.Arg428Trp) in KYNU (NM_003937.2), whereas homozygous missense variants in KYNU were found in the other two patients (c.989G > A (p.Arg330Gln) and c.326G > C (p.Trp109Ser)). Plasma and urine metabolomic analysis of two patients indicated a block along the tryptophan catabolic pathway and urine organic acid analysis showed excretion of xanthurenic acid. Biallelic loss-of-function mutations in KYNU were recently described as a cause of NAD deficiency with vertebral, cardiac, renal and limb defects; however, no hand hyperphalangism was described in those patients, and Catel-Manzke syndrome was not discussed as a differential diagnosis. In conclusion, we present unrelated patients identified with biallelic variants in KYNU leading to kynureninase deficiency and xanthurenic aciduria as a very likely cause of their hyperphalangism, heart defect, short stature, and developmental delay. We suggest performance of urine organic acid analysis in patients with suspected Catel-Manzke syndrome, particularly in those with cardiac or vertebral defects or without mutations in TGDS.