Cluster analysis of transcriptomic datasets to identify endotypes of idiopathic pulmonary fibrosis.

BACKGROUND: Considerable clinical heterogeneity in idiopathic pulmonary fibrosis (IPF) suggests the existence of multiple disease endotypes. Identifying these endotypes would improve our understanding of the pathogenesis of IPF and could allow for a biomarker-driven personalised medicine approach. We aimed to identify clinically distinct groups of patients with IPF that could represent distinct disease endotypes. METHODS: We co-normalised, pooled and clustered three publicly available blood transcriptomic datasets (total 220 IPF cases). We compared clinical traits across clusters and used gene enrichment analysis to identify biological pathways and processes that were over-represented among the genes that were differentially expressed across clusters. A gene-based classifier was developed and validated using three additional independent datasets (total 194 IPF cases). FINDINGS: We identified three clusters of patients with IPF with statistically significant differences in lung function (p=0.009) and mortality (p=0.009) between groups. Gene enrichment analysis implicated mitochondrial homeostasis, apoptosis, cell cycle and innate and adaptive immunity in the pathogenesis underlying these groups. We developed and validated a 13-gene cluster classifier that predicted mortality in IPF (high-risk clusters vs low-risk cluster: HR 4.25, 95% CI 2.14 to 8.46, p=3.7×10-5). INTERPRETATION: We have identified blood gene expression signatures capable of discerning groups of patients with IPF with significant differences in survival. These clusters could be representative of distinct pathophysiological states, which would support the theory of multiple endotypes of IPF. Although more work must be done to confirm the existence of these endotypes, our classifier could be a useful tool in patient stratification and outcome prediction in IPF.

Genome-wide association study of chronic sputum production implicates loci involved in mucus production and infection.

BACKGROUND: Chronic sputum production impacts on quality of life and is a feature of many respiratory diseases. Identification of the genetic variants associated with chronic sputum production in a disease agnostic sample could improve understanding of its causes and identify new molecular targets for treatment. METHODS: We conducted a genome-wide association study (GWAS) of chronic sputum production in UK Biobank. Signals meeting genome-wide significance (p<5×10-8) were investigated in additional independent studies, were fine-mapped and putative causal genes identified by gene expression analysis. GWASs of respiratory traits were interrogated to identify whether the signals were driven by existing respiratory disease among the cases and variants were further investigated for wider pleiotropic effects using phenome-wide association studies (PheWASs). RESULTS: From a GWAS of 9714 cases and 48 471 controls, we identified six novel genome-wide significant signals for chronic sputum production including signals in the human leukocyte antigen (HLA) locus, chromosome 11 mucin locus (containing MUC2, MUC5AC and MUC5B) and FUT2 locus. The four common variant associations were supported by independent studies with a combined sample size of up to 2203 cases and 17 627 controls. The mucin locus signal had previously been reported for association with moderate-to-severe asthma. The HLA signal was fine-mapped to an amino acid change of threonine to arginine (frequency 36.8%) in HLA-DRB1 (HLA-DRB1*03:147). The signal near FUT2 was associated with expression of several genes including FUT2, for which the direction of effect was tissue dependent. Our PheWAS identified a wide range of associations including blood cell traits, liver biomarkers, infections, gastrointestinal and thyroid-associated diseases, and respiratory disease. CONCLUSIONS: Novel signals at the FUT2 and mucin loci suggest that mucin fucosylation may be a driver of chronic sputum production even in the absence of diagnosed respiratory disease and provide genetic support for this pathway as a target for therapeutic intervention.

Genetics plays a limited role in predicting chronic obstructive pulmonary disease treatment response and exacerbation.

