Publisher Correction: Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Whole-genome sequencing of a sporadic primary immunodeficiency cohort.

Primary immunodeficiency (PID) is characterized by recurrent and often life-threatening infections, autoimmunity and cancer, and it poses major diagnostic and therapeutic challenges. Although the most severe forms of PID are identified in early childhood, most patients present in adulthood, typically with no apparent family history and a variable clinical phenotype of widespread immune dysregulation: about 25% of patients have autoimmune disease, allergy is prevalent and up to 10% develop lymphoid malignancies1-3. Consequently, in sporadic (or non-familial) PID genetic diagnosis is difficult and the role of genetics is not well defined. Here we address these challenges by performing whole-genome sequencing in a large PID cohort of 1,318 participants. An analysis of the coding regions of the genome in 886 index cases of PID found that disease-causing mutations in known genes that are implicated in monogenic PID occurred in 10.3% of these patients, and a Bayesian approach (BeviMed4) identified multiple new candidate PID-associated genes, including IVNS1ABP. We also examined the noncoding genome, and found deletions in regulatory regions that contribute to disease causation. In addition, we used a genome-wide association study to identify loci that are associated with PID, and found evidence for the colocalization of-and interplay between-novel high-penetrance monogenic variants and common variants (at the PTPN2 and SOCS1 loci). This begins to explain the contribution of common variants to the variable penetrance and phenotypic complexity that are observed in PID. Thus, using a cohort-based whole-genome-sequencing approach in the diagnosis of PID can increase diagnostic yield and further our understanding of the key pathways that influence immune responsiveness in humans.

REEV: review, evaluate and explain variants.

