Experimental and Computational Analysis of Newly Identified Pathogenic Mutations in the Creatine Transporter SLC6A8

Creatine is an essential metabolite for the storage and rapid supply of energy in muscle and nerve cells. In humans, impaired metabolism, transport, and distribution of creatine throughout tissues can cause varying forms of mental disability, also known as creatine deficiency syndrome (CDS). So far, 80 mutations in the creatine transporter (SLC6A8) have been associated to CDS. To better understand the effect of human genetic variants on the physiology of SLC6A8 and their possible impact on CDS, we studied 30 missense variants including 15 variants of unknown significance, two of which are reported here for the first time. We expressed these variants in HEK293 cells and explored their subcellular localization and transport activity. We also applied computational methods to predict variant effect and estimate site-specific changes in thermodynamic stability. To explore variants that might have a differential effect on the transporter's conformers along the transport cycle, we constructed homology models of the inward facing, and outward facing conformations. In addition, we used mass-spectrometry to study proteins that interact with wild type SLC6A8 and five selected variants in HEK293 cells. In silico models of the protein complexes revealed how two variants impact the interaction interface of SLC6A8 with other proteins and how pathogenic variants lead to an enrichment of ER protein partners. Overall, our integrated analysis disambiguates the pathogenicity of 15 variants of unknown significance revealing diverse mechanisms of pathogenicity, including two previously unreported variants obtained from patients suffering from the creatine deficiency syndrome.

Proteomic analysis of 92 circulating proteins and their effects in cardiometabolic diseases

BACKGROUND

Human plasma contains a wide variety of circulating proteins. These proteins can be important clinical biomarkers in disease and also possible drug targets. Large scale genomics studies of circulating proteins can identify genetic variants that lead to relative protein abundance.

METHODS

We conducted a meta-analysis on genome-wide association studies of autosomal chromosomes in 22,997 individuals of primarily European ancestry across 12 cohorts to identify protein quantitative trait loci (pQTL) for 92 cardiometabolic associated plasma proteins.

RESULTS

We identified 503 (337 cis and 166 trans) conditionally independent pQTLs, including several novel variants not reported in the literature. We conducted a sex-stratified analysis and found that 118 (23.5%) of pQTLs demonstrated heterogeneity between sexes. The direction of effect was preserved but there were differences in effect size and significance. Additionally, we annotate trans-pQTLs with nearest genes and report plausible biological relationships. Using Mendelian randomization, we identified causal associations for 18 proteins across 19 phenotypes, of which 10 have additional genetic colocalization evidence. We highlight proteins associated with a constellation of cardiometabolic traits including angiopoietin-related protein 7 (ANGPTL7) and Semaphorin 3F (SEMA3F).

CONCLUSION

Through large-scale analysis of protein quantitative trait loci, we provide a comprehensive overview of common variants associated with plasma proteins. We highlight possible biological relationships which may serve as a basis for further investigation into possible causal roles in cardiometabolic diseases.

Causal analysis of plasma IL-8 on carotid intima media thickness, a measure of subclinical atherosclerosis

BACKGROUND

We investigated the causality of IL-8 on carotid intima-media thickness (c-IMT), a measure of sub-clinical atherosclerosis.

METHODS

The IMPROVE is a multicenter European study (n = 3,711). The association of plasma IL-8 with c-IMT (mm) was estimated by quantile regression. Genotyping was performed using the Illumina CardioMetabo and Immuno chips. Replication was attempted in three independent studies and a meta-analysis was performed using a random model.

RESULTS

In IMPROVE, each unit increase in plasma IL-8 was associated with an increase in median c-IMT measures (all p<0·03) in multivariable analyses. Linear regression identified rs117518778 and rs8057084 as associated with IL-8 levels and with measures of c-IMT. The two SNPs were combined in an IL-8-increasing genetic risk that showed causality of IL-8 on c-IMT in IMPROVE and in the UK Biobank (n = 22,179). The effect of IL-8 on c-IMT measures was confirmed in PIVUS (n = 1,016) and MDCCC (n = 6,103). The association of rs8057084 with c-IMT was confirmed in PIVUS and UK Biobank with a pooled estimate effect (β) of -0·006 with 95%CI (-0·008- -0·003).

CONCLUSION

Our results indicate that genetic variants associated with plasma IL-8 also associate with c-IMT. However, we cannot infer causality of this association, as these variants lie outside of the IL8 locus.

FinnGen provides genetic insights from a well-phenotyped isolated population

Population isolates such as those in Finland benefit genetic research because deleterious alleles are often concentrated on a small number of low-frequency variants (0.1% ≤ minor allele frequency < 5%). These variants survived the founding bottleneck rather than being distributed over a large number of ultrarare variants. Although this effect is well established in Mendelian genetics, its value in common disease genetics is less explored^1,2. FinnGen aims to study the genome and national health register data of 500,000 Finnish individuals. Given the relatively high median age of participants (63 years) and the substantial fraction of hospital-based recruitment, FinnGen is enriched for disease end points. Here we analyse data from 224,737 participants from FinnGen and study 15 diseases that have previously been investigated in large genome-wide association studies (GWASs). We also include meta-analyses of biobank data from Estonia and the United Kingdom. We identified 30 new associations, primarily low-frequency variants, enriched in the Finnish population. A GWAS of 1,932 diseases also identified 2,733 genome-wide significant associations (893 phenome-wide significant (PWS), P < 2.6 × 10^-11) at 2,496 (771 PWS) independent loci with 807 (247 PWS) end points. Among these, fine-mapping implicated 148 (73 PWS) coding variants associated with 83 (42 PWS) end points. Moreover, 91 (47 PWS) had an allele frequency of <5% in non-Finnish European individuals, of which 62 (32 PWS) were enriched by more than twofold in Finland. These findings demonstrate the power of bottlenecked populations to find entry points into the biology of common diseases through low-frequency, high impact variants.

