GWAS meta-analysis of psoriasis identifies new susceptibility alleles impacting disease mechanisms and therapeutic targets

Psoriasis is a common, debilitating immune-mediated skin disease. Genetic studies have identified biological mechanisms of psoriasis risk, including those targeted by effective therapies. However, the genetic liability to psoriasis is not fully explained by variation at robustly identified risk loci. To move towards a saturation map of psoriasis susceptibility we meta-analysed 18 GWAS comprising 36,466 cases and 458,078 controls and identified 109 distinct psoriasis susceptibility loci, including 45 that have not been previously reported. These include susceptibility variants at loci in which the therapeutic targets IL17RA and AHR are encoded, and deleterious coding variants supporting potential new drug targets (including in STAP2, CPVL and POU2F3). We conducted a transcriptome-wide association study to identify regulatory effects of psoriasis susceptibility variants and cross-referenced these against single cell expression profiles in psoriasis-affected skin, highlighting roles for the transcriptional regulation of haematopoietic cell development and epigenetic modulation of interferon signalling in psoriasis pathobiology.

Development of antidrug antibodies against adalimumab maps to variation within the HLA-DR peptide-binding groove

Targeted biologic therapies can elicit an undesirable host immune response characterized by the development of antidrug antibodies (ADA), an important cause of treatment failure. The most widely used biologic across immune-mediated diseases is adalimumab, a tumor necrosis factor inhibitor. This study aimed to identify genetic variants that contribute to the development of ADA against adalimumab, thereby influencing treatment failure. In patients with psoriasis on their first course of adalimumab, in whom serum ADA had been evaluated 6-36 months after starting treatment, we observed a genome-wide association with ADA against adalimumab within the major histocompatibility complex (MHC). The association signal mapped to the presence of tryptophan at position 9 and lysine at position 71 of the HLA-DR peptide-binding groove, with both residues conferring protection against ADA. Underscoring their clinical relevance, these residues were also protective against treatment failure. Our findings highlight antigenic peptide presentation via MHC class II as a critical mechanism in the development of ADA against biologic therapies and downstream treatment response.

Pathogenic variants in SLF2 and SMC5 cause segmented chromosomes and mosaic variegated hyperploidy

Embryonic development is dictated by tight regulation of DNA replication, cell division and differentiation. Mutations in DNA repair and replication genes disrupt this equilibrium, giving rise to neurodevelopmental disease characterized by microcephaly, short stature and chromosomal breakage. Here, we identify biallelic variants in two components of the RAD18-SLF1/2-SMC5/6 genome stability pathway, SLF2 and SMC5, in 11 patients with microcephaly, short stature, cardiac abnormalities and anemia. Patient-derived cells exhibit a unique chromosomal instability phenotype consisting of segmented and dicentric chromosomes with mosaic variegated hyperploidy. To signify the importance of these segmented chromosomes, we have named this disorder Atelís (meaning - incomplete) Syndrome. Analysis of Atelís Syndrome cells reveals elevated levels of replication stress, partly due to a reduced ability to replicate through G-quadruplex DNA structures, and also loss of sister chromatid cohesion. Together, these data strengthen the functional link between SLF2 and the SMC5/6 complex, highlighting a distinct role for this pathway in maintaining genome stability.

MED27, SLC6A7, and MPPE1 Variants in a Complex Neurodevelopmental Disorder with Severe Dystonia

BACKGROUND

Despite advances in next generation sequencing technologies, the identification of variants of uncertain significance (VUS) can often hinder definitive diagnosis in patients with complex neurodevelopmental disorders.

OBJECTIVE

The objective of this study was to identify and characterize the underlying cause of disease in a family with two children with severe developmental delay associated with generalized dystonia and episodic status dystonicus, chorea, epilepsy, and cataracts.

METHODS

Candidate genes identified by autozygosity mapping and whole-exome sequencing were characterized using cellular and vertebrate model systems.

RESULTS

Homozygous variants were found in three candidate genes: MED27, SLC6A7, and MPPE1. Although the patients had features of MED27-related disorder, the SLC6A7 and MPPE1 variants were functionally investigated. SLC6A7 variant in vitro overexpression caused decreased proline transport as a result of reduced cell-surface expression, and zebrafish knockdown of slc6a7 exhibited developmental delay and fragile motor neuron morphology that could not be rescued by L-proline transporter-G396S RNA. Lastly, patient fibroblasts displayed reduced cell-surface expression of glycophosphatidylinositol-anchored proteins linked to MPPE1 dysfunction.

CONCLUSIONS

We report a family harboring a homozygous MED27 variant with additional loss-of-function SLC6A7 and MPPE1 gene variants, which potentially contribute to a blended phenotype caused by multilocus pathogenic variants. © 2022 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Differences in Clinical Features and Comorbid Burden between HLA-C∗06:02 Carrier Groups in >9,000 People with Psoriasis

The identification of robust endotypes-disease subgroups of clinical relevance-is fundamental to stratified medicine. We hypothesized that HLA-C∗06:02 status, the major genetic determinant of psoriasis, defines a psoriasis endotype of clinical relevance. Using two United Kingdom-based cross-sectional datasets-an observational severe-psoriasis study (Biomarkers of Systemic Treatment Outcomes in Psoriasis; n = 3,767) and a large population-based bioresource (UK Biobank, including n = 5,519 individuals with psoriasis)-we compared demographic, environmental, and clinical variables of interest in HLA-C∗06:02-positive (one or two copies of the HLA-C∗06:02 allele) with those in HLA-C∗06:02‒negative (no copies) individuals of European ancestry. We used multivariable regression analyses to account for mediation effects established a priori. We confirm previous observations that HLA-C∗06:02-positive status is associated with earlier age of psoriasis onset and extend findings to reveal an association with disease expressivity in females (Biomarkers of Systemic Treatment Outcomes in Psoriasis: P = 2.7 × 10-14, UK Biobank: P = 1.0 × 10-8). We also show HLA-C∗06:02-negative status to be associated with characteristic clinical features (large plaque disease, OR for HLA-C∗06:02 = 0.73, P = 7.4 × 10-4; nail involvement, OR = 0.70, P = 2.4 × 10-6); higher central adiposity (Biomarkers of Systemic Treatment Outcomes in Psoriasis: waist circumference difference of 2.0 cm, P = 8.4 × 10-4; UK Biobank: waist circumference difference of 1.4 cm, P = 1.5 × 10-4), especially in women; and a higher prevalence of other cardiometabolic comorbidities. These findings extend the clinical phenotype delineated by HLA-C∗06:02 and highlight its potential as an important biomarker to consider in future multimarker stratified medicine approaches.

Comparative genetic analysis of psoriatic arthritis and psoriasis for the discovery of genetic risk factors and risk prediction modelling

Objectives

Psoriatic arthritis (PsA) has a strong genetic component, and the identification of genetic risk factors could help identify the ~30% of psoriasis patients at high risk of developing PsA. Our objectives were to identify genetic risk factors and pathways that differentiate PsA from cutaneous‐only psoriasis (PsC) and to evaluate the performance of PsA risk prediction models.

Methods

