Abstract 4489: Using functional data from Roadmap Epigenomics to inform analysis of rare variants linked to gene expression in a large colorectal cancer study

To investigate the role of low frequency and rare genetic variation in colorectal cancer (CRC) susceptibility, the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO) and the Colorectal Cancer Family Registry (CCFR) conducted whole genome sequencing and imputed into genome-wide association studies (GWAS) of 14,718 CRC cases and 12,186 controls. These data provide a unique opportunity to investigate rare variants, which contribute to the majority of the variation in the genome. To improve power for discovering rare CRC susceptibility variants (<1% MAF), Roadmap Epigenomics data were used to construct biologically relevant testing sets of enhancers, promoters and exons for gene-based association testing across the genome. Since enhancers exert their effects by impacting expression of target genes, we defined enhancer-gene networks by linking enhancer(s) to target gene expression using Roadmap chromatin state maps and gene expression. Variants in linked enhancers from digestive and immune tissues were aggregated together with variants in the promoter and non-synonymous coding variants in the target gene. We tested 9,884 variant sets for association with CRC risk using the Mixed effects Score Test (MiST). Our most significant findings are for acyl-Coenzyme A dehydrogenase, C-2 to C-3 short chain precursor-ACADS (p = 1×10−4), AlkB homologs, including AlkB homolog 1-ALKBH1 (p = 2×10−4), and SRA stem-loop interacting RNA binding protein-SLIRP (p = 2×10−4). We will replicate these findings within the Colorectal Cancer Transdisciplinary Study (CORECT), as well as additional samples currently genotyped in CCFR and GECCO (over 25,000 CRC cases and controls). Although the top findings are statistically non-significant in this initial dataset, each of these genes linked to molecular pathways implicated in CRC carcinogenesis (fatty acid metabolism, DNA/RNA repair, and Nuclear Receptor signaling pathway, which interacts with the Wnt, beta-catenin pathways to result in a diverse array of cellular effects including altered cellular adhesion, tissue morphogenesis, and oncogenesis). Our current findings suggest that although functional insight can improve power for novel discovery, even larger sample sizes and/or pathway-based analyses are necessary to understand the role of rare variants in CRC carcinogenesis.

Citation Format: Stephanie A. Bien, Tabitha A. Harrison, Paul L. Auer, Flora Qu, Jeroen Huyghe, Barbara Banbury, Peyton Greenside, Goncalo R. Abecasis, Sonja I. Berndt, Stephane Bézieau, Hermann Brenner, Graham Casey, Andrew T. Chan, Jenny Chang-Claude, Sai Chen, Joshua D. Smith, Loic Le Marchand, Christopher Carlson, Polly A. Newcomb, Christian Fuchsberger, Marty L. Slattery, Hyun M. Kang, Emily White, John Potter, Steven J. Gallinger, Michael Hoffmeister, Stephen B. Gruber, Deborah A. Nickerson, Ulrike Peters, Anshul Kundaje, Li Hsu. Using functional data from Roadmap Epigenomics to inform analysis of rare variants linked to gene expression in a large colorectal cancer study. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 4489.

Large-Scale Exome-wide Association Analysis Identifies Loci for White Blood Cell Traits and Pleiotropy with Immune-Mediated Diseases

White blood cells play diverse roles in innate and adaptive immunity. Genetic association analyses of phenotypic variation in circulating white blood cell (WBC) counts from large samples of otherwise healthy individuals can provide insights into genes and biologic pathways involved in production, differentiation, or clearance of particular WBC lineages (myeloid, lymphoid) and also potentially inform the genetic basis of autoimmune, allergic, and blood diseases. We performed an exome array-based meta-analysis of total WBC and subtype counts (neutrophils, monocytes, lymphocytes, basophils, and eosinophils) in a multi-ancestry discovery and replication sample of ∼157,622 individuals from 25 studies. We identified 16 common variants (8 of which were coding variants) associated with one or more WBC traits, the majority of which are pleiotropically associated with autoimmune diseases. Based on functional annotation, these loci included genes encoding surface markers of myeloid, lymphoid, or hematopoietic stem cell differentiation (CD69, CD33, CD87), transcription factors regulating lineage specification during hematopoiesis (ASXL1, IRF8, IKZF1, JMJD1C, ETS2-PSMG1), and molecules involved in neutrophil clearance/apoptosis (C10orf54, LTA), adhesion (TNXB), or centrosome and microtubule structure/function (KIF9, TUBD1). Together with recent reports of somatic ASXL1 mutations among individuals with idiopathic cytopenias or clonal hematopoiesis of undetermined significance, the identification of a common regulatory 3' UTR variant of ASXL1 suggests that both germline and somatic ASXL1 mutations contribute to lower blood counts in otherwise asymptomatic individuals. These association results shed light on genetic mechanisms that regulate circulating WBC counts and suggest a prominent shared genetic architecture with inflammatory and autoimmune diseases.