Abstract 105: Proprotein Convertase Subtilisin/kexin Type 6 (Pcsk6) Plays A Role In Angiogenesis Regulation

Background: PCSK6 cleaves and activates growth factors involved in cell differentiation and Pcsk6-/- mice exhibit 25% embryonic lethality, while surviving offspring have blindness and cyclopia. We have shown that PCSK6 localizes to vascular smooth muscle cells of intra-plaque neovessels. We therefore hypothesize that the processes of angiogenesis, atherosclerosis and plaque neovascularisation share common molecular regulators and that PCSK6 plays a role in these processes.

Methods: Association between SNPs near the PCSK6 gene and levels of 90 plasma proteins was tested using data from 30,931 individuals in the SCALLOP consortium. Associations of angiogenic factors with PCSK6 (transcriptomic and proteomic) and patient symptomatology were explored in a biobank of plaque endarterectomies (BiKE). Morpholino technology was used for Pcsk6 knockdown in zebrafish, coupled with vascular phenotype studies. Pcsk6-/- mice were used to study the expression of vasculogenesis related genes and vascular morphology at baseline (12 weeks).

Results: Rs7178801 and rs11639051 in the PCSK6 gene were found to relate to plasma PDGFB and VEGFD levels. TGFB2, VEGFA, VEGFC, PDGFA and PLAUR associated with patients symptomatology in BiKE and correlated with PCSK6 levels. Gross phenotypic and histological examination of zebrafish embryos with ablated PCSK6 showed improper peripheral vascular patterning of intersegmental vessels along with cerebral and myocardial haemorrhage hemorrhage. Expression of angiogenic factors; Vegfa, Vegfb, Angpt1 and Tgfb2 were altered in the heart, liver and adipose tissue in Pcsk6-/- mice. Vessel density, length and branch numbers were reduced in the Pcsk6 -/- compared to the controls in the superficial but more significantly in the deep peripheral vessels. The smooth muscle cell coverage of the retinal vasculature was also reduced in the Pcsk6 -/- mice vs controls and Pdgfa was significantly downregulated compared to controls in mouse retinas.

Conclusions: PCSK6 associates with levels of angiogenic factors on the genetic, transcriptomic and proteomic level. Lack of Pcsk6 in mice and zebrafish led to vascular patterning defects. Further planned studies will examine PCSK6 in pathological neovascularization.

Abstract 217: PCSK6 Is A Key Regulator Of Immune Status In Mice And Its Ablation Increases Atherosclerotic Plaque Burden In A Bone Marrow Transplant Model

Introduction: Proprotein convertase subtilisins/kexins (PCSKs) activate cytokines and growth factors and have been implicated in various cancers. We have previously shown that PCSK6 is a key protease modulating smooth muscle cell response in atherosclerosis, but its expression also correlated positively with typical markers of T lymphocytes and macrophages in plaques and it was localised in the proximity of these cells. Here, we hypothesized that PCSK6 may be involved in modulating inflammatory responses and aimed to elucidate its role in a hyperlipidemic mice model.

Methods: Detailed immunophenotyping using histology, FACS-, OLINK- and ELISA-based analyses and primary cell cultures was used to compare Pcsk6 -/- and littermate controls. Atherosclerosis was evaluated in Ldlr -/- mice upon bone marrow transplant with Pcsk6 -/- or wild-type bone-marrow.

Results: At baseline Pcks6 -/- mice showed an enrichment of pro-inflammatory cytokines Ccl2, Ccl3, Ccl20, Cxcl1 and in particular Il17a and Il17f in plasma compared controls. Spleens of Pcsk6 -/- mice had an increased number of germinal centres and FACS analysis showed that they contained significantly more CD8+ T cells. Microarray analysis of spleens from Pcsk6 -/- vs. controls confirmed that T cell markers CD4, CD3E and CD3G were upregulated. In vitro , splenocytes isolated from Pcsk6 -/- mice secreted more IFN-γ, IL-2 and IL-10 than controls upon stimulation with α-CD3 and α-CD28 antibodies. Moreover, peritoneal macrophages from Pcsk6 -/- mice secreted more TNF-α, MCP-1, IL-6 and IL-10 compared to control mice upon LPS stimulation. Interestingly, bone marrow derived macrophages from Pcsk6 -/- mice were also more prone to lipid uptake. Finally, in vivo transplantation of Pcsk6 -/- bone marrow to Ldlr -/- mice led to increased atherosclerotic plaque burden compared to controls, as quantified in the aortic root.

