MED resulting from recessively inherited mutations in the gene encoding calcium-activated nucleotidase CANT1

Multiple Epiphyseal Dysplasia (MED) is a relatively mild skeletal dysplasia characterized by mild short stature, joint pain, and early-onset osteoarthropathy. Dominantly inherited mutations in COMP, MATN3, COL9A1, COL9A2, and COL9A3, and recessively inherited mutations in SLC26A2, account for the molecular basis of disease in about 80-85% of the cases. In two families with recurrent MED of an unknown molecular basis, we used exome sequencing and candidate gene analysis to identify homozygosity for recessively inherited missense mutations in CANT1, which encodes calcium-activated nucleotidase 1. The MED phenotype is thus allelic to the more severe Desbuquois dysplasia phenotype and the results identify CANT1 as a second locus for recessively inherited MED.

Encephalopathy caused by novel mutations in the CMP-sialic acid transporter, SLC35A1

Transport of activated nucleotide-sugars into the Golgi is critical for proper glycosylation and mutations in these transporters cause a group of rare genetic disorders termed congenital disorders of glycosylation. We performed exome sequencing on an individual with a profound neurological presentation and identified rare compound heterozygous mutations, p.Thr156Arg and p.Glu196Lys, in the CMP-sialic acid transporter, SLC35A1. Patient primary fibroblasts and serum showed a considerable decrease in the amount of N- and O-glycans terminating in sialic acid. Direct measurement of CMP-sialic acid transport into the Golgi showed a substantial decrease in overall rate of transport. Here we report the identification of the third patient with CMP-sialic acid transporter deficiency, who presented with severe neurological phenotype, but without hematological abnormalities.

Monoallelic and biallelic CREB3L1 variant causes mild and severe osteogenesis imperfecta, respectively

Purpose

Osteogenesis imperfecta (OI) is a heritable skeletal dysplasia. Dominant pathogenic variants in COL1A1 and COL1A2 explain the majority of OI cases. At least fifteen additional genes have been identified, but still do not account for all OI phenotypes that present. We sought the genetic cause of mild and lethal OI phenotypes in an unsolved family.

Methods

We performed exome sequencing on seven members of the family, both affected and unaffected.

Results

We identified a variant in Cyclic AMP Responsive Element Binding Protein 3-Like 1 (CREB3L1) in a consanguineous family. The variant caused a prenatal/perinatal lethal OI in homozygotes, similar to that seen in OI type II as a result of mutations in type I collagen genes, and a mild phenotype (fractures, blue sclerae) in multiple heterozygous family members. CREB3L1 encodes Old Astrocyte Specifically-Induced Substance (OASIS), an ER stress transducer. The variant disrupts a DNA-binding site and prevents OASIS from acting on its transcriptional targets including SEC24D, which encodes a component of the coat protein II (COPII) complex.

Conclusion

This report confirms that CREB3L1 is an OI-related gene and suggests the pathogenic mechanism of CREB3L1-associated OI involves the altered regulation of proteins involved in cellular secretion.

Influence of common and rare genetic variation on warfarin dose among African-Americans and European-Americans using the exome array

AIM

We conducted a genome-wide association study using the Illumina Exome Array to identify coding SNPs that may explain additional warfarin dose variability.

PATIENTS & METHODS

Analysis was performed after adjustment for clinical variables and genetic factors known to influence warfarin dose among 1680 warfarin users (838 European-Americans and 842 African-Americans). Replication was performed in an independent sample.

RESULTS

We confirmed the influence of known genetic variants on warfarin dose variability. Our study is the first to show the association between rs12772169 and warfarin dose in African-Americans. In addition, genes COX15 and FGF5 showed significant association in European-Americans.

CONCLUSION

We identified some novel genes/SNPs that underpin warfarin dose response. Further replication is needed to confirm our findings.

Mutations in ARMC9, which Encodes a Basal Body Protein, Cause Joubert Syndrome in Humans and Ciliopathy Phenotypes in Zebrafish

