VarFish: comprehensive DNA variant analysis for diagnostics and research

Whole genome sequencing in families with oligodontia

BACKGROUND/OBJECTIVES

Tooth agenesis (TA) is among the most common malformations in humans. Although several causative mutations have been described, the genetic cause often remains elusive. Here, we test whether whole genome sequencing (WGS) could bridge this diagnostic gap.

METHODS

In four families with TA, we assessed the dental phenotype using the Tooth Agenesis Code after intraoral examination and radiographic and photographic documentation. We performed WGS of index patients and subsequent segregation analysis.

RESULTS

We identified two variants of uncertain significance (a potential splice variant in PTH1R, and a 2.1 kb deletion abrogating a non-coding element in FGF7) and three pathogenic variants: a novel frameshift in the final exon of PITX2, a novel deletion in PAX9, and a known nonsense variant in WNT10A. Notably, the FGF7 variant was found in the patient, also featuring the WNT10A variant. While mutations in PITX2 are known to cause Axenfeld-Rieger syndrome 1 (ARS1) predominantly featuring ocular findings, accompanied by dental malformations, we found the PITX2 frameshift in a family with predominantly dental and varying ocular findings.

CONCLUSION

Severe TA predicts a genetic cause identifiable by WGS. Final exon PITX2 frameshifts can cause a predominantly dental form of ARS1.

Next-generation phenotyping in Nigerian children with Cornelia de Lange syndrome

Next-generation phenotyping (NGP) can be used to compute the similarity of dysmorphic patients to known syndromic diseases. So far, the technology has been evaluated in variant prioritization and classification, providing evidence for pathogenicity if the phenotype matched with other patients with a confirmed molecular diagnosis. In a Nigerian cohort of individuals with facial dysmorphism, we used the NGP tool GestaltMatcher to screen portraits prior to genetic testing and subjected individuals with high similarity scores to exome sequencing (ES). Here, we report on two individuals with global developmental delay, pulmonary artery stenosis, and genital and limb malformations for whom GestaltMatcher yielded Cornelia de Lange syndrome (CdLS) as the top hit. ES revealed a known pathogenic nonsense variant, NM_133433.4: c.598C>T; p.(Gln200*), as well as a novel frameshift variant c.7948dup; p.(Ile2650Asnfs*11) in NIPBL. Our results suggest that NGP can be used as a screening tool and thresholds could be defined for achieving high diagnostic yields in ES. Training the artificial intelligence (AI) with additional cases of the same ethnicity might further increase the positive predictive value of GestaltMatcher.

Next-generation phenotyping integrated in a national framework for patients with ultrarare disorders improves genetic diagnostics and yields new molecular findings

Individuals with ultrarare disorders pose a structural challenge for healthcare systems since expert clinical knowledge is required to establish diagnoses. In TRANSLATE NAMSE, a 3-year prospective study, we evaluated a novel diagnostic concept based on multidisciplinary expertise in Germany. Here we present the systematic investigation of the phenotypic and molecular genetic data of 1,577 patients who had undergone exome sequencing and were partially analyzed with next-generation phenotyping approaches. Molecular genetic diagnoses were established in 32% of the patients totaling 370 distinct molecular genetic causes, most with prevalence below 1:50,000. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were identified, mainly in individuals with neurodevelopmental disorders. Sequencing data of the subcohort that consented to computer-assisted analysis of their facial images with GestaltMatcher could be prioritized more efficiently compared with approaches based solely on clinical features and molecular scores. Our study demonstrates the synergy of using next-generation sequencing and phenotyping for diagnosing ultrarare diseases in routine healthcare and discovering novel etiologies by multidisciplinary teams.

LINE1-mediated epigenetic repression of androgen receptor transcription causes androgen insensitivity syndrome

