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.

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.

The Atlantic salmon genome provides insights into rediploidization.

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

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.

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.

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.

FILTUS: a desktop GUI for fast and efficient detection of disease-causing variants, including a novel autozygosity detector.

UNLABELLED: FILTUS is a stand-alone tool for working with annotated variant files, e.g. when searching for variants causing Mendelian disease. Very flexible in terms of input file formats, FILTUS offers efficient filtering and a range of downstream utilities, including statistical analysis of gene sharing patterns, detection of de novo mutations in trios, quality control plots and autozygosity mapping. The autozygosity mapping is based on a hidden Markov model and enables accurate detection of autozygous regions directly from exome-scale variant files. AVAILABILITY AND IMPLEMENTATION: FILTUS is written in Python and runs on Windows, Mac and Linux. Binaries and source code are freely available at http://folk.uio.no/magnusv/filtus.html and on GitHub: https://github.com/magnusdv/filtus Automatic installation is available via PyPI (e.g. pip install filtus). CONTACT: magnusdv@medisin.uio.no SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Kaufman oculocerebrofacial syndrome in sisters with novel compound heterozygous mutation in UBE3B.

A pair of sisters was ascertained for multiple congenital defects, including marked craniofacial dysmorphisms with blepharophimosis, and severe psychomotor delay. Two novel compound heterozygous mutations in UBE3B were identified in both the sisters by exome sequencing. These mutations include c.1A>G, which predicts p.Met1?, and a c.1773delC variant, predicted to cause a frameshift at p.Phe591fs. UBE3B encodes a widely expressed protein ubiquitin ligase E3B, which, when mutated in both alleles, causes Kaufman oculocerebrofacial syndrome. We report on the thorough clinical examination of the patients and review the state of art knowledge of this disorder.

DNA methylation and gene expression changes in monozygotic twins discordant for psoriasis: identification of epigenetically dysregulated genes.

Monozygotic (MZ) twins do not show complete concordance for many complex diseases; for example, discordance rates for autoimmune diseases are 20%-80%. MZ discordance indicates a role for epigenetic or environmental factors in disease. We used MZ twins discordant for psoriasis to search for genome-wide differences in DNA methylation and gene expression in CD4(+) and CD8(+) cells using Illumina's HumanMethylation27 and HT-12 expression assays, respectively. Analysis of these data revealed no differentially methylated or expressed genes between co-twins when analyzed separately, although we observed a substantial amount of small differences. However, combined analysis of DNA methylation and gene expression identified genes where differences in DNA methylation between unaffected and affected twins were correlated with differences in gene expression. Several of the top-ranked genes according to significance of the correlation in CD4(+) cells are known to be associated with psoriasis. Further, gene ontology (GO) analysis revealed enrichment of biological processes associated with the immune response and clustering of genes in a biological pathway comprising cytokines and chemokines. These data suggest that DNA methylation is involved in an epigenetic dysregulation of biological pathways involved in the pathogenesis of psoriasis. This is the first study based on data from MZ twins discordant for psoriasis to detect epigenetic alterations that potentially contribute to development of the disease.

Evidence for adaptive evolution of low-temperature stress response genes in a Pooideae grass ancestor.

Adaptation to temperate environments is common in the grass subfamily Pooideae, suggesting an ancestral origin of cold climate adaptation. Here, we investigated substitution rates of genes involved in low-temperature-induced (LTI) stress responses to test the hypothesis that adaptive molecular evolution of LTI pathway genes was important for Pooideae evolution. Substitution rates and signatures of positive selection were analyzed using 4330 gene trees including three warm climate-adapted species (maize (Zea mays), sorghum (Sorghum bicolor), and rice (Oryza sativa)) and five temperate Pooideae species (Brachypodium distachyon, wheat (Triticum aestivum), barley (Hordeum vulgare), Lolium perenne and Festuca pratensis). Nonsynonymous substitution rate differences between Pooideae and warm habitat-adapted species were elevated in LTI trees compared with all trees. Furthermore, signatures of positive selection were significantly stronger in LTI trees after the rice and Pooideae split but before the Brachypodium divergence (P < 0.05). Genome-wide heterogeneity in substitution rates was also observed, reflecting divergent genome evolution processes within these grasses. Our results provide evidence for a link between adaptation to cold habitats and adaptive evolution of LTI stress responses in early Pooideae evolution and shed light on a poorly understood chapter in the evolutionary history of some of the world's most important temperate crops.

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.

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.

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.

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.

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.

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.