shinyseg: a web application for flexible cosegregation and sensitivity analysis.

MOTIVATION: Cosegregation analysis is a powerful tool for identifying pathogenic genetic variants, but its implementation remains challenging. Existing software is either limited in scope or too demanding for many end users. Moreover, current solutions lack methods for assessing the robustness of cosegregation evidence, which is important due to its reliance on uncertain estimates. RESULTS: We present shinyseg, a comprehensive web application for clinical cosegregation analysis. Our app streamlines penetrance specification based on either liability classes or epidemiological data such as risks, hazard ratios, and age of onset distribution. In addition, it incorporates sensitivity analyses to assess the robustness of cosegregation evidence, and offers support in clinical interpretation. AVAILABILITY AND IMPLEMENTATION: The shinyseg app is freely available at https://chrcarrizosa.shinyapps.io/shinyseg, with documentation and complete R source code on https://chrcarrizosa.github.io/shinyseg and https://github.com/chrcarrizosa/shinyseg.

A 39 kb structural variant causing Lynch Syndrome detected by optical genome mapping and nanopore sequencing.

Lynch Syndrome (LS) is a hereditary cancer syndrome caused by pathogenic germline variants in one of the four mismatch repair (MMR) genes MLH1, MSH2, MSH6 and PMS2. It is characterized by a significantly increased risk of multiple cancer types, particularly colorectal and endometrial cancer, with autosomal dominant inheritance. Access to precise and sensitive methods for genetic testing is important, as early detection and prevention of cancer is possible when the variant is known. We present here two unrelated Norwegian families with family histories strongly suggestive of LS, where immunohistochemical and microsatellite instability analyses indicated presence of a pathogenic variant in MSH2, but targeted exon sequencing and multiplex ligation-dependent probe amplification (MLPA) were negative. Using Bionano optical genome mapping, we detected a 39 kb insertion in the MSH2 gene. Precise mapping of the insertion breakpoints and inserted sequence was performed by low-coverage whole-genome sequencing with an Oxford Nanopore MinION. The same variant was present in both families, and later found in other families from the same region of Norway, indicative of a founder event. To our knowledge, this is the first diagnosis of LS caused by a structural variant using these technologies. We suggest that structural variant detection be performed when LS is suspected but not confirmed with first-tier standard genetic testing.

Whole-exome sequencing in moyamoya patients of Northern-European origin identifies gene variants involved in Nitric Oxide metabolism: A pilot study.

Introduction: Moyamoya disease (MMD) is a chronic cerebrovascular steno-occlusive disease of largely unknown etiology. Variants in the RNF213 gene are strongly associated with MMD in East-Asia. In MMD patients of Northern-European origin, no predominant susceptibility variants have been identified so far. Research question: Are there specific candidate genes associated with MMD of Northern-European origin, including the known RNF213 gene? Can we establish a hypothesis for MMD phenotype and associated genetic variants identified for further research? Material and methods: Adult patients of Northern-European origin, treated surgically for MMD at Oslo University Hospital between October 2018 to January 2019 were asked to participate. WES was performed, with subsequent bioinformatic analysis and variant filtering. The selected candidate genes were either previously reported in MMD or known to be involved in angiogenesis. The variant filtering was based on variant type, location, population frequency, and predicted impact on protein function. Results: Analysis of WES data revealed nine variants of interest in eight genes. Five of those encode proteins involved in nitric oxide (NO) metabolism: NOS3, NR4A3, ITGAV, GRB7 and AGXT2. In the AGXT2 gene, a de novo variant was detected, not previously described in MMD. None harboured the p.R4810K missense variant in the RNF213 gene known to be associated with MMD in East-Asian patients. Discussion and conclusion: Our findings suggest a role for NO regulation pathways in Northern-European MMD and introduce AGXT2 as a new susceptibility gene. This pilot study warrants replication in larger patient cohorts and further functional investigations.

The use of segregation analysis in interpretation of sequence variants in SMAD3: A case report.