Genome‐wide meta‐analyses were conducted separately for 5,065 patients with PsA and 21,286 healthy controls and separately for 4,340 patients with PsA and 6,431 patients with PsC. The heritability of PsA was calculated as a single‐nucleotide polymorphism (SNP)–based heritability estimate (h²_SNP) and biologic pathways that differentiate PsA from PsC were identified using Priority Index software. The generalizability of previously published PsA risk prediction pipelines was explored, and a risk prediction model was developed with external validation.

Results

We identified a novel genome‐wide significant susceptibility locus for the development of PsA on chromosome 22q11 (rs5754467; P = 1.61 × 10⁻⁹), and key pathways that differentiate PsA from PsC, including NF‐κB signaling (adjusted P = 1.4 × 10⁻⁴⁵) and Wnt signaling (adjusted P = 9.5 × 10⁻⁵⁸). The heritability of PsA in this cohort was found to be moderate (h²_SNP = 0.63), which was similar to the heritability of PsC (h²_SNP = 0.61). We observed modest performance of published classification pipelines (maximum area under the curve 0.61), with similar performance of a risk model derived using the current data.

Conclusion

Key biologic pathways associated with the development of PsA were identified, but the investigation of risk classification revealed modest utility in the available data sets, possibly because many of the PsC patients included in the present study were receiving treatments that are also effective in PsA. Future predictive models of PsA should be tested in PsC patients recruited from primary care.

Genome-wide association meta-analysis identifies 29 new acne susceptibility loci

Acne vulgaris is a highly heritable skin disorder that primarily impacts facial skin. Severely inflamed lesions may leave permanent scars that have been associated with long-term psychosocial consequences. Here, we perform a GWAS meta-analysis comprising 20,165 individuals with acne from nine independent European ancestry cohorts. We identify 29 novel genome-wide significant loci and replicate 14 of the 17 previously identified risk loci, bringing the total number of reported acne risk loci to 46. Using fine-mapping and eQTL colocalisation approaches, we identify putative causal genes at several acne susceptibility loci that have previously been implicated in Mendelian hair and skin disorders, including pustular psoriasis. We identify shared genetic aetiology between acne, hormone levels, hormone-sensitive cancers and psychiatric traits. Finally, we show that a polygenic risk score calculated from our results explains up to 5.6% of the variance in acne liability in an independent cohort.

Exome Sequencing Highlights a Potential Role for Concealed Cardiomyopathies in Youthful Sudden Cardiac Death

BACKGROUND

Sudden cardiac arrest (SCA) and sudden unexplained death (SUD) are feared sequelae of many genetic heart diseases. In rare circumstances, pathogenic variants in cardiomyopathy-susceptibility genes may result in electrical instability leading to SCA/SUD before any structural manifestations of underlying cardiomyopathy are evident.

METHODS

Collectively, 38 unexplained SCA survivors (21 males; mean age at SCA 26.4±13.1 years), 68 autopsy-inconclusive SUD cases (46 males; mean age at death 20.4±9.0 years) without disease-causative variants in the channelopathy genes, and 973 ostensibly healthy controls were included. Following exome sequencing, ultrarare (minor allele frequency ≤0.00005 in any ethnic group within Genome Aggregation Database [gnomAD, N=141 456 individuals]) nonsynonymous variants identified in 24 Clinical Genome Resource adjudicated definitive/strong evidence cardiomyopathy-susceptibility genes were analyzed. Eligible variants were adjudicated as pathogenic, likely pathogenic, or variant of uncertain significance in accordance with current American College of Medical Genetics and Genomics guidelines.

RESULTS

Overall, 7 out of 38 (18.4%) SCA survivors and 14 out of 68 (20.5%) autopsy-inconclusive, channelopathic-negative SUD cases had at least one pathogenic/likely pathogenic or a variant of uncertain significance nonsynonymous variant within a strong evidence, cardiomyopathy-susceptibility gene. Following American College of Medical Genetics and Genomics criterion variant adjudication, a pathogenic or likely pathogenic variant was identified in 3 out of 38 (7.9%; P=0.05) SCA survivors and 8 out of 68 (11.8%; P=0.0002) autopsy-inconclusive SUD cases compared to 20 out of 973 (2.1%) European controls. Interestingly, the yield of pathogenic/likely pathogenic variants was significantly greater in autopsy-inconclusive SUD cases with documented interstitial fibrosis (4/11, 36%) compared with only 4 out of 57 (7%, P<0.02) SUD cases without ventricular fibrosis.

CONCLUSIONS

Our data further supports the inclusion of strong evidence cardiomyopathy-susceptibility genes on the genetic testing panels used to evaluate unexplained SCA survivors and autopsy-inconclusive/negative SUD decedents. However, to avoid diagnostic miscues, the careful interpretation of genetic test results in patients without overt phenotypes is vital.

Confirmation of Cause of Death Via Comprehensive Autopsy and Whole Exome Molecular Sequencing in People With Epilepsy and Sudden Unexpected Death

Background

Sudden cardiac arrest is the leading mode of death in the United States. Epilepsy affects 1% of Americans; yet epidemiological data show a prevalence of 4% in cases of sudden cardiac arrest. Sudden unexpected death in epilepsy (SUDEP) may share features with sudden cardiac arrest. The objective of this study was to report autopsy and genomic findings in a large cohort of SUDEP cases.

Methods and Results

Mayo Clinic Sudden Death Registry containing cases (ages 0–90 years) of sudden unexpected and unexplained deaths 1960 to present was queried. Exome sequencing performed on decedent cases. From 13 687 cases of sudden death, 656 (4.8%) had a history of seizures, including 368 confirmed by electroencephalography, 96 classified as SUDEP, 58 as non‐SUDEP, and 214 as unknown (insufficient records). Mean age of death in SUDEP was 37 (±19.7) years; 56 (58.3%) were male; 65% of deaths occurred at night; 54% were found in bed; and 80.6% were prone. Autopsies were obtained in 83 cases; bystander coronary artery disease was frequently reported as cause of death; nonspecific fibrosis was seen in 32.6% of cases, in structurally normal hearts. There were 4 cases of Dravet syndrome with pathogenic variants in SCN1A gene. Using whole exome sequencing in 11 cases, 18 ultrarare nonsynonymous variants were identified in 6 cases including CACNB2, RYR2, CLNB, CACNA1H, and CLCN2.

Conclusions

This study examined one of the largest single‐center US series of SUDEP cases. Several cases were reclassified as SUDEP, 15% had an ECG when alive, and 11 (11.4%) had blood for whole exome sequencing analysis. The most frequent antemortem genetic finding was pathogenic variants in SCN1A; postmortem whole exome sequencing identified 18 ultrarare variants.

Application of information theoretic feature selection and machine learning methods for the development of genetic risk prediction models

