Search

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.

Most individuals with rare diseases initially consult their primary care physician. For a subset of rare diseases, efficient diagnostic pathways are available. However, ultra-rare diseases often require both expert clinical knowledge and comprehensive genetic diagnostics, which poses structural challenges for public healthcare systems. To address these challenges within Germany, a novel structured diagnostic concept, based on multidisciplinary expertise at established university hospital centers for rare diseases (CRDs), was evaluated in the three year prospective study TRANSLATE NAMSE. A key goal of TRANSLATE NAMSE was to assess the clinical value of exome sequencing (ES) in the ultra-rare disease population. The aims of the present study were to perform a systematic investigation of the phenotypic and molecular genetic data of TRANSLATE NAMSE patients who had undergone ES in order to determine the yield of both ultra-rare diagnoses and novel gene-disease associations; and determine whether the complementary use of machine learning and artificial intelligence (AI) tools improved diagnostic effectiveness and efficiency.ES was performed for 1,577 patients (268 adult and 1,309 pediatric). Molecular genetic diagnoses were established in 499 patients (74 adult and 425 pediatric). A total of 370 distinct molecular genetic causes were established. The majority of these concerned known disorders, most of which were ultra-rare. During the diagnostic process, 34 novel and 23 candidate genotype-phenotype associations were delineated, mainly in individuals with neurodevelopmental disorders.To determine the likelihood that ES will lead to a molecular diagnosis in a given patient, based on the respective clinical features only, we developed a statistical framework called YieldPred. The genetic data of a subcohort of 224 individuals that also gave consent to the computer-assisted analysis of their facial images were processed with the AI tool Prioritization of Exome Data by Image Analysis (PEDIA) and showed superior performance in variant prioritization.The present analyses demonstrated that the novel structured diagnostic concept facilitated the identification of ultra-rare genetic disorders and novel gene-disease associations on a national level and that the machine learning and AI tools improved diagnostic effectiveness and efficiency for ultra-rare genetic disorders.

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.

Copy Number Variants (CNVs) are deletions, duplications or insertions larger than 50 base pairs. They account for a large percentage of the normal genome variation and play major roles in human pathology. While array-based approaches have long been used to detect them in clinical practice, whole-genome sequencing (WGS) bears the promise to allow concomitant exploration of CNVs and smaller variants. However, accurately calling CNVs from WGS remains a difficult computational task, for which a consensus is still lacking. In this paper, we explore practical calling options to reach the best compromise between sensitivity and sensibility. We show that callers based on different signal (paired-end reads, split reads, coverage depth) yield complementary results. We suggest approaches combining four selected callers (Manta, Delly, ERDS, CNVnator) and a regenotyping tool (SV2), and show that this is applicable in everyday practice in terms of computation time and further interpretation. We demonstrate the superiority of these approaches over array-based Comparative Genomic Hybridization (aCGH), specifically regarding the lack of resolution in breakpoint definition and the detection of potentially relevant CNVs. Finally, we confirm our results on the NA12878 benchmark genome, as well as one clinically validated sample. In conclusion, we suggest that WGS constitutes a timely and economically valid alternative to the combination of aCGH and whole-exome sequencing.

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.

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. Polish National Science Centre [UMO-2016/22/E/NZ5/00270]; Deutsche Forschungsgemeinschaft (DFG)German Research Foundation (DFG) [SP1532/3-1, SP1532/4-1, SP1532/5-1]; Max Planck FoundationFoundation CELLEX; Deutsches Zentrum fur Luft-und Raumfahrt (DLR)Helmholtz AssociationGerman Aerospace Centre (DLR) [01GM1925]; Projekt DEAL Open Access funding enabled and organized by Projekt DEAL. A.J. was supported by the grant from the Polish National Science Centre UMO-2016/22/E/NZ5/00270. M.S. is supported by grants from the Deutsche Forschungsgemeinschaft (DFG) (SP1532/31, SP1532/4-1 and SP1532/5-1), the Max Planck Foundation and the Deutsches Zentrum fur Luft-und Raumfahrt (DLR 01GM1925).

BackgroundGenes implicated in the Golgi and endosomal trafficking machinery are crucial for brain development, and mutations in them are particularly associated with postnatal microcephaly (POM).MethodsExome 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.ResultsWe 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.ConclusionIn 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.