Conclusions: PCSK6 ablation led to increased number of CD8+ T cells, as well as macrophage and cytokine activation. Transplantation of Pcsk6 -/- bone marrow resulted in an increase in atherosclerotic plaque burden compared to controls. Taken together, these results indicate that PCSK6 is a key regulator of the immune system, though the exact mechanisms involved require further investigation.

Technological readiness and implementation of genomic-driven precision medicine for complex diseases

The fields of human genetics and genomics have generated considerable knowledge about the mechanistic basis of many diseases. Genomic approaches to diagnosis, prognostication, prevention and treatment - genomic-driven precision medicine (GDPM) - may help optimize medical practice. Here, we provide a comprehensive review of GDPM of complex diseases across major medical specialties. We focus on technological readiness: how rapidly a test can be implemented into health care. Although these areas of medicine are diverse, key similarities exist across almost all areas. Many medical areas have, within their standards of care, at least one GDPM test for a genetic variant of strong effect that aids the identification/diagnosis of a more homogeneous subset within a larger disease group or identifies a subset with different therapeutic requirements. However, for almost all complex diseases, the majority of patients do not carry established single-gene mutations with large effects. Thus, research is underway that seeks to determine the polygenic basis of many complex diseases. Nevertheless, most complex diseases are caused by the interplay of genetic, behavioural and environmental risk factors, which will likely necessitate models for prediction and diagnosis that incorporate genetic and non-genetic data.

Abstract MP02: Proprotein Convertase Subtilisin/kexin Type 6 (PCSK6) Is Involved In Regulation Of Vasculogenesis

Background: PCSK6 cleaves and activates growth factors involved in cell differentiation and Pcsk6-/- mice exhibit 25% embryonic lethality, while surviving offspring have blindness and cyclopia. We have shown that PCSK6 localizes to vascular smooth muscle cells of intra-plaque neovessels.

Hypothesis: We hypothesize that a) the processes of vasculogenesis, atherosclerosis and plaque neovascularisation share common molecular regulators and b) that PCSK6 plays a common role in vascular remodeling.

Methods: Association between SNPs near the PCSK6 gene and levels of 90 plasma proteins was tested using data from 30,931 individuals in the SCALLOP consortium. Data from the FinnGen consortium (Release 3) were queried to explore the PCSK6 SNPs associated with specific disease phenotypes. Morpholino technology was used for Pcsk6 knockdown in zebrafish, coupled with vascular phenotype studies. Pcsk6-/- mice were used to study the expression of vasculogenesis related genes as well as vascular morphology.

Results: Rs7178801 and rs11639051 in the PCSK6 gene were found to relate to plasma PDGFB and VEGFD levels, while rs45482895 associated with diabetic retinopathy and disorders of the choroid and retina in humans. Gross phenotypic and histological examination of zebrafish embryos with ablated PCSK6 showed improper peripheral vascular patterning of intersegmental vessels, cerebral and myocardial haemorrhage, atrial dilatation and pericardial oedema. Gene expression levels of vasculogenesis related genes, i.e. Vegfa, Vegfb, Angpt1 and Tgfb2 were altered in the heart, liver and adipose tissue in Pcsk6-/- mice. In the retinas of Pcsk6-/- mice, Pdgfb was significantly downregulated compared to controls. Staining of the retinal vasculature showed that superficial retinal vasculature had a lower number of branching points and vessel width, while deep retinal vasculature covered a significantly smaller area in Pcsk6-/-

Conclusions: Here, we show that variants in the PCSK6 gene associate with growth factors implicated in vasculogenesis and retinopathy. Lack of Pcsk6 in mice and zebrafish led to vascular patterning defects. Further studies will aim to elucidate the mechanisms by which PCSK6 regulates vasculogenesis in homeostatic and pathological conditions

CDH6 and HAGH protein levels in plasma associate with Alzheimer's disease in APOE ε4 carriers

Many Alzheimer’s disease (AD) genes including Apolipoprotein E (APOE) are found to be expressed in blood-derived macrophages and thus may alter blood protein levels. We measured 91 neuro-proteins in plasma from 316 participants of the Rotterdam Study (incident AD = 161) using Proximity Extension Ligation assay. We studied the association of plasma proteins with AD in the overall sample and stratified by APOE. Findings from the Rotterdam study were replicated in 186 AD patients of the BioFINDER study. We further evaluated the correlation of these protein biomarkers with total tau (t-tau), phosphorylated tau (p-tau) and amyloid-beta (Aβ) 42 levels in cerebrospinal fluid (CSF) in the Amsterdam Dementia Cohort (N = 441). Finally, we conducted a genome-wide association study (GWAS) to identify the genetic variants determining the blood levels of AD-associated proteins. Plasma levels of the proteins, CDH6 (β = 0.638, P = 3.33 × 10⁻⁴) and HAGH (β = 0.481, P = 7.20 × 10⁻⁴), were significantly elevated in APOE ε4 carrier AD patients. The findings in the Rotterdam Study were replicated in the BioFINDER study for both CDH6 (β = 1.365, P = 3.97 × 10⁻³) and HAGH proteins (β = 0.506, P = 9.31 × 10⁻⁷) when comparing cases and controls in APOE ε4 carriers. In the CSF, CDH6 levels were positively correlated with t-tau and p-tau in the total sample as well as in APOE ε4 stratum (P < 1 × 10⁻³). The HAGH protein was not detected in CSF. GWAS of plasma CDH6 protein levels showed significant association with a cis-regulatory locus (rs111283466, P = 1.92 × 10⁻⁹). CDH6 protein is implicated in cell adhesion and synaptogenesis while HAGH protein is related to the oxidative stress pathway. Our findings suggest that these pathways may be altered during presymptomatic AD and that CDH6 and HAGH may be new blood-based biomarkers.