Androgen insensitivity syndrome (AIS) is a difference of sex development (DSD) characterized by different degrees of undervirilization in individuals with a 46,XY karyotype despite normal to high gonadal testosterone production. Classically, AIS is explained by hemizygous mutations in the X-chromosomal androgen receptor (AR) gene. Nevertheless, the majority of individuals with clinically diagnosed AIS do not carry an AR gene mutation. Here, we present a patient with a 46,XY karyotype, born with undervirilized genitalia, age-appropriate testosterone levels and no uterus, characteristic for AIS. Diagnostic whole exome sequencing (WES) showed a maternally inherited LINE1 (L1) retrotransposon insertion in the 5' untranslated region (5'UTR) of the AR gene. Long-read nanopore sequencing confirmed this as an insertion of a truncated L1 element of ≈ 2.7 kb and showed an increased DNA methylation at the L1 insertion site in patient-derived genital skin fibroblasts (GSFs) compared to healthy controls. The insertion coincided with reduced AR transcript and protein levels in patient-derived GSFs confirming the clinical diagnosis AIS. Our results underline the relevance of retrotransposons in human disease, and expand the growing list of human diseases associated with them.

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.

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.

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.

Enhancing Variant Prioritization in VarFish through On-Premise Computational Facial Analysis

Genomic variant prioritization is crucial for identifying disease-associated genetic variations. Integrating facial and clinical feature analyses into this process enhances performance. This study demonstrates the integration of facial analysis (GestaltMatcher) and Human Phenotype Ontology analysis (CADA) within VarFish, an open-source variant analysis framework. Challenges related to non-open-source components were addressed by providing an open-source version of GestaltMatcher, facilitating on-premise facial analysis to address data privacy concerns. Performance evaluation on 163 patients recruited from a German multi-center study of rare diseases showed PEDIA's superior accuracy in variant prioritization compared to individual scores. This study highlights the importance of further benchmarking and future integration of advanced facial analysis approaches aligned with ACMG guidelines to enhance variant classification.

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.

UNC93B1 variants underlie TLR7-dependent autoimmunity

UNC93B1 is critical for trafficking and function of nucleic acid-sensing Toll-like receptors (TLRs) TLR3, TLR7, TLR8, and TLR9, which are essential for antiviral immunity. Overactive TLR7 signaling induced by recognition of self-nucleic acids has been implicated in systemic lupus erythematosus (SLE). Here, we report UNC93B1 variants (E92G and R336L) in four patients with early-onset SLE. Patient cells or mouse macrophages carrying the UNC93B1 variants produced high amounts of TNF-α and IL-6 and upon stimulation with TLR7/TLR8 agonist, but not with TLR3 or TLR9 agonists. E92G causes UNC93B1 protein instability and reduced interaction with TLR7, leading to selective TLR7 hyperactivation with constitutive type I IFN signaling. Thus, UNC93B1 regulates TLR subtype-specific mechanisms of ligand recognition. Our findings establish a pivotal role for UNC93B1 in TLR7-dependent autoimmunity and highlight the therapeutic potential of targeting TLR7 in SLE.

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.

Case report: KPTN gene-related syndrome associated with a spectrum of neurodevelopmental anomalies including severe epilepsy

Biallelic variants in the kaptin gene KPTN were identified recently in individuals with a novel syndrome referred to as autosomal recessive intellectual developmental disorder 41 (MRT41). MRT41 is characterized by developmental delay, predominantly in language development, behavioral abnormalities, and epilepsy. Only about 15 affected individuals have been described in the literature, all with primary or secondary macrocephaly. Using exome sequencing, we identified three different biallelic variants in KPTN in five affected individuals from three unrelated families. In total, two KPTN variants were already reported as a loss of function variants. A novel splice site variant in KPTN was detected in two unrelated families of this study. The core phenotype with neurodevelopment delay was present in all patients. However, macrocephaly was not present in at least one patient. In total, two patients exhibited developmental and epileptic encephalopathies with generalized tonic-clonic seizures that were drug-resistant in one of them. Thus, we further delineate the KPTN-related syndrome, especially emphasizing the severity of epilepsy phenotypes and difficulties in treatment in patients of our cohort.

SODAR: managing multiomics study data and metadata