In the era of high throughput sequencing, special software is required for the clinical evaluation of genetic variants. We developed REEV (Review, Evaluate and Explain Variants), a user-friendly platform for clinicians and researchers in the field of rare disease genetics. Supporting data was aggregated from public data sources. We compared REEV with seven other tools for clinical variant evaluation. REEV (semi-)automatically fills individual ACMG criteria facilitating variant interpretation. REEV can store disease and phenotype data related to a case to use these for phenotype similarity measures. Users can create public permanent links for individual variants that can be saved as browser bookmarks and shared. REEV may help in the fast diagnostic assessment of genetic variants in a clinical as well as in a research context. REEV (https://reev.bihealth.org/) is free and open to all users and there is no login requirement.

Zebrafish as a model to investigate a biallelic gain-of-function variant in MSGN1, associated with a novel skeletal dysplasia syndrome.

BACKGROUND/OBJECTIVES: Rare genetic disorders causing specific congenital developmental abnormalities often manifest in single families. Investigation of disease-causing molecular features are most times lacking, although these investigations may open novel therapeutic options for patients. In this study, we aimed to identify the genetic cause in an Iranian patient with severe skeletal dysplasia and to model its molecular function in zebrafish embryos. RESULTS: The proband displays short stature and multiple skeletal abnormalities, including mesomelic dysplasia of the arms with complete humero-radio-ulna synostosis, arched clavicles, pelvic dysplasia, short and thin fibulae, proportionally short vertebrae, hyperlordosis and mild kyphosis. Exome sequencing of the patient revealed a novel homozygous c.374G > T, p.(Arg125Leu) missense variant in MSGN1 (NM_001105569). MSGN1, a basic-Helix-Loop-Helix transcription factor, plays a crucial role in formation of presomitic mesoderm progenitor cells/mesodermal stem cells during early developmental processes in vertebrates. Initial in vitro experiments show protein stability and correct intracellular localization of the novel variant in the nucleus and imply retained transcription factor function. To test the pathogenicity of the detected variant, we overexpressed wild-type and mutant msgn1 mRNA in zebrafish embryos and analyzed tbxta (T/brachyury/ntl). Overexpression of wild-type or mutant msgn1 mRNA significantly reduces tbxta expression in the tailbud compared to control embryos. Mutant msgn1 mRNA injected embryos depict a more severe effect, implying a gain-of-function mechanism. In vivo analysis on embryonic development was performed by clonal msgn1 overexpression in zebrafish embryos further demonstrated altered cell compartments in the presomitic mesoderm, notochord and pectoral fin buds. Detection of ectopic tbx6 and bmp2 expression in these embryos hint to affected downstream signals due to Msgn1 gain-of-function. CONCLUSION: In contrast to loss-of-function effects described in animal knockdown models, gain-of-function of MSGN1 explains the only mildly affected axial skeleton of the proband and rather normal vertebrae. In this context we observed notochord bending and potentially disruption of pectoral fin buds/upper extremity after overexpression of msgn1 in zebrafish embryos. The latter might result from Msgn1 function on mesenchymal stem cells or on chondrogenesis in these regions. In addition, we detected ectopic tbx6 and bmp2a expression after gain of Msgn1 function in zebrafish, which are interconnected to short stature, congenital scoliosis, limb shortening and prominent skeletal malformations in patients. Our findings highlight a rare, so far undescribed skeletal dysplasia syndrome associated with a gain-of-function mutation in MSGN1 and hint to its molecular downstream effectors.

RNA analysis and computer-aided facial phenotyping help to classify a novel TRIO splice site variant.

Pathogenic variants in TRIO, encoding the guanine nucleotide exchange factor, are associated with two distinct neurodevelopmental delay phenotypes: gain-of-function missense mutations within the spectrin repeats are causative for a severe developmental delay with macrocephaly (MIM: 618825), whereas loss-of-function missense variants in the GEF1 domain and truncating variants throughout the gene lead to a milder developmental delay and microcephaly (MIM: 617061). In three affected family members with mild intellectual disability/NDD and microcephaly, we detected a novel heterozygous TRIO variant at the last coding base of exon 31 (NM_007118.4:c.4716G>A). RNA analysis from patient-derived lymphoblastoid cells confirmed aberrant splicing resulting in the skipping of exon 31 (r.4615_4716del), leading to an in-frame deletion in the first Pleckstrin homology subdomain of the GEF1 domain: p.(Thr1539_Lys1572del). To test for a distinct gestalt, facial characteristics of the family members and 41 previously published TRIO cases were systematically evaluated via GestaltMatcher. Computational analysis of the facial gestalt suggests a distinguishable facial TRIO-phenotype not outlined in the existing literature.

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/.

Clinical and functional heterogeneity associated with the disruption of retinoic acid receptor beta.

PURPOSE: Dominant variants in the retinoic acid receptor beta (RARB) gene underlie a syndromic form of microphthalmia, known as MCOPS12, which is associated with other birth anomalies and global developmental delay with spasticity and/or dystonia. Here, we report 25 affected individuals with 17 novel pathogenic or likely pathogenic variants in RARB. This study aims to characterize the functional impact of these variants and describe the clinical spectrum of MCOPS12. METHODS: We used in vitro transcriptional assays and in silico structural analysis to assess the functional relevance of RARB variants in affecting the normal response to retinoids. RESULTS: We found that all RARB variants tested in our assays exhibited either a gain-of-function or a loss-of-function activity. Loss-of-function variants disrupted RARB function through a dominant-negative effect, possibly by disrupting ligand binding and/or coactivators' recruitment. By reviewing clinical data from 52 affected individuals, we found that disruption of RARB is associated with a more variable phenotype than initially suspected, with the absence in some individuals of cardinal features of MCOPS12, such as developmental eye anomaly or motor impairment. CONCLUSION: Our study indicates that pathogenic variants in RARB are functionally heterogeneous and associated with extensive clinical heterogeneity.

Novel noncanonical splice site variant causes mild CHD7-related disorder with variable intrafamilial expressivity.

The clinical diagnosis criteria for CHARGE syndrome have been revised several times in the last 25 years. Variable expressivity and reduced penetrance are known, particularly in mild and familial cases. Therefore, it has been proposed to include the detection of a pathogenic CHD7 variant as a major diagnostic criterion. However, intronic variants not located at the canonical splice site are still underrepresented in mutation databases, often because functional analysis is not performed in the diagnostic setting. Here, we report a two-generation family that did not meet the criteria for CHARGE syndrome, until the molecular findings were taken into account. By exome sequencing, we detected an intronic variant in a male individual, who presented with unilateral external ear malformation, bilateral semicircular canal aplasia, polydactyly, vertebral body fusion and a heart defect. The variant was inherited by his mother, who also had bilateral semicircular canal aplasia additionally to unilateral sensorineural hearing impairment, unilateral mandibular palpebral synkinesia, orofacial cleft, and dysphagia. Using RNA studies, we were able to demonstrate that aberrant splicing occurs at an upstream cryptic splice acceptor site, resulting in a frameshift and premature stop of translation. Our data show causality of the noncanonical intronic CHD7 variant and end the diagnostic odyssey of this unsolved phenotype of the family.

Primidone improves symptoms in TRPM3-linked developmental and epileptic encephalopathy with spike-and-wave activation in sleep.

Developmental and epileptic encephalopathy with continuous spike-and-wave activation in sleep (CSWS) or DEE-SWAS is an age-dependent disease, often accompanied by a decline in cognitive abilities. Early successful treatment of CSWS is associated with a better cognitive outcome. We retrospectively analyzed the clinical, electrophysiological, radiological, and genetic data of children with DEE-SWAS associated with melastatin-related transient receptor type 3 gene (TRPM3) missense variants. We report two unrelated children with pharmacoresistant DEE-SWAS and developmental delay/regression and different heterozygous de novo missense variants in the TRPM3 gene (NM_001366145.2; c.3397 T > C/p.Ser1133Pro, c.2004G > A/p.Val1002Met). The variant p.Val1002Met (previously known as p.Val990Met or p.Val837Met) and p.Ser1133Pro were recently shown to result in a gain-of-function effect. Based on this finding, previous drug resistance, and the experimentally demonstrated inhibitory effect of primidone on TRPM3, we initiated an individualized therapy with this drug. In both children, developmental regression was stopped, psychomotor development improved, and CSWS was no longer detectable. To our knowledge, this is the first report of a treatment with primidone in TRPM3-associated CSWS. Our results highlight the importance of early genetic diagnosis in patients with epilepsy and the possibility of precision medicine, which should be considered in the future in individuals with a TRPM3-linked DEE-SWAS.

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.

Recurrent de novo missense variants in GNB2 can cause syndromic intellectual disability.

PURPOSE: Binding proteins (G-proteins) mediate signalling pathways involved in diverse cellular functions and comprise Gα and Gßγ units. Human diseases have been reported for all five Gß proteins. A de novo missense variant in GNB2 was recently reported in one individual with developmental delay/intellectual disability (DD/ID) and dysmorphism. We aim to confirm GNB2 as a neurodevelopmental disease gene, and elucidate the GNB2-associated neurodevelopmental phenotype in a patient cohort. METHODS: We discovered a GNB2 variant in the index case via exome sequencing and sought individuals with GNB2 variants via international data-sharing initiatives. In silico modelling of the variants was assessed, along with multiple lines of evidence in keeping with American College of Medical Genetics and Genomics guidelines for interpretation of sequence variants. RESULTS: We identified 12 unrelated individuals with five de novo missense variants in GNB2, four of which are recurrent: p.(Ala73Thr), p.(Gly77Arg), p.(Lys89Glu) and p.(Lys89Thr). All individuals have DD/ID with variable dysmorphism and extraneurologic features. The variants are located at the universally conserved shared interface with the Gα subunit, which modelling suggests weaken this interaction. CONCLUSION: Missense variants in GNB2 cause a congenital neurodevelopmental disorder with variable syndromic features, broadening the spectrum of multisystem phenotypes associated with variants in genes encoding G-proteins.

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.

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.

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 truncating variants in MAPKAPK5 cause a new developmental disorder involving neurological, cardiac, and facial anomalies combined with synpolydactyly.

PURPOSE: This study aimed to identify the genetic cause of a new multiple congenital anomalies syndrome observed in three individuals from two unrelated families. METHODS: Clinical assessment was conducted prenatally and at different postnatal stages. Genetic studies included exome sequencing (ES) combined with single-nucleotide polymorphism (SNP) array based homozygosity mapping and trio ES. Dermal fibroblasts were used for functional assays. RESULTS: A clinically recognizable syndrome characterized by severe developmental delay, variable brain anomalies, congenital heart defects, dysmorphic facial features, and a distinctive type of synpolydactyly with an additional hypoplastic digit between the fourth and fifth digits of hands and/or feet was identified. Additional features included eye abnormalities, hearing impairment, and electroencephalogram anomalies. ES detected different homozygous truncating variants in MAPKAPK5 in both families. Patient-derived cells showed no expression of MAPKAPK5 protein isoforms and reduced levels of the MAPKAPK5-interacting protein ERK3. F-actin recovery after latrunculin B treatment was found to be less efficient in patient-derived fibroblasts than in control cells, supporting a role of MAPKAPK5 in F-actin polymerization. CONCLUSION: Our data indicate that loss-of-function variants in MAPKAPK5 result in a severe developmental disorder and reveal a major role of this gene in human brain, heart, and limb development.

Clinical and molecular delineation of PUS3-associated neurodevelopmental disorders.

Biallelic variants in PUS3 have recently been recognized as a rare cause of neurodevelopmental disorders. Pseudouridine synthase-3 encoded by PUS3 is an enzyme important for modification of various RNAs, including transfer RNA (tRNA). Here we present the clinical and genetic features of 21 individuals with biallelic PUS3 variants: seven new and 14 previously reported individuals, where clinical features of two were updated. The clinical and genetic information were collected through collaborations or by literature search. All individuals were characterized by the local clinicians and the gene variants were identified by next generation sequencing (NGS) based methodologies. The clinical picture was dominated by global developmental delay, epilepsy, hypotonia and microcephaly. Gray sclera, which has previously been suggested to be a characteristic feature of PUS3-associated phenotypes, was reported in only seven individuals. The patients had some dysmorphic facial features, but a recognizable gestalt was not observed. In conclusion, homozygous and compound heterozygous PUS3 variants lead to a rare neurodevelopmental disorder. Further functional studies are necessary to understand involvement of PUS3 and tRNA biogenesis in normal and abnormal brain development.

Phenotypic expansion of the BPTF-related neurodevelopmental disorder with dysmorphic facies and distal limb anomalies.

Neurodevelopmental disorder with dysmorphic facies and distal limb anomalies (NEDDFL), defined primarily by developmental delay/intellectual disability, speech delay, postnatal microcephaly, and dysmorphic features, is a syndrome resulting from heterozygous variants in the dosage-sensitive bromodomain PHD finger chromatin remodeler transcription factor BPTF gene. To date, only 11 individuals with NEDDFL due to de novo BPTF variants have been described. To expand the NEDDFL phenotypic spectrum, we describe the clinical features in 25 novel individuals with 20 distinct, clinically relevant variants in BPTF, including four individuals with inherited changes in BPTF. In addition to the previously described features, individuals in this cohort exhibited mild brain abnormalities, seizures, scoliosis, and a variety of ophthalmologic complications. These results further support the broad and multi-faceted complications due to haploinsufficiency of BPTF.

TGDS pathogenic variants cause Catel-Manzke syndrome without hyperphalangy.

Catel-Manzke syndrome, also known as micrognathia-digital-syndrome, is a rare autosomal recessive disorder characterized by the combination of the two cardinal features Pierre-Robin sequence and bilateral hyperphalangy leading to ulnar clinodactyly (ulnar curvature of the phalanges) and radial deviation (radial angulation at the metacarpophalangeal joint) of the index fingers. Individuals without one of these major hallmarks or with additional hand malformations have been described as atypical or Catel-Manzke-like syndrome. Biallelic TGDS pathogenic variants have thus far been detected in eight individuals with typical Catel-Manzke syndrome and in one fetus with additional features. Here we report on two individuals with TGDS pathogenic variants who presented with mild radial deviation and ulnar clinodactyly of the index fingers but without radiologic signs of hyperphalangy. Furthermore, both individuals have disproportionate short stature, a feature that has not yet been associated with Catel-Manzke syndrome. Our data broaden the phenotypic spectrum of TGDS-associated Catel-Manzke syndrome and expand the indication for diagnostic testing.

Plasma Kynurenine Predicts Severity and Complications of Heart Failure and Associates with Established Biochemical and Clinical Markers of Disease.

BACKGROUND: Kynurenine, a metabolite of the L-tryptophan pathway, plays a pivotal role in neuro-inflammation, cancer immunology, and cardiovascular inflammation, and has been shown to predict cardiovascular events. OBJECTIVES: It was our objective to increase the body of data regarding the value of kynurenine as a biomarker in chronic heart failure (CHF). METHODS: We investigated the predictive value of plasma kynurenine in a CHF cohort (CHF, n = 114); in a second cohort of defibrillator carriers with CHF (AICD, n = 156), we determined clinical and biochemical determinants of the marker which was measured by enzyme immunoassay. RESULTS: In the CHF cohort, both kynurenine and NT-proBNP increased with NYHA class. Univariate binary logistic regression showed kynurenine to predict death within a 6-month follow-up (OR 1.43, 95% CI 1.03-2.00, p = 0.033) whereas NT-proBNP did not contribute significantly. Kynurenine, like NT-proBNP, was able to discriminate at a 30% threshold of left ventricular ejection fraction (LVEF; AUC-ROC, both 0.74). Kynurenine correlated inversely with LVEF (ϱ = -0.394), glomerular filtration fraction (GFR; ϱ = -0.615), and peak VO2 (ϱ = -0.626). Moreover, there was a strong correlation of kynurenine with NT-proBNP (ϱ = 0.615). In the AICD cohort, multiple linear regression analysis demonstrated highly significant associations of kynurenine with GFR, hsCRP, and tryptophan, as well as a significant impact of age. CONCLUSIONS: This work speaks in favor of kynurenine being a new and valuable biomarker of CHF, with particular attention placed on its ability to predict mortality and reflect exercise capacity.