Joubert syndrome (JS) is a recessive neurodevelopmental disorder characterized by hypotonia, ataxia, abnormal eye movements, and variable cognitive impairment. It is defined by a distinctive brain malformation known as the "molar tooth sign" on axial MRI. Subsets of affected individuals have malformations such as coloboma, polydactyly, and encephalocele, as well as progressive retinal dystrophy, fibrocystic kidney disease, and liver fibrosis. More than 35 genes have been associated with JS, but in a subset of families the genetic cause remains unknown. All of the gene products localize in and around the primary cilium, making JS a canonical ciliopathy. Ciliopathies are unified by their overlapping clinical features and underlying mechanisms involving ciliary dysfunction. In this work, we identify biallelic rare, predicted-deleterious ARMC9 variants (stop-gain, missense, splice-site, and single-exon deletion) in 11 individuals with JS from 8 families, accounting for approximately 1% of the disorder. The associated phenotypes range from isolated neurological involvement to JS with retinal dystrophy, additional brain abnormalities (e.g., heterotopia, Dandy-Walker malformation), pituitary insufficiency, and/or synpolydactyly. We show that ARMC9 localizes to the basal body of the cilium and is upregulated during ciliogenesis. Typical ciliopathy phenotypes (curved body shape, retinal dystrophy, coloboma, and decreased cilia) in a CRISPR/Cas9-engineered zebrafish mutant model provide additional support for ARMC9 as a ciliopathy-associated gene. Identifying ARMC9 mutations as a cause of JS takes us one step closer to a full genetic understanding of this important disorder and enables future functional work to define the central biological mechanisms underlying JS and other ciliopathies.

Genome sequencing and carrier testing: decisions on categorization and whether to disclose results of carrier testing

PURPOSE

We investigated the use of genome sequencing for preconception carrier testing. Genome sequencing could identify one or more of thousands of X-linked or autosomal recessive conditions that could be disclosed during preconception or prenatal counseling. Therefore, a framework that helps both clinicians and patients understand the possible range of findings is needed to respect patient preferences by ensuring that information about only the desired types of genetic conditions are provided to a given patient.

METHODS

We categorized gene-condition pairs into groups using a previously developed taxonomy of genetic conditions. Patients could elect to receive results from these categories. A Return of Results Committee (RORC) developed inclusion and exclusion criteria for each category.

RESULTS

To date, the RORC has categorized 728 gene-condition pairs: 177 are categorized as life span-limiting, 406 are categorized as serious, 93 are categorized as mild, 41 are categorized as unpredictable, and 11 are categorized as adult-onset. An additional 64 gene-condition pairs were excluded from reporting to patients or put on a watch list, generally because evidence that a gene and condition were associated was limited.

CONCLUSION

Categorization of gene-condition pairs using our taxonomy simplifies communication regarding patient preferences for carrier information from a genomic test.Genet Med advance online publication 12 January 2017.

Abstract 1300: Genetic predictors of gene expression associated with risk of colorectal cancer

To date, genome-wide association studies (GWAS) have reported common variants in over 50 loci

with weak to moderate effects on CRC risk. These genetic factors in aggregate explain only a small

fraction of familial risk of CRC. To aid in the discovery of novel CRC loci, we integrated large

transcriptome data, including those generated in the Genotype-Tissue Expression (GTEx) Project in

genetic association analyses of CRC. The computational method, PrediXcan, was used to predict

transcript levels in relevant tissues and perform gene-level association tests with CRC. Prediction

models were developed using whole blood transcriptomes (n=922) from the depression genes and

networks (DGN), as well as colon transcriptomes (transverse n=169 and sigmoid n=124) from

GTEx datasets, along with high-density genotyping data from the same subjects. Genetically

determined expression levels were tested for association with CRC in 12,186 cases and 14,718

controls from GECCO-CCFR and suggestive associations (false discovery rate = 0.2) were evaluated

in 7,481 cases and 17,796 controls from the Asia Colorectal Cancer Consortium (ACCC) and 22,974

cases and 14,392 controls from the Colorectal Transdisciplinary (CORECT) study. We attempted to

replicate novel associations for eight genes and found statistically significant associations with

CXCR1 (OR=1.21 (1.10-1.33), p-value=7.8x10-5) and CXCR2 (OR=1.24 (1.11-1.38), p-value=9.9x10-

5). We also recovered previous associations at six known GWAS loci, thereby providing additional

support for putative target genes. CXCR1 and CXCR2 are therapeutic targets for the anticancer agent

Reparixin, which is currently being investigated in a stage II clinical trial for triple negative breast

cancer. As such, these findings provide preliminary support for new molecular targets that could