BACKGROUND: While representing a significant improvement, the introduction of next-generation sequencing in genetic diagnosis also prompted new challenges. Despite widely recognized consensus guidelines for the interpretation of sequence variants, many variants remain unclassified or are discordantly interpreted. In heritable thoracic aortic aneurysms with dissection (HTAAD), most cases are caused by a heterozygous, private missense mutation, possibly contributing to the relatively common reports of variants with uncertain significance in this group. Segregation analysis necessitates advanced likelihood-based methods typically inaccessible to non-experts and is hampered by reduced penetrance, possible phenocopies, and non-availability of DNA from deceased relatives. METHODS: In this report, challenges in variant interpretation and the use of segregation analyses were illustrated in two families with a suspected HTAAD disorder. The R package segregatr, a novel implementation of full-likelihood Bayes factor (FLB), was performed to explore the cosegregation of the variants in these families. CONCLUSION: Using the R package segregatr, cosegregation in the reported families concluded with strong and supporting evidence for pathogenicity. Surveillance of families in a multidisciplinary team enabling systematic phenotype description for standardized segregation analysis with a robust calculation method may be imperative for reliable variant interpretation in HTAAD.

QuickPed: an online tool for drawing pedigrees and analysing relatedness.

BACKGROUND: The ubiquity of pedigrees in many scientific areas calls for versatile and user-friendly software. Previously published online pedigree tools have limited support for complex pedigrees and do not provide analysis of relatedness between pedigree members. RESULTS: We introduce QuickPed, a web application for interactive pedigree creation and analysis. It supports complex inbreeding and comes with a rich built-in library of common and interesting pedigrees. The program calculates all standard coefficients of relatedness, including inbreeding, kinship and identity coefficients, and offers specialised plots for visualising relatedness. It also implements a novel algorithm for describing pairwise relationships in words. CONCLUSION: QuickPed is a user-friendly pedigree tool aimed at researchers, case workers and teachers. It contains a number of features not found in other similar tools, and represents a significant addition to the body of pedigree software by making advanced relatedness analyses available for non-bioinformaticians.

Joint DNA-based disaster victim identification.

We address computational and statistical aspects of DNA-based identification of victims in the aftermath of disasters. Current methods and software for such identification typically consider each victim individually, leading to suboptimal power of identification and potential inconsistencies in the statistical summary of the evidence. We resolve these problems by performing joint identification of all victims, using the complete genetic data set. Individual identification probabilities, conditional on all available information, are derived from the joint solution in the form of posterior pairing probabilities. A closed formula is obtained for the a priori number of possible joint solutions to a given DVI problem. This number increases quickly with the number of victims and missing persons, posing computational challenges for brute force approaches. We address this complexity with a preparatory sequential step aiming to reduce the search space. The examples show that realistic cases are handled efficiently. User-friendly implementations of all methods are provided in the R package dvir, freely available on all platforms.

Making decisions in missing person identification cases with low statistical power.

The present work proposes a general strategy for dealing with missing person identification cases through DNA-database search. Our main example is the identification of abducted children in the last civic-dictatorship of Argentina, known as the "Missing Grandchildren of Argentina". Particularly we focus on those pedigrees where few, or only distant relatives of the missing person are available, resulting in low statistical power. For such complex cases we provide a statistical method for selecting a likelihood ratio (LR) threshold for each pedigree based on error rates. Furthermore, we provide an open-source user friendly software for computing LR thresholds and error rates. The strategy described in the paper could be applied to other large-scale cases of DNA-based identification hampered by low statistical power.

Differential Glial Activation in Early Epileptogenesis-Insights From Cell-Specific Analysis of DNA Methylation and Gene Expression in the Contralateral Hippocampus.

Background and Aims: Morphological changes in mesial temporal lobe epilepsy with hippocampal sclerosis (mTLE-HS) are well-characterized. Yet, it remains elusive whether these are a consequence of seizures or originate from a hitherto unknown underlying pathology. We recently published data on changes in gene expression and DNA methylation in the ipsilateral hippocampus (ILH) using the intracortical kainate mouse model of mTLE-HS. In order to explore the effects of epileptic activity alone and also to further disentangle what triggers morphological alterations, we investigated glial and neuronal changes in gene expression and DNA methylation in the contralateral hippocampus (CLH). Methods: The intracortical kainic acid mouse model of mTLE-HS was used to elicit status epilepticus. Hippocampi contralateral to the injection site from eight kainate-injected and eight sham mice were extracted and shock frozen at 24 h post-injection. Glial and neuronal nuclei were sorted by flow cytometry. Alterations in gene expression and DNA methylation were assessed using reduced representation bisulfite sequencing and RNA sequencing. The R package edgeR was used for statistical analysis. Results: The CLH featured substantial, mostly cell-specific changes in both gene expression and DNA methylation in glia and neurons. While changes in gene expression overlapped to a great degree between CLH and ILH, alterations in DNA methylation did not. In the CLH, we found a significantly lower number of glial genes up- and downregulated compared to previous results from the ILH. Furthermore, several genes and pathways potentially involved in anti-epileptogenic effects were upregulated in the CLH. By comparing gene expression data from the CLH to previous results from the ILH (featuring hippocampal sclerosis), we derive potential upstream targets for epileptogenesis, including glial Cox2 and Cxcl10. Conclusion: Despite the absence of morphological changes, the CLH displays substantial changes in gene expression and DNA methylation. We find that gene expression changes related to potential anti-epileptogenic effects seem to dominate compared to the pro-epileptogenic effects in the CLH and speculate whether this imbalance contributes to prevent morphological alterations like neuronal death and reactive gliosis.

