Andrew C. Edmondson

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a “cosmopolitan” tagging approach to capture the genetic diversity across ∼2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Loss-of-function variants in endothelial lipase are a cause of elevated HDL cholesterol in humans

Elevated plasma concentrations of HDL cholesterol (HDL-C) are associated with protection from atherosclerotic cardiovascular disease. Animal models indicate that decreased expression of endothelial lipase (LIPG) is inversely associated with HDL-C levels, and genome-wide association studies have identified LIPG variants as being associated with HDL-C levels in humans. We hypothesized that loss-of-function mutations in LIPG may result in elevated HDL-C and therefore performed deep resequencing of LIPG exons in cases with elevated HDL-C levels and controls with decreased HDL-C levels. We identified a significant excess of nonsynonymous LIPG variants unique to cases with elevated HDL-C. In vitro lipase activity assays demonstrated that these variants significantly decreased endothelial lipase activity. In addition, a meta-analysis across 5 cohorts demonstrated that the low-frequency Asn396Ser variant is significantly associated with increased HDL-C, while the common Thr111Ile variant is not. Functional analysis confirmed that the Asn396Ser variant has significantly decreased lipase activity both in vitro and in vivo, while the Thr111Ile variant has normal lipase activity. Our results establish that loss-of-function mutations in LIPG lead to increased HDL-C levels and support the idea that inhibition of endothelial lipase may be an effective mechanism to raise HDL-C.

Hepatic sortilin regulates both apolipoprotein B secretion and LDL catabolism

Genome-wide association studies (GWAS) have identified a genetic variant at a locus on chromosome 1p13 that is associated with reduced risk of myocardial infarction, reduced plasma levels of LDL cholesterol (LDL-C), and markedly increased expression of the gene sortilin-1 (SORT1) in liver. Sortilin is a lysosomal sorting protein that binds ligands both in the Golgi apparatus and at the plasma membrane and traffics them to the lysosome. We previously reported that increased hepatic sortilin expression in mice reduced plasma LDL-C levels. Here we show that increased hepatic sortilin not only reduced hepatic apolipoprotein B (APOB) secretion, but also increased LDL catabolism, and that both effects were dependent on intact lysosomal targeting. Loss-of-function studies demonstrated that sortilin serves as a bona fide receptor for LDL in vivo in mice. Our data are consistent with a model in which increased hepatic sortilin binds intracellular APOB-containing particles in the Golgi apparatus as well as extracellular LDL at the plasma membrane and traffics them to the lysosome for degradation. We thus provide functional evidence that genetically increased hepatic sortilin expression both reduces hepatic APOB secretion and increases LDL catabolism, providing dual mechanisms for the very strong association between increased hepatic sortilin expression and reduced plasma LDL-C levels in humans.

Dense Genotyping of Candidate Gene Loci Identifies Variants Associated With High-Density Lipoprotein Cholesterol

Background—Plasma levels of high-density lipoprotein cholesterol (HDL-C) are known to be heritable, but only a fraction of the heritability is explained. We used a high-density genotyping array containing single-nucleotide polymorphisms (SNPs) from HDL-C candidate genes selected on known biology of HDL-C metabolism, mouse genetic studies, and human genetic association studies. SNP selection was based on tagging SNPs and included low-frequency nonsynonymous SNPs. Methods and Results—Association analysis in a cohort containing extremes of HDL-C (case-control, n=1733) provided a discovery phase, with replication in 3 additional populations for a total meta-analysis in 7857 individuals. We replicated the majority of loci identified through genome-wide association studies and present on the array (including ABCA1, APOA1/C3/A4/A5, APOB, APOE/C1/C2, CETP, CTCF-PRMT8, FADS1/2/3, GALNT2, LCAT, LILRA3, LIPC, LIPG, LPL, LRP4, SCARB1, TRIB1, ZNF664) and provide evidence that suggests an association in several previously unreported candidate gene loci (including ABCG1, GPR109A/B/81, NFKB1, PON1/2/3/4). There was evidence for multiple, independent association signals in 5 loci, including association with low-frequency nonsynonymous variants. Conclusions—Genetic loci associated with HDL-C are likely to harbor multiple, independent causative variants, frequently with opposite effects on the HDL-C phenotype. Cohorts comprising subjects at the extremes of the HDL-C distribution may be efficiently used in a case-control discovery of quantitative traits.