potentially repurpose a putative cancer therapeutic. These findings highlight the utility integrating

transcriptome data for novel discovery and biological insight of risk loci.

Citation Format: Stephanie A. Bien, Xingyi Guo, Yu-Ru Su, Tabitha A. Harrison, Conghui Qu, Yingchang Lu, Jiron Long, Sai Chen, Andrew T. Chan, David V. Conti, Hyun M. Kang, Michael Hoffmeister, Thomas J. Hudson, Mark A. Jenkins, Loic Le Marchand, Polly A. Newcomb, Martha L. Slattery, Emily White, Goncalo R. Abeçasis, Stephen B. Gruber, Deborah A. Nickerson, Stephanie L. Schmit, Graham Casey, Li Hsu, Wei Zheng, Ulrike Peters, GECCO-CCFR-AAAC-CORECT. Genetic predictors of gene expression associated with risk of colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1300. doi:10.1158/1538-7445.AM2017-1300

Sequencing of sporadic Attention-Deficit Hyperactivity Disorder (ADHD) identifies novel and potentially pathogenic de novo variants and excludes overlap with genes associated with autism spectrum disorder

Attention-Deficit Hyperactivity Disorder (ADHD) has high heritability; however, studies of common variation account for <5% of ADHD variance. Using data from affected participants without a family history of ADHD, we sought to identify de novo variants that could account for sporadic ADHD. Considering a total of 128 families, two analyses were conducted in parallel: first, in 11 unaffected parent/affected proband trios (or quads with the addition of an unaffected sibling) we completed exome sequencing. Six de novo missense variants at highly conserved bases were identified and validated from four of the 11 families: the brain-expressed genes TBC1D9, DAGLA, QARS, CSMD2, TRPM2, and WDR83. Separately, in 117 unrelated probands with sporadic ADHD, we sequenced a panel of 26 genes implicated in intellectual disability (ID) and autism spectrum disorder (ASD) to evaluate whether variation in ASD/ID-associated genes were also present in participants with ADHD. Only one putative deleterious variant (Gln600STOP) in CHD1L was identified; this was found in a single proband. Notably, no other nonsense, splice, frameshift, or highly conserved missense variants in the 26 gene panel were identified and validated. These data suggest that de novo variant analysis in families with independently adjudicated sporadic ADHD diagnosis can identify novel genes implicated in ADHD pathogenesis. Moreover, that only one of the 128 cases (0.8%, 11 exome, and 117 MIP sequenced participants) had putative deleterious variants within our data in 26 genes related to ID and ASD suggests significant independence in the genetic pathogenesis of ADHD as compared to ASD and ID phenotypes. © 2017 Wiley Periodicals, Inc.

Presynaptic congenital myasthenic syndrome with a homozygous sequence variant in LAMA5 combines myopia, facial tics, and failure of neuromuscular transmission

Defects in genes encoding the isoforms of the laminin alpha subunit have been linked to various phenotypic manifestations, including brain malformations, muscular dystrophy, ocular defects, cardiomyopathy, and skin abnormalities. We report here a severe defect of neuromuscular transmission in a consanguineous patient with a homozygous variant in the laminin alpha-5 subunit gene (LAMA5). The variant c.8046C>T (p.Arg2659Trp) is rare and has a predicted deleterious effect. The affected individual, who also carries a rare homozygous sequence variant in LAMA1, had muscle weakness, myopia, and facial tics. Magnetic resonance imaging of brain showed mild volume loss and periventricular T2 prolongation. Repetitive nerve stimulation revealed 50% decrement of compound muscle action potential amplitudes and 250% facilitation immediately after exercise, Endplate studies identified a profound reduction of the endplate potential quantal content and endplates with normal postsynaptic folding that were denuded or partially occupied by small nerve terminals. Expression studies revealed that p.Arg2659Trp caused decreased binding of laminin alpha-5 to SV2A and impaired laminin-521 cell-adhesion and cell projection support in primary neuronal cultures. In summary, this report describing severe neuromuscular transmission failure in a patient with a LAMA5 mutation expands the list of phenotypes associated with defects in genes encoding alpha-laminins.

Mutations in IFT-A satellite core component genes IFT43 and IFT121 produce short rib polydactyly syndrome with distinctive campomelia

Background

