Christian M. Shaffer

Inactivating mutations in NPC1L1 and protection from coronary heart disease

BACKGROUND Ezetimibe lowers plasma levels of low-density lipoprotein (LDL) cholesterol by inhibiting the activity of the Niemann-Pick C1-like 1 (NPC1L1) protein. However, whether such inhibition reduces the risk of coronary heart disease is not known. Human mutations that inactivate a gene encoding a drug target can mimic the action of an inhibitory drug and thus can be used to infer potential effects of that drug. METHODS We sequenced the exons of NPC1L1 in 7364 patients with coronary heart disease and in 14,728 controls without such disease who were of European, African, or South Asian ancestry. We identified carriers of inactivating mutations (nonsense, splice-site, or frameshift mutations). In addition, we genotyped a specific inactivating mutation (p.Arg406X) in 22,590 patients with coronary heart disease and in 68,412 controls. We tested the association between the presence of an inactivating mutation and both plasma lipid levels and the risk of coronary heart disease. RESULTS With sequencing, we identified 15 distinct NPC1L1 inactivating mutations; approximately 1 in every 650 persons was a heterozygous carrier for 1 of these mutations. Heterozygous carriers of NPC1L1 inactivating mutations had a mean LDL cholesterol level that was 12 mg per deciliter (0.31 mmol per liter) lower than that in noncarriers (P=0.04). Carrier status was associated with a relative reduction of 53% in the risk of coronary heart disease (odds ratio for carriers, 0.47; 95% confidence interval, 0.25 to 0.87; P=0.008). In total, only 11 of 29,954 patients with coronary heart disease had an inactivating mutation (carrier frequency, 0.04%) in contrast to 71 of 83,140 controls (carrier frequency, 0.09%). CONCLUSIONS Naturally occurring mutations that disrupt NPC1L1 function were found to be associated with reduced plasma LDL cholesterol levels and a reduced risk of coronary heart disease. (Funded by the National Institutes of Health and others.).

Calmodulin Mutations Associated With Recurrent Cardiac Arrest in Infants

Background— Life-threatening disorders of heart rhythm may arise during infancy and can result in the sudden and tragic death of a child. We performed exome sequencing on 2 unrelated infants presenting with recurrent cardiac arrest to discover a genetic cause. Methods and Results— We ascertained 2 unrelated infants (probands) with recurrent cardiac arrest and dramatically prolonged QTc interval who were both born to healthy parents. The 2 parent-child trios were investigated with the use of exome sequencing to search for de novo genetic variants. We then performed follow-up candidate gene screening on an independent cohort of 82 subjects with congenital long-QT syndrome without an identified genetic cause. Biochemical studies were performed to determine the functional consequences of mutations discovered in 2 genes encoding calmodulin. We discovered 3 heterozygous de novo mutations in either CALM1 or CALM2, 2 of the 3 human genes encoding calmodulin, in the 2 probands and in 2 additional subjects with recurrent cardiac arrest. All mutation carriers were infants who exhibited life-threatening ventricular arrhythmias combined variably with epilepsy and delayed neurodevelopment. Mutations altered residues in or adjacent to critical calcium binding loops in the calmodulin carboxyl-terminal domain. Recombinant mutant calmodulins exhibited several-fold reductions in calcium binding affinity. Conclusions— Human calmodulin mutations disrupt calcium ion binding to the protein and are associated with a life-threatening condition in early infancy. Defects in calmodulin function will disrupt important calcium signaling events in heart, affecting membrane ion channels, a plausible molecular mechanism for potentially deadly disturbances in heart rhythm during infancy.

A gene expression fingerprint of C. elegans embryonic motor neurons

BackgroundDifferential gene expression specifies the highly diverse cell types that constitute the nervous system. With its sequenced genome and simple, well-defined neuroanatomy, the nematode C. elegans is a useful model system in which to correlate gene expression with neuron identity. The UNC-4 transcription factor is expressed in thirteen embryonic motor neurons where it specifies axonal morphology and synaptic function. These cells can be marked with an unc-4::GFP reporter transgene. Here we describe a powerful strategy, Micro-Array Profiling of C. elegans cells (MAPCeL), and confirm that this approach provides a comprehensive gene expression profile of unc-4::GFP motor neurons in vivo.ResultsFluorescence Activated Cell Sorting (FACS) was used to isolate unc-4::GFP neurons from primary cultures of C. elegans embryonic cells. Microarray experiments detected 6,217 unique transcripts of which ~1,000 are enriched in unc-4::GFP neurons relative to the average nematode embryonic cell. The reliability of these data was validated by the detection of known cell-specific transcripts and by expression in UNC-4 motor neurons of GFP reporters derived from the enriched data set. In addition to genes involved in neurotransmitter packaging and release, the microarray data include transcripts for receptors to a remarkably wide variety of signaling molecules. The added presence of a robust array of G-protein pathway components is indicative of complex and highly integrated mechanisms for modulating motor neuron activity. Over half of the enriched genes (537) have human homologs, a finding that could reflect substantial overlap with the gene expression repertoire of mammalian motor neurons.ConclusionWe have described a microarray-based method, MAPCeL, for profiling gene expression in specific C. elegans motor neurons and provide evidence that this approach can reveal candidate genes for key roles in the differentiation and function of these cells. These methods can now be applied to generate a gene expression map of the C. elegans nervous system.