Genome-wide association and Mendelian randomisation analysis provide insights into the pathogenesis of heart failure

Heart failure (HF) is a leading cause of morbidity and mortality worldwide. A small proportion of HF cases are attributable to monogenic cardiomyopathies and existing genome-wide association studies (GWAS) have yielded only limited insights, leaving the observed heritability of HF largely unexplained. We report results from a GWAS meta-analysis of HF comprising 47,309 cases and 930,014 controls. Twelve independent variants at 11 genomic loci are associated with HF, all of which demonstrate one or more associations with coronary artery disease (CAD), atrial fibrillation, or reduced left ventricular function, suggesting shared genetic aetiology. Functional analysis of non-CAD-associated loci implicate genes involved in cardiac development (MYOZ1, SYNPO2L), protein homoeostasis (BAG3), and cellular senescence (CDKN1A). Mendelian randomisation analysis supports causal roles for several HF risk factors, and demonstrates CAD-independent effects for atrial fibrillation, body mass index, and hypertension. These findings extend our knowledge of the pathways underlying HF and may inform new therapeutic strategies.

Multiplex proteomics identifies novel CSF and plasma biomarkers of early Alzheimer's disease

To date, the development of disease-modifying therapies for Alzheimer’s disease (AD) has largely focused on the removal of amyloid beta Aβ fragments from the CNS. Proteomic profiling of patient fluids may help identify novel therapeutic targets and biomarkers associated with AD pathology. Here, we applied the Olink™ ProSeek immunoassay to measure 270 CSF and plasma proteins across 415 Aβ- negative cognitively normal individuals (Aβ- CN), 142 Aβ-positive CN (Aβ+ CN), 50 Aβ- mild cognitive impairment (MCI) patients, 75 Aβ+ MCI patients, and 161 Aβ+ AD patients from the Swedish BioFINDER study. A validation cohort included 59 Aβ- CN, 23 Aβ- + CN, 44 Aβ- MCI and 53 Aβ+ MCI. To compare protein concentrations in patients versus controls, we applied multiple linear regressions adjusting for age, gender, medications, smoking and mean subject-level protein concentration, and corrected findings for false discovery rate (FDR, q < 0.05). We identified, and replicated, altered levels of ten CSF proteins in Aβ+ individuals, including CHIT1, SMOC2, MMP-10, LDLR, CD200, EIF4EBP1, ALCAM, RGMB, tPA and STAMBP (− 0.14 < d < 1.16; q < 0.05). We also identified and replicated alterations of six plasma proteins in Aβ+ individuals OSM, MMP-9, HAGH, CD200, AXIN1, and uPA (− 0.77 < d < 1.28; q < 0.05). Multiple analytes associated with cognitive performance and cortical thickness (q < 0.05). Plasma biomarkers could distinguish AD dementia (AUC = 0.94, 95% CI = 0.87–0.98) and prodromal AD (AUC = 0.78, 95% CI = 0.68–0.87) from CN. These findings reemphasize the contributions of immune markers, phospholipids, angiogenic proteins and other biomarkers downstream of, and potentially orthogonal to, Aβ- and tau in AD, and identify candidate biomarkers for earlier detection of neurodegeneration.

Application of Machine Learning Methods to Predict Non-Alcoholic Steatohepatitis (NASH) in Non-Alcoholic Fatty Liver (NAFL) Patients

Non-alcoholic fatty liver disease (NAFLD) is the leading cause of chronic liver disease worldwide. NAFLD patients have excessive liver fat (steatosis), without other liver diseases and without excessive alcohol consumption. NAFLD consists of a spectrum of conditions: benign steatosis or non-alcoholic fatty liver (NAFL), steatosis accompanied by inflammation and fibrosis or nonalcoholic steatohepatitis (NASH), and cirrhosis. Given a lack of clinical biomarkers and its asymptomatic nature, NASH is under-diagnosed. We use electronic health records from the Optum Analytics to (1) identify patients diagnosed with benign steatosis and NASH, and (2) train machine learning classifiers for NASH and healthy (non-NASH) populations to (3) predict NASH disease status on patients diagnosed with NAFL. Summarized temporal lab data for alanine aminotransferase, aspartate aminotransferase, and platelet counts, with basic demographic information and type 2 diabetes status were included in the models.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Exome-wide association study of plasma lipids in >300,000 individuals