Prioritising family members for genotyping in missing person cases: A general approach combining the statistical power of exclusion and inclusion.

Missing person identification typically involves genetic matching of a person of interest against relatives of the missing person. In cases with few available relatives, exhumations or other substantial efforts may be necessary in order to secure adequate statistical power. We propose a simulation approach for solving prioritisation problems arising in such cases. Conditioning on the already typed individuals we estimate the power of each alternative, both to detect the true person, and to exclude false candidates. Graphical summaries of the simulations are given in complementary power plots, facilitating interpretation and decision making. Through a series of examples originating from the well-known Missing grandchildren of Argentina we demonstrate that our method may untangle complex prioritisation problems and other power-related questions. In particular we offer novel insights in recent cases where only children of the potential match are available for testing. We also show that X-chromosomal markers may give high statistical power in missing person identification, but that this requires careful selection of relatives for genotyping. All simulations, power calculations and plots are done with the R package forrel.

Phenotypic spectrum and transcriptomic profile associated with germline variants in TRAF7.

PURPOSE: Somatic variants in tumor necrosis factor receptor-associated factor 7 (TRAF7) cause meningioma, while germline variants have recently been identified in seven patients with developmental delay and cardiac, facial, and digital anomalies. We aimed to define the clinical and mutational spectrum associated with TRAF7 germline variants in a large series of patients, and to determine the molecular effects of the variants through transcriptomic analysis of patient fibroblasts. METHODS: We performed exome, targeted capture, and Sanger sequencing of patients with undiagnosed developmental disorders, in multiple independent diagnostic or research centers. Phenotypic and mutational comparisons were facilitated through data exchange platforms. Whole-transcriptome sequencing was performed on RNA from patient- and control-derived fibroblasts. RESULTS: We identified heterozygous missense variants in TRAF7 as the cause of a developmental delay-malformation syndrome in 45 patients. Major features include a recognizable facial gestalt (characterized in particular by blepharophimosis), short neck, pectus carinatum, digital deviations, and patent ductus arteriosus. Almost all variants occur in the WD40 repeats and most are recurrent. Several differentially expressed genes were identified in patient fibroblasts. CONCLUSION: We provide the first large-scale analysis of the clinical and mutational spectrum associated with the TRAF7 developmental syndrome, and we shed light on its molecular etiology through transcriptome studies.

Mitochondrial genome-wide association study of migraine - the HUNT Study.

BACKGROUND: Variation in mitochondrial DNA (mtDNA) has been indicated in migraine pathogenesis, but genetic studies to date have focused on candidate variants, with sparse findings. We aimed to perform the first mitochondrial genome-wide association study of migraine, examining both single variants and mitochondrial haplogroups. METHODS: In total, 71,860 participants from the population-based Nord-Trøndelag Health Study were genotyped. We excluded samples not passing quality control for nuclear genotypes, in addition to samples with low call rate and closely maternally related. We analysed 775 mitochondrial DNA variants in 4021 migraine cases and 14,288 headache-free controls, using logistic regression. In addition, we analysed 3831 cases and 13,584 controls who could be reliably assigned to a mitochondrial haplogroup. Lastly, we attempted to replicate previously reported mitochondrial DNA candidate variants. RESULTS: Neither of the mitochondrial variants or haplogroups were associated with migraine. In addition, none of the previously reported mtDNA candidate variants replicated in our data. CONCLUSIONS: Our findings do not support a major role of mitochondrial genetic variation in migraine pathophysiology, but a larger sample is needed to detect rare variants and future studies should also examine heteroplasmic variation, epigenetic changes and copy-number variation.

Neuronal and glial DNA methylation and gene expression changes in early epileptogenesis.