In view of the growth of clinical risk prediction models using genetic data, there is an increasing need for studies that use appropriate methods to select the optimum number of features from a large number of genetic variants with a high degree of redundancy between features due to linkage disequilibrium (LD). Filter feature selection methods based on information theoretic criteria, are well suited to this challenge and will identify a subset of the original variables that should result in more accurate prediction. However, data collected from cohort studies are often high-dimensional genetic data with potential confounders presenting challenges to feature selection and risk prediction machine learning models. Patients with psoriasis are at high risk of developing a chronic arthritis known as psoriatic arthritis (PsA). The prevalence of PsA in this patient group can be up to 30% and the identification of high risk patients represents an important clinical research which would allow early intervention and a reduction of disability. This also provides us with an ideal scenario for the development of clinical risk prediction models and an opportunity to explore the application of information theoretic criteria methods. In this study, we developed the feature selection and psoriatic arthritis (PsA) risk prediction models that were applied to a cross-sectional genetic dataset of 1462 PsA cases and 1132 cutaneous-only psoriasis (PsC) cases using 2-digit HLA alleles imputed using the SNP2HLA algorithm. We also developed stratification method to mitigate the impact of potential confounder features and illustrate that confounding features impact the feature selection. The mitigated dataset was used in training of seven supervised algorithms. 80% of data was randomly used for training of seven supervised machine learning methods using stratified nested cross validation and 20% was selected randomly as a holdout set for internal validation. The risk prediction models were then further validated in UK Biobank dataset containing data on 1187 participants and a set of features overlapping with the training dataset.Performance of these methods has been evaluated using the area under the curve (AUC), accuracy, precision, recall, F1 score and decision curve analysis(net benefit). The best model is selected based on three criteria: the ‘lowest number of feature subset’ with the ‘maximal average AUC over the nested cross validation’ and good generalisability to the UK Biobank dataset. In the original dataset, with over 100 different bootstraps and seven feature selection (FS) methods, HLA_C_*06 was selected as the most informative genetic variant. When the dataset is mitigated the single most important genetic features based on rank was identified as HLA_B_*27 by the seven different feature selection methods, consistent with previous analyses of this data using regression based methods. However, the predictive accuracy of these single features in post mitigation was found to be moderate (AUC= 0.54 (internal cross validation), AUC=0.53 (internal hold out set), AUC=0.55(external data set)). Sequentially adding additional HLA features based on rank improved the performance of the Random Forest classification model where 20 2-digit features selected by Interaction Capping (ICAP) demonstrated (AUC= 0.61 (internal cross validation), AUC=0.57 (internal hold out set), AUC=0.58 (external dataset)). The stratification method for mitigation of confounding features and filter information theoretic feature selection can be applied to a high dimensional dataset with the potential confounders.

100,000 Genomes Pilot on Rare-Disease Diagnosis in Health Care - Preliminary Report

BACKGROUND

The U.K. 100,000 Genomes Project is in the process of investigating the role of genome sequencing in patients with undiagnosed rare diseases after usual care and the alignment of this research with health care implementation in the U.K. National Health Service. Other parts of this project focus on patients with cancer and infection.

METHODS

We conducted a pilot study involving 4660 participants from 2183 families, among whom 161 disorders covering a broad spectrum of rare diseases were present. We collected data on clinical features with the use of Human Phenotype Ontology terms, undertook genome sequencing, applied automated variant prioritization on the basis of applied virtual gene panels and phenotypes, and identified novel pathogenic variants through research analysis.

RESULTS

Diagnostic yields varied among family structures and were highest in family trios (both parents and a proband) and families with larger pedigrees. Diagnostic yields were much higher for disorders likely to have a monogenic cause (35%) than for disorders likely to have a complex cause (11%). Diagnostic yields for intellectual disability, hearing disorders, and vision disorders ranged from 40 to 55%. We made genetic diagnoses in 25% of the probands. A total of 14% of the diagnoses were made by means of the combination of research and automated approaches, which was critical for cases in which we found etiologic noncoding, structural, and mitochondrial genome variants and coding variants poorly covered by exome sequencing. Cohortwide burden testing across 57,000 genomes enabled the discovery of three new disease genes and 19 new associations. Of the genetic diagnoses that we made, 25% had immediate ramifications for clinical decision making for the patients or their relatives.

CONCLUSIONS

Our pilot study of genome sequencing in a national health care system showed an increase in diagnostic yield across a range of rare diseases. (Funded by the National Institute for Health Research and others.).

New Homozygous Missense MSMO1 Mutation in Two Siblings with SC4MOL Deficiency Presenting with Psoriasiform Dermatitis

Sterol-C4-methyl oxidase (SC4MOL) deficiency was recently described as an autosomal recessive cholesterol biosynthesis disorder caused by mutations in the MSMO1 (sometimes also referred to as SC4MOL) gene. To date, 5 patients from 4 unrelated families with SC4MOL deficiency have been reported. Diagnosis can be challenging as the biochemical accumulation of methylsterols can affect global development and cause skin and ocular pathology. Herein, we describe 2 siblings from a consanguineous Turkish family with SC4MOL deficiency presenting with psoriasiform dermatitis, ocular abnormalities (nystagmus, optic hypoplasia, myopia, and strabismus), severe intellectual disability, and growth and motor delay. We undertook whole-exome sequencing and identified a new homozygous missense mutation c.81A>C; p.Asn27Thr in MSMO1. Segregation analysis in all available family members confirmed recessive inheritance of the mutation. The siblings were treated with a combination of oral and topical statin and cholesterol which resulted in clinical improvement. This study demonstrates how genomics-based diagnosis and therapy can be helpful in clinical practice.

Autosomal recessive hypotrichosis with loose anagen hairs associated with TKFC mutations

BACKGROUND

Loose anagen hair is a rare form of impaired hair anchorage in which anagen hairs that lack inner and outer root sheaths can be gently and painlessly plucked from the scalp. This condition usually occurs in children and is often self-limiting. A genetic basis for the disorder has been suggested but not proven. A better understanding the aetiology of loose anagen hair may improve prevention and treatment strategies.

OBJECTIVES

To identify a possible genetic basis of loose anagen hair using next-generation DNA sequencing and functional analysis of variants identified.

METHODS

In this case study, whole-exome sequencing analysis of a pedigree with one affected individual with features of loose anagen hair was performed.

RESULTS

The patient was found to be compound heterozygous for two single-nucleotide substitutions in TKFC resulting in the following missense mutations: c.574G> C (p.Gly192Arg) and c.682C> T (p.Arg228Trp). Structural analysis of human TKFC showed that both mutations are located near the active site cavity. Kinetic assays of recombinant proteins bearing either of these amino acid substitutions showed almost no dihydroxyacetone kinase or D-glyceraldehyde kinase activity, and FMN cyclase activity reduced to just 10% of wildtype catalytic activity.

CONCLUSIONS

TKFC missense mutations may predispose to the development of loose anagen hairs. Identification of this new biochemical pathobiology expands the metabolic and genetic basis of hypotrichosis.

Frequency of Pathogenic Germline Variants in CDH1, BRCA2, CHEK2, PALB2, BRCA1, and TP53 in Sporadic Lobular Breast Cancer

BACKGROUND

Invasive lobular breast cancer (ILC) accounts for approximately 15% of invasive breast carcinomas and is commonly associated with lobular carcinoma in situ (LCIS). Both have been shown to have higher familial risks than the more common ductal cancers. However, there are little data on the prevalence of the known high and moderate penetrance breast cancer predisposition genes in ILC. The aim of this study was to assess the frequency of germline variants in CDH1, BRCA2, BRCA1, CHEK2, PALB2, and TP53 in sporadic ILC and LCIS diagnosed in women ages ≤60 years.

METHODS

Access Array technology (Fluidigm) was used to amplify all exons of CDH1, BRCA2, BRCA1, TP53, CHEK2, and PALB2 using a custom-made targeted sequencing panel in 1,434 cases of ILC and 368 cases of pure LCIS together with 1,611 controls.