Skeletal ciliopathies comprise a spectrum of ciliary malfunction disorders that have a profound effect on the skeleton. Most common among these disorders is short rib polydactyly syndrome (SRPS), a recessively inherited perinatal lethal condition characterized by a long narrow chest, markedly shortened long bones, polydactyly and, often, multi-organ system involvement. SRPS shows extensive locus heterogeneity with mutations in genes encoding proteins that participate in cilia formation and/or function.

Results

Herein we describe mutations in IFT43, a satellite member of the retrograde IFT-A complex, that produce a form of SRPS with unusual bending of the ribs and appendicular bones. These newly described IFT43 mutations disrupted cilia formation, produced abnormalities in cartilage growth plate architecture thus contributing to altered endochondral ossification. We further show that the IFT43 SRPS phenotype is similar to SRPS resulting from mutations in the gene encoding IFT121 (WDR35), a direct interactor with IFT43.

Conclusions

This study defines a new IFT43-associated phenotype, identifying an additional locus for SRPS. The data demonstrate that IFT43 is essential for ciliogenesis and that the mutations disrupted the orderly proliferation and differentiation of growth plate chondrocytes, resulting in a severe effect on endochondral ossification and mineralization. Phenotypic similarities with SRPS cases resulting from mutations in the gene encoding the IFT43 direct interacting protein IFT121 suggests that similar mechanisms may be disrupted by defects in these two IFT-A satellite interactors.

Electronic supplementary material

The online version of this article (doi:10.1186/s13630-017-0051-y) contains supplementary material, which is available to authorized users.

Activation of a cryptic splice site in the mitochondrial elongation factor GFM1 causes combined OXPHOS deficiency

We report the clinical, biochemical, and molecular findings in two brothers with encephalopathy and multi-systemic disease. Abnormal transferrin glycoforms were suggestive of a type I congenital disorder of glycosylation (CDG). While exome sequencing was negative for CDG related candidate genes, the testing revealed compound heterozygous mutations in the mitochondrial elongation factor G gene (GFM1). One of the mutations had been reported previously while the second, novel variant was found deep in intron 6, activating a cryptic splice site. Functional studies demonstrated decreased GFM1 protein levels, suggested disrupted assembly of mitochondrial complexes III and V and decreased activities of mitochondrial complexes I and IV, all indicating combined OXPHOS deficiency.

Loss-of-Function Mutations in YY1AP1 Lead to Grange Syndrome and a Fibromuscular Dysplasia-Like Vascular Disease

Fibromuscular dysplasia (FMD) is a heterogeneous group of non-atherosclerotic and non-inflammatory arterial diseases that primarily involves the renal and cerebrovascular arteries. Grange syndrome is an autosomal-recessive condition characterized by severe and early-onset vascular disease similar to FMD and variable penetrance of brachydactyly, syndactyly, bone fragility, and learning disabilities. Exome-sequencing analysis of DNA from three affected siblings with Grange syndrome identified compound heterozygous nonsense variants in YY1AP1, and homozygous nonsense or frameshift YY1AP1 variants were subsequently identified in additional unrelated probands with Grange syndrome. YY1AP1 encodes yin yang 1 (YY1)-associated protein 1 and is an activator of the YY1 transcription factor. We determined that YY1AP1 localizes to the nucleus and is a component of the INO80 chromatin remodeling complex, which is responsible for transcriptional regulation, DNA repair, and replication. Molecular studies revealed that loss of YY1AP1 in vascular smooth muscle cells leads to cell cycle arrest with decreased proliferation and increased levels of the cell cycle regulator p21/WAF/CDKN1A and disrupts TGF-β-driven differentiation of smooth muscle cells. Identification of YY1AP1 mutations as a cause of FMD indicates that this condition can result from underlying genetic variants that significantly alter the phenotype of vascular smooth muscle cells.

DUOX2 Gene Mutation Manifesting as Resistance to Thyrotropin Phenotype

Resistance to thyrotropin (TSH) (RTSH; defined by elevated TSH and a normal or hypoplastic thyroid gland) can be caused by mutations in genes encoding the TSH receptor and PAX8, and it has been linked to a locus on chromosome 15. In two nonconsanguineous families with nongoitrous euthyroid hyperthyrotropinemia, typical of the RTSH phenotype, exome analysis identified five rare DUOX2 gene variants (p.A649E, p.P1391A, p.R885L, p.G488R, and p.SF965-6SfsX29) found to be pathogenic. This form of nongoitrous dyshormonogenesis masquerades both clinically and biochemically as RTSH. Accordingly, mutations in DUOX2 should be added to those of SLC26A4 as causes of RTSH.

