Rare variant contribution to human disease in 281,104 UK Biobank exomes.

Genome-wide association studies have uncovered thousands of common variants associated with human disease, but the contribution of rare variants to common disease remains relatively unexplored. The UK Biobank contains detailed phenotypic data linked to medical records for approximately 500,000 participants, offering an unprecedented opportunity to evaluate the effect of rare variation on a broad collection of traits1,2. Here we study the relationships between rare protein-coding variants and 17,361 binary and 1,419 quantitative phenotypes using exome sequencing data from 269,171 UK Biobank participants of European ancestry. Gene-based collapsing analyses revealed 1,703 statistically significant gene-phenotype associations for binary traits, with a median odds ratio of 12.4. Furthermore, 83% of these associations were undetectable via single-variant association tests, emphasizing the power of gene-based collapsing analysis in the setting of high allelic heterogeneity. Gene-phenotype associations were also significantly enriched for loss-of-function-mediated traits and approved drug targets. Finally, we performed ancestry-specific and pan-ancestry collapsing analyses using exome sequencing data from 11,933 UK Biobank participants of African, East Asian or South Asian ancestry. Our results highlight a significant contribution of rare variants to common disease. Summary statistics are publicly available through an interactive portal ( http://azphewas.com/ )

Prevalence and clinical challenges among adults with primary immunodeficiency and recombination-activating gene deficiency.

Letter to the Edito

Phenotypic Characterization of EIF2AK4 Mutation Carriers in a Large Cohort of Patients Diagnosed Clinically With Pulmonary Arterial Hypertension.

BACKGROUND: Pulmonary arterial hypertension (PAH) is a rare disease with an emerging genetic basis. Heterozygous mutations in the gene encoding the bone morphogenetic protein receptor type 2 (BMPR2) are the commonest genetic cause of PAH, whereas biallelic mutations in the eukaryotic translation initiation factor 2 alpha kinase 4 gene (EIF2AK4) are described in pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Here, we determine the frequency of these mutations and define the genotype-phenotype characteristics in a large cohort of patients diagnosed clinically with PAH. METHODS: Whole-genome sequencing was performed on DNA from patients with idiopathic and heritable PAH and with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis recruited to the National Institute of Health Research BioResource-Rare Diseases study. Heterozygous variants in BMPR2 and biallelic EIF2AK4 variants with a minor allele frequency of <1:10 000 in control data sets and predicted to be deleterious (by combined annotation-dependent depletion, PolyPhen-2, and sorting intolerant from tolerant predictions) were identified as potentially causal. Phenotype data from the time of diagnosis were also captured. RESULTS: Eight hundred sixty-four patients with idiopathic or heritable PAH and 16 with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis were recruited. Mutations in BMPR2 were identified in 130 patients (14.8%). Biallelic mutations in EIF2AK4 were identified in 5 patients with a clinical diagnosis of pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Furthermore, 9 patients with a clinical diagnosis of PAH carried biallelic EIF2AK4 mutations. These patients had a reduced transfer coefficient for carbon monoxide (Kco; 33% [interquartile range, 30%-35%] predicted) and younger age at diagnosis (29 years; interquartile range, 23-38 years) and more interlobular septal thickening and mediastinal lymphadenopathy on computed tomography of the chest compared with patients with PAH without EIF2AK4 mutations. However, radiological assessment alone could not accurately identify biallelic EIF2AK4 mutation carriers. Patients with PAH with biallelic EIF2AK4 mutations had a shorter survival. CONCLUSIONS: Biallelic EIF2AK4 mutations are found in patients classified clinically as having idiopathic and heritable PAH. These patients cannot be identified reliably by computed tomography, but a low Kco and a young age at diagnosis suggests the underlying molecular diagnosis. Genetic testing can identify these misclassified patients, allowing appropriate management and early referral for lung transplantation

Phenotypic Characterization of <i>EIF2AK4</i> Mutation Carriers in a Large Cohort of Patients Diagnosed Clinically With Pulmonary Arterial Hypertension

Background: Pulmonary arterial hypertension (PAH) is a rare disease with an emerging genetic basis. Heterozygous mutations in the gene encoding the bone morphogenetic protein receptor type 2 ( BMPR2 ) are the commonest genetic cause of PAH, whereas biallelic mutations in the eukaryotic translation initiation factor 2 alpha kinase 4 gene ( EIF2AK4 ) are described in pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Here, we determine the frequency of these mutations and define the genotype-phenotype characteristics in a large cohort of patients diagnosed clinically with PAH. Methods: Whole-genome sequencing was performed on DNA from patients with idiopathic and heritable PAH and with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis recruited to the National Institute of Health Research BioResource–Rare Diseases study. Heterozygous variants in BMPR2 and biallelic EIF2AK4 variants with a minor allele frequency of &lt;1:10 000 in control data sets and predicted to be deleterious (by combined annotation-dependent depletion, PolyPhen-2, and sorting intolerant from tolerant predictions) were identified as potentially causal. Phenotype data from the time of diagnosis were also captured. Results: Eight hundred sixty-four patients with idiopathic or heritable PAH and 16 with pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis were recruited. Mutations in BMPR2 were identified in 130 patients (14.8%). Biallelic mutations in EIF2AK4 were identified in 5 patients with a clinical diagnosis of pulmonary veno-occlusive disease/pulmonary capillary hemangiomatosis. Furthermore, 9 patients with a clinical diagnosis of PAH carried biallelic EIF2AK4 mutations. These patients had a reduced transfer coefficient for carbon monoxide (K co ; 33% [interquartile range, 30%–35%] predicted) and younger age at diagnosis (29 years; interquartile range, 23–38 years) and more interlobular septal thickening and mediastinal lymphadenopathy on computed tomography of the chest compared with patients with PAH without EIF2AK4 mutations. However, radiological assessment alone could not accurately identify biallelic EIF2AK4 mutation carriers. Patients with PAH with biallelic EIF2AK4 mutations had a shorter survival. Conclusions: Biallelic EIF2AK4 mutations are found in patients classified clinically as having idiopathic and heritable PAH. These patients cannot be identified reliably by computed tomography, but a low K co and a young age at diagnosis suggests the underlying molecular diagnosis. Genetic testing can identify these misclassified patients, allowing appropriate management and early referral for lung transplantation. </jats:sec

Publisher Correction: Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper

How common are single gene mutations as a cause for lacunar stroke? A targeted gene panel study.

OBJECTIVES: To determine the frequency of rare and pertinent disease-causing variants in small vessel disease (SVD)-associated genes (such as NOTCH3, HTRA1, COL4A1, COL4A2, FOXC1, TREX1, and GLA) in cerebral SVD, we performed targeted gene sequencing in 950 patients with younger-onset apparently sporadic SVD stroke using a targeted sequencing panel. METHODS: We designed a high-throughput sequencing panel to identify variants in 15 genes (7 known SVD genes, 8 SVD-related disorder genes). The panel was used to screen a population of 950 patients with younger-onset (≤70 years) MRI-confirmed SVD stroke, recruited from stroke centers across the United Kingdom. Variants were filtered according to their frequency in control databases, predicted effect, presence in curated variant lists, and combined annotation dependent depletion scores. Whole genome sequencing and genotyping were performed on a subset of patients to provide a direct comparison of techniques. The frequency of known disease-causing and pertinent variants of uncertain significance was calculated. RESULTS: We identified previously reported variants in 14 patients (8 cysteine-changing NOTCH3 variants in 11 patients, 2 HTRA1 variants in 2 patients, and 1 missense COL4A1 variant in 1 patient). In addition, we identified 29 variants of uncertain significance in 32 patients. CONCLUSION: Rare monogenic variants account for about 1.5% of younger onset lacunar stroke. Most are cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy variants, but the second most common gene affected is HTRA1. A high-throughput sequencing technology platform is an efficient, reliable method to screen for such mutations.Funding for the project was provided by the National Institute for Health Research (NIHR, grant number RG65966). This work was supported by a BHF Programme Grant (RG/16/4/32218). Infrastructural support was provided by the Cambridge University Hospitals Trust Biomedical Research Centre. RYYT is supported by the Agency for Science, Technology and Research (Singapore). KD is a HSST trainee supported by NHS Health Education England. HM is supported by a NIHR Senior Investigator award

ADA2 deficiency complicated by EBV-driven lymphoproliferative disease.

A 29-year old male with recurrent respiratory and skin infections, anaemia and neutropaenia during childhood required immunoglobulin replacement for antibody deficiency from age 16. He remained relatively well until age 28 when he presented with a two-week history of fatigue, sore throat, fever and productive cough. He was found to have EBV viraemia and splenomegaly and a diagnosis of EBV-driven lymphoproliferative disease was made following bone marrow trephine. Family history was notable with three siblings: a healthy sister and two brothers with anaemia and neutropaenia; one who succumbed to septicaemia secondary to neutropaenic enterocolitis age 5 and another who developed intestinal vasculitis and antibody deficiency and had a successful haemopoetic stem cell transplant. The proband's DNA underwent targeted sequencing of 279 genes associated with immunodeficiency (GRID panel). The best candidates were two ADA2 variants, p.Arg169Gln (R169Q) and p.Asn370Lys (N370K). Sanger sequencing and co-segregation of variants in the parents, unaffected sister and all three affected brothers was fully consistent with compound heterozygous inheritance. Subsequent whole genome sequencing of the proband identified no other potential causal variants. ADA2 activity was consistent with a diagnosis of ADA2 deficiency in affected family members. This is the first description of EBV-driven lymphoproliferative disease in ADA2 deficiency. ADA2 deficiency may cause susceptibility to severe EBV-induced disease and we would recommend that EBV status and viral load is monitored in patients with this diagnosis and allogeneic SCT is considered at an early stage for patients whose ADA2 deficiency is associated with significant complications

Next-generation sequencing for the diagnosis of MYH9-RD: Predicting pathogenic variants

The heterogeneous manifestations of MYH9-related disorder (MYH9-RD), characterized by macrothrombocytopenia, Döhle-like inclusion bodies in leukocytes, bleeding of variable severity with, in some cases, ear, eye, kidney, and liver involvement, make the diagnosis for these patients still challenging in clinical practice. We collected phenotypic data and analyzed the genetic variants in more than 3,000 patients with a bleeding or platelet disorder. Patients were enrolled in the BRIDGE-BPD and ThromboGenomics Projects and their samples processed by high throughput sequencing (HTS). We identified 50 patients with a rare variant in MYH9. All patients had macrothrombocytes and all except two had thrombocytopenia. Some degree of bleeding diathesis was reported in 41 of the 50 patients. Eleven patients presented hearing impairment, three renal failure and two elevated liver enzymes. Among the 28 rare variants identified in MYH9, 12 were novel. HTS was instrumental in diagnosing 23 patients (46%). Our results confirm the clinical heterogeneity of MYH9-RD and show that, in the presence of an unclassified platelet disorder with macrothrombocytes, MYH9-RD should always be considered. A HTS-based strategy is a reliable method to reach a conclusive diagnosis of MYH9-RD in clinical practice.status: publishe

Spontaneous Coronary Artery Dissection: Insights on Rare Genetic Variation from Genome Sequencing

Background - Spontaneous coronary artery dissection (SCAD) occurs when an epicardial coronary artery is narrowed or occluded by an intramural hematoma. SCAD mainly affects women and is associated with pregnancy and systemic arteriopathies, particularly fibromuscular dysplasia. Variants in several genes, such as those causing connective tissue disorders, have been implicated; however, the genetic architecture is poorly understood. Here, we aim to better understand the diagnostic yield of rare variant genetic testing among a cohort of SCAD survivors and to identify genes or gene-sets that have a significant enrichment of rare variants.Methods - We sequenced a cohort of 384 SCAD survivors from the UK, alongside 13,722 UK Biobank controls and a validation cohort of 92 SCAD survivors. We performed a research diagnostic screen for pathogenic variants, and exome-wide and gene-set rare variant collapsing analyses.Results - The majority of patients within both cohorts are female, 29% of the study cohort and 14% validation cohort have a remote arteriopathy. Four cases across the two cohorts had a diagnosed connective tissue disorder. We identified pathogenic or likely pathogenic variants in seven genes (PKD1, COL3A1, SMAD3, TGFB2, LOX, MYLK, and YY1AP1) in 14/384 cases in the study cohort and in 1/92 cases in the validation cohort. In our rare variant collapsing analysis, PKD1 was the highest ranked gene and several functionally plausible genes were enriched for rare variants, although no gene achieved study-wide statistical significance. Gene-set enrichment analysis suggested a role for additional genes involved in renal function.Conclusions - By studying the largest sequenced cohort of SCAD survivors we demonstrate that, based on current knowledge, only a small proportion have a pathogenic variant that could explain their disease. Our findings strengthen the overlap between SCAD and renal and connective tissue disorders and we highlight several new genes for future validation.</div

A dominant gain-of-function mutation in universal tyrosine kinase SRC causes thrombocytopenia, myelofibrosis, bleeding, and bone pathologies

The Src family kinase (SFK) member SRC is a major target in drug development because it is activated in many human cancers, yet deleterious SRC germline mutations have not been reported. We used genome sequencing and Human Phenotype Ontology patient coding to identify a gain-of-function mutation in SRC causing thrombocytopenia, myelofibrosis, bleeding, and bone pathologies in nine cases. Modeling of the E527K substitution predicts loss of SRC's self-inhibitory capacity, which we confirmed with in vitro studies showing increased SRC kinase activity and enhanced Tyr(419) phosphorylation in COS-7 cells overexpressing E527K SRC. The active form of SRC predominates in patients' platelets, resulting in enhanced overall tyrosine phosphorylation. Patients with myelofibrosis have hypercellular bone marrow with trilineage dysplasia, and their stem cells grown in vitro form more myeloid and megakaryocyte (MK) colonies than control cells. These MKs generate platelets that are dysmorphic, low in number, highly variable in size, and have a paucity of α-granules. Overactive SRC in patient-derived MKs causes a reduction in proplatelet formation, which can be rescued by SRC kinase inhibition. Stem cells transduced with lentiviral E527K SRC form MKs with a similar defect and enhanced tyrosine phosphorylation levels. Patient-derived and E527K-transduced MKs show Y419 SRC-positive stained podosomes that induce altered actin organization. Expression of mutated src in zebrafish recapitulates patients' blood and bone phenotypes. Similar studies of platelets and MKs may reveal the mechanism underlying the severe bleeding frequently observed in cancer patients treated with next-generation SFK inhibitors.status: publishe