RESULTS

Case-control analysis revealed an excess of pathogenic variants in BRCA2, CHEK2, PALB2, and CDH1 in women with ILC. CHEK2 was the only gene that showed an association with pure LCIS [OR = 9.90; 95% confidence interval (CI), 3.42-28.66, P = 1.4 × 10^-5] with a larger effect size seen in LCIS compared with ILC (OR = 4.31; 95% CI, 1.61-11.58, P = 1.7 × 10^-3).

CONCLUSIONS

Eleven percent of patients with ILC ages ≤40 years carried germline variants in known breast cancer susceptibility genes.

IMPACT

Women with ILC ages ≤40 years should be offered genetic screening using a panel of genes that includes BRCA2, CHEK2, PALB2, and CDH1.

Phenotypic and genetic spectrum of epilepsy with myoclonic atonic seizures

OBJECTIVE

We aimed to describe the extent of neurodevelopmental impairments and identify the genetic etiologies in a large cohort of patients with epilepsy with myoclonic atonic seizures (MAE).

METHODS

We deeply phenotyped MAE patients for epilepsy features, intellectual disability, autism spectrum disorder, and attention-deficit/hyperactivity disorder using standardized neuropsychological instruments. We performed exome analysis (whole exome sequencing) filtered on epilepsy and neuropsychiatric gene sets to identify genetic etiologies.

RESULTS

We analyzed 101 patients with MAE (70% male). The median age of seizure onset was 34 months (range = 6-72 months). The main seizure types were myoclonic atonic or atonic in 100%, generalized tonic-clonic in 72%, myoclonic in 69%, absence in 60%, and tonic seizures in 19% of patients. We observed intellectual disability in 62% of patients, with extremely low adaptive behavioral scores in 69%. In addition, 24% exhibited symptoms of autism and 37% exhibited attention-deficit/hyperactivity symptoms. We discovered pathogenic variants in 12 (14%) of 85 patients, including five previously published patients. These were pathogenic genetic variants in SYNGAP1 (n = 3), KIAA2022 (n = 2), and SLC6A1 (n = 2), as well as KCNA2, SCN2A, STX1B, KCNB1, and MECP2 (n = 1 each). We also identified three new candidate genes, ASH1L, CHD4, and SMARCA2 in one patient each.

SIGNIFICANCE

MAE is associated with significant neurodevelopmental impairment. MAE is genetically heterogeneous, and we identified a pathogenic genetic etiology in 14% of this cohort by exome analysis. These findings suggest that MAE is a manifestation of several etiologies rather than a discrete syndromic entity.

Psoriasis and Genetics

Psoriasis is a common inflammatory skin disease caused by the interplay between multiple genetic and environmental risk factors. This review summarises recent progress in elucidating the genetic basis of psoriasis, particularly through large genome-wide association studies. We illustrate the power of genetic analyses for disease stratification. Psoriasis can be stratified by phenotype (common plaque versus rare pustular variants), or by outcome (prognosis, comorbidities, response to treatment); recent progress has been made in delineating the genetic contribution in each of these areas. We also highlight how genetic data can directly inform the development of effective psoriasis treatments.

Coagulation Factor XIII-A Subunit Missense Mutation in the Pathobiology of Autosomal Dominant Multiple Dermatofibromas

Dermatofibromas are common benign skin lesions, the etiology of which is poorly understood. We identified two unrelated pedigrees in which there was autosomal dominant transmission of multiple dermatofibromas. Whole exome sequencing revealed a rare shared heterozygous missense variant in the F13A1 gene encoding factor XIII subunit A (FXIII-A), a transglutaminase involved in hemostasis, wound healing, tumor growth, and apoptosis. The variant (p.Lys679Met) has an allele frequency of 0.0002 and is predicted to be a damaging mutation. Recombinant human Lys679Met FXIII-A demonstrated reduced fibrin crosslinking activity in vitro. Of note, the treatment of fibroblasts with media containing Lys679Met FXIII-A led to enhanced adhesion, proliferation, and type I collagen synthesis. Immunostaining revealed co-localization between FXIII-A and α4β1 integrins, more prominently for Lys679Met FXIII-A than the wild type. In addition, both the α4β1 inhibitors and the mutation of the FXIII-A Isoleucine-Leucine-Aspartate-Threonine (ILDT) motif prevented Lys679Met FXIII-A-dependent proliferation and collagen synthesis of fibroblasts. Our data suggest that the Lys679Met mutation may lead to a conformational change in the FXIII-A protein that enhances α4-integrin binding and provides insight into an unexpected role for FXIII-A in the pathobiology of familial dermatofibroma.

Importance of Variant Interpretation in Whole-Exome Molecular Autopsy: Population-Based Case Series

BACKGROUND

Potentially lethal cardiac channelopathies/cardiomyopathies may underlie a substantial portion of sudden unexplained death in the young (SUDY). The whole-exome molecular autopsy represents the latest approach to postmortem genetic testing for SUDY. However, proper variant adjudication in the setting of SUDY can be challenging.

METHODS

From January 2012 through December 2013, 25 consecutive cases of SUDY from 1 to 40 years of age (average age at death 27±5.7 years; 13 white, 12 black) from Cook County, Illinois, were referred after a negative (n=16) or equivocal (n=9) conventional autopsy. A whole-exome molecular autopsy with analysis of 99 sudden death-susceptibility genes was performed. The predicted pathogenicity of ultrarare, nonsynonymous variants was determined using the American College of Medical Genetics guidelines.

RESULTS

Overall, 27 ultrarare nonsynonymous variants were seen in 16/25 (64%) victims of SUDY. Among black individuals, 9/12 (75%) had an ultrarare nonsynonymous variant compared with 7/13 (54%) white individuals. Of the 27 variants, 10 were considered pathogenic or likely pathogenic in 7/25 (28%) individuals in accordance with the American College of Medical Genetics guidelines. Pathogenic/likely pathogenic variants were identified in 5/16 (31%) of autopsy-negative cases and in 2/6 (33%) victims of SUDY with equivocal findings of cardiomyopathy. Overall, 6 pathogenic/likely pathogenic variants in 4/25 (16%) cases were congruent with the phenotypic findings at autopsy and therefore considered clinically actionable.

CONCLUSIONS

Whole-exome molecular autopsy with gene-specific surveillance is an effective approach for the detection of potential pathogenic variants in SUDY cases. However, systematic variant adjudication is crucial to ensure accurate and proper care for surviving family members.

Cardiac Genetic Predisposition in Sudden Infant Death Syndrome

BACKGROUND

Sudden infant death syndrome (SIDS) is a leading cause of postneonatal mortality. Genetic heart diseases (GHDs) underlie some cases of SIDS.

OBJECTIVES

This study aimed to determine the spectrum and prevalence of GHD-associated mutations as a potential monogenic basis for SIDS.

METHODS

A cohort of 419 unrelated SIDS cases (257 male; average age 2.7 ± 1.9 months) underwent whole exome sequencing and a targeted analysis of 90 GHD-susceptibility genes. The yield of "potentially informative," ultra-rare variants (minor allele frequency <0.00005) in GHD-associated genes was assessed.

RESULTS