Platelet-Related Variants Identified by Exomechip Meta-analysis in 157,293 Individuals

Platelet production, maintenance, and clearance are tightly controlled processes indicative of platelets' important roles in hemostasis and thrombosis. Platelets are common targets for primary and secondary prevention of several conditions. They are monitored clinically by complete blood counts, specifically with measurements of platelet count (PLT) and mean platelet volume (MPV). Identifying genetic effects on PLT and MPV can provide mechanistic insights into platelet biology and their role in disease. Therefore, we formed the Blood Cell Consortium (BCX) to perform a large-scale meta-analysis of Exomechip association results for PLT and MPV in 157,293 and 57,617 individuals, respectively. Using the low-frequency/rare coding variant-enriched Exomechip genotyping array, we sought to identify genetic variants associated with PLT and MPV. In addition to confirming 47 known PLT and 20 known MPV associations, we identified 32 PLT and 18 MPV associations not previously observed in the literature across the allele frequency spectrum, including rare large effect (FCER1A), low-frequency (IQGAP2, MAP1A, LY75), and common (ZMIZ2, SMG6, PEAR1, ARFGAP3/PACSIN2) variants. Several variants associated with PLT/MPV (PEAR1, MRVI1, PTGES3) were also associated with platelet reactivity. In concurrent BCX analyses, there was overlap of platelet-associated variants with red (MAP1A, TMPRSS6, ZMIZ2) and white (PEAR1, ZMIZ2, LY75) blood cell traits, suggesting common regulatory pathways with shared genetic architecture among these hematopoietic lineages. Our large-scale Exomechip analyses identified previously undocumented associations with platelet traits and further indicate that several complex quantitative hematological, lipid, and cardiovascular traits share genetic factors.

Exome Genotyping Identifies Pleiotropic Variants Associated with Red Blood Cell Traits

Red blood cell (RBC) traits are important heritable clinical biomarkers and modifiers of disease severity. To identify coding genetic variants associated with these traits, we conducted meta-analyses of seven RBC phenotypes in 130,273 multi-ethnic individuals from studies genotyped on an exome array. After conditional analyses and replication in 27,480 independent individuals, we identified 16 new RBC variants. We found low-frequency missense variants in MAP1A (rs55707100, minor allele frequency [MAF] = 3.3%, p = 2 × 10(-10) for hemoglobin [HGB]) and HNF4A (rs1800961, MAF = 2.4%, p < 3 × 10(-8) for hematocrit [HCT] and HGB). In African Americans, we identified a nonsense variant in CD36 associated with higher RBC distribution width (rs3211938, MAF = 8.7%, p = 7 × 10(-11)) and showed that it is associated with lower CD36 expression and strong allelic imbalance in ex vivo differentiated human erythroblasts. We also identified a rare missense variant in ALAS2 (rs201062903, MAF = 0.2%) associated with lower mean corpuscular volume and mean corpuscular hemoglobin (p < 8 × 10(-9)). Mendelian mutations in ALAS2 are a cause of sideroblastic anemia and erythropoietic protoporphyria. Gene-based testing highlighted three rare missense variants in PKLR, a gene mutated in Mendelian non-spherocytic hemolytic anemia, associated with HGB and HCT (SKAT p < 8 × 10(-7)). These rare, low-frequency, and common RBC variants showed pleiotropy, being also associated with platelet, white blood cell, and lipid traits. Our association results and functional annotation suggest the involvement of new genes in human erythropoiesis. We also confirm that rare and low-frequency variants play a role in the architecture of complex human traits, although their phenotypic effect is generally smaller than originally anticipated.

Fine-Mapping of Common Genetic Variants Associated with Colorectal Tumor Risk Identified Potential Functional Variants

Genome-wide association studies (GWAS) have identified many common single nucleotide polymorphisms (SNPs) associated with colorectal cancer risk. These SNPs may tag correlated variants with biological importance. Fine-mapping around GWAS loci can facilitate detection of functional candidates and additional independent risk variants. We analyzed 11,900 cases and 14,311 controls in the Genetics and Epidemiology of Colorectal Cancer Consortium and the Colon Cancer Family Registry. To fine-map genomic regions containing all known common risk variants, we imputed high-density genetic data from the 1000 Genomes Project. We tested single-variant associations with colorectal tumor risk for all variants spanning genomic regions 250-kb upstream or downstream of 31 GWAS-identified SNPs (index SNPs). We queried the University of California, Santa Cruz Genome Browser to examine evidence for biological function. Index SNPs did not show the strongest association signals with colorectal tumor risk in their respective genomic regions. Bioinformatics analysis of SNPs showing smaller P-values in each region revealed 21 functional candidates in 12 loci (5q31.1, 8q24, 11q13.4, 11q23, 12p13.32, 12q24.21, 14q22.2, 15q13, 18q21, 19q13.1, 20p12.3, and 20q13.33). We did not observe evidence of additional independent association signals in GWAS-identified regions. Our results support the utility of integrating data from comprehensive fine-mapping with expanding publicly available genomic databases to help clarify GWAS associations and identify functional candidates that warrant more onerous laboratory follow-up. Such efforts may aid the eventual discovery of disease-causing variant(s).