Scientists employing omics in life science studies face challenges such as the modeling of multiassay studies, recording of all relevant parameters, and managing many samples with their metadata. They must manage many large files that are the results of the assays or subsequent computation. Users with diverse backgrounds, ranging from computational scientists to wet-lab scientists, have dissimilar needs when it comes to data access, with programmatic interfaces being favored by the former and graphical ones by the latter. We introduce SODAR, the system for omics data access and retrieval. SODAR is a software package that addresses these challenges by providing a web-based graphical user interface for managing multiassay studies and describing them using the ISA (Investigation, Study, Assay) data model and the ISA-Tab file format. Data storage is handled using the iRODS data management system, which handles large quantities of files and substantial amounts of data. SODAR also offers programmable APIs and command-line access for metadata and file storage. SODAR supports complex omics integration studies and can be easily installed. The software is written in Python 3 and freely available at https://github.com/bihealth/sodar-server under the MIT license.

Phenotype-aware prioritisation of rare Mendelian disease variants

A molecular diagnosis from the analysis of sequencing data in rare Mendelian diseases has a huge impact on the management of patients and their families. Numerous patient phenotype-aware variant prioritisation (VP) tools have been developed to help automate this process, and shorten the diagnostic odyssey, but performance statistics on real patient data are limited. Here we identify, assess, and compare the performance of all up-to-date, freely available, and programmatically accessible tools using a whole-exome, retinal disease dataset from 134 individuals with a molecular diagnosis. All tools were able to identify around two-thirds of the genetic diagnoses as the top-ranked candidate, with LIRICAL performing best overall. Finally, we discuss the challenges to overcome most cases remaining undiagnosed after current, state-of-the-art practices.

WARS1 and SARS1: Two tRNA synthetases implicated in autosomal recessive microcephaly

Aminoacylation of transfer RNA (tRNA) is a key step in protein biosynthesis, carried out by highly specific aminoacyl-tRNA synthetases (ARSs). ARSs have been implicated in autosomal dominant and autosomal recessive human disorders. Autosomal dominant variants in tryptophanyl-tRNA synthetase 1 (WARS1) are known to cause distal hereditary motor neuropathy and Charcot-Marie-Tooth disease, but a recessively inherited phenotype is yet to be clearly defined. Seryl-tRNA synthetase 1 (SARS1) has rarely been implicated in an autosomal recessive developmental disorder. Here, we report five individuals with biallelic missense variants in WARS1 or SARS1, who presented with an overlapping phenotype of microcephaly, developmental delay, intellectual disability, and brain anomalies. Structural mapping showed that the SARS1 variant is located directly within the enzyme's active site, most likely diminishing activity, while the WARS1 variant is located in the N-terminal domain. We further characterize the identified WARS1 variant by showing that it negatively impacts protein abundance and is unable to rescue the phenotype of a CRISPR/Cas9 wars1 knockout zebrafish model. In summary, we describe two overlapping autosomal recessive syndromes caused by variants in WARS1 and SARS1, present functional insights into the pathogenesis of the WARS1-related syndrome and define an emerging disease spectrum: ARS-related developmental disorders with or without microcephaly.

Computational approaches for predicting variant impact: An overview from resources, principles to applications

One objective of human genetics is to unveil the variants that contribute to human diseases. With the rapid development and wide use of next-generation sequencing (NGS), massive genomic sequence data have been created, making personal genetic information available. Conventional experimental evidence is critical in establishing the relationship between sequence variants and phenotype but with low efficiency. Due to the lack of comprehensive databases and resources which present clinical and experimental evidence on genotype-phenotype relationship, as well as accumulating variants found from NGS, different computational tools that can predict the impact of the variants on phenotype have been greatly developed to bridge the gap. In this review, we present a brief introduction and discussion about the computational approaches for variant impact prediction. Following an innovative manner, we mainly focus on approaches for non-synonymous variants (nsSNVs) impact prediction and categorize them into six classes. Their underlying rationale and constraints, together with the concerns and remedies raised from comparative studies are discussed. We also present how the predictive approaches employed in different research. Although diverse constraints exist, the computational predictive approaches are indispensable in exploring genotype-phenotype relationship.

Pathogenic Variants in Cardiomyopathy Disorder Genes Underlie Pediatric Myocarditis—Further Impact of Heterozygous Immune Disorder Gene Variants?