Meta-analysis of genome-wide association studies of HDL cholesterol response to statins

BACKGROUND

In addition to lowering low density lipoprotein cholesterol (LDL-C), statin therapy also raises high density lipoprotein cholesterol (HDL-C) levels. Inter-individual variation in HDL-C response to statins may be partially explained by genetic variation.

METHODS AND RESULTS

We performed a meta-analysis of genome-wide association studies (GWAS) to identify variants with an effect on statin-induced high density lipoprotein cholesterol (HDL-C) changes. The 123 most promising signals with p<1×10-4 from the 16 769 statin-treated participants in the first analysis stage were followed up in an independent group of 10 951 statin-treated individuals, providing a total sample size of 27 720 individuals. The only associations of genome-wide significance (p<5×10-8) were between minor alleles at the CETP locus and greater HDL-C response to statin treatment.

CONCLUSIONS

Based on results from this study that included a relatively large sample size, we suggest that CETP may be the only detectable locus with common genetic variants that influence HDL-C response to statins substantially in individuals of European descent. Although CETP is known to be associated with HDL-C, we provide evidence that this pharmacogenetic effect is independent of its association with baseline HDL-C levels.

Guidelines for Large-Scale Sequence-Based Complex Trait Association Studies: Lessons Learned from the NHLBI Exome Sequencing Project

Massively parallel whole-genome sequencing (WGS) data have ushered in a new era in human genetics. These data are now being used to understand the role of rare variants in complex traits and to advance the goals of precision medicine. The technological and computing advances that have enabled us to generate WGS data on thousands of individuals have also outpaced our ability to perform analyses in scientifically and statistically rigorous and thoughtful ways. The past several years have witnessed the application of whole-exome sequencing (WES) to complex traits and diseases. From our analysis of NHLBI Exome Sequencing Project (ESP) data, not only have a number of important disease and complex trait association findings emerged, but our collective experience offers some valuable lessons for WGS initiatives. These include caveats associated with generating automated pipelines for quality control and analysis of rare variants; the importance of studying minority populations; sample size requirements and efficient study designs for identifying rare-variant associations; and the significance of incidental findings in population-based genetic research. With the ESP as an example, we offer guidance and a framework on how to conduct a large-scale association study in the era of WGS.

denovo-db: a compendium of human de novo variants