BACKGROUND AND AIMS: Mesial Temporal Lobe Epilepsy is characterized by progressive changes of both neurons and glia, also referred to as epileptogenesis. No curative treatment options, apart from surgery, are available. DNA methylation (DNAm) is a potential upstream mechanism in epileptogenesis and may serve as a novel therapeutic target. To our knowledge, this is the first study to investigate epilepsy-related DNAm, gene expression (GE) and their relationship, in neurons and glia. METHODS: We used the intracortical kainic acid injection model to elicit status epilepticus. At 24 hours post injection, hippocampi from eight kainic acid- (KA) and eight saline-injected (SH) mice were extracted and shock frozen. Separation into neurons and glial nuclei was performed by flow cytometry. Changes in DNAm and gene expression were measured with reduced representation bisulfite sequencing (RRBS) and mRNA-sequencing (mRNAseq). Statistical analyses were performed in R with the edgeR package. RESULTS: We observed fulminant DNAm- and GE changes in both neurons and glia at 24 hours after initiation of status epilepticus. The vast majority of these changes were specific for either neurons or glia. At several epilepsy-related genes, like HDAC11, SPP1, GAL, DRD1 and SV2C, significant differential methylation and differential gene expression coincided. CONCLUSION: We found neuron- and glia-specific changes in DNAm and gene expression in early epileptogenesis. We detected single genetic loci in several epilepsy-related genes, where DNAm and GE changes coincide, worth further investigation. Further, our results may serve as an information source for neuronal and glial alterations in both DNAm and GE in early epileptogenesis.

Coexistence of Congenital Adrenal Hyperplasia and Autoimmune Addison's Disease.

Background: Underlying causes of adrenal insufficiency include congenital adrenal hyperplasia (CAH) and autoimmune adrenocortical destruction leading to autoimmune Addison's disease (AAD). Here, we report a patient with a homozygous stop-gain mutation in 3ß-hydroxysteroid dehydrogenase type 2 (3ßHSD2), in addition to impaired steroidogenesis due to AAD. Case Report: Whole exome sequencing revealed an extremely rare homozygous nonsense mutation in exon 2 of the HSD3B2 gene, leading to a premature stop codon (NM_000198.3: c.15C>A, p.Cys5Ter) in a patient with AAD and premature ovarian insufficiency. Scrutiny of old medical records revealed that the patient was initially diagnosed with CAH with hyperandrogenism and severe salt-wasting shortly after birth. However, the current steroid profile show complete adrenal insufficiency including low production of pregnenolone, dehydroepiandrosterone (DHEA) and DHEA sulfate (DHEA-S), without signs of overtreatment with steroids. Conclusion: To the best of our knowledge, this is the first description of autoimmune adrenalitis in a patient with 3ßHSD2 deficiency and suggests a possible association between AAD and inborn errors of the steroidogenesis.

Identification and characterization of rare toll-like receptor 3 variants in patients with autoimmune Addison's disease.

Autoimmune Addison's disease (AAD) is a classic organ-specific autoimmune disease characterized by an immune-mediated attack on the adrenal cortex. As most autoimmune diseases, AAD is believed to be caused by a combination of genetic and environmental factors, and probably interactions between the two. Persistent viral infections have been suggested to play a triggering role, by invoking inflammation and autoimmune destruction. The inability of clearing infections can be due to aberrations in innate immunity, including mutations in genes involved in the recognition of conserved microbial patterns. In a whole exome sequencing study of anonymized AAD patients, we discovered several rare variants predicted to be damaging in the gene encoding Toll-like receptor 3 (TLR3). TLR3 recognizes double stranded RNAs, and is therefore a major factor in antiviral defense. We here report the occurrence and functional characterization of five rare missense variants in TLR3 of patients with AAD. Most of these variants occurred together with a common TLR3 variant that has been associated with a wide range of immunopathologies. The biological implications of these variants on TLR3 function were evaluated in a cell-based assay, revealing a partial loss-of-function effect of three of the rare variants. In addition, rare mutations in other members of the TLR3-interferon (IFN) signaling pathway were detected in the AAD patients. Together, these findings indicate a potential role for TLR3 and downstream signaling proteins in the pathogenesis in a subset of AAD patients.

Further delineation of an entity caused by CREBBP and EP300 mutations but not resembling Rubinstein-Taybi syndrome.