Myocarditis is an inflammatory disease of the heart. Pediatric myocarditis with the dilated cardiomyopathy (DCM) phenotype may be caused by likely pathogenic or pathogenic genetic variants [(L)P] in cardiomyopathy (CMP) genes. Systematic analysis of immune disorder gene defects has not been performed so far. We analyzed 12 patients with biopsy-proven myocarditis and the DCM phenotype together with their parents using whole-exome sequencing (WES). The WES data were filtered for rare pathogenic variants in CMP (n = 89) and immune disorder genes (n = 631). Twelve children with a median age of 2.9 (1.0–6.8) years had a mean left ventricular ejection fraction of 28% (22–32%) and myocarditis was confirmed by endomyocardial biopsy. Patients with primary immunodeficiency were excluded from the study. Four patients underwent implantation of a ventricular assist device and subsequent heart transplantation. Genetic analysis of the 12 families revealed an (L)P variant in the CMP gene in 8/12 index patients explaining DCM. Screening of recessive immune disorder genes identified a heterozygous (L)P variant in 3/12 index patients. This study supports the genetic impact of CMP genes for pediatric myocarditis with the DCM phenotype. Piloting the idea that additional immune-related genetic defects promote myocarditis suggests that the presence of heterozygous variants in these genes needs further investigation. Altered cilium function might play an additional role in inducing inflammation in the context of CMP.

Genetic pain loss disorders

Genetic pain loss includes congenital insensitivity to pain (CIP), hereditary sensory neuropathies and, if autonomic nerves are involved, hereditary sensory and autonomic neuropathy (HSAN). This heterogeneous group of disorders highlights the essential role of nociception in protecting against tissue damage. Patients with genetic pain loss have recurrent injuries, burns and poorly healing wounds as disease hallmarks. CIP and HSAN are caused by pathogenic genetic variants in >20 genes that lead to developmental defects, neurodegeneration or altered neuronal excitability of peripheral damage-sensing neurons. These genetic variants lead to hyperactivity of sodium channels, disturbed haem metabolism, altered clathrin-mediated transport and impaired gene regulatory mechanisms affecting epigenetic marks, long non-coding RNAs and repetitive elements. Therapies for pain loss disorders are mainly symptomatic but the first targeted therapies are being tested. Conversely, chronic pain remains one of the greatest unresolved medical challenges, and the genes and mechanisms associated with pain loss offer new targets for analgesics. Given the progress that has been made, the coming years are promising both in terms of targeted treatments for pain loss disorders and the development of innovative pain medicines based on knowledge of these genetic diseases.

Hyper-IgE and Carcinoma in CADINS Disease

Background

Atopic dermatitis (AD) affects up to 25% of children and 10% of adults in Western countries. When severe or recurrent infections and exceedingly elevated serum IgE levels occur in AD patients, an inborn error of immunity (IEI) may be suspected. The International Union of Immunological Societies classification lists variants in different genes responsible for so-called Hyper-IgE syndromes. Diagnosing an underlying IEI may influence treatment strategies.

Methods

Clinical and diagnostic workup of family members are presented including a detailed immunological description and histology of the carcinoma. Functional testing of the novel variant in CARD11 underlying ‘CARD11-associated atopy with dominant interference of NF-kB signaling’ (CADINS) was performed.

Results

We report on an 18-year-old patient with a long-standing history of infections, accompanied by hypogammaglobulinemia, intermittent agranulocytosis, atopy, eosinophilia and colitis. The working diagnosis of common variable immunodeficiency was revised when a novel heterozygous CARD11 variant [c.223C>T; p.(Arg75Trp)] was identified. Functional studies confirmed this variant to have a dominant negative (DN) effect, as previously described in patients with CADINS. Five other family members were affected by severe atopy associated with the above variant, but not hypogammaglobulinemia. Malignancies occurred in two generations: an HPV-positive squamous cell carcinoma and a cutaneous T-cell lymphoma. So far, one patient is under treatment with dupilumab, which has shown marked benefit in controlling severe eczema.

Conclusion

The phenotypic spectrum associated with heterozygous CARD11 DN mutations is broad. Partial T-cell deficiency, diminished IFN-γ cytokine and increased IL-4 production, were identified as disease-causing mechanisms. Malignant disease associated with germline CARD11 DN variants has only been reported sporadically. HPV vaccination in teenage years, and cytology screening analogous with routine cervical swabs may be recommended. Treatment with dupilumab, a monoclonal antibody blocking interleukin-4- and interleukin-13 signaling, may be of benefit in controlling severe and extended AD for some patients as reported for STAT3 loss-of-function.