Purpose To expand the recent description of a new neurodevelopmental syndrome related to alterations in CDK19. Methods Individuals were identified through international collaboration. Functional studies included autophosphorylation assays for CDK19 Gly28Arg and Tyr32His variants and in vivo zebrafish assays of the CDK19(G28R) and CDK19(Y32H). Results We describe 11 unrelated individuals (age range: 9 months to 14 years) with de novo missense variants mapped to the kinase domain of CDK19, including two recurrent changes at residues Tyr32 and Gly28. In vitro autophosphorylation and substrate phosphorylation assays revealed that kinase activity of protein was lower for p.Gly28Arg and higher for p.Tyr32His substitutions compared with that of the wild-type protein. Injection of CDK19 messenger RNA (mRNA) with either the Tyr32His or the Gly28Arg variants using in vivo zebrafish model significantly increased fraction of embryos with morphological abnormalities. Overall, the phenotype of the now 14 individuals with CDK19-related disorder includes universal developmental delay and facial dysmorphism, hypotonia (79%), seizures (64%), ophthalmologic anomalies (64%), and autism/autistic traits (56%). Conclusion CDK19 de novo missense variants are responsible for a novel neurodevelopmental disorder. Both kinase assay and zebrafish experiments showed that the pathogenetic mechanism may be more diverse than previously thought.

Scientists employing omics in life science studies face challenges such as the modeling of multi assay 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 multi assay 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 athttps://github.com/bihealth/sodar-serverunder the MIT license.

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is a highly heritable facial disorder. To date, systematic investigations of the contribution of rare variants in non-coding regions to nsCL/P etiology are sparse. Here, we re-analyzed available whole-genome sequence (WGS) data from 211 European case-parent trios with nsCL/P and identified 13,522

KCNN2 encodes the small conductance calcium-activated potassium channel 2 (SK2). Rodent models with spontaneous Kcnn2 mutations show abnormal gait and locomotor activity, tremor and memory deficits, but human disorders related to KCNN2 variants are largely unknown. Using exome sequencing, we identified a de novo KCNN2 frameshift deletion in a patient with learning disabilities, cerebellar ataxia and white matter abnormalities on brain MRI. This discovery prompted us to collect data from nine additional patients with de novo KCNN2 variants (one nonsense, one splice site, six missense variants and one in-frame deletion) and one family with a missense variant inherited from the affected mother. We investigated the functional impact of six selected variants on SK2 channel function using the patch-clamp technique. All variants tested but one, which was reclassified to uncertain significance, led to a loss-of-function of SK2 channels. Patients with KCNN2 variants had motor and language developmental delay, intellectual disability often associated with early-onset movement disorders comprising cerebellar ataxia and/or extrapyramidal symptoms. Altogether, our findings provide evidence that heterozygous variants, likely causing a haploinsufficiency of the KCNN2 gene, lead to novel autosomal dominant neurodevelopmental movement disorders mirroring phenotypes previously described in rodents.

Human Phenotype Ontology (HPO)-based analysis has become standard for genomic diagnostics of rare diseases. Current algorithms use a variety of semantic and statistical approaches to prioritize the typically long lists of genes with candidate pathogenic variants. These algorithms do not provide robust estimates of the strength of the predictions beyond the placement in a ranked list, nor do they provide measures of how much any individual phenotypic observation has contributed to the prioritization result. However, given that the overall success rate of genomic diagnostics is only around 25%–50% or less in many cohorts, a good ranking cannot be taken to imply that the gene or disease at rank one is necessarily a good candidate. Here, we present an approach to genomic diagnostics that exploits the likelihood ratio (LR) framework to provide an estimate of (1) the posttest probability of candidate diagnoses, (2) the LR for each observed HPO phenotype, and (3) the predicted pathogenicity of observed genotypes. LIkelihood Ratio Interpretation of Clinical AbnormaLities (LIRICAL) placed the correct diagnosis within the first three ranks in 92.9% of 384 case reports comprising 262 Mendelian diseases, and the correct diagnosis had a mean posttest probability of 67.3%. Simulations show that LIRICAL is robust to many typically encountered forms of genomic and phenomic noise. In summary, LIRICAL provides accurate, clinically interpretable results for phenotype-driven genomic diagnostics.