BACKGROUND: Combination treatments, targeting multiple disease processes, benefit subjects with acute exacerbation of chronic obstructive pulmonary disease (AECOPD). However, predicting treatment response and exacerbation risk remain challenging. OBJECTIVE: To identify genetic associations with AECOPD risk and response to combination therapy (fluticasone furoate, umeclidinium bromide and vilanterol). METHODS: The genetic basis of AECOPD disease was investigated in 19,841 subjects from 23 clinical studies and 2 disease cohorts to identify exacerbation disease targets. AECOPD pharmacogenetic effects were examined in 8439 moderate to severe COPD patients with exacerbation rate, lung function and quality of life endpoints; results were followed up in an additional 2201 subjects. RESULTS: We did not identify significant associations in the AECOPD disease analysis. In the AECOPD pharmacogenetics analysis, rs56195836 (MAPK8) was significantly associated with moderate to severe exacerbation rate in subjects on fluticasone furoate with baseline blood eosinophils ≥150 cells/µl (P = 1.8 × 10-8). Post-hoc, one variant was associated with on-treatment moderate to severe exacerbation rate stratifying by exacerbation history. AZU1 rs1962343 was significantly associated in subjects with frequent moderate exacerbation history when treated with fluticasone furoate/vilanterol (P = 1.1 × 10-8). Neither of these signals was supported in independent follow-up. CONCLUSION: Common genetic variants do not play major roles in AECOPD disease nor predict response to triple therapy or its components in moderate to very severe COPD.

Pharmacogenetic meta-analysis of baseline risk factors, pharmacodynamic, efficacy and tolerability endpoints from two large global cardiovascular outcomes trials for darapladib.

Darapladib, a lipoprotein-associated phospholipase A2 (Lp-PLA2) inhibitor, failed to demonstrate efficacy for the primary endpoints in two large phase III cardiovascular outcomes trials, one in stable coronary heart disease patients (STABILITY) and one in acute coronary syndrome (SOLID-TIMI 52). No major safety signals were observed but tolerability issues of diarrhea and odor were common (up to 13%). We hypothesized that genetic variants associated with Lp-PLA2 activity may influence efficacy and tolerability and therefore performed a comprehensive pharmacogenetic analysis of both trials. We genotyped patients within the STABILITY and SOLID-TIMI 52 trials who provided a DNA sample and consent (n = 13,577 and 10,404 respectively, representing 86% and 82% of the trial participants) using genome-wide arrays with exome content and performed imputation using a 1000 Genomes reference panel. We investigated baseline and change from baseline in Lp-PLA2 activity, two efficacy endpoints (major coronary events and myocardial infarction) as well as tolerability parameters at genome-wide and candidate gene level using a meta-analytic approach. We replicated associations of published loci on baseline Lp-PLA2 activity (APOE, CELSR2, LPA, PLA2G7, LDLR and SCARB1) and identified three novel loci (TOMM5, FRMD5 and LPL) using the GWAS-significance threshold P≤5E-08. Review of the PLA2G7 gene (encoding Lp-PLA2) within these datasets identified V279F null allele carriers as well as three other rare exonic null alleles within various ethnic groups, however none of these variants nor any other loci associated with Lp-PLA2 activity at baseline were associated with any of the drug response endpoints. The analysis of darapladib efficacy endpoints, despite low power, identified six low frequency loci with main genotype effect (though with borderline imputation scores) and one common locus (minor allele frequency 0.24) with genotype by treatment interaction effect passing the GWAS-significance threshold. This locus conferred risk in placebo subjects, hazard ratio (HR) 1.22 with 95% confidence interval (CI) 1.11-1.33, but was protective in darapladib subjects, HR 0.79 (95% CI 0.71-0.88). No major loci for tolerability were found. Thus, genetic analysis confirmed and extended the influence of lipoprotein loci on Lp-PLA2 levels, identified some novel null alleles in the PLA2G7 gene, and only identified one potentially efficacious subgroup within these two large clinical trials.

A phenome-wide association study of a lipoprotein-associated phospholipase A2 loss-of-function variant in 90 000 Chinese adults.