FAIR Genomes metadata schema promoting Next Generation Sequencing data reuse in Dutch healthcare and research

The genomes of thousands of individuals are profiled within Dutch healthcare and research each year. However, this valuable genomic data, associated clinical data and consent are captured in different ways and stored across many systems and organizations. This makes it difficult to discover rare disease patients, reuse data for personalized medicine and establish research cohorts based on specific parameters. FAIR Genomes aims to enable NGS data reuse by developing metadata standards for the data descriptions needed to FAIRify genomic data while also addressing ELSI issues. We developed a semantic schema of essential data elements harmonized with international FAIR initiatives. The FAIR Genomes schema v1.1 contains 110 elements in 9 modules. It reuses common ontologies such as NCIT, DUO and EDAM, only introducing new terms when necessary. The schema is represented by a YAML file that can be transformed into templates for data entry software (EDC) and programmatic interfaces (JSON, RDF) to ease genomic data sharing in research and healthcare. The schema, documentation and MOLGENIS reference implementation are available at https://fairgenomes.org .

GeneTerpret: a customizable multilayer approach to genomic variant prioritization and interpretation

Background

Variant interpretation is the main bottleneck in medical genomic sequencing efforts. This usually involves genome analysts manually searching through a multitude of independent databases, often with the aid of several, mostly independent, computational tools. To streamline variant interpretation, we developed the GeneTerpret platform which collates data from current interpretation tools and databases, and applies a phenotype-driven query to categorize the variants identified in the genome(s). The platform assigns quantitative validity scores to genes by query and assembly of the genotype–phenotype data, sequence homology, molecular interactions, expression data, and animal models. It also uses the American College of Medical Genetics and Genomics (ACMG) criteria to categorize variants into five tiers of pathogenicity. The final output is a prioritized list of potentially causal variants/genes.

Results

We tested GeneTerpret by comparing its performance to expert-curated genes (ClinGen’s gene-validity database) and variant pathogenicity reports (DECIPHER database). Output from GeneTerpret was 97.2% and 83.5% concordant with the expert-curated sources, respectively. Additionally, similar concordance was observed when GeneTerpret’s performance was compared with our internal expert-interpreted clinical datasets.

Conclusions

GeneTerpret is a flexible platform designed to streamline the genome interpretation process, through a unique interface, with improved ease, speed and accuracy. This modular and customizable system allows the user to tailor the component-programs in the analysis process to their preference. GeneTerpret is available online at https://geneterpret.com.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01166-3.

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.

Cutevariant: a standalone GUI-based desktop application to explore genetic variations from an annotated VCF file

Summary

Cutevariant is a graphical user interface (GUI)-based desktop application designed to filter variations from annotated VCF file. The application imports data into a local SQLite database where complex filter queries can be built either from GUI controllers or using a domain-specific language called Variant Query Language. Cutevariant provides more features than existing applications and is fully customizable thanks to a complete plugins architecture.

Availability and implementation

Cutevariant is distributed as a multiplatform client-side software under an open source license and is available at https://github.com/labsquare/cutevariant.

An Introductory Overview of Open-Source and Commercial Software Options for the Analysis of Forensic Sequencing Data

The top challenges of adopting new methods to forensic DNA analysis in routine laboratories are often the capital investment and the expertise required to implement and validate such methods locally. In the case of next-generation sequencing, in the last decade, several specifically forensic commercial options became available, offering reliable and validated solutions. Despite this, the readily available expertise to analyze, interpret and understand such data is still perceived to be lagging behind. This review gives an introductory overview for the forensic scientists who are at the beginning of their journey with implementing next-generation sequencing locally and because most in the field do not have a bioinformatics background may find it difficult to navigate the new terms and analysis options available. The currently available open-source and commercial software for forensic sequencing data analysis are summarized here to provide an accessible starting point for those fairly new to the forensic application of massively parallel sequencing.

Powerful use of automated prioritization of candidate variants in genetic hearing loss with extreme etiologic heterogeneity