Clinical Sequencing Exploratory Research Consortium: Accelerating Evidence-Based Practice of Genomic Medicine

Despite rapid technical progress and demonstrable effectiveness for some types of diagnosis and therapy, much remains to be learned about clinical genome and exome sequencing (CGES) and its role within the practice of medicine. The Clinical Sequencing Exploratory Research (CSER) consortium includes 18 extramural research projects, one National Human Genome Research Institute (NHGRI) intramural project, and a coordinating center funded by the NHGRI and National Cancer Institute. The consortium is exploring analytic and clinical validity and utility, as well as the ethical, legal, and social implications of sequencing via multidisciplinary approaches; it has thus far recruited 5,577 participants across a spectrum of symptomatic and healthy children and adults by utilizing both germline and cancer sequencing. The CSER consortium is analyzing data and creating publically available procedures and tools related to participant preferences and consent, variant classification, disclosure and management of primary and secondary findings, health outcomes, and integration with electronic health records. Future research directions will refine measures of clinical utility of CGES in both germline and somatic testing, evaluate the use of CGES for screening in healthy individuals, explore the penetrance of pathogenic variants through extensive phenotyping, reduce discordances in public databases of genes and variants, examine social and ethnic disparities in the provision of genomics services, explore regulatory issues, and estimate the value and downstream costs of sequencing. The CSER consortium has established a shared community of research sites by using diverse approaches to pursue the evidence-based development of best practices in genomic medicine.

Association Between Absolute Neutrophil Count and Variation at TCIRG1: The NHLBI Exome Sequencing Project

Neutrophils are a key component of innate immunity. Individuals with low neutrophil count are susceptible to frequent infections. Linkage and association between congenital neutropenia and a single rare missense variant in TCIRG1 have been reported in a single family. Here, we report on nine rare missense variants at evolutionarily conserved sites in TCIRG1 that are associated with lower absolute neutrophil count (ANC; p = 0.005) in 1,058 participants from three cohorts: Atherosclerosis Risk in Communities (ARIC), Coronary Artery Risk Development in Young Adults (CARDIA), and Jackson Heart Study (JHS) of the NHLBI Grand Opportunity Exome Sequencing Project (GO ESP). These results validate the effects of TCIRG1 coding variation on ANC and suggest that this gene may be associated with a spectrum of mild to severe effects on ANC.

Identification of Rare Variants in ATP8B4 as a Risk Factor for Systemic Sclerosis by Whole-Exome Sequencing

OBJECTIVE

To determine the contribution of rare variants as genetic modifiers of the expressivity, penetrance, and severity of systemic sclerosis (SSc).

METHODS

We performed whole-exome sequencing of 78 European American patients with SSc, including 35 patients without pulmonary arterial hypertension (PAH) and 43 patients with PAH. Association testing of case-control probability for rare variants was performed using the unified sequence kernel association test with optimal kernel weighting and small sample adjustment by comparing all SSc patients with a reference population of 3,179 controls from the Exome Sequencing Project 5,500 exome data set. Replication genotyping was performed in an independent sample of 3,263 patients (415 patients with SSc and 2,848 controls). We conducted expression profiling of messenger RNA from 61 SSc patients (19 without PAH and 42 with PAH) and 41 corresponding controls.

RESULTS

The ATP8B4 gene was associated with a significant increase in the risk of SSc (P = 2.77 × 10(-7)). Among the 64 ATP8B4 variants tested, a single missense variant, c.1308C>G (F436L, rs55687265), provided the most compelling evidence of association (P = 9.35 × 10(-10), odds ratio [OR] 6.11), which was confirmed in the replication cohort (P = 0.012, OR 1.86) and meta-analysis (P = 1.92 × 10(-7), OR 2.5). Genes involved in E3 ubiquitin-protein ligase complex (ASB10) and cyclic nucleotide gated channelopathies (CNGB3) as well as HLA-DRB5 and HSPB2 (heat-shock protein 27) provided additional evidence of association (P < 10(-5)). Differential ATP8B4 expression was observed among the SSc patients compared to the controls (P = 0.0005).

CONCLUSION

ATP8B4 may represent a putative genetic risk factor for SSc and pulmonary vascular complications.

ALG1-CDG: Clinical and Molecular Characterization of 39 Unreported Patients