Biobanks and Electronic Medical Records: Enabling Cost-Effective Research

Linking of data from electronic medical records to biological specimens enables cost-effective and rapid genomic analyses. The use of electronic medical record data linked to biological specimens in health care settings is expected to enable cost-effective and rapid genomic analyses. Here, we present a model that highlights potential advantages for genomic discovery and describe the operational infrastructure that facilitated multiple simultaneous discovery efforts.

Novel rare variants in congenital cardiac arrhythmia genes are frequent in drug-induced torsades de pointes.

Marked prolongation of the QT interval and polymorphic ventricular tachycardia following medication (drug-induced long QT syndrome, diLQTS) is a severe adverse drug reaction (ADR) that phenocopies congenital long QT syndrome (cLQTS) and is one of the leading causes for drug withdrawal and relabeling. We evaluated the frequency of rare non-synonymous variants in genes contributing to the maintenance of heart rhythm in cases of diLQTS using targeted capture coupled to next-generation sequencing. Eleven of 31 diLQTS subjects (36%) carried a novel missense mutation in genes with known congenital arrhythmia associations or with a known cLQTS mutation. In the 26 Caucasian subjects, 23% carried a highly conserved rare variant predicted to be deleterious to protein function in these genes compared with only 2–4% in public databases (P<0.003). We conclude that the rare variation in genes responsible for congenital arrhythmia syndromes is frequent in diLQTS. Our findings demonstrate that diLQTS is a pharmacogenomic syndrome predisposed by rare genetic variants.

Exome sequencing implicates an increased burden of rare potassium channel variants in the risk of drug-induced long QT interval syndrome.

OBJECTIVES The aim of this study was to test the hypothesis that rare variants are associated with drug-induced long QT interval syndrome (diLQTS) and torsades de pointes. BACKGROUND diLQTS is associated with the potentially fatal arrhythmia torsades de pointes. The contribution of rare genetic variants to the underlying genetic framework predisposing to diLQTS has not been systematically examined. METHODS We performed whole-exome sequencing on 65 diLQTS patients and 148 drug-exposed control subjects of European descent. We used rare variant analyses (variable threshold and sequence kernel association test) and gene-set analyses to identify genes enriched with rare amino acid coding (AAC) variants associated with diLQTS. Significant associations were reanalyzed by comparing diLQTS patients with 515 ethnically matched control subjects from the National Heart, Lung, and Blood Grand Opportunity Exome Sequencing Project. RESULTS Rare variants in 7 genes were enriched in the diLQTS patients according to the sequence kernel association test or variable threshold compared with drug-exposed controls (p < 0.001). Of these, we replicated the diLQTS associations for KCNE1 and ACN9 using 515 Exome Sequencing Project control subjects (p < 0.05). A total of 37% of the diLQTS patients also had 1 or more rare AAC variants compared with 21% of control subjects (p = 0.009), in a pre-defined set of 7 congenital long QT interval syndrome (cLQTS) genes encoding potassium channels or channel modulators (KCNE1, KCNE2, KCNH2, KCNJ2, KCNJ5, KCNQ1, AKAP9). CONCLUSIONS By combining whole-exome sequencing with aggregated rare variant analyses, we implicate rare variants in KCNE1 and ACN9 as risk factors for diLQTS. Moreover, diLQTS patients were more burdened by rare AAC variants in cLQTS genes encoding potassium channel modulators, supporting the idea that multiple rare variants, notably across cLQTS genes, predispose to diLQTS.

The use of next generation sequencing technology to study the effect of radiation therapy on mitochondrial DNA mutation

The human mitochondrial genome has an exclusively maternal mode of inheritance. Mitochondrial DNA (mtDNA) is particularly vulnerable to environmental insults due in part to an underdeveloped DNA repair system, limited to base excision and homologous recombination repair. Radiation exposure to the ovaries may cause mtDNA mutations in oocytes, which may in turn be transmitted to offspring. We hypothesized that the children of female cancer survivors who received radiation therapy may have an increased rate of mtDNA heteroplasmy mutations, which conceivably could increase their risk of developing cancer and other diseases. We evaluated 44 DNA blood samples from 17 Danish and 1 Finnish families (18 mothers and 26 children). All mothers had been treated for cancer as children and radiation doses to their ovaries were determined based on medical records and computational models. DNA samples were sequenced for the entire mitochondrial genome using the Illumina GAII system. Mothers age at sample collection was positively correlated with mtDNA heteroplasmy mutations. There was evidence of heteroplasmy inheritance in that 9 of the 18 families had at least one child who inherited at least one heteroplasmy site from his or her mother. No significant difference in single nucleotide polymorphisms between mother and offspring, however, was observed. Radiation therapy dose to ovaries also was not significantly associated with the heteroplasmy mutation rate among mothers and children. No evidence was found that radiotherapy for pediatric cancer is associated with the mitochondrial genome mutation rate in female cancer survivors and their children.