Hand hyperphalangism leading to shortened index fingers with ulnar deviation, hallux valgus, mild facial dysmorphism and respiratory compromise requiring assisted ventilation are the key features of Chitayat syndrome. This condition results from the recurrent heterozygous missense variant NM_006494.2:c.266A>G; p.(Tyr89Cys) in ERF on chromosome 19q13.2, encoding the ETS2 repressor factor (ERF) protein. The pathomechanism of Chitayat syndrome is unknown. To date, seven individuals with Chitayat syndrome and the recurrent pathogenic ERF variant have been reported in the literature. Here, we describe six additional individuals, among them only one presenting with a history of assisted ventilation, and the remaining presenting with variable pulmonary phenotypes, including one individual without any obvious pulmonary manifestations. Our findings widen the phenotype spectrum caused by the recurrent pathogenic variant in ERF, underline Chitayat syndrome as a cause of isolated skeletal malformations and therefore contribute to the improvement of diagnostic strategies in individuals with hand hyperphalangism.

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.

Background Paroxysmal tonic upgaze (PTU), defined as an involuntary upward movement of the eyes, has been considered as a benign phenomenon but may also be associated with ataxia and developmental delay. Methods We report eight children with PTU; six of them also exhibiting symptoms of ataxia and/or developmental delay. Treatment with carbonic anhydrase inhibition was offered to children with persisting and/or severe forms. Results Whole-exome sequencing and genome-wide array analysis (n = 7) did not reveal mutations in the three known genes associated with PTU (CACNA1A, GRID2, SEPSECS), whereas by MLPA a heterozygous deletion of exon 31 of the CACNA1A gene could be detected in one patient, her mother and two further family members. Further exome and array analysis showed no recurrent variants in potentially novel PTU-related genes in more than one patient. A de novo variant at a highly conserved position in the SIM1 gene was detected in one patient, for which a pathogenic effect could be speculated. Carbonic anhydrase inhibition was started in five children and proved at least partially effective in all of them. Conclusion Irrespective of the clinical background and the molecular basic mechanism of PTU, therapeutic carbonic anhydrase inhibition was effective in all five children (acetazolamide, n = 3; sultiame, n = 2) who received this treatment.

Supplemental Digital Content is available in the text., 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

Background 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. Affected individuals and methods 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. Results Molecular studies revealed 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 and a delayed Brefeldin A-induced retrograde transport was shown in two and glycosylation abnormalities were present in fibroblasts from three individuals. Conclusion 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. (words: 260) This article is protected by copyright. All rights reserved.

Purpose: Haploinsufficiency of USP7, located at chromosome 16p13.2, has recently been reported in seven individuals with neurodevelopmental phenotypes, including developmental delay/intellectual disability (DD/ID), autism spectrum disorder (ASD), seizures, and hypogonadism. Further, USP7 was identified to critically incorporate into the MAGEL2-USP7-TRIM27 (MUST), such that pathogenic variants in USP7 lead to altered endosomal F-actin polymerization and dysregulated protein recycling. Methods: We report 16 newly identified individuals with heterozygous USP7 variants, identified by genome or exome sequencing or by chromosome microarray analysis. Clinical features were evaluated by review of medical records. Additional clinical information was obtained on the seven previously reported individuals to fully elucidate the phenotypic expression associated with USP7 haploinsufficiency. Results: The clinical manifestations of these 23 individuals suggest a syndrome characterized by DD/ID, hypotonia, eye anomalies,feeding difficulties, GERD, behavioral anomalies, and ASD, and more specific phenotypes of speech delays including a nonverbal phenotype and abnormal brain magnetic resonance image findings including white matter changes based on neuroradiologic examination. Conclusion: The consistency of clinical features among all individuals presented regardless of de novo USP7 variant type supports haploinsufficiency as a mechanism for pathogenesis and refines the clinical impact faced by affected individuals and caregivers.

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).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.We detected homozygous truncating variants inIn aggregate, our findings show that pathogenic variants in

Loss of function variants of GLI3 are associated with a variety of forms of polysyndactyly: Pallister-Hall syndrome (PHS), Greig-Cephalopolysyndactyly syndrome (GCPS), and isolated polysyndactyly (IPD). Variants affecting the N-terminal and C-terminal thirds of the GLI3 protein have been associated with GCPS, those within the central third with PHS. Cases of IPD have been attributed to variants affecting the C-terminal third of the GLI3 protein. In this study, we further investigate these genotype-phenotype correlations. Sequencing of GLI3 was performed in patients with clinical findings suggestive of a GLI3-associated syndrome. Additionally, we searched the literature for reported cases of either manifestation with mutations in the GLI3 gene. Here, we report 48 novel cases from 16 families with polysyndactyly in whom we found causative variants in GLI3 and a review on 314 previously reported GLI3 variants. No differences in location of variants causing either GCPS or IPD were found. Review of published data confirmed the association of PHS and variants affecting the GLI3 protein's central third. We conclude that the observed manifestations of GLI3 variants as GCPS or IPD display different phenotypic severities of the same disorder and propose a binary division of GLI3-associated disorders in either PHS or GCPS/polysyndactyly.