An Arabidopsis homologue of bacterial RecA that complements an E. coli recA deletion is targeted to plant mitochondria

Homologous recombination results in the exchange and rearrangement of DNA, and thus generates genetic variation in living organisms. RecA is known to function in all bacteria as the central enzyme catalyzing strand transfer and has functional homologues in eukaryotes. Most of our knowledge of homologous recombination in eukaryotes is limited to processes in the nucleus. The mitochondrial genomes of higher plants contain repeated sequences that are known to undergo frequent rearrangements and recombination events. However, very little is known about the proteins involved or the biochemical mechanisms of DNA recombination in plant mitochondria. We provide here the first report of an Arabidopsis thaliana homologue of Escherichia coli RecA that is targeted to mitochondria. The mtrecA gene has a putative mitochondrial presequence identified from the A. thaliana genome database. This nuclear gene encodes a predicted product that shows highest sequence homology to chloroplast RecA and RecA proteins from proteobacteria. When fused to the GFP coding sequence, the predicted presequence was able to target the fusion protein to isolated mitochondria but not to chloroplasts. The mitochondrion-specific localization of the mtrecA gene product was confirmed by Western analysis using polyclonal antibodies raised against a synthetic peptide from a unique region of the mature mtRecA. The Arabidopsis mtrecA gene partially complemented a recA deletion in E. coli, enhancing survival after exposure to DNA-damaging agents. These results suggest a possible role for mtrecA in homologous recombination and/or repair in Arabidopsis mitochondria.

Germline gain-of-function mutations in AFF4 cause a developmental syndrome functionally linking the super elongation complex and cohesin

Transcriptional elongation is critical for gene expression regulation during embryogenesis. The super elongation complex (SEC) governs this process by mobilizing paused RNA polymerase II (RNAP2). Using exome sequencing, we discovered missense mutations in AFF4, a core component of the SEC, in three unrelated probands with a new syndrome that phenotypically overlaps Cornelia de Lange syndrome (CdLS) that we have named CHOPS syndrome (C for cognitive impairment and coarse facies, H for heart defects, O for obesity, P for pulmonary involvement and S for short stature and skeletal dysplasia). Transcriptome and chromatin immunoprecipitation sequencing (ChIP-seq) analyses demonstrated similar alterations of genome-wide binding of AFF4, cohesin and RNAP2 in CdLS and CHOPS syndrome. Direct molecular interaction of the SEC, cohesin and RNAP2 was demonstrated. These data support a common molecular pathogenesis for CHOPS syndrome and CdLS caused by disturbance of transcriptional elongation due to alterations in genome-wide binding of AFF4 and cohesin.

Characterization of a mitochondrially targeted single-stranded DNA-binding protein in Arabidopsis thaliana

A gene encoding a predicted mitochondrially targeted single-stranded DNA binding protein (mtSSB) was identified in the Arabidopsis thaliana genome sequence. This gene (At4g11060) codes for a protein of 201 amino acids, including a 28-residue putative mitochondrial targeting transit peptide. Protein sequence alignment shows high similarity between the mtSSB protein and single-stranded DNA binding proteins (SSB) from bacteria, including residues conserved for SSB function. Phylogenetic analysis indicates a close relationship between this protein and other mitochondrially targeted SSB proteins. The predicted targeting sequence was fused with the GFP coding region, and the organellar localization of the expressed fusion protein was determined. Specific targeting to mitochondria was observed in in-vitro import experiments and by transient expression of a GFP fusion construct in Arabidopsis leaves after microprojectile bombardment. The mature mtSSB coding region was overexpressed in Escherichia coli and the protein was purified for biochemical characterization. The purified protein binds single-stranded, but not double-stranded, DNA. MtSSB stimulates the homologous strand-exchange activity of E. coli RecA. These results indicate that mtSSB is a functional homologue of the E. coli SSB, and that it may play a role in mitochondrial DNA recombination.