We screened variants on an exome-focused genotyping array in >300,000 participants (replication in >280,000 participants) and identified 444 independent variants in 250 loci significantly associated with total cholesterol (TC), high-density-lipoprotein cholesterol (HDL-C), low-density-lipoprotein cholesterol (LDL-C), and/or triglycerides (TG). At two loci (JAK2 and A1CF), experimental analysis in mice showed lipid changes consistent with the human data. We also found that: (i) beta-thalassemia trait carriers displayed lower TC and were protected from coronary artery disease (CAD); (ii) excluding the CETP locus, there was not a predictable relationship between plasma HDL-C and risk for age-related macular degeneration; (iii) only some mechanisms of lowering LDL-C appeared to increase risk for type 2 diabetes (T2D); and (iv) TG-lowering alleles involved in hepatic production of TG-rich lipoproteins (TM6SF2 and PNPLA3) tracked with higher liver fat, higher risk for T2D, and lower risk for CAD, whereas TG-lowering alleles involved in peripheral lipolysis (LPL and ANGPTL4) had no effect on liver fat but decreased risks for both T2D and CAD.

Improving Assessment of Drug Safety Through Proteomics: Early Detection and Mechanistic Characterization of the Unforeseen Harmful Effects of Torcetrapib

Supplemental Digital Content is available in the text.

Background:

Early detection of adverse effects of novel therapies and understanding of their mechanisms could improve the safety and efficiency of drug development. We have retrospectively applied large-scale proteomics to blood samples from ILLUMINATE (Investigation of Lipid Level Management to Understand its Impact in Atherosclerotic Events), a trial of torcetrapib (a cholesterol ester transfer protein inhibitor), that involved 15 067 participants at high cardiovascular risk. ILLUMINATE was terminated at a median of 550 days because of significant absolute increases of 1.2% in cardiovascular events and 0.4% in mortality with torcetrapib. The aims of our analysis were to determine whether a proteomic analysis might reveal biological mechanisms responsible for these harmful effects and whether harmful effects of torcetrapib could have been detected early in the ILLUMINATE trial with proteomics.

Methods:

A nested case-control analysis of paired plasma samples at baseline and at 3 months was performed in 249 participants assigned to torcetrapib plus atorvastatin and 223 participants assigned to atorvastatin only. Within each treatment arm, cases with events were matched to controls 1:1. Main outcomes were a survey of 1129 proteins for discovery of biological pathways altered by torcetrapib and a 9-protein risk score validated to predict myocardial infarction, stroke, heart failure, or death.

Results:

Plasma concentrations of 200 proteins changed significantly with torcetrapib. Their pathway analysis revealed unexpected and widespread changes in immune and inflammatory functions, as well as changes in endocrine systems, including in aldosterone function and glycemic control. At baseline, 9-protein risk scores were similar in the 2 treatment arms and higher in participants with subsequent events. At 3 months, the absolute 9-protein derived risk increased in the torcetrapib plus atorvastatin arm compared with the atorvastatin-only arm by 1.08% (P=0.0004). Thirty-seven proteins changed in the direction of increased risk of 49 proteins previously associated with cardiovascular and mortality risk.

Conclusions:

Heretofore unknown effects of torcetrapib were revealed in immune and inflammatory functions. A protein-based risk score predicted harm from torcetrapib within just 3 months. A protein-based risk assessment embedded within a large proteomic survey may prove to be useful in the evaluation of therapies to prevent harm to patients.

Clinical Trial Registration:

URL: https://www.clinicaltrials.gov. Unique identifier: NCT00134264.

Genetic variants affecting cross-sectional lung function in adults show little or no effect on longitudinal lung function decline

BACKGROUND

Genome-wide association studies have identified numerous genetic regions that influence cross-sectional lung function. Longitudinal decline in lung function also includes a heritable component but the genetic determinants have yet to be defined.

OBJECTIVES

We aimed to determine whether regions associated with cross-sectional lung function were also associated with longitudinal decline and to seek novel variants which influence decline.

METHODS

We analysed genome-wide data from 4167 individuals from the Busselton Health Study cohort, who had undergone spirometry (12 695 observations across eight time points). A mixed model was fitted and weighted risk scores were calculated for the joint effect of 26 known regions on baseline and longitudinal changes in FEV₁ and FEV₁/FVC. Potential additional regions of interest were identified and followed up in two independent cohorts.

RESULTS

The 26 regions previously associated with cross-sectional lung function jointly showed a strong effect on baseline lung function (p=4.44×10^-16 for FEV₁/FVC) but no effect on longitudinal decline (p=0.160 for FEV₁/FVC). This was replicated in an independent cohort. 39 additional regions of interest (48 variants) were identified; these associations were not replicated in two further cohorts.

CONCLUSIONS

Previously identified genetic variants jointly have a strong effect on cross-sectional lung function in adults but little or no effect on the rate of decline of lung function. It is possible that they influence COPD risk through lung development. Although no genetic variants have yet been associated with lung function decline at stringent genome-wide significance, longitudinal change in lung function is heritable suggesting that there is scope for future discoveries.

Genetic invalidation of Lp-PLA2 as a therapeutic target: Large-scale study of five functional Lp-PLA2-lowering alleles