Congenital disorders of glycosylation (CDG) arise from pathogenic mutations in over 100 genes leading to impaired protein or lipid glycosylation. ALG1 encodes a β1,4 mannosyltransferase that catalyzes the addition of the first of nine mannose moieties to form a dolichol-lipid linked oligosaccharide intermediate required for proper N-linked glycosylation. ALG1 mutations cause a rare autosomal recessive disorder termed ALG1-CDG. To date 13 mutations in 18 patients from 14 families have been described with varying degrees of clinical severity. We identified and characterized 39 previously unreported cases of ALG1-CDG from 32 families and add 26 new mutations. Pathogenicity of each mutation was confirmed based on its inability to rescue impaired growth or hypoglycosylation of a standard biomarker in an alg1-deficient yeast strain. Using this approach we could not establish a rank order comparison of biomarker glycosylation and patient phenotype, but we identified mutations with a lethal outcome in the first two years of life. The recently identified protein-linked xeno-tetrasaccharide biomarker, NeuAc-Gal-GlcNAc2 , was seen in all 27 patients tested. Our study triples the number of known patients and expands the molecular and clinical correlates of this disorder.

Race-Specific Influence of CYP4F2 on Dose and Risk of Hemorrhage Among Warfarin Users

OBJECTIVE

The p.V433M in cytochrome P450 4F2 (rs2108622, CYP4F2*3) is associated with a higher warfarin dose and lower risk of hemorrhage among European Americans. We evaluate the influence of CYP4F2*3 on warfarin dose, time to target international normalized ratio (INR), and stable dose, proportion of time spent in target range (PTTR), as well as the risk of overanticoagulation and hemorrhage among European and African Americans.

DESIGN

CYP4F2*3 was genotyped in 1238 patients initiated on warfarin in a prospective inception cohort. Multivariable linear regression was used to assess warfarin dose and PTTR; proportional hazards analysis was performed to evaluate time to target INR and stable dose, overanticoagulation, and hemorrhage.

SETTING

Two outpatient anticoagulation clinics.

PARTICIPANTS

A total of 1238 anticoagulated patients.

OUTCOMES

Warfarin dose (mg/day), time to target INR and stable dose, PTTR, overanticoagulation (INR more than 4), and major hemorrhage.

RESULTS