BackgroundSince the beginning of the coronavirus disease 2019 (COVID-19) pandemic, there has been increasing demand to identify predictors of severe clinical course in patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Human leukocyte antigen alleles (HLA) have been suggested as potential genetic host factors. We sought to evaluate this hypothesis by conducting an international multicenter study using HLA sequencing with subsequent independent validation.MethodsWe analyzed a total of 332 samples. First, we enrolled 233 patients in Germany, Spain, and Switzerland for HLA and whole exome sequencing. Furthermore, we validated our results in a public data set (United States, n=99). Patients older than 18 years presenting with COVID-19 were included, representing the full spectrum of the disease. HLA candidate alleles were identified in the derivation cohort (n=92) and tested in two independent validation cohorts (n=240).ResultsWe identified HLA-C* 04:01 as a novel genetic predictor for severe clinical course in COVID-19. Carriers of HLA-C* 04:01 had twice the risk of intubation when infected with SARS-CoV-2 (hazard ratio 2.1, adjusted p-value=0.0036). Importantly, these findings were successfully replicated in an independent data set. Furthermore, our findings are biologically plausible, as HLA-C* 04:01 has fewer predicted bindings sites with relevant SARS-CoV-2 peptides as compared to other HLA alleles. Exome sequencing confirmed findings from HLA analysis.ConclusionsHLA-C* 04:01 carriage is associated with a twofold increased risk of intubation in patients infected with SARS-CoV-2. Testing for HLA-C* 04:01 could have clinical implications to identify high-risk patients and individualize management.

Huntington's disease and subcortical vascular dementia display similar dementing features, shaped by different degrees of striatal atrophy, deep white matter degeneration and tau pathology. To investigate the hypothesis that Huntington's disease transcriptomic hallmarks may provide a window into potential protective genes upregulated during brain acute and subacute ischemia, we compared RNA sequencing signatures in the most affected brain areas of 2 widely used experimental mouse models: Huntington's disease, (R6/2, striatum and cortex and Q175, hippocampus) and brain ischemia-subcortical vascular dementia (BCCAS, striatum, cortex and hippocampus). We identified a cluster of 55 shared genes significantly differentially regulated in both models and we screened these in 2 different mouse models of Alzheimer's disease, and 96 early-onset familial and apparently sporadic small vessel ischemic disease patients. Our data support the prevalent role of transcriptional regulation upon genetic coding variability of known neuroprotective genes (Egr2, Fos, Ptgs2, Itga5, Cdkn1a, Gsn, Npas4, Btg2, Cebpb) and provide a list of potential additional ones likely implicated in different dementing disorders and worth further investigation.

VarFish is a user-friendly web application for the quality control, filtering, prioritization, analysis, and user-based annotation of panel and exome variant data for 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. Filtered variants 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 user 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.

Human Phenotype Ontology (HPO)-based analysis has become standard for genomic diagnostics of rare diseases. Current algorithms use a variety of semantic and statistical approaches to prioritize the typically long lists of genes with candidate pathogenic variants. These algorithms do not provide robust estimates of the strength of the predictions beyond the placement in a ranked list, nor do they provide measures of how much any individual phenotypic observation has contributed to the prioritization result. However, given that the overall success rate of genomic diagnostics is only around 25–50% or less in many cohorts, a good ranking cannot be taken to imply that the gene or disease at rank one is necessarily a good candidate. Likelihood ratios (LR) are statistics for summarizing diagnostic accuracy, providing a measure of how much more (or less) a patient with a disease has a particular test result compared to patients without the disease. Here, we present an approach to genomic diagnostics that exploits the LR framework to provide an estimate of (1) the posttest probability of candidate diagnoses; (2) the LR for each observed HPO phenotype, and (3) the predicted pathogenicity of observed genotypes. LIkelihood Ratio Interpretation of Clinical AbnormaLities (LIRICAL) placed the correct diagnosis within the first three ranks in 92.9% of 384 cases reports comprising 262 Mendelian diseases, with the correct diagnosis having a mean posttest probability of 67.3%. Simulations show that LIRICAL is robust to many typically encountered forms of genomic and phenomic noise. In summary, LIRICAL provides accurate, clinically interpretable results for phenotype-driven genomic diagnostics.