A naturally occurring variant of endothelial lipase associated with elevated HDL exhibits impaired synthesis

Human endothelial lipase (EL) is a member of a family of lipases and phospholipases that are involved in the metabolism of plasma lipoproteins. EL displays a preference to hydrolyze lipids in HDL. We report here that a naturally occurring low frequency coding variant in the EL gene (LIPG), glycine-26 to serine (G26S), is significantly more common in African-American individuals with elevated HDL cholesterol (HDL-C) levels. To test the hypothesis that this variant results in reduced EL function, we extensively characterized and compared the catalytic and noncatalytic functions of the G26S variant and wild-type (WT) EL. While the catalytic-specific activity of G26S EL is similar to WT EL, its secretion is markedly reduced. Consistent with this observation, we found that carriers of the G26S variant had significantly reduced plasma levels of EL protein. Thus, this N-terminal variant results in reduced secretion of EL protein, plausibly leading to increased HDL-C levels.

Mining the LIPG Allelic Spectrum Reveals the Contribution of Rare and Common Regulatory Variants to HDL Cholesterol

Genome-wide association studies (GWAS) have successfully identified loci associated with quantitative traits, such as blood lipids. Deep resequencing studies are being utilized to catalogue the allelic spectrum at GWAS loci. The goal of these studies is to identify causative variants and missing heritability, including heritability due to low frequency and rare alleles with large phenotypic impact. Whereas rare variant efforts have primarily focused on nonsynonymous coding variants, we hypothesized that noncoding variants in these loci are also functionally important. Using the HDL-C gene LIPG as an example, we explored the effect of regulatory variants identified through resequencing of subjects at HDL-C extremes on gene expression, protein levels, and phenotype. Resequencing a portion of the LIPG promoter and 5′ UTR in human subjects with extreme HDL-C, we identified several rare variants in individuals from both extremes. Luciferase reporter assays were used to measure the effect of these rare variants on LIPG expression. Variants conferring opposing effects on gene expression were enriched in opposite extremes of the phenotypic distribution. Minor alleles of a common regulatory haplotype and noncoding GWAS SNPs were associated with reduced plasma levels of the LIPG gene product endothelial lipase (EL), consistent with its role in HDL-C catabolism. Additionally, we found that a common nonfunctional coding variant associated with HDL-C (rs2000813) is in linkage disequilibrium with a 5′ UTR variant (rs34474737) that decreases LIPG promoter activity. We attribute the gene regulatory role of rs34474737 to the observed association of the coding variant with plasma EL levels and HDL-C. Taken together, the findings show that both rare and common noncoding regulatory variants are important contributors to the allelic spectrum in complex trait loci.

Gene-based interaction analysis by incorporating external linkage disequilibrium information

Gene–gene interactions have an important role in complex human diseases. Detection of gene–gene interactions has long been a challenge due to their complexity. The standard method aiming at detecting SNP–SNP interactions may be inadequate as it does not model linkage disequilibrium (LD) among SNPs in each gene and may lose power due to a large number of comparisons. To improve power, we propose a principal component (PC)-based framework for gene-based interaction analysis. We analytically derive the optimal weight for both quantitative and binary traits based on pairwise LD information. We then use PCs to summarize the information in each gene and test for interactions between the PCs. We further extend this gene-based interaction analysis procedure to allow the use of imputation dosage scores obtained from a popular imputation software package, MACH, which incorporates multilocus LD information. To evaluate the performance of the gene-based interaction tests, we conducted extensive simulations under various settings. We demonstrate that gene-based interaction tests are more powerful than SNP-based tests when more than two variants interact with each other; moreover, tests that incorporate external LD information are generally more powerful than those that use genotyped markers only. We also apply the proposed gene-based interaction tests to a candidate gene study on high-density lipoprotein. As our method operates at the gene level, it can be applied to a genome-wide association setting and used as a screening tool to detect gene–gene interactions.