Whole-exome and whole-genome sequencing have facilitated the large-scale discovery of de novo variants in human disease. To date, most de novo discovery through next-generation sequencing focused on congenital heart disease and neurodevelopmental disorders (NDDs). Currently, de novo variants are one of the most significant risk factors for NDDs with a substantial overlap of genes involved in more than one NDD. To facilitate better usage of published data, provide standardization of annotation, and improve accessibility, we created denovo-db (http://denovo-db.gs.washington.edu), a database for human de novo variants. As of July 2016, denovo-db contained 40 different studies and 32,991 de novo variants from 23,098 trios. Database features include basic variant information (chromosome location, change, type); detailed annotation at the transcript and protein levels; severity scores; frequency; validation status; and, most importantly, the phenotype of the individual with the variant. We included a feature on our browsable website to download any query result, including a downloadable file of the full database with additional variant details. denovo-db provides necessary information for researchers to compare their data to other individuals with the same phenotype and also to controls allowing for a better understanding of the biology of de novo variants and their contribution to disease.

Destabilization of the IFT-B cilia core complex due to mutations in IFT81 causes a Spectrum of Short-Rib Polydactyly Syndrome

Short-rib polydactyly syndromes (SRPS) and Asphyxiating thoracic dystrophy (ATD) or Jeune Syndrome are recessively inherited skeletal ciliopathies characterized by profound skeletal abnormalities and are frequently associated with polydactyly and multiorgan system involvement. SRPS are produced by mutations in genes that participate in the formation and function of primary cilia and usually result from disruption of retrograde intraflagellar (IFT) transport of the cilium. Herein we describe a new spectrum of SRPS caused by mutations in the gene IFT81, a key component of the IFT-B complex essential for anterograde transport. In mutant chondrocytes, the mutations led to low levels of IFT81 and mutant cells produced elongated cilia, had altered hedgehog signaling, had increased post-translation modification of tubulin, and showed evidence of destabilization of additional anterograde transport complex components. These findings demonstrate the importance of IFT81 in the skeleton, its role in the anterograde transport complex, and expand the number of loci associated with SRPS.

Somatic mosaicism for a lethal TRPV4 mutation results in non-lethal metatropic dysplasia

Dominant mutations in TRPV4, which encodes the Transient Receptor Potential Cation Channel Subfamily V Member 4 calcium channel, result in a series of musculoskeletal disorders that include a set of peripheral neuropathies and a broad phenotypic spectrum of skeletal dysplasias. The skeletal phenotypes range from brachyolmia, in which there is scoliosis with mild short stature, through perinatal lethal metatropic dysplasia. We describe a case with phenotypic findings consistent with metatropic dysplasia, but in whom no TRPV4 mutation was detected by Sanger sequence analysis. Exome sequence analysis identified a known lethal metatropic dysplasia mutation, TRPV4^L618P , which was present at lower frequency than would be expected for a heterozygous change. The affected individual was shown to be a somatic mosaic for the mutation, providing an explanation for the milder than expected phenotype. The data illustrate that high-throughput sequencing of genomic DNA can facilitate detection of mosaicism with higher sensitivity than Sanger sequence analysis and identify a new genetic mechanism for metatropic dysplasia. © 2016 Wiley Periodicals, Inc.

SRD5A3-CDG: Expanding the phenotype of a congenital disorder of glycosylation with emphasis on adult onset features

Increasing numbers of congenital disorders of glycosylation (CDG) have been reported recently resulting in an expansion of the phenotypes associated with this group of disorders. SRD5A3 codes for polyprenol reductase which converts polyprenol to dolichol. This is a major pathway for dolichol biosynthesis for N-glycosylation, O-mannosylation, C-mannosylation, and GPI anchor synthesis. We present the features of five individuals (three children and two adults) with mutations in SRD5A3 focusing on the variable eye and skin involvement. We compare that to 13 affected individuals from the literature including five adults allowing us to delineate the features that may develop over time with this disorder including kyphosis, retinitis pigmentosa, and cataracts. © 2016 Wiley Periodicals, Inc.

IFT52 mutations destabilize anterograde complex assembly, disrupt ciliogenesis and result in short rib polydactyly syndrome

The short-rib polydactyly syndromes (SRPS) encompass a radiographically and genetically heterogeneous group of skeletal ciliopathies that are characterized by a long narrow chest, short extremities, and variable occurrence of polydactyly. Radiographic abnormalities include undermineralization of the calvarium, shortened and bowed appendicular bones, trident shaped acetabula and polydactyly. In a case of SRPS we identified compound heterozygosity for mutations in IFT52, which encodes a component of the anterograde intraflagellar transport complex. The IFT52 mutant cells synthesized a significantly reduced amount of IFT52 protein, leading to reduced synthesis of IFT74, IFT81, IFT88 and ARL13B, other key anterograde complex members. Ciliogenesis was also disrupted in the mutant cells, with a 60% reduction in the presence of cilia on mutant cells and loss of cilia length regulation for the cells with cilia. These data demonstrate that IFT52 is essential for anterograde complex integrity and for the biosynthesis and maintenance of cilia. The data identify a new locus for SRPS and show that IFT52 mutations result in a ciliopathy with primary effects on the skeleton.

Association of Exome Sequences With Cardiovascular Traits Among Blacks in the Jackson Heart Study

BACKGROUND

The correlation of null alleles with human phenotypes can provide insight into gene function in humans. In individuals of African ancestry, we set out to identify null and damaging missense variants, and test these variants for association with a range of cardiovascular phenotypes.

METHODS AND RESULTS

We performed whole-exome sequencing in 3223 black individuals from the Jackson Heart Study and found a total of 729 666 variant sites with minor allele frequency <5%, including 17 263 null variants and 49 929 missense variants predicted to be damaging by in silico algorithms. We tested null and damaging missense variants within each gene for association with 36 cardiovascular traits. We found 3 associations that met our prespecified level of significance (α=1.1×10(-7)). Null and damaging missense variants in PCSK9 were associated with 36 mg/dL lower low-density lipoprotein cholesterol (P=3×10(-21)). Three individuals in their 50s with complete PCSK9 deficiency (each compound heterozygote for PCSK9 p.Y142X and p.C679X) were identified, with one having a coronary artery calcification score in the 83rd percentile despite a low-density lipoprotein cholesterol of 32 mg/dL. A damaging missense variant in HBQ1 (p.G52A) was associated with a 2 pg/cell lower mean corpuscular hemoglobin (P=9×10(-13)) and rare damaging missense variants in VPS13A with higher red blood cell distribution width (P=9.9×10(-8)).

CONCLUSIONS

A limited number of null/damaging alleles with a large effect on cardiovascular traits were detectable in ≈3000 black individuals.

Abstract 5230: Large scale whole genome sequencing with imputation into GWAS improves our understanding of the genetic architecture of colorectal cancer

Whole-genome sequencing (WGS) has started a new era in human genetics in which data can be used to more fully understand the role of genetic variation in common complex diseases, including the role of less frequent and rare variants and structural variation. To explore the impact of these variants on colorectal cancer risk we conducted the first large scale WGS study for colorectal cancer (CRC) including 1,961 CRC cases and 981 controls. These WGS data as well as those from the Haplotype Reference Consortium were imputed in 13,104 CRC cases and 15,521 controls with genome-wide association study (GWAS) data that are part of the Colorectal Cancer Family Registry (CCFR) and the Genetics and Epidemiology of Colorectal Cancer Consortium (GECCO). Focusing on rare and less frequent variants, insertions and deletions we observed potentially novel variants: a less frequent variant (MAF = 0.026) on chromosome 5 located in NREP/STARD4-AS1 (p = value 4E-08); and a novel rare multi-allelic variant (MAF = 0.003) on chromosome 9 near KLF9 and TRPM3 (p-value 2E-09; the other allele of this multi-allelic variant had a MAF of 0.0003 and p-value of 0.55). Furthermore, we observed an independent locus close to the known region 8q24 that was located upstream of GSDMC (MAF = 0.16, p-value 5E-08). Within the known region 8q23/EIF3H we identified several low frequency variants with similar MAF (0.0181 to 0.0204) including a 6bp deletion with p-values between 4E-08 and 1E-09 that were independent of the common variant signal in this region. In addition, we identified statistically significant (p&lt;5E-08) deletions, insertions, and an essential splice site within known GWAS loci that present interesting candidates for functional studies. We will follow up these findings in independent samples from the Colorectal Cancer Transdisciplinary Study (CORECT) and CCFR, as well as additional samples currently genotyped in GECCO. In conclusion, next generation sequencing combined with imputation in large GWAS data sets has the potential to identify novel low frequency and rare genetic variants, aid fine-mapping of known CRC susceptibility loci and point to interesting functional candidates.

Citation Format: Jeroen Huyghe, Sai Chen, Hyun M. Kang, Tabitha A. Harrison, Sonja I. Berndt, Stephane Bézieau, Hermann Brenner, Graham Casey, Andrew T. Chan, Jenny Chang-Claude, Gallinger J. Steven, Stephen B. Gruber, Andrea Gsur, Michael Hoffmeister, Thomas J. Hudson, Loic Le Marchand, Polly A. Newcomb, John D. Potter, Conghui Qu, Martha L. Slattery, Joshua D. Smith, Emily White, Li Hsu, Goncalo R. Abecasis, Deborah A. Nickerson, Ulrike Peters. Large scale whole genome sequencing with imputation into GWAS improves our understanding of the genetic architecture of colorectal cancer. [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 5230.

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 (&lt;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.

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.

Exome sequencing in HFE C282Y homozygous men with extreme phenotypes identifies a GNPAT variant associated with severe iron overload

To identify polymorphisms associated with variability of iron overload severity in HFE-associated hemochromatosis, we performed exome sequencing of DNA from 35 male HFE C282Y homozygotes with either markedly increased iron stores (n = 22; cases) or with normal or mildly increased iron stores (n = 13; controls). The 35 participants, residents of the United States, Canada, and Australia, reported no or light alcohol consumption. Sequencing data included 82,068 single-nucleotide variants, and 10,337 genes were tested for a difference between cases and controls. A variant in the GNPAT gene showed the most significant association with severe iron overload (P = 3 × 10(-6) ; P = 0.033 by the likelihood ratio test after correction for multiple comparisons). Sixteen of twenty-two participants with severe iron overload had glyceronephosphate O-acyltransferase (GNPAT) polymorphism p.D519G (rs11558492; 15 heterozygotes, one homozygote). No control participant had this polymorphism. To examine functional consequences of GNPAT deficiency, we performed small interfering RNA-based knockdown of GNPAT in the human liver-derived cell line, HepG2/C3A. This knockdown resulted in a >17-fold decrease in expression of the messenger RNA encoding the iron-regulatory hormone, hepcidin.

CONCLUSION

GNPAT p.D519G is associated with a high-iron phenotype in HFE C282Y homozygotes and may participate in hepcidin regulation.

A homozygous missense variant in type I keratin KRT25 causes autosomal recessive woolly hair

BACKGROUND

Woolly hair (WH) is a hair abnormality that is primarily characterised by tightly curled hair with abnormal growth.

METHODS

In two unrelated consanguineous Pakistani families with non-syndromic autosomal recessive (AR) WH, homozygosity mapping and linkage analysis identified a locus within 17q21.1-q22, which contains the type I keratin gene cluster. A DNA sample from an affected individual from each family underwent exome sequencing.

RESULTS

A homozygous missense variant c.950T>C (p.(Leu317Pro)) within KRT25 segregated with ARWH in both families, and has a combined maximum two-point LOD score of 7.9 at ϴ=0. The KRT25 variant is predicted to result in disruption of the second α-helical rod domain and the entire protein structure, thus possibly interfering with heterodimerisation of K25 with type II keratins within the inner root sheath (IRS) of the hair follicle and the medulla of the hair shaft.

CONCLUSIONS

Our findings implicate a novel gene involved in human hair abnormality, and are consistent with the curled, fragile hair found in mice with Krt25 mutations, and further support the role of IRS-specific type I keratins in hair follicle development and maintenance of hair texture.

Prospective participant selection and ranking to maximize actionable pharmacogenetic variants and discovery in the eMERGE Network

Background

In an effort to return actionable results from variant data to electronic health records (EHRs), participants in the Electronic Medical Records and Genomics (eMERGE) Network are being sequenced with the targeted Pharmacogenomics Research Network sequence platform (PGRNseq). This cost-effective, highly-scalable, and highly-accurate platform was created to explore rare variation in 84 key pharmacogenetic genes with strong drug phenotype associations.

Methods

To return Clinical Laboratory Improvement Amendments (CLIA) results to our participants at the Group Health Cooperative, we sequenced the DNA of 900 participants (61 % female) with non-CLIA biobanked samples. We then selected 450 of those to be re-consented, to redraw blood, and ultimately to validate CLIA variants in anticipation of returning the results to the participant and EHR. These 450 were selected using an algorithm we designed to harness data from self-reported race, diagnosis and procedure codes, medical notes, laboratory results, and variant-level bioinformatics to ensure selection of an informative sample. We annotated the multi-sample variant call format by a combination of SeattleSeq and SnpEff tools, with additional custom variables including evidence from ClinVar, OMIM, HGMD, and prior clinical associations.

Results

We focused our analyses on 27 actionable genes, largely driven by the Clinical Pharmacogenetics Implementation Consortium. We derived a ranking system based on the total number of coding variants per participant (75.2±14.7), and the number of coding variants with high or moderate impact (11.5±3.9). Notably, we identified 11 stop-gained (1 %) and 519 missense (20 %) variants out of a total of 1785 in these 27 genes. Finally, we prioritized variants to be returned to the EHR with prior clinical evidence of pathogenicity or annotated as stop-gain for the following genes: CACNA1S and RYR1 (malignant hyperthermia); SCN5A, KCNH2, and RYR2 (arrhythmia); and LDLR (high cholesterol).

Conclusions

The incorporation of genetics into the EHR for clinical decision support is a complex undertaking for many reasons including lack of prior consent for return of results, lack of biospecimens collected in a CLIA environment, and EHR integration. Our study design accounts for these hurdles and is an example of a pilot system that can be utilized before expanding to an entire health system.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-015-0181-z) contains supplementary material, which is available to authorized users.