A genome-wide evaluation of the effects of ionizing radiation on mutation induction in the mouse germline has identified multisite de novo mutations (MSDNs) as marker for previous exposure. Here we present the results of a small pilot study of whole genome sequencing in offspring of soldiers who served in radar units on weapon systems that were emitting high-frequency radiation. We found cases of exceptionally high MSDN rates as well as an increased mean in our cohort: While a MSDN mutation is detected in average in 1 out of 5 offspring of unexposed controls, we observed 12 MSDNs in altogether 18 offspring, including a family with 6 MSDNs in 3 offspring. Moreover, we found two translocations, also resulting from neighboring mutations. Our findings indicate that MSDNs might be suited in principle for the assessment of DNA damage from ionizing radiation also in humans. However, as exact person-related dose values in risk groups are usually not available, the interpretation of MSDNs in single families would benefit from larger molecular epidemiologic studies on this new biomarker.

BackgroundSingle cell omics technologies present unique opportunities for biomedical and life sciences from lab to clinic, but the high dimensional nature of such data poses challenges for computational analysis and interpretation. Furthermore, FAIR data management as well as data privacy and security become crucial when working with clinical data, especially in cross-institutional and translational settings. Existing solutions are either bound to the desktop of one researcher or come with dependencies on vendor-specific technology for cloud storage or user authentication.ResultsTo facilitate analysis and interpretation of single-cell data by users without bioinformatics expertise, we present SCelVis, a flexible, interactive and user-friendly app for web-based visualization of pre-processed single-cell data. Users can survey multiple interactive visualizations of their single cell expression data and cell annotation, and download raw or processed data for further offline analysis. SCelVis can be run both on the desktop and cloud systems, accepts input from local and various remote sources using standard and open protocols, and allows for hosting data in the cloud and locally.MethodsSCelVis is implemented in Python using Dash by Plotly. It is available as a standalone application as a Python package, via Conda/Bioconda and as a Docker image. All components are available as open source under the permissive MIT license and are based on open standards and interfaces, enabling further development and integration with third party pipelines and analysis components. The GitHub repository is https://github.com/bihealth/scelvis.

Summary. Managing raw sequencing data and conversion into sequences (demultiplexing) remains a challenging topic for groups running sequencing devices. They face many challenges in such efforts and solutions range from manual management of spreadsheets to very complex and customized LIMS systems handling much more than just sequencing raw data. In this manuscript, we describe the software package DigestiFlow that focuses on the management of Illumina flow cell sample sheets and raw data. Namely, it allows for automated extraction of flow cell raw data information, management of sample sheets, and the automated (and thus reproducible) demultiplexing of Illumina base calls data. Availability and Implementation. The software is available under the MIT license at https://github.com/bihealth/digestiflow-server. The client and demux software components are available via Bioconda.