In 2016, we described that missense variants in parts of exons 30 and 31 of CREBBP can cause a phenotype that differs from Rubinstein-Taybi syndrome (RSTS). Here we report on another 11 patients with variants in this region of CREBBP (between bp 5,128 and 5,614) and two with variants in the homologous region of EP300. None of the patients show characteristics typical for RSTS. The variants were detected by exome sequencing using a panel for intellectual disability in all but one individual, in whom Sanger sequencing was performed upon clinical recognition of the entity. The main characteristics of the patients are developmental delay (90%), autistic behavior (65%), short stature (42%), and microcephaly (43%). Medical problems include feeding problems (75%), vision (50%), and hearing (54%) impairments, recurrent upper airway infections (42%), and epilepsy (21%). Major malformations are less common except for cryptorchidism (46% of males), and cerebral anomalies (70%). Individuals with variants between bp 5,595 and 5,614 of CREBBP show a specific phenotype (ptosis, telecanthi, short and upslanted palpebral fissures, depressed nasal ridge, short nose, anteverted nares, short columella, and long philtrum). 3D face shape demonstrated resemblance to individuals with a duplication of 16p13.3 (the region that includes CREBBP), possibly indicating a gain of function. The other affected individuals show a less specific phenotype. We conclude that there is now more firm evidence that variants in these specific regions of CREBBP and EP300 result in a phenotype that differs from RSTS, and that this phenotype may be heterogeneous.

PUF60 variants cause a syndrome of ID, short stature, microcephaly, coloboma, craniofacial, cardiac, renal and spinal features.

PUF60 encodes a nucleic acid-binding protein, a component of multimeric complexes regulating RNA splicing and transcription. In 2013, patients with microdeletions of chromosome 8q24.3 including PUF60 were found to have developmental delay, microcephaly, craniofacial, renal and cardiac defects. Very similar phenotypes have been described in six patients with variants in PUF60, suggesting that it underlies the syndrome. We report 12 additional patients with PUF60 variants who were ascertained using exome sequencing: six through the Deciphering Developmental Disorders Study and six through similar projects. Detailed phenotypic analysis of all patients was undertaken. All 12 patients had de novo heterozygous PUF60 variants on exome analysis, each confirmed by Sanger sequencing: four frameshift variants resulting in premature stop codons, three missense variants that clustered within the RNA recognition motif of PUF60 and five essential splice-site (ESS) variant. Analysis of cDNA from a fibroblast cell line derived from one of the patients with an ESS variants revealed aberrant splicing. The consistent feature was developmental delay and most patients had short stature. The phenotypic variability was striking; however, we observed similarities including spinal segmentation anomalies, congenital heart disease, ocular colobomata, hand anomalies and (in two patients) unilateral renal agenesis/horseshoe kidney. Characteristic facial features included micrognathia, a thin upper lip and long philtrum, narrow almond-shaped palpebral fissures, synophrys, flared eyebrows and facial hypertrichosis. Heterozygote loss-of-function variants in PUF60 cause a phenotype comprising growth/developmental delay and craniofacial, cardiac, renal, ocular and spinal anomalies, adding to disorders of human development resulting from aberrant RNA processing/spliceosomal function.

Novel UCHL1 mutations reveal new insights into ubiquitin processing.

Recessive loss of function of the neuronal ubiquitin hydrolase UCHL1 has been implicated in early-onset progressive neurodegeneration (MIM no. 615491), so far only in one family. In this study a second family is characterized, and the functional consequences of the identified mutations in UCHL1 are explored. Three siblings developed childhood-onset optic atrophy, followed by spasticity and ataxia. Whole exome sequencing identified compound heterozygous variants in UCHL1, c.533G > A (p.Arg178Gln) and c.647C > A (p.Ala216Asp), cosegregating with the phenotype. Enzymatic activity of purified recombinant proteins analysed by ubiquitin hydrolase assays showed a 4-fold increased hydrolytic activity of the recombinant UCHL1 mutant Arg178Gln compared to wild type, whereas the Ala216Asp protein was insoluble. Structural 3D analysis of UCHL1 by computer modelling suggests that Arg178 is a rate-controlling residue in catalysis which is partly abolished in the Arg178Gln mutant and, consequently, the Arg178Gln mutant increases the enzymatic turnover. UCHL1 protein levels in fibroblasts measured by targeted mass spectrometry showed a total amount of UCHL1 in control fibroblasts about 4-fold higher than in the patients. Hence, studies of the identified missense variants reveal surprisingly different functional consequences as the insoluble Ala216Asp variant leads to loss of function, whereas the Arg178Gln leads to increased enzyme activity. The reported patients have remarkably preserved cognition, and we propose that the increased enzyme activity of the Arg178Gln variant offers a protective effect on cognitive function. This study establishes the importance of UCHL1 in neurodegeneration, provides new mechanistic insight about ubiquitin processing, and underlines the complexity of the different roles of UCHL1.