Variant prioritization of exome sequencing (ES) data for molecular diagnosis of sensorineural hearing loss (SNHL) with extreme etiologic heterogeneity poses a significant challenge. This study used an automated variant prioritization system (“EVIDENCE”) to analyze SNHL patient data and assess its diagnostic accuracy. We performed ES of 263 probands manifesting mild to moderate or higher degrees of SNHL. Candidate variants were classified according to the 2015 American College of Medical Genetics guidelines, and we compared the accuracy, call rates, and efficiency of variant prioritizations performed manually by humans or using EVIDENCE. In our in silico panel, 21 synthetic cases were successfully analyzed by EVIDENCE. In our cohort, the ES diagnostic yield for SNHL by manual analysis was 50.19% (132/263) and 50.95% (134/263) by EVIDENCE. EVIDENCE processed ES data 24-fold faster than humans, and the concordant call rate between humans and EVIDENCE was 97.72% (257/263). Additionally, EVIDENCE outperformed human accuracy, especially at discovering causative variants of rare syndromic deafness, whereas flexible interpretations that required predefined specific genotype–phenotype correlations were possible only by manual prioritization. The automated variant prioritization system remarkably facilitated the molecular diagnosis of hearing loss with high accuracy and efficiency, fostering the popularization of molecular genetic diagnosis of SNHL.

Exome Sequencing of 5 Families with Severe Early-Onset Periodontitis

Periodontitis is characterized by alveolar bone loss leading to tooth loss. A small proportion of patients develop severe periodontitis at the juvenile or adolescent age without exposure to the main risk factors of the disease. It is considered that these cases carry rare variants with large causal effects, but the specific variants are largely unknown. In this study, we performed exome sequencing of 5 families with children who developed stage IV, grade C, periodontitis between 3 and 18 y of age. In 1 family, we found compound heterozygous variants in the gene CTSC (p.R272H, p.G139R), 1 of which was previously identified in a family with prepubertal periodontitis. Subsequent targeted resequencing of the CTSC gene in 24 patients <25 y of age (stage IV, grade C) identified the known mutation p.I453V (odds ratio = 4.06, 95% CI = 1.6 to 10.3, P = 0.001), which was previously reported to increase the risk for adolescent periodontitis. An affected sibling of another family carried a homozygous deleterious mutation in the gene TUT7 (p.R560Q, CADD score >30 [Combined Annotation Dependent Depletion]), which is implicated in regulation of interleukin 6 expression. Two other affected siblings shared heterozygous deleterious mutations in the interacting genes PADI1 and FLG (both CADD = 36), which contribute to the integrity of the environment-tissue barrier interface. Additionally, we found predicted deleterious mutations in the periodontitis risk genes ABCA1, GLT6D1, and SIGLEC5. We conclude that the CTSC variants p.R272H and p.I453V have different expressivity and diagnostic relevance for prepubertal and adolescent periodontitis, respectively. We propose additional causal variants for early-onset periodontitis, which also locate within genes that carry known susceptibility variants for common forms. However, the genetic architecture of juvenile periodontitis is complex and differs among the affected siblings of the sequenced families.

Genome sequencing data analysis for rare disease gene discovery

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

One Gene, Many Facets: Multiple Immune Pathway Dysregulation in SOCS1 Haploinsufficiency

Background

Inborn errors of immunity (IEI) present with a large phenotypic spectrum of disease, which can pose diagnostic and therapeutic challenges. Suppressor of cytokine signaling 1 (SOCS1) is a key negative regulator of cytokine signaling, and has recently been associated with a novel IEI. Of patients described to date, it is apparent that SOCS1 haploinsufficiency has a pleiotropic effect in humans.

Objective

We sought to investigate whether dysregulation of immune pathways, in addition to STAT1, play a role in the broad clinical manifestations of SOCS1 haploinsufficiency.

Methods

We assessed impacts of reduced SOCS1 expression across multiple immune cell pathways utilizing patient cells and CRISPR/Cas9 edited primary human T cells.

Results