Minor allele frequency for the CYP4F2*3 variant was 30.3% among European Americans and 8.4% among African Americans. CYP4F2*3 was associated with higher dose among European Americans but not African Americans. Compared to CYP4F2*1/*1, *1/*3 was associated with a statistically nonsignificant increase in dose (4.5%, p=0.22) and *3/*3 was associated with a statistically significant increase in dose (13.2%, p=0.02). CYP4F2 genotype did not influence time to target INR, time to stable dose, or PTTR in either race group. CYP4F2*3/*3 was associated with a 31% lower risk of over anticoagulation (p=0.06). Incidence of hemorrhage was lower among participants with CYP4F2 *3/*3 compared with *1/*3 or *1/*1 (incidence rate ratio = 0.45, 95% confidence interval 0.14-1.11, p=0.09). After controlling for covariates, CYP4F2 *3/*3 was associated with a 52% lower risk of hemorrhage, although this was not statistically significant (p=0.24).

CONCLUSION

Possession of CYP4F2*3 variant influences warfarin dose among European Americans but not African Americans. The CYP4F2-dose, CYP4F2-overanticoagulation, and CYP4F2-hemorrhage association follows a recessive pattern with possession of CYP4F2*3/*3 genotype likely demonstrating a protective effect. These findings need further confirmation.

Statin-induced expression change of INSIG1 in lymphoblastoid cell lines correlates with plasma triglyceride statin response in a sex-specific manner

Statins are widely prescribed to lower plasma low-density lipoprotein (LDL) cholesterol levels. They also modestly reduce plasma triglyceride (TG), an independent cardiovascular disease risk factor, in most people. The mechanism and inter-individual variability of TG statin response is poorly understood. We measured statin-induced gene expression changes in lymphoblastoid cell lines derived from 150 participants of a simvastatin clinical trial and identified 23 genes (false discovery rate, FDR=15%) with expression changes correlated with plasma TG response. The correlation of insulin-induced gene 1 (INSIG1) expression changes with TG response (rho=0.32, q=0.11) was driven by men (interaction P=0.0055). rs73161338 was associated with INSIG1 expression changes (P=5.4 × 10⁻⁵) and TG response in two statin clinical trials (P=0.0048), predominantly in men. A combined model including INSIG1 expression level and splicing changes accounted for 29.5% of plasma TG statin response variance in men (P=5.6 × 10⁻⁶). Our results suggest that INSIG1 variation may contribute to statin-induced changes in plasma TG in a sex-specific manner.

Disrupted nitric oxide signaling due to GUCY1A3 mutations increases risk for moyamoya disease, achalasia and hypertension

Moyamoya disease (MMD) is a progressive vasculopathy characterized by occlusion of the terminal portion of the internal carotid arteries and its branches, and the formation of compensatory moyamoya collateral vessels. Homozygous mutations in GUCY1A3 have been reported as a cause of MMD and achalasia. Probands (n = 96) from unrelated families underwent sequencing of GUCY1A3. Functional studies were performed to confirm the pathogenicity of identified GUCY1A3 variants. Two affected individuals from the unrelated families were found to have compound heterozygous mutations in GUCY1A3. MM041 was diagnosed with achalasia at 4 years of age, hypertension and MMD at 18 years of age. MM149 was diagnosed with MMD and hypertension at the age of 20 months. Both individuals carry one allele that is predicted to lead to haploinsufficiency and a second allele that is predicted to produce a mutated protein. Biochemical studies of one of these alleles, GUCY1A3 Cys517Tyr, showed that the mutant protein (a subunit of soluble guanylate cyclase) has a significantly blunted signaling response with exposure to nitric oxide (NO). GUCY1A3 missense and haploinsufficiency mutations disrupt NO signaling leading to MMD and hypertension, with or without achalasia.

Mitral regurgitation as a phenotypic manifestation of nonphotosensitive trichothiodystrophy due to a splice variant in MPLKIP

Background

Nonphotosensitive trichothiodystrophy (TTDN) is a rare autosomal recessive disorder of neuroectodermal origin. The condition is marked by hair abnormalities, intellectual impairment, nail dystrophies and susceptibility to infections but with no UV sensitivity.

Methods

We identified three consanguineous Pakistani families with varied TTDN features and used homozygosity mapping, linkage analysis, and Sanger and exome sequencing in order to identify pathogenic variants. Haplotype analysis was performed and haplotype age estimated. A splicing assay was used to validate the effect of the MPLKIP splice variant on expression.

Results

Affected individuals from all families exhibit several TTDN features along with a heart-specific feature, i.e. mitral regurgitation. Exome sequencing in the probands from families ED168 and ED241 identified a homozygous splice mutation c.339 + 1G > A within MPLKIP. The same splice variant co-segregates with TTDN in a third family ED210. The MPLKIP splice variant was not found in public databases, e.g. the Exome Aggregation Consortium, and in unrelated Pakistani controls. Functional analysis of the splice variant confirmed intron retention, which leads to protein truncation and loss of a phosphorylation site. Haplotype analysis identified a 585.1-kb haplotype which includes the MPLKIP variant, supporting the existence of a founder haplotype that is estimated to be 25,900 years old.

Conclusion

This study extends the allelic and phenotypic spectra of MPLKIP-related TTDN, to include a splice variant that causes cardiomyopathy as part of the TTDN phenotype.

Autosomal-Recessive Hearing Impairment Due to Rare Missense Variants within S1PR2

The sphingosine-1-phosphate receptors (S1PRs) are a well-studied class of transmembrane G protein-coupled sphingolipid receptors that mediate multiple cellular processes. However, S1PRs have not been previously reported to be involved in the genetic etiology of human traits. S1PR2 lies within the autosomal-recessive nonsyndromic hearing impairment (ARNSHI) locus DFNB68 on 19p13.2. From exome sequence data we identified two pathogenic S1PR2 variants, c.323G>C (p.Arg108Pro) and c.419A>G (p.Tyr140Cys). Each of these variants co-segregates with congenital profound hearing impairment in consanguineous Pakistani families with maximum LOD scores of 6.4 for family DEM4154 and 3.3 for family PKDF1400. Neither S1PR2 missense variant was reported among ∼120,000 chromosomes in the Exome Aggregation Consortium database, in 76 unrelated Pakistani exomes, or in 720 Pakistani control chromosomes. Both DNA variants affect highly conserved residues of S1PR2 and are predicted to be damaging by multiple bioinformatics tools. Molecular modeling predicts that these variants affect binding of sphingosine-1-phosphate (p.Arg108Pro) and G protein docking (p.Tyr140Cys). In the previously reported S1pr2(-/-) mice, stria vascularis abnormalities, organ of Corti degeneration, and profound hearing loss were observed. Additionally, hair cell defects were seen in both knockout mice and morphant zebrafish. Family PKDF1400 presents with ARNSHI, which is consistent with the lack of gross malformations in S1pr2(-/-) mice, whereas family DEM4154 has lower limb malformations in addition to hearing loss. Our findings suggest the possibility of developing therapies against hair cell damage (e.g., from ototoxic drugs) through targeted stimulation of S1PR2.

Pharmacogenetic allele nomenclature: International workgroup recommendations for test result reporting

This article provides nomenclature recommendations developed by an international workgroup to increase transparency and standardization of pharmacogenetic (PGx) result reporting. Presently, sequence variants identified by PGx tests are described using different nomenclature systems. In addition, PGx analysis may detect different sets of variants for each gene, which can affect interpretation of results. This practice has caused confusion and may thereby impede the adoption of clinical PGx testing. Standardization is critical to move PGx forward.

Pathogenic FBN1 variants in familial thoracic aortic aneurysms and dissections

Marfan syndrome (MFS) due to mutations in FBN1 is a known cause of thoracic aortic aneurysms and acute aortic dissections (TAAD) associated with pleiotropic manifestations. Genetic predisposition to TAAD can also be inherited in families in the absence of syndromic features, termed familial TAAD (FTAAD), and several causative genes have been identified to date. FBN1 mutations can also be identified in FTAAD families, but the frequency of these mutations has not been established. We performed exome sequencing of 183 FTAAD families and identified pathogenic FBN1 variants in five (2.7%) of these families. We also identified eight additional FBN1 rare variants that could not be unequivocally classified as disease-causing in six families. FBN1 sequencing should be considered in individuals with FTAAD even without significant systemic features of MFS.

LOX Mutations Predispose to Thoracic Aortic Aneurysms and Dissections

RATIONALE

Mutations in several genes have been identified that are responsible for 25% of families with familial thoracic aortic aneurysms and dissections. However, the causative gene remains unknown in 75% of families.

OBJECTIVES

To identify the causative mutation in families with autosomal dominant inheritance of thoracic aortic aneurysms and dissections.

METHODS AND RESULTS

Exome sequencing was used to identify the mutation responsible for a large family with thoracic aortic aneurysms and dissections. A heterozygous rare variant, c.839G>T (p.Ser280Arg), was identified in LOX, encoding a lysyl oxidase, that segregated with disease in the family. Sanger and exome sequencing was used to investigate mutations in LOX in an additional 410 probands from unrelated families. Additional LOX rare variants that segregated with disease in families were identified, including c.125G>A (p.Trp42*), c.604G>T (p.Gly202*), c.743C>T (p.Thr248Ile), c.800A>C (p.Gln267Pro), and c.1044T>A (p.Ser348Arg). The altered amino acids cause haploinsufficiency for LOX or are located at a highly conserved LOX catalytic domain, which is relatively invariant in the population. Expression of the LOX variants p.Ser280Arg and p.Ser348Arg resulted in significantly lower lysyl oxidase activity when compared with the wild-type protein. Individuals with LOX variants had fusiform enlargement of the root and ascending thoracic aorta, leading to ascending aortic dissections.

CONCLUSIONS

These data, along with previous studies showing that the deficiency of LOX in mice or inhibition of lysyl oxidases in turkeys and rats causes aortic dissections, support the conclusion that rare genetic variants in LOX predispose to thoracic aortic disease.

Association of Arrhythmia-Related Genetic Variants With Phenotypes Documented in Electronic Medical Records

IMPORTANCE

Large-scale DNA sequencing identifies incidental rare variants in established Mendelian disease genes, but the frequency of related clinical phenotypes in unselected patient populations is not well established. Phenotype data from electronic medical records (EMRs) may provide a resource to assess the clinical relevance of rare variants.

OBJECTIVE

To determine the clinical phenotypes from EMRs for individuals with variants designated as pathogenic by expert review in arrhythmia susceptibility genes.

DESIGN, SETTING, AND PARTICIPANTS

This prospective cohort study included 2022 individuals recruited for nonantiarrhythmic drug exposure phenotypes from October 5, 2012, to September 30, 2013, for the Electronic Medical Records and Genomics Network Pharmacogenomics project from 7 US academic medical centers. Variants in SCN5A and KCNH2, disease genes for long QT and Brugada syndromes, were assessed for potential pathogenicity by 3 laboratories with ion channel expertise and by comparison with the ClinVar database. Relevant phenotypes were determined from EMRs, with data available from 2002 (or earlier for some sites) through September 10, 2014.

EXPOSURES

One or more variants designated as pathogenic in SCN5A or KCNH2.

MAIN OUTCOMES AND MEASURES

Arrhythmia or electrocardiographic (ECG) phenotypes defined by International Classification of Diseases, Ninth Revision (ICD-9) codes, ECG data, and manual EMR review.

RESULTS

Among 2022 study participants (median age, 61 years [interquartile range, 56-65 years]; 1118 [55%] female; 1491 [74%] white), a total of 122 rare (minor allele frequency <0.5%) nonsynonymous and splice-site variants in 2 arrhythmia susceptibility genes were identified in 223 individuals (11% of the study cohort). Forty-two variants in 63 participants were designated potentially pathogenic by at least 1 laboratory or ClinVar, with low concordance across laboratories (Cohen κ = 0.26). An ICD-9 code for arrhythmia was found in 11 of 63 (17%) variant carriers vs 264 of 1959 (13%) of those without variants (difference, +4%; 95% CI, -5% to +13%; P = .35). In the 1270 (63%) with ECGs, corrected QT intervals were not different in variant carriers vs those without (median, 429 vs 439 milliseconds; difference, -10 milliseconds; 95% CI, -16 to +3 milliseconds; P = .17). After manual review, 22 of 63 participants (35%) with designated variants had any ECG or arrhythmia phenotype, and only 2 had corrected QT interval longer than 500 milliseconds.

CONCLUSIONS AND RELEVANCE

Among laboratories experienced in genetic testing for cardiac arrhythmia disorders, there was low concordance in designating SCN5A and KCNH2 variants as pathogenic. In an unselected population, the putatively pathogenic genetic variants were not associated with an abnormal phenotype. These findings raise questions about the implications of notifying patients of incidental genetic findings.

Novel mutations in the genes TGM1 and ALOXE3 underlying autosomal recessive congenital ichthyosis

BACKGROUND

Ichthyoses are clinically characterized by scaling or hyperkeratosis of the skin or both. It can be an isolated condition limited to the skin or appear secondarily with involvement of other cutaneous or systemic abnormalities.

METHODS

The present study investigated clinical and molecular characterization of three consanguineous families (A, B, C) segregating two different forms of autosomal recessive congenital ichthyosis (ARCI). Linkage in three consanguineous families (A, B, C) segregating two different forms of ARCI was searched by typing microsatellite and single nucleotide polymorphism marker analysis. Sequencing of the two genes TGM1 and ALOXE3 was performed by the dideoxy chain termination method.

RESULTS

Genome-wide linkage analysis established linkage in family A to TGM1 gene on chromosome 14q11 and in families B and C to ALOXE3 gene on chromosome 17p13. Subsequently, sequencing of these genes using samples from affected family members led to the identification of three novel mutations: a missense variant p.Trp455Arg in TGM1 (family A); a nonsense variant p.Arg140* in ALOXE3 (family B); and a complex rearrangement in ALOXE3 (family C).

CONCLUSION

The present study further extends the spectrum of mutations in the two genes involved in causing ARCI. Characterizing the clinical spectrum resulting from mutations in the TGM1 and ALOXE3 genes will improve diagnosis and may direct clinical care of the family members.

Burden of potentially pathologic copy number variants is higher in children with isolated congenital heart disease and significantly impairs covariate-adjusted transplant-free survival

OBJECTIVES

Copy number variants (CNVs) are duplications or deletions of genomic regions. Large CNVs are potentially pathogenic and are overrepresented in children with congenital heart disease (CHD). We sought to determine the frequency of large CNVs in children with isolated CHD, and to evaluate the relationship of these potentially pathogenic CNVs with transplant-free survival.

METHODS

These cases are derived from a prospective cohort of patients with nonsyndromic CHD (n = 422) identified before first surgery. Healthy pediatric controls (n = 500) were obtained from the electronic Medical Records and Genetic Epidemiology Network, and CNV frequency was contrasted for CHD cases and controls. CNVs were determined algorithmically; subsequently screened for >95% overlap between 2 methods, size (>300 kb), quality score, overlap with a gene, and novelty (absent from databases of known, benign CNVs); and separately validated by quantitative polymerase chain reaction. Survival likelihoods for cases were calculated using Cox proportional hazards modeling to evaluate the joint effect of CNV burden and known confounders on transplant-free survival.

RESULTS

Children with nonsyndromic CHD had a higher burden of potentially pathogenic CNVs compared with pediatric controls (12.1% vs 5.0%; P = .00016). Presence of a CNV was associated with significantly decreased transplant-free survival after surgery (hazard ratio, 3.42; 95% confidence interval, 1.66-7.09; P = .00090) with confounder adjustment.

CONCLUSIONS

We confirm that children with isolated CHD have a greater burden of rare/large CNVs. We report a novel finding that these CNVs are associated with an adjusted 2.55-fold increased risk of death or transplant. These data suggest that CNV burden is an important modifier of survival after surgery for CHD.

Mandibular hypoplasia, deafness, progeroid features and lipodystrophy (MDPL) syndrome in the context of inherited lipodystrophies

BACKGROUND

Lipodystrophies are a large heterogeneous group of genetic or acquired disorders characterized by generalized or partial fat loss, usually associated with metabolic complications such as diabetes mellitus, hypertriglyceridemia and hepatic steatosis. Many efforts have been made in the last years in identifying the genetic etiologies of several lipodystrophy forms, although some remain to be elucidated.

METHODS

We report here the clinical description of a woman with a rare severe lipodystrophic and progeroid syndrome associated with hypertriglyceridemia and diabetes whose genetic bases have been clarified through whole-exome sequencing (WES) analysis.

RESULTS

This article reports the 5th MDPL (Mandibular hypoplasia, deafness, progeroid features, and lipodystrophy syndrome) patient with the same de novo p.S605del mutation in POLD1. We provided further genetic evidence that this is a disease-causing mutation along with a plausible molecular mechanism responsible for this recurring event. Moreover we overviewed the current classification of the inherited forms of lipodystrophy, along with their underlying molecular basis.

CONCLUSIONS

Progress in the identification of lipodystrophy genes will help in better understanding the role of the pathways involved in the complex physiology of fat. This will lead to new targets towards develop innovative therapeutic strategies for treating the disorder and its metabolic complications, as well as more common forms of adipose tissue redistribution as observed in the metabolic syndrome and type 2 diabetes.

PRIMUS: improving pedigree reconstruction using mitochondrial and Y haplotypes

PRIMUS is a pedigree reconstruction algorithm that uses estimates of genome-wide identity by descent to reconstruct pedigrees consistent with observed genetic data. However, when genetic data for individuals within a pedigree are missing, often multiple pedigrees can be reconstructed that fit the data. We report a major expansion of PRIMUS that uses mitochondrial (mtDNA) and non-recombining Y chromosome (NRY) haplotypes to eliminate many pedigree structures that are inconsistent with the genetic data. We demonstrate that discordances in mtDNA and NRY haplotypes substantially reduce the number of potential pedigrees, and often lead to the identification of the correct pedigree.

AVAILABILITY AND IMPLEMENTATION

We have implemented PRIMUS updates in PERL and it is available at primus.gs.washington.edu.

Rare and Coding Region Genetic Variants Associated With Risk of Ischemic Stroke: The NHLBI Exome Sequence Project

IMPORTANCE

Stroke is the second leading cause of death and the third leading cause of years of life lost. Genetic factors contribute to stroke prevalence, and candidate gene and genome-wide association studies (GWAS) have identified variants associated with ischemic stroke risk. These variants often have small effects without obvious biological significance. Exome sequencing may discover predicted protein-altering variants with a potentially large effect on ischemic stroke risk.

OBJECTIVE

To investigate the contribution of rare and common genetic variants to ischemic stroke risk by targeting the protein-coding regions of the human genome.

DESIGN, SETTING, AND PARTICIPANTS

The National Heart, Lung, and Blood Institute (NHLBI) Exome Sequencing Project (ESP) analyzed approximately 6000 participants from numerous cohorts of European and African ancestry. For discovery, 365 cases of ischemic stroke (small-vessel and large-vessel subtypes) and 809 European ancestry controls were sequenced; for replication, 47 affected sibpairs concordant for stroke subtype and an African American case-control series were sequenced, with 1672 cases and 4509 European ancestry controls genotyped. The ESP's exome sequencing and genotyping started on January 1, 2010, and continued through June 30, 2012. Analyses were conducted on the full data set between July 12, 2012, and July 13, 2013.

MAIN OUTCOMES AND MEASURES

Discovery of new variants or genes contributing to ischemic stroke risk and subtype (primary analysis) and determination of support for protein-coding variants contributing to risk in previously published candidate genes (secondary analysis).

RESULTS

We identified 2 novel genes associated with an increased risk of ischemic stroke: a protein-coding variant in PDE4DIP (rs1778155; odds ratio, 2.15; P = 2.63 × 10(-8)) with an intracellular signal transduction mechanism and in ACOT4 (rs35724886; odds ratio, 2.04; P = 1.24 × 10(-7)) with a fatty acid metabolism; confirmation of PDE4DIP was observed in affected sibpair families with large-vessel stroke subtype and in African Americans. Replication of protein-coding variants in candidate genes was observed for 2 previously reported GWAS associations: ZFHX3 (cardioembolic stroke) and ABCA1 (large-vessel stroke).

CONCLUSIONS AND RELEVANCE

Exome sequencing discovered 2 novel genes and mechanisms, PDE4DIP and ACOT4, associated with increased risk for ischemic stroke. In addition, ZFHX3 and ABCA1 were discovered to have protein-coding variants associated with ischemic stroke. These results suggest that genetic variation in novel pathways contributes to ischemic stroke risk and serves as a target for prediction, prevention, and therapy.

Expanding the Molecular and Clinical Phenotype of SSR4-CDG

Congenital disorders of glycosylation (CDG) are a group of mostly autosomal recessive disorders primarily characterized by neurological abnormalities. Recently, we described a single CDG patient with a de novo mutation in the X-linked gene, Signal Sequence Receptor 4 (SSR4). We performed whole-exome sequencing to identify causal variants in several affected individuals who had either an undifferentiated neurological disorder or unsolved CDG of unknown etiology based on abnormal transferrin glycosylation. We now report eight affected males with either de novo (4) or inherited (4) loss of function mutations in SSR4. Western blot analysis revealed that the mutations caused a complete loss of SSR4 protein. In nearly all cases, the abnormal glycosylation of serum transferrin was only slightly above the accepted normal cutoff range.

Rare A2ML1 variants confer susceptibility to otitis media

A duplication variant within the middle ear-specific gene A2ML1 cosegregates with otitis media in an indigenous Filipino pedigree (LOD score = 7.5 at reduced penetrance) and lies within a founder haplotype that is also shared by 3 otitis-prone European-American and Hispanic-American children but is absent in non-otitis-prone children and >62,000 next-generation sequences. We identified seven additional A2ML1 variants in six otitis-prone children. Collectively, our studies support a role for A2ML1 in the pathophysiology of otitis media.