Genotype and risk of major bleeding during warfarin treatment

AIM To determine whether genetic variants associated with warfarin dose variability were associated with increased risk of major bleeding during warfarin therapy. MATERIALS & METHODS Using Vanderbilts DNA biobank we compared the prevalence of CYP2C9, VKORC1 and CYP4F2 variants in 250 cases with major bleeding and 259 controls during warfarin therapy. RESULTS CYP2C9*3 was the only allele that differed significantly among cases (14.2%) and controls (7.8%; p = 0.022). In the 214 (85.6%) cases with a major bleed 30 or more days after warfarin initiation, CYP2C9*3 was the only variant associated with bleeding (adjusted odds ratio: 2.05; 95% CI: 1.04, 4.04). CONCLUSION The CYP2C9*3 allele may double the risk of major bleeding among patients taking warfarin for 30 or more days.

Phenome Wide Association Studies demonstrating pleiotropy of genetic variants within FTO with and without adjustment for body mass index

Phenome-wide association studies (PheWAS) have demonstrated utility in validating genetic associations derived from traditional genetic studies as well as identifying novel genetic associations. Here we used an electronic health record (EHR)-based PheWAS to explore pleiotropy of genetic variants in the fat mass and obesity associated gene (FTO), some of which have been previously associated with obesity and type 2 diabetes (T2D). We used a population of 10,487 individuals of European ancestry with genome-wide genotyping from the Electronic Medical Records and Genomics (eMERGE) Network and another population of 13,711 individuals of European ancestry from the BioVU DNA biobank at Vanderbilt genotyped using Illumina HumanExome BeadChip. A meta-analysis of the two study populations replicated the well-described associations between FTO variants and obesity (odds ratio [OR] = 1.25, 95% Confidence Interval = 1.11–1.24, p = 2.10 × 10−9) and FTO variants and T2D (OR = 1.14, 95% CI = 1.08–1.21, p = 2.34 × 10−6). The meta-analysis also demonstrated that FTO variant rs8050136 was significantly associated with sleep apnea (OR = 1.14, 95% CI = 1.07–1.22, p = 3.33 × 10−5); however, the association was attenuated after adjustment for body mass index (BMI). Novel phenotype associations with obesity-associated FTO variants included fibrocystic breast disease (rs9941349, OR = 0.81, 95% CI = 0.74–0.91, p = 5.41 × 10−5) and trends toward associations with non-alcoholic liver disease and gram-positive bacterial infections. FTO variants not associated with obesity demonstrated other potential disease associations including non-inflammatory disorders of the cervix and chronic periodontitis. These results suggest that genetic variants in FTO may have pleiotropic associations, some of which are not mediated by obesity.

A genome-wide association study of heparin-induced thrombocytopenia using an electronic medical record

Heparin-induced thrombocytopenia (HIT) is an unpredictable, potentially catastrophic adverse effect of heparin treatment resulting from an immune response to platelet factor 4 (PF4)/heparin complexes. No genome-wide evaluations have been performed to identify potential genetic influences on HIT. Here, we performed a genome-wide association study (GWAS) and candidate gene study using HIT cases and controls identified using electronic medical records (EMRs) coupled to a DNA biobank and attempted to replicate GWAS associations in an independent cohort. We subsequently investigated influences of GWAS-associated single nucleotide polymorphisms (SNPs) on PF4/heparin antibodies in non-heparin treated individuals. In a recessive model, we observed significant SNP associations (odds ratio [OR] 18.52; 95% confidence interval [CI] 6.33-54.23; p=3.18×10(-9)) with HIT near the T-Cell Death-Associated Gene 8 (TDAG8). These SNPs are in linkage disequilibrium with a missense TDAG8 SNP. TDAG8 SNPs trended toward an association with HIT in replication analysis (OR 5.71; 0.47-69.22; p=0.17), and the missense SNP was associated with PF4/heparin antibody levels and positive PF4/heparin antibodies in non-heparin treated patients (OR 3.09; 1.14-8.13; p=0.02). In the candidate gene study, SNPs at HLA-DRA were nominally associated with HIT (OR 0.25; 0.15-0.44; p=2.06×10(-6)). Further study of TDAG8 and HLA-DRA SNPs is warranted to assess their influence on the risk of developing HIT.