Overall, 53 of 419 (12.6%) SIDS cases had ≥1 "potentially informative," GHD-associated variant. The yield was 14.9% (21 of 141) for mixed-European ancestry cases and 11.5% (32 of 278) for European ancestry SIDS cases. Infants older than 4 months were more likely to host a "potentially informative" GHD-associated variant. There was significant overrepresentation of ultra-rare nonsynonymous variants in European SIDS cases (18 of 278 [6.5%]) versus European control subjects (30 of 973 [3.1%]; p = 0.013) when combining all 4 major cardiac channelopathy genes (KCNQ1, KCNH2, SCN5A, and RYR2). According to the American College of Medical Genetics guidelines, only 18 of 419 (4.3%) SIDS cases hosted a "pathogenic" or "likely pathogenic" variant.

CONCLUSIONS

Less than 15% of more than 400 SIDS cases had a "potentially informative" variant in a GHD-susceptibility gene, predominantly in the 4- to 12-month age group. Only 4.3% of cases possessed immediately clinically actionable variants. Consistent with previous studies, ultra-rare, nonsynonymous variants within the major cardiac channelopathy-associated genes were overrepresented in SIDS cases in infants of European ethnicity. These findings have major implications for the investigation of SIDS cases and families.

Mitochondria function associated genes contribute to Parkinson's Disease risk and later age at onset

Mitochondrial dysfunction has been implicated in the etiology of monogenic Parkinson's disease (PD). Yet the role that mitochondrial processes play in the most common form of the disease; sporadic PD, is yet to be fully established. Here, we comprehensively assessed the role of mitochondrial function-associated genes in sporadic PD by leveraging improvements in the scale and analysis of PD GWAS data with recent advances in our understanding of the genetics of mitochondrial disease. We calculated a mitochondrial-specific polygenic risk score (PRS) and showed that cumulative small effect variants within both our primary and secondary gene lists are significantly associated with increased PD risk. We further reported that the PRS of the secondary mitochondrial gene list was significantly associated with later age at onset. Finally, to identify possible functional genomic associations we implemented Mendelian randomization, which showed that 14 of these mitochondrial function-associated genes showed functional consequence associated with PD risk. Further analysis suggested that the 14 identified genes are not only involved in mitophagy, but implicate new mitochondrial processes. Our data suggests that therapeutics targeting mitochondrial bioenergetics and proteostasis pathways distinct from mitophagy could be beneficial to treating the early stage of PD.

Frequency of pathogenic germline variants in BRCA1, BRCA2, PALB2, CHEK2 and TP53 in ductal carcinoma in situ diagnosed in women under the age of 50 years

INTRODUCTION

Ductal carcinoma in situ (DCIS) is a non-obligate precursor of invasive ductal breast cancer, and approximately 20% of screen-detected tumours are pure DCIS. Most risk factors for breast cancer have similar associations with DCIS and IDC; however, there is limited data on the prevalence of the known high and moderate penetrance breast cancer predisposition genes in DCIS and which women with DCIS should be referred for genetic screening. The aim of this study was to assess the frequency of germline variants in BRCA2, BRCA1, CHEK2, PALB2 and TP53 in DCIS in women aged less than 50 years of age.

METHODS

After DNA extraction from the peripheral blood, Access Array technology (Fluidigm) was used to amplify all exons of these five known breast cancer predisposition genes using a custom made targeted sequencing panel in 655 cases of pure DCIS presenting in women under the age of 50 years together with 1611 controls.

RESULTS

Case-control analysis revealed an excess of pathogenic variants in BRCA2 (OR = 27.96, 95%CI 6.56-119.26, P = 2.0 × 10^-10) and CHEK2 (OR = 8.04, 95%CI 2.93-22.05, P = 9.0 × 10^-6), with weaker associations with PALB2 (P = 0.003), BRCA1 (P = 0.007) and TP53 (P = 0.02). For oestrogen receptor (ER)-positive DCIS the frequency of pathogenic variants was 9% under the age of 50 (14% with a family history of breast cancer) and 29% under the age of 40 (42% with a family history of breast cancer). For ER-negative DCIS, the frequency was 9% (16% with a family history of breast cancer) and 8% (11% with a family history of breast cancer) under the ages of 50 and 40, respectively.

CONCLUSIONS

This study has shown that breast tumourigenesis in women with pathogenic variants in BRCA2, CHEK2, PALB2, BRCA1 and TP53 can involve a DCIS precursor stage and that the focus of genetic testing in DCIS should be on women under the age of 40 with ER-positive DCIS.

Bi-allelic Loss-of-Function CACNA1B Mutations in Progressive Epilepsy-Dyskinesia

The occurrence of non-epileptic hyperkinetic movements in the context of developmental epileptic encephalopathies is an increasingly recognized phenomenon. Identification of causative mutations provides an important insight into common pathogenic mechanisms that cause both seizures and abnormal motor control. We report bi-allelic loss-of-function CACNA1B variants in six children from three unrelated families whose affected members present with a complex and progressive neurological syndrome. All affected individuals presented with epileptic encephalopathy, severe neurodevelopmental delay (often with regression), and a hyperkinetic movement disorder. Additional neurological features included postnatal microcephaly and hypotonia. Five children died in childhood or adolescence (mean age of death: 9 years), mainly as a result of secondary respiratory complications. CACNA1B encodes the pore-forming subunit of the pre-synaptic neuronal voltage-gated calcium channel Ca_v2.2/N-type, crucial for SNARE-mediated neurotransmission, particularly in the early postnatal period. Bi-allelic loss-of-function variants in CACNA1B are predicted to cause disruption of Ca²⁺ influx, leading to impaired synaptic neurotransmission. The resultant effect on neuronal function is likely to be important in the development of involuntary movements and epilepsy. Overall, our findings provide further evidence for the key role of Ca_v2.2 in normal human neurodevelopment.

Author Correction: Novel mutations in PIEZO1 cause an autosomal recessive generalized lymphatic dysplasia with non-immune hydrops fetalis

This Article contains an error in the last sentence of the 'Variant analysis suggests they are pathogenic' section of the Results, which incorrectly reads 'No truncated PIEZO1 protein products were identified in western blot analysis in GLD1:II.3 and GLD2:II.2 (Fig. 2, Supplementary Fig. 6), suggesting that the truncated protein is not stable and therefore degraded.' This should read 'No full-size PIEZO1 protein products were identified in western blot analysis in GLD1:II.3 and GLD2:II.2 (Fig. 2, Supplementary Fig. 6); the three nonsense mutations are predicted to lead to premature termination of the protein, hence it is possible that those truncated proteins will be non-functional or even unstable and degraded.' The error has not been fixed in the PDF or HTML versions of the Article.

Sequencing of human genomes with nanopore technology