BACKGROUND: Lipoprotein-associated phospholipase A2 (Lp-PLA2) has been implicated in development of atherosclerosis; however, recent randomized trials of Lp-PLA2 inhibition reported no beneficial effects on vascular diseases. In East Asians, a loss-of-function variant in the PLA2G7 gene can be used to assess the effects of genetically determined lower Lp-PLA2 METHODS: PLA2G7 V279F (rs76863441) was genotyped in 91 428 individuals randomly selected from the China Kadoorie Biobank of 0.5 M participants recruited in 2004-08 from 10 regions of China, with 7 years' follow-up. Linear regression was used to assess effects of V279F on baseline traits. Logistic regression was conducted for a range of vascular and non-vascular diseases, including 41 ICD-10 coded disease categories. RESULTS: PLA2G7 V279F frequency was 5% overall (range 3-7% by region), and 9691 (11%) participants had at least one loss-of-function variant. V279F was not associated with baseline blood pressure, adiposity, blood glucose or lung function. V279F was not associated with major vascular events [7141 events; odds ratio (OR) = 0.98 per F variant, 95% confidence interval (CI) 0.90-1.06] or other vascular outcomes, including major coronary events (922 events; 0.96, 0.79-1.18) and stroke (5967 events; 1.00, 0.92-1.09). Individuals with V279F had lower risks of diabetes (7031 events; 0.91, 0.84-0.98) and asthma (182 events; 0.53, 0.28-0.98), but there was no association after adjustment for multiple testing. CONCLUSIONS: Lifelong lower Lp-PLA2 activity was not associated with major risks of vascular or non-vascular diseases in Chinese adults. Using functional genetic variants in large-scale prospective studies with linkage to a range of health outcomes is a valuable approach to inform drug development and repositioning.

First administration of the Fc-attenuated anti-ß amyloid antibody GSK933776 to patients with mild Alzheimer's disease: a randomized, placebo-controlled study.

OBJECTIVE: To assess the safety, tolerability, pharmacokinetics, and pharmacodynamics of the Fc-inactivated anti-ß amyloid (Aß) monoclonal antibody (mAb) GSK933776 in patients with mild Alzheimer's disease (AD) or mild cognitive impairment (MCI). METHODS: This was a two-part, single blind, placebo-controlled, first-time-in-human (FTIH) study of single (n = 18) and repeat dose (n = 32) intravenous GSK933776 0.001-6 mg/kg (ClinicalTrials.gov: NCT00459550). Additional safety data from an open-label, uncontrolled, single dose study of intravenous GSK933776 1-6 mg/kg (n = 18) are included (ClinicalTrials.gov: NCT01424436). RESULTS: There were no cases of amyloid-related imaging abnormalities-edema (ARIA-E) or -hemorrhage (ARIA-H) after GSK933776 administration in both studies. Three patients across the two studies developed anti-GSK933776 antibodies. Plasma GSK933776 half-life (t1/2) was 10-15 days after repeat dosing. After each of three administrations of GSK933776, plasma levels of total Aß42 and Aß increased whereas plasma levels of free Aß decreased dose dependently; no changes were observed for placebo. For total Aß42 the peak:trough ratio was ≤2 at doses ≥3 mg/kg; for total Aß the ratio was ≤2 at 6 mg/kg. CSF concentrations of Aß showed increases from baseline to week 12 for Aß X-38 (week 12:baseline ratio: 1.65; 95%CI: 1.38, 1.93) and Aß X-42 (week 12:baseline ratio: 1.18; 95%CI: 1.06, 1.30) for values pooled across doses. CONCLUSION: In this FTIH study the Fc-inactivated anti-Aß mAb GSK933776 engaged its target in plasma and CSF without causing brain ARIA-E/H in patients with mild AD or MCI. TRIAL REGISTRATION: ClinicalTrials.gov NCT00459550.

Detection of genotyping errors by Hardy-Weinberg equilibrium testing.

Genotyping data sets may contain errors that, in some instances, lead to false conclusions. Deviation from Hardy-Weinberg equilibrium (HWE) in random samples may be indicative of problematic assays. This study has analysed 107,000 genotypes generated by TaqMan, RFLP, sequencing or mass spectrometric methods from 443 single-nucleotide polymorphisms (SNPs). These SNPs are distributed both within genes and in intergenic regions. Genotype distributions for 36 out of 313 assays (11.5%) whose minor allele frequencies were >0.05 deviated from HWE (P<0.05). Some of the possible reasons for this deviation were explored: assays for five SNPs proved nonspecific, and genotyping errors were identified in 21 SNPs. For the remaining 10 SNPs, no reasons for deviation from HWE were identified. We demonstrate the successful identification of a proportion of nonspecific assays, and assays harbouring genotyping error. Consequently, our current high-throughput genotyping system incorporates tests for both assay specificity and deviation from HWE, to minimise the genotype error rate and therefore improve data quality.