SOCS1 haploinsufficiency phenotypes straddled across the International Union of Immunological Societies classifications of IEI. We found that reduced SOCS1 expression led to dysregulation of multiple intracellular pathways in immune cells. STAT1 phosphorylation is enhanced, comparably with STAT1 gain-of-function mutations, and STAT3 phosphorylation is similarly reduced with concurrent reduction of Th17 cells. Furthermore, reduced SOCS1 E3 ligase function was associated with increased FAK1 in immune cells, and increased AKT and p70 ribosomal protein S6 kinase phosphorylation. We also found Toll-like receptor responses are increased in SOCS1 haploinsufficiency patients.

Conclusions

SOCS1 haploinsufficiency is a pleiotropic monogenic IEI. Dysregulation of multiple immune cell pathways may explain the variable clinical phenotype associated with this new condition. Knowledge of these additional dysregulated immune pathways is important when considering the optimum management for SOCS1 haploinsufficient patients.

MutationTaster2021

Here we present an update to MutationTaster, our DNA variant effect prediction tool. The new version uses a different prediction model and attains higher accuracy than its predecessor, especially for rare benign variants. In addition, we have integrated many sources of data that only became available after the last release (such as gnomAD and ExAC pLI scores) and changed the splice site prediction model. To more easily assess the relevance of detected known disease mutations to the clinical phenotype of the patient, MutationTaster now provides information on the diseases they cause. Further changes represent a major overhaul of the interfaces to increase user-friendliness whilst many changes under the hood have been designed to accelerate the processing of uploaded VCF files. We also offer an API for the rapid automated query of smaller numbers of variants from within other software. MutationTaster2021 integrates our disease mutation search engine, MutationDistiller, to prioritise variants from VCF files using the patient's clinical phenotype. The novel version is available at https://www.genecascade.org/MutationTaster2021/. This website is free and open to all users and there is no login requirement.

Pathogenic Variants Associated With Dilated Cardiomyopathy Predict Outcome in Pediatric Myocarditis

BACKGROUND

Myocarditis is one of the most common causes leading to heart failure in children and a possible genetic background has been postulated. We sought to characterize the clinical and genetic characteristics in patients with myocarditis ≤18 years of age to predict outcome.

METHODS

A cohort of 42 patients (Genetics in Pediatric Myocarditis) with biopsy-proven myocarditis underwent genetic testing with targeted panel sequencing of cardiomyopathy-associated genes. Genetics in Pediatric Myocarditis patients were divided into subgroups according to the phenotype of dilated cardiomyopathy (DCM) at presentation, resulting in 22 patients without DCM (myocarditis without phenotype of DCM) and 20 patients with DCM (myocarditis with phenotype of DCM).

RESULTS

Myocarditis with phenotype of DCM patients (median age 1.4 years) were younger than myocarditis without phenotype of DCM patients (median age 16.1 years; P<0.001) and were corresponding to heart failure-like and coronary syndrome-like phenotypes, respectively. At least one likely pathogenic/pathogenic variant was identified in 9 out of 42 patients (22%), 8 of them were heterozygous, and 7 out of 9 were in myocarditis with phenotype of DCM. Likely pathogenic/pathogenic variants were found in genes validated for primary DCM (BAG3, DSP, LMNA, MYH7, TNNI3, TNNT2, and TTN). Rare variant enrichment analysis revealed significant accumulation of high-impact disease variants in myocarditis with phenotype of DCM versus healthy individuals (P=0.0003). Event-free survival was lower (P=0.008) in myocarditis with phenotype of DCM patients compared with myocarditis without phenotype of DCM and primary DCM.

CONCLUSIONS

We report heterozygous likely pathogenic/pathogenic variants in biopsy-proven pediatric myocarditis. Myocarditis patients with DCM phenotype were characterized by early-onset heart failure, significant enrichment of likely pathogenic/pathogenic variants, and poor outcome. These phenotype-specific and age group-specific findings will be useful for personalized management of these patients. Genetic evaluation in children newly diagnosed with myocarditis and DCM phenotype is warranted.

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.

Genome sequencing in families with congenital limb malformations

