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.

Longitudinal changes of serum cytokines in patients with chronic low back pain and Modic changes

OBJECTIVES

To explore serum cytokine levels over time in patients with chronic low back pain (cLBP) and Modic changes (MCs), difference in change between treatment groups in the Antibiotics in Modic Changes (AIM) study and associations between change in cytokines and low back pain.

METHODS

Serum concentrations of 39 cytokines were measured at baseline and 1 year from 73 participants in the AIM study; 30 randomized to placebo, 43 to Amoxicillin. Low back pain intensity was measured by numeric rating scale. Change in cytokine levels over time were assessed by paired t-tests. Difference in change in cytokine levels between treatment groups and associations between changes in LBP and cytokine levels were assessed by linear regression models. Networks of cytokine changes in each treatment groups were explored by Pearson's correlations.

RESULTS

Five cytokines changed from baseline to 1 year, (mean change, log transformed values with CI) C-X-C motif chemokine ligand (CXCL) 10 (IP-10) (0.11 (0.01-0.20)), CXCL13 (0.61 (0.00-0.12)), C-C motif chemokine ligand (CCL)26 (0.05 (0.01-0.1)), granulocyte macrophage-colony stimulating factor (GM-CSF) (-0.12 (-0.23 to 0.00)) and CXCL11 (0.12 (0.03-0.22)). Treatment group only influenced change in CCL21 (β 0.07 (0.01-0.12)), and IL-6 (β -0.17 (-0.30 to -0.03)). Change in CXCL13 (β 2.43 (0.49-4.38)), CCL27 (β 3.07 (0.46-5.69)), IL-8 (β 1.83 (0.08-3.58)) and CCL19 (β 3.10 (0.86-5.43)) were associated with change in LBP. The correlation networks of cytokine changes demonstrate small differences between treatment groups.

CONCLUSIONS

Cytokine levels are relatively stable over time in our sample, with little difference between treatment groups. Some cytokines may be associated with LBP intensity. The differences between the correlation networks suggest that long-term Amoxicillin-treatment may have longstanding effects to be further explored.

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.

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.

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.

Homozygous KIDINS220 loss-of-function variants in fetuses with cerebral ventriculomegaly and limb contractures

Heterozygous mutations in KIDINS220 were recently suggested a cause of spastic paraplegia, intellectual disability, nystagmus and obesity. All patients carried terminal nonsense de novo mutations that seemed to escape nonsense-mediated mRNA decay. The mechanism for pathogenicity is yet unexplained, as it seems that heterozygous loss-of-function variants of KIDINS220 are generally well tolerated. We present a consanguineous couple who experienced four pregnancy terminations due to repeated findings in the fetuses comprising enlarged cerebral ventricles and limb contractures. Exome sequencing in two of the aborted fetuses revealed a shared homozygous frameshift variant in exon 24 in KIDINS220. Sanger sequencing of the variant in available family members showed complete segregation with the affection status, resulting in a LOD score of 2.5 under an autozygous inheritance model. mRNA studies revealed destruction of the original splice site, resulting in an out-of-frame transcript and introduction of a premature termination codon in exon 25. Premature termination codons in this position are likely to cause activation of nonsense-mediated mRNA decay and result in complete absence of KIDINS220 protein in individuals homozygous for the variant. The phenotype of the presented fetuses overlaps with findings in functional studies of knockout Kidins220 mice embryos that are non-viable with enlarged cerebral ventricles. The human fetuses also exhibit several similarities to the milder phenotype described in patients with heterozygous KIDINS220 mutations. We hence propose that the identified homozygous loss-of-function variant in KIDINS220 causes the phenotype in the presented fetuses, and that this represents a hitherto undescribed severe autosomal recessive neurodevelopmental disorder.

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.

A potential founder variant in CARMIL2/RLTPR in three Norwegian families with warts, molluscum contagiosum, and T-cell dysfunction

Background