Aims Darapladib, a potent inhibitor of lipoprotein-associated phospholipase A2 (Lp-PLA2), has not reduced risk of cardiovascular disease outcomes in recent randomized trials. We aimed to test whether Lp-PLA2 enzyme activity is causally relevant to coronary heart disease. Methods In 72,657 patients with coronary heart disease and 110,218 controls in 23 epidemiological studies, we genotyped five functional variants: four rare loss-of-function mutations (c.109+2T > C (rs142974898), Arg82His (rs144983904), Val279Phe (rs76863441), Gln287Ter (rs140020965)) and one common modest-impact variant (Val379Ala (rs1051931)) in PLA2G7, the gene encoding Lp-PLA2. We supplemented de-novo genotyping with information on a further 45,823 coronary heart disease patients and 88,680 controls in publicly available databases and other previous studies. We conducted a systematic review of randomized trials to compare effects of darapladib treatment on soluble Lp-PLA2 activity, conventional cardiovascular risk factors, and coronary heart disease risk with corresponding effects of Lp-PLA2-lowering alleles. Results Lp-PLA2 activity was decreased by 64% ( p = 2.4 × 10-25) with carriage of any of the four loss-of-function variants, by 45% ( p < 10-300) for every allele inherited at Val279Phe, and by 2.7% ( p = 1.9 × 10-12) for every allele inherited at Val379Ala. Darapladib 160 mg once-daily reduced Lp-PLA2 activity by 65% ( p < 10-300). Causal risk ratios for coronary heart disease per 65% lower Lp-PLA2 activity were: 0.95 (0.88-1.03) with Val279Phe; 0.92 (0.74-1.16) with carriage of any loss-of-function variant; 1.01 (0.68-1.51) with Val379Ala; and 0.95 (0.89-1.02) with darapladib treatment. Conclusions In a large-scale human genetic study, none of a series of Lp-PLA2-lowering alleles was related to coronary heart disease risk, suggesting that Lp-PLA2 is unlikely to be a causal risk factor.

The Ser82 RAGE Variant Affects Lung Function and Serum RAGE in Smokers and sRAGE Production In Vitro

INTRODUCTION

Genome-Wide Association Studies have identified associations between lung function measures and Chronic Obstructive Pulmonary Disease (COPD) and chromosome region 6p21 containing the gene for the Advanced Glycation End Product Receptor (AGER, encoding RAGE). We aimed to (i) characterise RAGE expression in the lung, (ii) identify AGER transcripts, (iii) ascertain if SNP rs2070600 (Gly82Ser C/T) is associated with lung function and serum sRAGE levels and (iv) identify whether the Gly82Ser variant is functionally important in altering sRAGE levels in an airway epithelial cell model.

METHODS

Immunohistochemistry was used to identify RAGE protein expression in 26 human tissues and qPCR was used to quantify AGER mRNA in lung cells. Gene expression array data was used to identify AGER expression during lung development in 38 fetal lung samples. RNA-Seq was used to identify AGER transcripts in lung cells. sRAGE levels were assessed in cells and patient serum by ELISA. BEAS2B-R1 cells were transfected to overexpress RAGE protein with either the Gly82 or Ser82 variant and sRAGE levels identified.

RESULTS

Immunohistochemical assessment of 6 adult lung samples identified high RAGE expression in the alveoli of healthy adults and individuals with COPD. AGER/RAGE expression increased across developmental stages in human fetal lung at both the mRNA (38 samples) and protein levels (20 samples). Extensive AGER splicing was identified. The rs2070600T (Ser82) allele is associated with higher FEV1, FEV1/FVC and lower serum sRAGE levels in UK smokers. Using an airway epithelium model overexpressing the Gly82 or Ser82 variants we found that HMGB1 activation of the RAGE-Ser82 receptor results in lower sRAGE production.

CONCLUSIONS

This study provides new information regarding the expression profile and potential role of RAGE in the human lung and shows a functional role of the Gly82Ser variant. These findings advance our understanding of the potential mechanisms underlying COPD particularly for carriers of this AGER polymorphism.

Trans-ancestry meta-analyses identify rare and common variants associated with blood pressure and hypertension

High blood pressure is a major risk factor for cardiovascular disease and premature death. However, there is limited knowledge on specific causal genes and pathways. To better understand the genetics of blood pressure, we genotyped 242,296 rare, low-frequency and common genetic variants in up to 192,763 individuals and used ∼155,063 samples for independent replication. We identified 30 new blood pressure- or hypertension-associated genetic regions in the general population, including 3 rare missense variants in RBM47, COL21A1 and RRAS with larger effects (>1.5 mm Hg/allele) than common variants. Multiple rare nonsense and missense variant associations were found in A2ML1, and a low-frequency nonsense variant in ENPEP was identified. Our data extend the spectrum of allelic variation underlying blood pressure traits and hypertension, provide new insights into the pathophysiology of hypertension and indicate new targets for clinical intervention.

Genome-wide association study of copy number variation with lung function identifies a novel signal of association near BANP for forced vital capacity

BACKGROUND

Genome-wide association studies of Single Nucleotide Polymorphisms (SNPs) have identified 55 SNPs associated with lung function. However, little is known about the effect of copy number variants (CNVs) on lung function, although CNVs represent a significant proportion of human genetic polymorphism. To assess the effect of CNVs on lung function quantitative traits, we measured copy number at 2788 previously characterised, common copy number variable regions in 6 independent cohorts (n = 24,237) using intensity data from SNP genotyping experiments. We developed a pipeline for genome-wide association analysis and meta-analysis of CNV genotypes measured across multiple studies using SNP genotype array intensity data from different platform technologies. We then undertook cohort-level genome-wide association studies of CNV with lung function in a subset of 4 cohorts (n < =12,403) with lung function measurements and meta-analysed the results. Follow-up was undertaken for CNVs which were well tagged by SNPs, in up to 146,871 individuals.