Evaluation of the effect of UGT1A1 polymorphisms on dolutegravir pharmacokinetics.

AIM: To evaluate potential pharmacogenetic effects of UGT1A1 polymorphisms on the pharmacokinetics (PK) of dolutegravir (Tivicay®; ViiV Healthcare, NC, USA), an HIV-1 integrase inhibitor. PATIENTS & METHODS: Analysis of pooled data from nine Phase I and II clinical studies was undertaken for 89 subjects receiving repeat dolutegravir 50 mg once daily (tablet formulation) who were genotyped for known UGT1A1 functional variants. RESULTS: Geometric mean ratio (92% CI) for subjects carrying low (*28/*28 and *28/*37) and reduced activity (*1/*6, *1/*28, *1/*37, *28/*36 and *36/*37) polymorphisms compared with subjects with normal activity (*1/*1 and *1/*36) showed decreased oral clearance (CL/F; 0.765 [92% CI: 0.659-0.889]), increased area under the concentration-time curve (AUC(0-τ); 1.307 [1.125-1.518]) and C(max) (1.221 [1.063-1.402]), respectively. CONCLUSION: Increased dolutegravir exposure in carriers of UGT1A1 reduced function polymorphisms is not clinically significant based on accumulated safety data so dose adjustment in these individuals is not required.

Modulation of ß-amyloid by a single dose of GSK933776 in patients with mild Alzheimer's disease: a phase I study.

INTRODUCTION: In this study, we evaluated the safety and pharmacodynamic effects of the Fc-inactivated anti-ß-amyloid (anti-Aß) monoclonal antibody GSK933776 in patients with mild Alzheimer's disease and mild cognitive impairment. Aß and tau levels were investigated in cerebrospinal fluid (CSF), and the relationship between Aß levels and Aß modulation in plasma was explored. The feasibility of a continuous sampling method using a lumbar catheter was assessed. METHODS: This trial was a phase I, open-label, uncontrolled, single-dose, exploratory experimental medicine study of intravenous GSK933776 at doses of 1 mg/kg, 3 mg/kg or 6 mg/kg (n = 6/group). The time course of plasma and CSF concentrations of GSK933776 and Aß was assessed. Sample size was based on feasibility, and no formal statistical analyses were performed. RESULTS: Following administration of GSK933776 at doses of 1 mg/kg, 3 mg/kg and 6 mg/kg, there were decreases from baseline in CSF Aß1-42 (from 0 to 12 hours) by 22.8 pg/ml (6.2%), 43.5 pg/ml (9.2%) and 60.5 pg/ml (12.5%), respectively. Plasma concentrations of total Aß18-35 and Aß4228-42 increased immediately after infusion and CSF tau concentration increased slightly, but did not significantly change, following administration of all doses of GSK933776. Pharmacokinetics confirmed the presence of GSK933776 in the CSF, which exhibited a dose-response relationship. One patient underwent minor surgery without sequelae following a ruptured lumbar catheter. CONCLUSION: GSK933776 demonstrated pharmacological activity and target engagement in CSF and plasma, and the continuous sampling method via a catheter successfully assessed the Aß changes following single-dose administration of GSK933776. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01424436. Registered 4 August 2011.

An 18-kDa translocator protein (TSPO) polymorphism explains differences in binding affinity of the PET radioligand PBR28.

[(11)C]PBR28 binds the 18-kDa Translocator Protein (TSPO) and is used in positron emission tomography (PET) to detect microglial activation. However, quantitative interpretations of signal are confounded by large interindividual variability in binding affinity, which displays a trimodal distribution compatible with a codominant genetic trait. Here, we tested directly for an underlying genetic mechanism to explain this. Binding affinity of PBR28 was measured in platelets isolated from 41 human subjects and tested for association with polymorphisms in TSPO and genes encoding other proteins in the TSPO complex. Complete agreement was observed between the TSPO Ala147Thr genotype and PBR28 binding affinity phenotype (P value=3.1 × 10(-13)). The TSPO Ala147Thr polymorphism predicts PBR28 binding affinity in human platelets. As all second-generation TSPO PET radioligands tested hitherto display a trimodal distribution in binding affinity analogous to PBR28, testing for this polymorphism may allow quantitative interpretation of TSPO PET studies with these radioligands.