Four patients from three Norwegian families presented with a common skin phenotype of warts, molluscum contagiosum, and dermatitis since early childhood, and various other immunological features. Warts are a common manifestation of human papilloma virus (HPV), but when they are overwhelming, disseminated and/or persistent, and presenting together with other immunological features, a primary immunodeficiency disease (PIDD) may be suspected.

Methods and results

The four patients were exome sequenced as part of a larger study for detecting genetic causes of primary immunodeficiencies. No disease‐causing variants were identified in known primary immunodeficiency genes or in other disease‐related OMIM genes. However, the same homozygous missense variant in CARMIL2 (also known as RLTPR) was identified in all four patients. In each family, the variant was located within a narrow region of homozygosity, representing a potential region of autozygosity. CARMIL2 is a protein of undetermined function. A role in T‐cell activation has been suggested and the mouse protein homolog (Rltpr) is essential for costimulation of T‐cell activation via CD28, and for the development of regulatory T cells. Immunophenotyping demonstrated reduced regulatory, CD4+ memory, and CD4+ follicular T cells in all four patients. In addition, they all seem to have a deficiency in IFNγ ‐synthesis in CD4+ T cells and NK cells.

Conclusions

We report a novel primary immunodeficiency, and a differential molecular diagnosis to CXCR4‐,DOCK8‐,GATA2‐,MAGT1‐,MCM4‐,STK4‐,RHOH‐,TMC6‐, and TMC8‐related diseases. The specific variant may represent a Norwegian founder variant segregating on a population‐specific haplotype.

Exome Sequencing Fails to Identify the Genetic Cause of Aicardi Syndrome

Aicardi syndrome (AS) is a well-characterized neurodevelopmental disorder with an unknown etiology. In this study, we performed whole-exome sequencing in 11 female patients with the diagnosis of AS, in order to identify the disease-causing gene. In particular, we focused on detecting variants in the X chromosome, including the analysis of variants with a low number of sequencing reads, in case of somatic mosaicism. For 2 of the patients, we also sequenced the exome of the parents to search for de novo mutations. We did not identify any genetic variants likely to be damaging. Only one single missense variant was identified by the de novo analyses of the 2 trios, and this was considered benign. The failure to identify a disease gene in this study may be due to technical limitations of our study design, including the possibility that the genetic aberration leading to AS is situated in a non-exonic region or that the mutation is somatic and not detectable by our approach. Alternatively, it is possible that AS is genetically heterogeneous and that 11 patients are not sufficient to reveal the causative genes. Future studies of AS should consider designs where also non-exonic regions are explored and apply a sequencing depth so that also low-grade somatic mosaicism can be detected.

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.

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

Summary: 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.

A founder mutation p.H701P identified as a major cause of SPG7 in Norway

BACKGROUND AND PURPOSE

SPG7 is one of the most common forms of autosomal recessive hereditary spastic paraplegia. The phenotype has been shown to be heterogeneous, varying from a complex spastic ataxia to pure spastic paraplegia or pure ataxia. The aim of this study was to clinically and genetically characterize patients with SPG7 in Norway.

METHODS

Six Norwegian families with a clinical diagnosis of hereditary spastic paraplegia were diagnosed with SPG7 through Sanger sequencing and whole-exome sequencing. Haplotypes were established to identify a possible founder mutation. All patients were thoroughly examined and the clinical and molecular findings are described.

RESULTS

The core phenotype was spastic paraparesis with ataxia, bladder disturbances and progressive external ophthalmoplegia. The variant p.H701P was identified in homozygous state in one family and in compound heterozygous state in three families. Haplotype analysis of seven surrounding single nucleotide polymorphisms supports that this variant resides on a founder haplotype. Four of the families were compound heterozygous for the previously well-described p.A510V variant.

CONCLUSION

SPG7 is a common subgroup of hereditary spinocerebellar disorders in Norway. The broad phenotype in the Norwegian SPG7 population illustrates the challenges with the traditional dichotomous classification of hereditary spinocerebellar disorders into hereditary spastic paraplegia or hereditary ataxia. A Norwegian founder mutation p.H701P was identified in four out of six families, making it a major cause of SPG7 in Norway.

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.

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.