RESULTS

We generated robust copy number calls for 1962 out of 2788 (70 %) known CNV regions genome-wide, with 1103 measured with compatible class frequencies in at least 2 cohorts. We report a novel CNV association (discovery P = 0.0007) with Forced Vital Capacity (FVC) downstream of BANP on chromosome 16 that shows evidence of replication by a tag SNP in two independent studies (replication P = 0.004). In addition, we provide suggestive evidence (discovery P = 0.0002) for a role of complex copy number variation at a previously reported lung function locus, containing the rootletin gene CROCC, that is not tagged by SNPs.

CONCLUSIONS

We demonstrate how common CNV regions can be reliably and consistently called across cohorts, using an existing calling algorithm and rigorous quality control steps, using SNP genotyping array intensity data. Although many common biallelic CNV regions were well-tagged by common SNPs, we also identified associations with untagged mulitallelic CNV regions thereby illustrating the potential of our approach to identify some of the missing heritability of complex traits.

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

HTR4 gene structure and altered expression in the developing lung

BACKGROUND

Meta-analyses of genome-wide association studies (GWAS) have identified single nucleotide polymorphisms (SNPs) spanning the 5-hydroxytryptamine receptor 4 (5-HT₄R) gene (HTR4) associated with lung function. The aims of this study were to i) investigate the expression profile of HTR4 in adult and fetal lung tissue and cultured airway cells, ii) further define HTR4 gene structure and iii) explore the potential functional implications of key SNPs using a bioinformatic approach.

METHODS

Following reverse transcription (RT)-PCR in human brain, 5' rapid amplification of cDNA ends (5' RACE) was used to examine the exonic structure of HTR4 at the 5' end. Quantitative (Q)-PCR was used to quantify HTR4 mRNA expression in total RNA from cultured airway cells and whole lung tissue. Publically available gene microarray data on fetal samples of estimated gestational age 7-22 weeks were mined for HTR4 expression. Immunohistochemistry (IHC; in adult and fetal lung tissue) and a radioligand binding assay (in cultured airway cells) were used to analyze 5-HT₄R protein expression.

RESULTS

IHC in adult lung, irrespective of the presence of chronic obstructive pulmonary disease (COPD), suggested low level expression of 5-HT₄R protein, which was most prominent in alveolar pneumocytes. There was evidence of differential 5-HT₄R protein levels during gestation in fetal lung, which was also evident in gene expression microarray data. HTR4 mRNA expression, assessed by Q-PCR, was <0.5% relative to brain in total adult lung tissue and in human airway smooth muscle (HASM) and bronchial epithelial cells (HBEC) derived from adult donors. Radioligand binding experiments also indicated that HBEC and HASM cells did not express a significant 5-HT₄R population. 5' RACE in brain identified a novel N-terminal variant, containing an extended N-terminal sequence. The functional significance of key HTR4 SNPs was investigated using the encyclopedia of DNA elements consortium (ENCODE) dataset. These analyses identified multiple alterations in regulatory motifs for transcription factors implicated in lung development, including Foxp1.

CONCLUSIONS

Taken together, these data suggest a role for HTR4 in lung development, which may at least in part explain the genetic association with lung function.

Functional IL6R 358Ala allele impairs classical IL-6 receptor signaling and influences risk of diverse inflammatory diseases

Inflammation, which is directly regulated by interleukin-6 (IL-6) signaling, is implicated in the etiology of several chronic diseases. Although a common, non-synonymous variant in the IL-6 receptor gene (IL6R Asp358Ala; rs2228145 A>C) is associated with the risk of several common diseases, with the 358Ala allele conferring protection from coronary heart disease (CHD), rheumatoid arthritis (RA), atrial fibrillation (AF), abdominal aortic aneurysm (AAA), and increased susceptibility to asthma, the variant's effect on IL-6 signaling is not known. Here we provide evidence for the association of this non-synonymous variant with the risk of type 1 diabetes (T1D) in two independent populations and confirm that rs2228145 is the major determinant of the concentration of circulating soluble IL-6R (sIL-6R) levels (34.6% increase in sIL-6R per copy of the minor allele 358Ala; rs2228145 [C]). To further investigate the molecular mechanism of this variant, we analyzed expression of IL-6R in peripheral blood mononuclear cells (PBMCs) in 128 volunteers from the Cambridge BioResource. We demonstrate that, although 358Ala increases transcription of the soluble IL6R isoform (P = 8.3×10⁻²²) and not the membrane-bound isoform, 358Ala reduces surface expression of IL-6R on CD4+ T cells and monocytes (up to 28% reduction per allele; P≤5.6×10⁻²²). Importantly, reduced expression of membrane-bound IL-6R resulted in impaired IL-6 responsiveness, as measured by decreased phosphorylation of the transcription factors STAT3 and STAT1 following stimulation with IL-6 (P≤5.2×10⁻⁷). Our findings elucidate the regulation of IL-6 signaling by IL-6R, which is causally relevant to several complex diseases, identify mechanisms for new approaches to target the IL-6/IL-6R axis, and anticipate differences in treatment response to IL-6 therapies based on this common IL6R variant.