The underlying genetic mechanisms and early pathological events of children with primary cardiomyopathy (CMP) are insufficiently characterized. In this study, we aimed to characterize the mutational spectrum of primary CMP in a large cohort of patients ���18 years referred to a tertiary center. Eighty unrelated index patients with pediatric primary CMP underwent genetic testing with a panel-based next-generation sequencing approach of 89 genes. At least one pathogenic or probably pathogenic variant was identified in 30/80 (38%) index patients. In all CMP subgroups, patients carried most frequently variants of interest in sarcomere genes suggesting them as a major contributor in pediatric primary CMP. In MYH7, MYBPC3, and TNNI3, we identified 18 pathogenic/probably pathogenic variants (MYH7 n = 7, MYBPC3 n = 6, TNNI3 n = 5, including one homozygous (TNNI3 c.24+2T>A) truncating variant. Protein and transcript level analysis on heart biopsies from individuals with homozygous mutation of TNNI3 revealed that the TNNI3 protein is absent and associated with upregulation of the fetal isoform TNNI1. The present study further supports the clinical importance of sarcomeric mutation-not only in adult-but also in pediatric primary CMP. TNNI3 is the third most important disease gene in this cohort and complete loss of TNNI3 leads to severe pediatric CMP.

Data analysis often entails a multitude of heterogeneous steps, from the application of various command line tools to the usage of scripting languages like R or Python for the generation of plots and tables. It is widely recognized that data analyses should ideally be conducted in a reproducible way. Reproducibility enables technical validation and regeneration of results on the original or even new data. However, reproducibility alone is by no means sufficient to deliver an analysis that is of lasting impact (i.e. sustainable) for the field, or even just one research group. We postulate that it is equally important to ensure adaptability and transparency. The former describes the ability to modify the analysis to answer extended or slightly different research questions. The latter describes the ability to understand the analysis in order to judge whether it is not only technically, but methodologically valid. Here, we analyze the properties needed for a data analysis to become reproducible, adaptable, and transparent, and show how the popular workflow management system Snakemake can be used to fulfill all these needs., This work was supported by the Netherlands Organisation for Scientific Research (NWO) (VENI grant 016.Veni.173.076, Johannes Köster), the German Research Foundation (SFB 876, Johannes Köster), and Google LLC (Vanessa Sochat and Johannes Köster).

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.

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.

Immune escape is one of the hallmarks of cancer and several new treatment approaches attempt to modulate and restore the immune system’s capability to target cancer cells. At the heart of the immune recognition process lies antigen presentation from somatic mutations. These neo-epitopes are emerging as attractive targets for cancer immunotherapy and new strategies for rapid identification of relevant candidates have become a priority. We carefully screen TCGA data sets for recurrent somatic amino acid exchanges and apply MHC class I binding predictions. We propose a method for in silico selection and prioritization of candidates which have a high potential for neo-antigen generation and are likely to appear in multiple patients. While the percentage of patients carrying a specific neo-epitope and HLA-type combination is relatively small, the sheer number of new patients leads to surprisingly high reoccurence numbers. We identify 769 epitopes which are expected to occur in 77629 patients per year. While our candidate list will definitely contain false positives, the results provide an objective order for wet-lab testing of reusable neo-epitopes. Thus recurrent neo-epitopes may be suitable to supplement existing personalized T cell treatment approaches with precision treatment options.

Exomiser is an application that prioritizes genes and variants in next-generation sequencing (NGS) projects for novel disease-gene discovery or differential diagnostics of Mendelian disease. Exomiser comprises a suite of algorithms for prioritizing exome sequences using random-walk analysis of protein interaction networks, clinical relevance and cross-species phenotype comparisons, as well as a wide range of other computational filters for variant frequency, predicted pathogenicity and pedigree analysis. In this protocol, we provide a detailed explanation of how to install Exomiser and use it to prioritize exome sequences in a number of scenarios. Exomiser requires ~3 GB of RAM and roughly 15–90 s of computing time on a standard desktop computer to analyze a variant call format (VCF) file. Exomiser is freely available for academic use from http://www.sanger.ac.uk/science/tools/exomiser.

Motivation: Large insertions of novel sequence are an important type of structural variants. Previous studies used traditional de novo assemblers for assembling non-mapping high-throughput sequencing (HTS) or capillary reads and then tried to anchor them in the reference using paired read information. Results: We present approaches for detecting insertion breakpoints and targeted assembly of large insertions from HTS paired data: BASIL and ANISE. On near identity repeats that are hard for assemblers, ANISE employs a repeat resolution step. This results in far better reconstructions than obtained by the compared methods. On simulated data, we found our insert assembler to be competitive with the de novo assemblers ABYSS and SGA while yielding already anchored inserted sequence as opposed to unanchored contigs as from ABYSS/SGA. On real-world data, we detected novel sequence in a human individual and thoroughly validated the assembled sequence. ANISE was found to be superior to the competing tool MindTheGap on both simulated and real-world data. Availability and implementation: ANISE and BASIL are available for download at http://www.seqan.de/projects/herbarium under a permissive open source license. Contact: manuel.holtgrewe@fu-berlin.de or knut.reinert@fu-berlin.de Supplementary information: Supplementary data are available at Bioinformatics online.

We present Bioconda (https://bioconda.github.io), a distribution of bioinformatics software for the lightweight, multiplatform and language-agnostic package manager Conda. Currently, Bioconda offers a collection of over 3000 software packages, which is continuously maintained, updated, and extended by a growing global community of more than 200 contributors. Bioconda improves analysis reproducibility by allowing users to define isolated environments with defined software versions, all of which are easily installed and managed without administrative privileges.