Whole-genome sequencing (WGS) is becoming widely used in clinical medicine in diagnostic contexts and to inform treatment choice. Here we evaluate the potential of the Oxford Nanopore Technologies (ONT) MinION long-read sequencer for routine WGS by sequencing the reference sample NA12878 and the genome of an individual with ataxia-pancytopenia syndrome and severe immune dysregulation. We develop and apply a novel reference panel-free analytical method to infer and then exploit phase information which improves single-nucleotide variant (SNV) calling performance from otherwise modest levels. In the clinical sample, we identify and directly phase two non-synonymous de novo variants in SAMD9L, (OMIM #159550) inferring that they lie on the same paternal haplotype. Whilst consensus SNV-calling error rates from ONT data remain substantially higher than those from short-read methods, we demonstrate the substantial benefits of analytical innovation. Ongoing improvements to base-calling and SNV-calling methodology must continue for nanopore sequencing to establish itself as a primary method for clinical WGS.

The ADAMTS13-VWF axis is dysregulated in chronic thromboembolic pulmonary hypertension

Chronic thromboembolic pulmonary hypertension (CTEPH) is an important consequence of pulmonary embolism that is associated with abnormalities in haemostasis. We investigated the ADAMTS13–von Willebrand factor (VWF) axis in CTEPH, including its relationship with disease severity, inflammation, ABO groups and ADAMTS13 genetic variants.

ADAMTS13 and VWF plasma antigen levels were measured in patients with CTEPH (n=208), chronic thromboembolic disease without pulmonary hypertension (CTED) (n=35), resolved pulmonary embolism (n=28), idiopathic pulmonary arterial hypertension (n=30) and healthy controls (n=68). CTEPH genetic ABO associations and protein quantitative trait loci were investigated. ADAMTS13–VWF axis abnormalities were assessed in CTEPH and healthy control subsets by measuring ADAMTS13 activity, D-dimers and VWF multimeric size.

Patients with CTEPH had decreased ADAMTS13 (adjusted β −23.4%, 95% CI −30.9– −15.1%, p<0.001) and increased VWF levels (β +75.5%, 95% CI 44.8–113%, p<0.001) compared to healthy controls. ADAMTS13 levels remained low after reversal of pulmonary hypertension by pulmonary endarterectomy surgery and were equally reduced in CTED. We identified a genetic variant near the ADAMTS13 gene associated with ADAMTS13 protein that accounted for ∼8% of the variation in levels.

The ADAMTS13–VWF axis is dysregulated in CTEPH. This is unrelated to pulmonary hypertension, disease severity or markers of systemic inflammation and implicates the ADAMTS13–VWF axis in CTEPH pathobiology.

The ADAMTS-13–VWF axis is dysregulated in chronic thromboembolism with and without pulmonary hypertension and is implicated in the pathogenesishttp://ow.ly/J9SC30nh5T0

Genome-wide meta-analysis implicates mediators of hair follicle development and morphogenesis in risk for severe acne

Acne vulgaris is a highly heritable common, chronic inflammatory disease of the skin for which five genetic risk loci have so far been identified. Here, we perform a genome-wide association study of 3823 cases and 16,144 controls followed by meta-analysis with summary statistics from a previous study, with a total sample size of 26,722. We identify 20 independent association signals at 15 risk loci, 12 of which have not been previously implicated in the disease. Likely causal variants disrupt the coding region of WNT10A and a P63 transcription factor binding site in SEMA4B. Risk alleles at the 1q25 locus are associated with increased expression of LAMC2, in which biallelic loss-of-function mutations cause the blistering skin disease epidermolysis bullosa. These findings indicate that variation affecting the structure and maintenance of the skin, in particular the pilosebaceous unit, is a critical aspect of the genetic predisposition to severe acne.

Acne vulgaris is a chronic inflammation of the skin, the genetic basis of which is incompletely understood. Here, Petridis et al. perform GWAS and meta-analysis for acne in 26,722 individuals and identify 12 novel risk loci that implicate structure and maintenance of the skin in severe acne risk.

Exome-Wide Rare Variant Analyses in Sudden Infant Death Syndrome

OBJECTIVE

To determine whether a monogenic basis explains sudden infant death syndrome (SIDS) using an exome-wide focus.

STUDY DESIGN

A cohort of 427 unrelated cases of SIDS (257 male; average age = 2.7 ± 1.9 months) underwent whole-exome sequencing. Exome-wide rare variant analyses were carried out with 278 SIDS cases of European ancestry (173 male; average age = 2.7 ± 1.98 months) and 973 ethnic-matched controls based on 6 genetic models. Ingenuity Pathway Analysis also was performed. The cohort was collected in collaboration with coroners, medical examiners, and pathologists by St George's University of London, United Kingdom, and Mayo Clinic, Rochester, Minnesota. Whole-exome sequencing was performed at the Genomic Laboratory, Kings College London, United Kingdom, or Mayo Clinic's Medical Genome Facility, Rochester, Minnesota.

RESULTS

Although no exome-wide significant (P < 2.5 × 10-6) difference in burden of ultra-rare variants was detected for any gene, 405 genes had a greater prevalence (P < .05) of ultra-rare nonsynonymous variants among cases with 17 genes at P < .005. Some of these potentially overrepresented genes may represent biologically plausible novel candidate genes for a monogenic basis for a portion of patients with SIDS. The top canonical pathway identified was glucocorticoid biosynthesis (P = .01).

CONCLUSIONS

The lack of exome-wide significant genetic associations indicates an extreme heterogeneity of etiologies underlying SIDS. Our approach to understanding the genetic mechanisms of SIDS has far reaching implications for the SIDS research community as a whole and may catalyze new evidence-based SIDS research across multiple disciplines. Perturbations in glucocorticoid biosynthesis may represent a novel SIDS-associated biological pathway for future SIDS investigative research.

WNT10A mutation results in severe tooth agenesis in a family of three sisters

OBJECTIVES

To identify the genetic basis of severe tooth agenesis in a family of three affected sisters.

PATIENTS AND METHODS

A family of three sisters with severe tooth agenesis was recruited for whole-exome sequencing to identify potential genetic variation responsible for this penetrant phenotype. The unaffected father was tested for specific mutations using Sanger sequencing. Gene discovery was supplemented with in situ hybridization to localize gene expression during human tooth development.

RESULTS

We report a nonsense heterozygous mutation in exon 2 of WNT10A c.321C>A[p.Cys107*] likely to be responsible for the severe tooth agenesis identified in this family through the creation of a premature stop codon, resulting in truncation of the amino acid sequence and therefore loss of protein function. In situ hybridization showed expression of WNT10A in odontogenic epithelium during the early and late stages of human primary tooth development.

CONCLUSIONS

WNT10A has previously been associated with both syndromic and non-syndromic forms of tooth agenesis, and this report further expands our knowledge of genetic variation underlying non-syndromic forms of this condition. We also demonstrate expression of WNT10A in the epithelial compartment of human tooth germs during development.

A genome-wide association study for extremely high intelligence

We used a case-control genome-wide association (GWA) design with cases consisting of 1238 individuals from the top 0.0003 (~170 mean IQ) of the population distribution of intelligence and 8172 unselected population-based controls. The single-nucleotide polymorphism heritability for the extreme IQ trait was 0.33 (0.02), which is the highest so far for a cognitive phenotype, and significant genome-wide genetic correlations of 0.78 were observed with educational attainment and 0.86 with population IQ. Three variants in locus ADAM12 achieved genome-wide significance, although they did not replicate with published GWA analyses of normal-range IQ or educational attainment. A genome-wide polygenic score constructed from the GWA results accounted for 1.6% of the variance of intelligence in the normal range in an unselected sample of 3414 individuals, which is comparable to the variance explained by GWA studies of intelligence with substantially larger sample sizes. The gene family plexins, members of which are mutated in several monogenic neurodevelopmental disorders, was significantly enriched for associations with high IQ. This study shows the utility of extreme trait selection for genetic study of intelligence and suggests that extremely high intelligence is continuous genetically with normal-range intelligence in the population.

Report from the fifth international consensus meeting to harmonize core outcome measures for atopic eczema/dermatitis clinical trials (HOME initiative)

This is the report from the fifth meeting of the Harmonising Outcome Measures for Eczema initiative (HOME V). The meeting was held on 12-14 June 2017 in Nantes, France, with 81 participants. The main aims of the meeting were (i) to achieve consensus over the definition of the core domain of long-term control and how to measure it and (ii) to prioritize future areas of research for the measurement of the core domain of quality of life (QoL) in children. Moderated whole-group and small-group consensus discussions were informed by presentations of qualitative studies, systematic reviews and validation studies. Small-group allocations were performed a priori to ensure that each group included different stakeholders from a variety of geographical regions. Anonymous whole-group voting was carried out using handheld electronic voting pads according to predefined consensus rules. It was agreed by consensus that the long-term control domain should include signs, symptoms, quality of life and a patient global instrument. The group agreed that itch intensity should be measured when assessing long-term control of eczema in addition to the frequency of itch captured by the symptoms domain. There was no recommendation of an instrument for the core outcome domain of quality of life in children, but existing instruments were assessed for face validity and feasibility, and future work that will facilitate the recommendation of an instrument was agreed upon.

Dysfunction of NaV1.4, a skeletal muscle voltage-gated sodium channel, in sudden infant death syndrome: a case-control study

BACKGROUND

Sudden infant death syndrome (SIDS) is the leading cause of post-neonatal infant death in high-income countries. Central respiratory system dysfunction seems to contribute to these deaths. Excitation that drives contraction of skeletal respiratory muscles is controlled by the sodium channel NaV1.4, which is encoded by the gene SCN4A. Variants in NaV1.4 that directly alter skeletal muscle excitability can cause myotonia, periodic paralysis, congenital myopathy, and myasthenic syndrome. SCN4A variants have also been found in infants with life-threatening apnoea and laryngospasm. We therefore hypothesised that rare, functionally disruptive SCN4A variants might be over-represented in infants who died from SIDS.

METHODS

We did a case-control study, including two consecutive cohorts that included 278 SIDS cases of European ancestry and 729 ethnically matched controls without a history of cardiovascular, respiratory, or neurological disease. We compared the frequency of rare variants in SCN4A between groups (minor allele frequency <0·00005 in the Exome Aggregation Consortium). We assessed biophysical characterisation of the variant channels using a heterologous expression system.

FINDINGS

Four (1·4%) of the 278 infants in the SIDS cohort had a rare functionally disruptive SCN4A variant compared with none (0%) of 729 ethnically matched controls (p=0·0057).

INTERPRETATION

Rare SCN4A variants that directly alter NaV1.4 function occur in infants who had died from SIDS. These variants are predicted to significantly alter muscle membrane excitability and compromise respiratory and laryngeal function. These findings indicate that dysfunction of muscle sodium channels is a potentially modifiable risk factor in a subset of infant sudden deaths.

FUNDING

UK Medical Research Council, the Wellcome Trust, National Institute for Health Research, the British Heart Foundation, Biotronik, Cardiac Risk in the Young, Higher Education Funding Council for England, Dravet Syndrome UK, the Epilepsy Society, the Eunice Kennedy Shriver National Institute of Child Health & Human Development of the National Institutes of Health, and the Mayo Clinic Windland Smith Rice Comprehensive Sudden Cardiac Death Program.

Text-mined phenotype annotation and vector-based similarity to improve identification of similar phenotypes and causative genes in monogenic disease patients

The genetic diagnosis of rare monogenic diseases using exome/genome sequencing requires the true causal variant(s) to be identified from tens of thousands of observed variants. Typically a virtual gene panel approach is taken whereby only variants in genes known to cause phenotypes resembling the patient under investigation are considered. With the number of known monogenic gene-disease pairs exceeding 5,000, manual curation of personalized virtual panels using exhaustive knowledge of the genetic basis of the human monogenic phenotypic spectrum is challenging. We present improved probabilistic methods for estimating phenotypic similarity based on Human Phenotype Ontology annotation. A limitation of existing methods for evaluating a disease's similarity to a reference set is that reference diseases are typically represented as a series of binary (present/absent) observations of phenotypic terms. We evaluate a quantified disease reference set, using term frequency in phenotypic text descriptions to approximate term relevance. We demonstrate an improved ability to identify related diseases through the use of a quantified reference set, and that vector space similarity measures perform better than established information content-based measures. These improvements enable the generation of bespoke virtual gene panels, facilitating more accurate and efficient interpretation of genomic variant profiles from individuals with rare Mendelian disorders. These methods are available online at https://atlas.genetics.kcl.ac.uk/~jake/cgi-bin/patient_sim.py.

Biallelic Mutations in ATP5F1D, which Encodes a Subunit of ATP Synthase, Cause a Metabolic Disorder

ATP synthase, H+ transporting, mitochondrial F1 complex, δ subunit (ATP5F1D; formerly ATP5D) is a subunit of mitochondrial ATP synthase and plays an important role in coupling proton translocation and ATP production. Here, we describe two individuals, each with homozygous missense variants in ATP5F1D, who presented with episodic lethargy, metabolic acidosis, 3-methylglutaconic aciduria, and hyperammonemia. Subject 1, homozygous for c.245C>T (p.Pro82Leu), presented with recurrent metabolic decompensation starting in the neonatal period, and subject 2, homozygous for c.317T>G (p.Val106Gly), presented with acute encephalopathy in childhood. Cultured skin fibroblasts from these individuals exhibited impaired assembly of F1FO ATP synthase and subsequent reduced complex V activity. Cells from subject 1 also exhibited a significant decrease in mitochondrial cristae. Knockdown of Drosophila ATPsynδ, the ATP5F1D homolog, in developing eyes and brains caused a near complete loss of the fly head, a phenotype that was fully rescued by wild-type human ATP5F1D. In contrast, expression of the ATP5F1D c.245C>T and c.317T>G variants rescued the head-size phenotype but recapitulated the eye and antennae defects seen in other genetic models of mitochondrial oxidative phosphorylation deficiency. Our data establish c.245C>T (p.Pro82Leu) and c.317T>G (p.Val106Gly) in ATP5F1D as pathogenic variants leading to a Mendelian mitochondrial disease featuring episodic metabolic decompensation.

Exome Sequencing and Genotyping Identify a Rare Variant in NLRP7 Gene Associated With Ulcerative Colitis

BACKGROUND AND AIMS

Although genome-wide association studies [GWAS] in inflammatory bowel disease [IBD] have identified a large number of common disease susceptibility alleles for both Crohn's disease [CD] and ulcerative colitis [UC], a substantial fraction of IBD heritability remains unexplained, suggesting that rare coding genetic variants may also have a role in pathogenesis. We used high-throughput sequencing in families with multiple cases of IBD, followed by genotyping of cases and controls, to investigate whether rare protein-altering genetic variants are associated with susceptibility to IBD.

METHODS

Whole-exome sequencing was carried out in 10 families in whom three or more individuals were affected with IBD. A stepwise filtering approach was applied to exome variants, to identify potential causal variants. Follow-up genotyping was performed in 6025 IBD cases [2948 CD; 3077 UC] and 7238 controls.

RESULTS

Our exome variant analysis revealed coding variants in the NLRP7 gene that were present in affected individuals in two distinct families. Genotyping of the two variants, p.S361L and p.R801H, in IBD cases and controls showed that the p.S361L variant was significantly associated with an increased risk of ulcerative colitis [odds ratio 4.79, p = 0.0039] and IBD [odds ratio 3.17, p = 0.037]. A combined analysis of both variants showed suggestive association with an increased risk of IBD [odds ratio 2.77, p = 0.018].

CONCLUSIONS

The results suggest that NLRP7 signalling and inflammasome formation may be a significant component in the pathogenesis of IBD.

De novo mutations implicate novel genes in systemic lupus erythematosus

The omnigenic model of complex disease stipulates that the majority of the heritability will be explained by the effects of common variation on genes in the periphery of core disease pathways. Rare variant associations, expected to explain far less of the heritability, may be enriched in core disease genes and thus will be instrumental in the understanding of complex disease pathogenesis and their potential therapeutic targets. Here, using complementary whole-exome sequencing, high-density imputation, and in vitro cellular assays, we identify candidate core genes in the pathogenesis of systemic lupus erythematosus (SLE). Using extreme-phenotype sampling, we sequenced the exomes of 30 SLE parent-affected-offspring trios and identified 14 genes with missense de novo mutations (DNM), none of which are within the >80 SLE susceptibility loci implicated through genome-wide association studies. In a follow-up cohort of 10, 995 individuals of matched European ancestry, we imputed genotype data to the density of the combined UK10K-1000 genomes Phase III reference panel across the 14 candidate genes. Gene-level analyses indicate three functional candidates: DNMT3A, PRKCD, and C1QTNF4. We identify a burden of rare variants across PRKCD associated with SLE risk (P = 0.0028), and across DNMT3A associated with two severe disease prognosis sub-phenotypes (P = 0.0005 and P = 0.0033). We further characterise the TNF-dependent functions of the third candidate gene C1QTNF4 on NF-κB activation and apoptosis, which are inhibited by the p.His198Gln DNM. Our results identify three novel genes in SLE susceptibility and support extreme-phenotype sampling and DNM gene discovery to aid the search for core disease genes implicated through rare variation.

MED12, TERT promoter and RBM15 mutations in primary and recurrent phyllodes tumours

BACKGROUND

MED12 and TERT promoter mutations have been shown to be the most common somatic mutations in phyllodes tumours (PTs). The aims of this study were to determine the frequency of these mutations in recurrent PTs, assess whether TERT promoter mutations could be helpful in distinguishing fibroadenomas (FAs) from PTs and identify novel mutations that may be driving malignant progression.

METHODS

MED12 and the TERT promoter were Sanger sequenced in 75 primary PTs, 21 recurrences, 19 single FAs and 2 cases of multiple FAs with benign PTs. Whole-exome sequencing was performed on one borderline PT.

RESULTS

Recurrent PTs and multiple FAs showed temporal discordance in MED12 but not TERT. Recurrent samples did acquire TERT mutations, with recurrent benign PTs more likely to have mutations in both genes. TERT mutations were not helpful in differentiating between benign PTs and FAs in cases of multiple FAs/PTs. Exome sequencing revealed a nonsense mutation in RBM15 and Sanger sequencing revealed another three RBM15 mutations in malignant/borderline PTs.

CONCLUSIONS

This study has shown that MED12 mutations can be heterogeneous in both synchronous and recurrent PTs unlike TERT mutations. We have also shown that RBM15 mutations may be important in the pathogenesis of borderline/malignant PTs.

Neurologic Phenotypes Associated With Mutations in RTN4IP1 (OPA10) in Children and Young Adults

Importance

Neurologic disorders with isolated symptoms or complex syndromes are relatively frequent among mitochondrial inherited diseases. Recessive RTN4IP1 gene mutations have been shown to cause isolated and syndromic optic neuropathies.

Objective

To define the spectrum of clinical phenotypes associated with mutations in RTN4IP1 encoding a mitochondrial quinone oxidoreductase.

Design, Setting, and Participants

This study involved 12 individuals from 11 families with severe central nervous system diseases and optic atrophy. Targeted and whole-exome sequencing were performed-at Hospital Angers (France), Institute of Neurology Milan (Italy), Imagine Institute Paris (France), Helmoltz Zentrum of Munich (Germany), and Beijing Genomics Institute (China)-to clarify the molecular diagnosis of patients. Each patient's neurologic, ophthalmologic, magnetic resonance imaging, and biochemical features were investigated. This study was conducted from May 1, 2014, to June 30, 2016.

Main Outcomes and Measures

Recessive mutations in RTN4IP1 were identified. Clinical presentations ranged from isolated optic atrophy to severe encephalopathies.

Results

Of the 12 individuals in the study, 6 (50%) were male and 6 (50%) were female. They ranged in age from 5 months to 32 years. Of the 11 families, 6 (5 of whom were consanguineous) had a member or members who presented isolated optic atrophy with the already reported p.Arg103His or the novel p.Ile362Phe, p.Met43Ile, and p.Tyr51Cys amino acid changes. The 5 other families had a member or members who presented severe neurologic syndromes with a common core of symptoms, including optic atrophy, seizure, intellectual disability, growth retardation, and elevated lactate levels. Additional clinical features of those affected were deafness, abnormalities on magnetic resonance images of the brain, stridor, and abnormal electroencephalographic patterns, all of which eventually led to death before age 3 years. In these patients, novel and very rare homozygous and compound heterozygous mutations were identified that led to the absence of the protein and complex I disassembly as well as mild mitochondrial network fragmentation.

Conclusions and Relevance

A broad clinical spectrum of neurologic features, ranging from isolated optic atrophy to severe early-onset encephalopathies, is associated with RTN4IP1 biallelic mutations and should prompt RTN4IP1 screening in both syndromic neurologic presentations and nonsyndromic recessive optic neuropathies.

Novel GFM2 variants associated with early-onset neurological presentations of mitochondrial disease and impaired expression of OXPHOS subunits

Mitochondrial diseases are characterised by clinical, molecular and functional heterogeneity, reflecting their bi-genomic control. The nuclear gene GFM2 encodes mtEFG2, a protein with an essential role during the termination stage of mitochondrial translation. We present here two unrelated patients harbouring different and previously unreported compound heterozygous (c.569G>A, p.(Arg190Gln); c.636delA, p.(Glu213Argfs*3)) and homozygous (c.275A>C, p.(Tyr92Ser)) recessive variants in GFM2 identified by whole exome sequencing (WES) together with histochemical and biochemical findings to support the diagnoses of pathological GFM2 variants in each case. Both patients presented similarly in early childhood with global developmental delay, raised CSF lactate and abnormalities on cranial MRI. Sanger sequencing of familial samples confirmed the segregation of bi-allelic GFM2 variants with disease, while investigations into steady-state mitochondrial protein levels revealed respiratory chain subunit defects and loss of mtEFG2 protein in muscle. These data demonstrate the effects of defective mtEFG2 function, caused by previously unreported variants, confirming pathogenicity and expanding the clinical phenotypes associated with GFM2 variants.