Interleukin-6 (IL-6) is a complex cytokine, which plays a critical role in the regulation of inflammatory responses. Genetic variation in the IL-6 receptor gene is associated with the risk of several human diseases with an inflammatory component, including coronary heart disease, rheumatoid arthritis, and asthma. A common non-synonymous single nucleotide polymorphism in this gene (Asp358Ala) has been suggested to be the causal variant in this region by affecting the circulatory concentrations of soluble IL-6R (sIL-6R). In this study we extend the genetic association of this variant to type 1 diabetes and provide evidence that this variant exerts its functional mechanism by regulating the balance between sIL-6R (generated through cleavage of the surface receptor and by alternative splicing of a soluble IL6R isoform) and membrane-bound IL-6R. These data show for the first time that the minor allele of this non-synonymous variant (Ala358) directly controls the surface levels of IL-6R on individual immune cells and that these differences in protein levels translate into a functional impairment in IL-6R signaling. These findings may have implications for clinical trials targeting inflammatory mechanisms involving IL-6R signaling and may provide tools for identifying patients with specific benefit from therapeutic intervention in the IL-6R signaling pathway.

Genes contributing to pain sensitivity in the normal population: an exome sequencing study

Sensitivity to pain varies considerably between individuals and is known to be heritable. Increased sensitivity to experimental pain is a risk factor for developing chronic pain, a common and debilitating but poorly understood symptom. To understand mechanisms underlying pain sensitivity and to search for rare gene variants (MAF<5%) influencing pain sensitivity, we explored the genetic variation in individuals' responses to experimental pain. Quantitative sensory testing to heat pain was performed in 2,500 volunteers from TwinsUK (TUK): exome sequencing to a depth of 70× was carried out on DNA from singletons at the high and low ends of the heat pain sensitivity distribution in two separate subsamples. Thus in TUK1, 101 pain-sensitive and 102 pain-insensitive were examined, while in TUK2 there were 114 and 96 individuals respectively. A combination of methods was used to test the association between rare variants and pain sensitivity, and the function of the genes identified was explored using network analysis. Using causal reasoning analysis on the genes with different patterns of SNVs by pain sensitivity status, we observed a significant enrichment of variants in genes of the angiotensin pathway (Bonferroni corrected p = 3.8×10⁻⁴). This pathway is already implicated in animal models and human studies of pain, supporting the notion that it may provide fruitful new targets in pain management. The approach of sequencing extreme exome variation in normal individuals has provided important insights into gene networks mediating pain sensitivity in humans and will be applicable to other common complex traits.

Chronic widespread pain is a complex clinical problem. Identification of underlying genetic factors would shed light on the biology of pain and offer targets for novel therapies. We aimed to identify rare genetic variants in the normal population associated with pain sensation by performing exome sequencing on individuals who were more or less sensitive to heat pain. While we did not identify any single variants having large effect, we did observe major group differences between the sensitive and insensitive individuals. Network analysis suggested a role for the angiotensin pathway, which previous work in animal models has suggested is important in pain mediation. Our results cast light on the genetic factors underlying normal pain sensation in humans and the utility of exome analyses. It suggests that further exploration of the angiotensin pathway may reveal novel targets for the treatment of pain.

Genetic variants associated with Lp(a) lipoprotein level and coronary disease

BACKGROUND

An increased level of Lp(a) lipoprotein has been identified as a risk factor for coronary artery disease that is highly heritable. The genetic determinants of the Lp(a) lipoprotein level and their relevance for the risk of coronary disease are incompletely understood.

METHODS

We used a novel gene chip containing 48,742 single-nucleotide polymorphisms (SNPs) in 2100 candidate genes to test for associations in 3145 case subjects with coronary disease and 3352 control subjects. Replication was tested in three independent populations involving 4846 additional case subjects with coronary disease and 4594 control subjects.

RESULTS

Three chromosomal regions (6q26-27, 9p21, and 1p13) were strongly associated with the risk of coronary disease. The LPA locus on 6q26-27 encoding Lp(a) lipoprotein had the strongest association. We identified a common variant (rs10455872) at the LPA locus with an odds ratio for coronary disease of 1.70 (95% confidence interval [CI], 1.49 to 1.95) and another independent variant (rs3798220) with an odds ratio of 1.92 (95% CI, 1.48 to 2.49). Both variants were strongly associated with an increased level of Lp(a) lipoprotein, a reduced copy number in LPA (which determines the number of kringle IV-type 2 repeats), and a small Lp(a) lipoprotein size. Replication studies confirmed the effects of both variants on the Lp(a) lipoprotein level and the risk of coronary disease. A meta-analysis showed that with a genotype score involving both LPA SNPs, the odds ratios for coronary disease were 1.51 (95% CI, 1.38 to 1.66) for one variant and 2.57 (95% CI, 1.80 to 3.67) for two or more variants. After adjustment for the Lp(a) lipoprotein level, the association between the LPA genotype score and the risk of coronary disease was abolished.

CONCLUSIONS

We identified two LPA variants that were strongly associated with both an increased level of Lp(a) lipoprotein and an increased risk of coronary disease. Our findings provide support for a causal role of Lp(a) lipoprotein in coronary disease.