The extensive clinical and genetic heterogeneity of congenital limb malformation calls for comprehensive genome-wide analysis of genetic variation. Genome sequencing (GS) has the potential to identify all genetic variants. Here we aim to determine the diagnostic potential of GS as a comprehensive one-test-for-all strategy in a cohort of undiagnosed patients with congenital limb malformations. We collected 69 cases (64 trios, 1 duo, 5 singletons) with congenital limb malformations with no molecular diagnosis after standard clinical genetic testing and performed genome sequencing. We also developed a framework to identify potential noncoding pathogenic variants. We identified likely pathogenic/disease-associated variants in 12 cases (17.4%) including four in known disease genes, and one repeat expansion in HOXD13. In three unrelated cases with ectrodactyly, we identified likely pathogenic variants in UBA2, establishing it as a novel disease gene. In addition, we found two complex structural variants (3%). We also identified likely causative variants in three novel high confidence candidate genes. We were not able to identify any noncoding variants. GS is a powerful strategy to identify all types of genomic variants associated with congenital limb malformation, including repeat expansions and complex structural variants missed by standard diagnostic approaches. In this cohort, no causative noncoding SNVs could be identified.

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.

Novel sequencing technologies and bioinformatic tools for deciphering the non-coding genome

High-throughput sequencing techniques have significantly increased the molecular diagnosis rate for patients with monogenic disorders. This is primarily due to a substantially increased identification rate of disease mutations in the coding sequence, primarily SNVs and indels. Further progress is hampered by difficulties in the detection of structural variants and the interpretation of variants outside the coding sequence. In this review, we provide an overview about how novel sequencing techniques and state-of-the-art algorithms can be used to discover small and structural variants across the whole genome and introduce bioinformatic tools for the prediction of effects variants may have in the non-coding part of the genome.

Compound Heterozygous Frameshift Mutations in MESD Cause a Lethal Syndrome Suggestive of Osteogenesis Imperfecta Type XX

Multiple genes are known to be associated with osteogenesis imperfecta (OI), a phenotypically and genetically heterogenous bone disorder, marked predominantly by low bone mineral density and increased risk of fractures. Recently, mutations affecting MESD, which encodes for a chaperone required for trafficking of the low-density lipoprotein receptors LRP5 and LRP6 in the endoplasmic reticulum, were described to cause autosomal-recessive OI XX in homozygous children. In the present study, whole-exome sequencing of three stillbirths in one family was performed to evaluate the presence of a hereditary disorder. To further characterize the skeletal phenotype, fetal autopsy, bone histology, and quantitative backscattered electron imaging (qBEI) were performed, and the results were compared with those from an age-matched control with regular skeletal phenotype. In each of the affected individuals, compound heterozygous mutations in MESD exon 2 and exon 3 were detected. Based on the skeletal phenotype, which was characterized by multiple intrauterine fractures and severe skeletal deformity, OI XX was diagnosed in these individuals. Histological evaluation of MESD specimens revealed an impaired osseous development with an altered osteocyte morphology and reduced canalicular connectivity. Moreover, analysis of bone mineral density distribution by qBEI indicated an impaired and more heterogeneous matrix mineralization in individuals with MESD mutations than in controls. In contrast to the previously reported phenotypes of individuals with OI XX, the more severe phenotype in the present study is likely explained by a mutation in exon 2, located within the chaperone domain of MESD, that leads to a complete loss of function, which indicates the relevance of MESD in early skeletal development. © 2021 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR)..

Computational resources for identification of cancer biomarkers from omics data

Cancer is one of the most prevailing, deadly and challenging diseases worldwide. The advancement in technology led to the generation of different types of omics data at each genome level that may potentially improve the current status of cancer patients. These data have tremendous applications in managing cancer effectively with improved outcome in patients. This review summarizes the various computational resources and tools housing several types of omics data related to cancer. Major categorization of resources includes-cancer-associated multiomics data repositories, visualization/analysis tools for omics data, machine learning-based diagnostic, prognostic, and predictive biomarker tools, and data analysis algorithms employing the multiomics data. The review primarily focuses on providing comprehensive information on the open-source multiomics tools and data repositories, owing to their broader applicability, economic-benefit and usability. Sections including the comparative analysis, tools applicability and possible future directions have also been discussed in detail. We hope that this information will significantly benefit the researchers and clinicians, especially those with no sound background in bioinformatics and who lack sufficient data analysis skills to interpret something from the plethora of cancer-specific data generated nowadays.

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.