Pharmacogenomic study of heart failure and candesartan response from the CHARM programme.

AIMS: The Candesartan in Heart failure Assessment of Reduction in Mortality and morbidity (CHARM) programme consisted of three parallel, randomized, double-blind clinical trials comparing candesartan with placebo in patients with heart failure (HF) categorized according to left ventricular ejection fraction and tolerability to an angiotensin-converting enzyme inhibitor. We conducted a pharmacogenomic study of the CHARM trials with the objective of identifying genetic predictors of HF progression and of the efficacy and safety of treatment with candesartan. METHODS: We performed genome-wide association studies in 2727 patients of European ancestry from CHARM-Overall and stratified by CHARM study according to preserved and reduced ejection fraction and according to assignment to the interventional treatment with candesartan. We tested genetic association with the composite endpoint of cardiovascular death or hospitalization for heart failure for drug efficacy in candesartan-treated patients and for HF progression using patients from both candesartan and placebo arms. The safety endpoints for response to candesartan were hyperkalaemia, renal dysfunction, hypotension, and change in systolic blood pressure between baseline and 6 weeks of treatment. To support our observations, we conducted a genome-wide gene-level collapsing analysis from whole-exome sequencing data with the composite cardiovascular endpoint. RESULTS: We found that the A allele (14% allele frequency) of the genetic variant rs66886237 at 8p21.3 near the gene GFRA2 was associated with the composite cardiovascular endpoint in 1029 HF patients with preserved ejection fraction from the CHARM-Preserved study (hazard ratio: 1.91, 95% confidence interval: 1.55-2.35; P = 1.7 × 10-9 ). The association was independent of candesartan treatment, and the genetic variant was not associated with the cardiovascular endpoint in patients with reduced ejection fraction. None of the genome-wide association studies for candesartan safety or efficacy conducted in patients treated with candesartan passed the significance threshold. We found no significant association from the gene-level collapsing analysis. CONCLUSIONS: We have identified a candidate genetic variant potentially predictive of the progression of heart failure in patients with preserved ejection fraction. The findings require further replication, and we cannot exclude the possibility that the results may be chance findings.

The genomics of heart failure: design and rationale of the HERMES consortium.

AIMS: The HERMES (HEart failure Molecular Epidemiology for Therapeutic targetS) consortium aims to identify the genomic and molecular basis of heart failure. METHODS AND RESULTS: The consortium currently includes 51 studies from 11 countries, including 68 157 heart failure cases and 949 888 controls, with data on heart failure events and prognosis. All studies collected biological samples and performed genome-wide genotyping of common genetic variants. The enrolment of subjects into participating studies ranged from 1948 to the present day, and the median follow-up following heart failure diagnosis ranged from 2 to 116 months. Forty-nine of 51 individual studies enrolled participants of both sexes; in these studies, participants with heart failure were predominantly male (34-90%). The mean age at diagnosis or ascertainment across all studies ranged from 54 to 84 years. Based on the aggregate sample, we estimated 80% power to genetic variant associations with risk of heart failure with an odds ratio of ≥1.10 for common variants (allele frequency ≥ 0.05) and ≥1.20 for low-frequency variants (allele frequency 0.01-0.05) at P < 5 × 10-8 under an additive genetic model. CONCLUSIONS: HERMES is a global collaboration aiming to (i) identify the genetic determinants of heart failure; (ii) generate insights into the causal pathways leading to heart failure and enable genetic approaches to target prioritization; and (iii) develop genomic tools for disease stratification and risk prediction.

Genetics of symptom remission in outpatients with COVID-19.

We conducted a genome-wide association study of time to remission of COVID-19 symptoms in 1723 outpatients with at least one risk factor for disease severity from the COLCORONA clinical trial. We found a significant association at 5p13.3 (rs1173773; P = 4.94 × 10-8) near the natriuretic peptide receptor 3 gene (NPR3). By day 15 of the study, 44%, 54% and 59% of participants with 0, 1, or 2 copies of the effect allele respectively, had symptom remission. In 851 participants not treated with colchicine (placebo), there was a significant association at 9q33.1 (rs62575331; P = 2.95 × 10-8) in interaction with colchicine (P = 1.19 × 10-5) without impact on risk of hospitalisations, highlighting a possibly shared mechanistic pathway. By day 15 of the study, 46%, 62% and 64% of those with 0, 1, or 2 copies of the effect allele respectively, had symptom remission. The findings need to be replicated and could contribute to the biological understanding of COVID-19 symptom remission.

Genetic meta-analysis of cancer diagnosis following statin use identifies new associations and implicates human leukocyte antigen (HLA) in women.

We sought to perform a genomic evaluation of the risk of incident cancer in statin users, free of cancer at study entry. Patients who previously participated in two phase IV trials (TNT and IDEAL) with genetic data were used (npooled = 11,196). A GWAS meta-analysis using Cox modeling for the prediction of incident cancer was conducted in the pooled cohort and sex-stratified. rs13210472 (near HLA-DOA gene) was associated with higher risk of incident cancer amongst women with prevalent coronary artery disease (CAD) taking statins (hazard ratio [HR]: 2.66, 95% confidence interval [CI]: 1.88-3.76, P = 3.5 × 10-8). Using the UK Biobank and focusing exclusively on women statin users with CAD (nfemale = 2952), rs13210472 remained significantly associated with incident cancer (HR: 1.71, 95% CI: 1.14-2.56, P = 9.0 × 10-3). The association was not observed in non-statin users. In this genetic meta-analysis, we have identified a variant in women statin users with prevalent CAD that was associated with incident cancer, possibly implicating the human leukocyte antigen pathway.

Pharmacogenomics of the Efficacy and Safety of Colchicine in COLCOT.

BACKGROUND: The randomized, placebo-controlled COLCOT (Colchicine Cardiovascular Outcomes Trial) has shown the benefits of colchicine 0.5 mg daily to lower the rate of ischemic cardiovascular events in patients with a recent myocardial infarction. Here, we conducted a post hoc pharmacogenomic study of COLCOT with the aim to identify genetic predictors of the efficacy and safety of treatment with colchicine. METHODS: There were 1522 participants of European ancestry from the COLCOT trial available for the pharmacogenomic study of COLCOT trial. The pharmacogenomic study's primary cardiovascular end point was defined as for the main trial, as time to first occurrence of cardiovascular death, resuscitated cardiac arrest, myocardial infarction, stroke, or urgent hospitalization for angina requiring coronary revascularization. The safety end point was time to the first report of gastrointestinal events. Patients' DNA was genotyped using the Illumina Global Screening array followed by imputation. We performed a genome-wide association study in colchicine-treated patients. RESULTS: None of the genetic variants passed the genome-wide association study significance threshold for the primary cardiovascular end point conducted in 702 patients in the colchicine arm who were compliant to medication. The genome-wide association study for gastrointestinal events was conducted in all 767 patients in the colchicine arm and found 2 significant association signals, one with lead variant rs6916345 (hazard ratio, 1.89 [95% CI, 1.52-2.35], P=7.41×10-9) in a locus which colocalizes with Crohn disease, and one with lead variant rs74795203 (hazard ratio, 2.51 [95% CI, 1.82-3.47]; P=2.70×10-8), an intronic variant in gene SEPHS1. The interaction terms between the genetic variants and treatment with colchicine versus placebo were significant. CONCLUSIONS: We found 2 genomic regions associated with gastrointestinal events in patients treated with colchicine. Those findings will benefit from replication to confirm that some patients may have genetic predispositions to lower tolerability of treatment with colchicine.

A genetic association study of heart failure: more evidence for the role of BAG3 in idiopathic dilated cardiomyopathy.

AIMS: Few investigations have been conducted to identify genetic determinants of common, polygenetic forms of heart failure (HF), and only a limited number of these genetic associations have been validated by multiple groups. METHODS AND RESULTS: We performed a case-control study to further investigate the potential impact of 14 previously reported candidate genes on the risk of HF and specific HF sub-types. We also performed an exploratory genome-wide study. We included 799 patients with HF and 1529 controls. After adjusting for age, sex, and genetic ancestry, we found that the C allele of rs2234962 in BAG3 was associated with a decreased risk of idiopathic dilated cardiomyopathy (odds ratio 0.42, 95% confidence interval 0.25-0.68, P = 0.0005), consistent with a previous report. No association for the other primary variants or exploratory genome-wide study was found. CONCLUSIONS: Our findings provide independent replication for the association between a common coding variant (rs2234962) in BAG3 and the risk of idiopathic dilated cardiomyopathy.

Nuclear receptor gene polymorphisms and warfarin dose requirements in the Quebec Warfarin Cohort.

Warfarin is primarily metabolized by cytochrome 2C9, encoded by gene CYP2C9. Here, we investigated whether variants in nuclear receptor genes which regulate the expression of CYP2C9 are associated with warfarin response. We used data from 906 warfarin users from the Quebec Warfarin Cohort (QWC) and tested the association of warfarin dose requirement at 3 months following the initiation of therapy in nine nuclear receptor genes: NR1I3, NR1I2, NR3C1, ESR1, GATA4, RXRA, VDR, CEBPA, and HNF4A. Three correlated SNPs in the VDR gene (rs4760658, rs11168292, and rs11168293) were associated with dose requirements of warfarin (P = 2.68 × 10-5, P = 5.81 × 10-4, and P = 5.94 × 10-4, respectively). Required doses of warfarin were the highest for homozygotes of the minor allele at the VDR variants (P < 0.0026). Variants in the VDR gene were associated with the variability in response to warfarin, emphasizing the possible clinical relevance of nuclear receptor gene variants on the inter-individual variability in drug metabolism.

Identification of the genetic determinants responsible for retinal degeneration in families of Mexican descent.

PURPOSE: To investigate the clinical characteristics and genetic basis of inherited retinal degeneration (IRD) in six unrelated pedigrees from Mexico. METHODS: A complete ophthalmic evaluation including measurement of visual acuities, Goldman kinetic or Humphrey dynamic perimetry, Amsler test, fundus photography, and color vision testing was performed. Family history and blood samples were collected from available family members. DNA from members of two pedigrees was examined for known mutations using the APEX ARRP genotyping microarray and one pedigree using the APEX LCA genotyping microarray. The remaining three pedigrees were analyzed using a custom-designed targeted capture array covering the exons of 233 known retinal degeneration genes. Sequencing was performed on Illumina HiSeq. Reads were mapped against hg19, and variants were annotated using GATK and filtered by exomeSuite. Segregation and ethnicity-matched control sample analyses were performed by dideoxy sequencing. RESULTS: Six pedigrees with IRD were analyzed. Nine rare or novel, potentially pathogenic variants segregating with the phenotype were detected in IMPDH1, USH2A, RPE65, ABCA4, and FAM161A genes. Among these, six were known mutations while the remaining three changes in USH2A, RPE65, and FAM161A genes have not been previously reported to be associated with IRD. Analysis of 100 ethnicity-matched controls did not detect the presence of these three novel variants indicating, these are rare variants in the Mexican population. CONCLUSIONS: Screening patients diagnosed with IRD from Mexico identified six known mutations and three rare or novel potentially damaging variants in IMPDH1, USH2A, RPE65, ABCA4, and FAM161A genes that segregated with disease.

genipe: an automated genome-wide imputation pipeline with automatic reporting and statistical tools.

Genotype imputation is now commonly performed following genome-wide genotyping experiments. Imputation increases the density of analyzed genotypes in the dataset, enabling fine-mapping across the genome. However, the process of imputation using the most recent publicly available reference datasets can require considerable computation power and the management of hundreds of large intermediate files. We have developed genipe, a complete genome-wide imputation pipeline which includes automatic reporting, imputed data indexing and management, and a suite of statistical tests for imputed data commonly used in genetic epidemiology (Sequence Kernel Association Test, Cox proportional hazards for survival analysis, and linear mixed models for repeated measurements in longitudinal studies). AVAILABILITY AND IMPLEMENTATION: The genipe package is an open source Python software and is freely available for non-commercial use (CC BY-NC 4.0) at https://github.com/pgxcentre/genipe Documentation and tutorials are available at http://pgxcentre.github.io/genipe CONTACT: louis-philippe.lemieux.perreault@statgen.org or marie-pierre.dube@statgen.orgSupplementary information: Supplementary data are available at Bioinformatics online.

Pharmacogenomic determinants of the cardiovascular effects of dalcetrapib.

BACKGROUND: Dalcetrapib did not improve clinical outcomes, despite increasing high-density lipoprotein cholesterol by 30%. These results differ from other evidence supporting high-density lipoprotein as a therapeutic target. Responses to dalcetrapib may vary according to patients' genetic profile. METHODS AND RESULTS: We conducted a pharmacogenomic evaluation using a genome-wide approach in the dal-OUTCOMES study (discovery cohort, n=5749) and a targeted genotyping panel in the dal-PLAQUE-2 imaging trial (support cohort, n=386). The primary endpoint for the discovery cohort was a composite of cardiovascular events. The change from baseline in carotid intima-media thickness on ultrasonography at 6 and 12 months was evaluated as supporting evidence. A single-nucleotide polymorphism was found to be associated with cardiovascular events in the dalcetrapib arm, identifying the ADCY9 gene on chromosome 16 (rs1967309; P=2.41×10(-8)), with 8 polymorphisms providing P<10(-6) in this gene. Considering patients with genotype AA at rs1967309, there was a 39% reduction in the composite cardiovascular endpoint with dalcetrapib compared with placebo (hazard ratio, 0.61; 95% confidence interval, 0.41-0.92). In patients with genotype GG, there was a 27% increase in events with dalcetrapib versus placebo. Ten single-nucleotide polymorphism in the ADCY9 gene, the majority in linkage disequilibrium with rs1967309, were associated with the effect of dalcetrapib on intima-media thickness (P<0.05). Marker rs2238448 in ADCY9, in linkage disequilibrium with rs1967309 (r(2)=0.8), was associated with both the effects of dalcetrapib on intima-media thickness in dal-PLAQUE-2 (P=0.009) and events in dal-OUTCOMES (P=8.88×10(-8); hazard ratio, 0.67; 95% confidence interval, 0.58-0.78). CONCLUSIONS: The effects of dalcetrapib on atherosclerotic outcomes are determined by correlated polymorphisms in the ADCY9 gene. CLINICAL TRIAL INFORMATION: URL: http://www.clinicaltrials.gov. Unique identifiers: NCT00658515 and NCT01059682.

CKM and LILRB5 are associated with serum levels of creatine kinase.

BACKGROUND: Statins (HMG-CoA reductase inhibitors) are the most prescribed class of lipid-lowering drugs for the treatment and prevention of cardiovascular disease. Creatine kinase (CK) is a commonly used biomarker to assist in the diagnosis of statin-induced myotoxicity but the normal range of CK concentrations is wide, which limits its use as a diagnostic biomarker. METHODS AND RESULTS: We conducted a genome-wide association study of serum CK levels in 3412 statin users. Patients were recruited in Quebec, Canada, and genotyped on Illumina Human610-Quad and an iSelect panel enriched for lipid homeostasis, hypertension, and drug metabolism genes. We found a strong association signal between serum levels of CK and the muscle CK (CKM) gene (rs11559024: P=3.69×10(-16); R(2)=0.02) and with the leukocyte immunoglobulin-like receptor subfamily B member 5 (LILRB5) gene (rs2361797: P=1.96×10(-10); R(2)=0.01). Genetic variants in those 2 genes were independently associated with CK levels in statin users. Results were successfully replicated in 5330 participants from the Montreal Heart Institute Biobank in statin users for CKM (rs11559024: P=4.32×10(-16); R(2)=0.02) and LILRB5 (rs12975366 P=4.45×10(-10); R(2)=0.01) and statin nonusers (P=4.08×10(-7), R(2)=0.01; P=3.17×10(-9), R(2)=0.02, respectively). CONCLUSIONS: This is the first genome-wide study to report on the underlying genetic determinants of CK variation in a population of statin users. We found statistically significant association for variants in the CKM and LILRB5 genes.

Mutations in a novel serine protease PRSS56 in families with nanophthalmos.

PURPOSE: Nanophthalmos is a rare genetic ocular disorder in which the eyes of affected individuals are abnormally small. Patients suffer from severe hyperopia as a result of their markedly reduced axial lengths, but otherwise are capable of seeing well unlike other more general forms of microphthalmia. To date one gene for nanophthalmos has been identified, encoding the membrane-type frizzled related protein MFRP. Identification of additional genes for nanophthalmos will improve our understanding of normal developmental regulation of eye growth. METHODS: We ascertained a cohort of families from eastern Canada and Mexico with familial nanophthalmos. We performed high density microsatellite and high density single nucleotide polymorphism (SNP) genotyping to identify potential chromosomal regions of linkage. We sequenced coding regions of genes in the linked interval by traditional PCR-based Sanger capillary electrophoresis methods. We cloned and sequenced a novel cDNA from a putative causal gene to verify gene structure. RESULTS: We identified a linked locus on chromosome 2q37 with a peak logarithm (base 10) of odds (LOD) score of 4.7. Sequencing of coding exons of all genes in the region identified multiple segregating variants in one gene, recently annotated as serine protease gene (PRSS56), coding for a predicted trypsin serine protease-like protein. One of our families was homozygous for a predicted pathogenic missense mutation, one family was compound heterozygous for two predicted pathogenic missense mutations, and one family was compound heterozygous for a predicted pathogenic missense mutation plus a frameshift leading to obligatory truncation of the predicted protein. The PRSS56 gene structure in public databases is based on a virtual transcript assembled from overlapping incomplete cDNA clones; we have now validated the structure of a full-length transcript from embryonic mouse brain RNA. CONCLUSIONS: PRSS56 is a good candidate for the causal gene for nanophthalmos in our families.

Familial ventricular aneurysms and septal defects map to chromosome 10p15.

AIMS: Although ventricular septal defects (VSD) are the most common congenital heart lesion, familial clustering has been described only in rare instances. The aim of this study was to identify genetic factors and chromosomal regions contributing to VSD. METHODS AND RESULTS: A unique, large kindred segregating various forms of septal pathologies-including VSD, ventricular septal aneurysms, and atrial septal defects (ASD)-was ascertained and characterized clinically and genetically. Eighteen family members in three generations could be studied, out of whom 10 are affected (2 ASD, 3 septal aneurysm, 4 VSD, and 1 tetralogy of Fallot). Parametric multipoint LOD scores reach significance on chromosome 10p15.3-10p15.2 (max. 3.29). The LOD score support interval is in a gene-poor region where deletions have been reported to associate with septal defects, but that is distinct from the DiGeorge syndrome 2 region on 10p. Multiple linkage analysis scenarios suggest that tetralogy of Fallot is a phenocopy and genetically distinct from the autosomal dominant form of septal pathologies observed in this family. CONCLUSION: This study maps a rare familial form of VSD/septal aneurysms to chromosome 10p15 and extends the spectrum of the genetic heterogeneity of septal pathologies. Fine mapping, haplotype construction, and resequencing will provide a unique opportunity to study the pathogenesis of septal defects and shed light on molecular mechanisms of septal development.

Partitioning of copy-number genotypes in pedigrees.

BACKGROUND: Copy number variations (CNVs) and polymorphisms (CNPs) have only recently gained the genetic community's attention. Conservative estimates have shown that CNVs and CNPs might affect more than 10% of the genome and that they may be at least as important as single nucleotide polymorphisms in assessing human variability. Widely used tools for CNP analysis have been implemented in Birdsuite and PLINK for the purpose of conducting genetic association studies based on the unpartitioned total number of CNP copies provided by the intensities from Affymetrix's Genome-Wide Human SNP Array. Here, we are interested in partitioning copy number variations and polymorphisms in extended pedigrees for the purpose of linkage analysis on familial data. RESULTS: We have developed CNGen, a new software for the partitioning of copy number polymorphism using the integrated genotypes from Birdsuite with the Affymetrix platform. The algorithm applied to familial trios or extended pedigrees can produce partitioned copy number genotypes with distinct parental alleles. We have validated the algorithm using simulations on a complex pedigree structure using frequencies calculated from a real dataset of 300 genotyped samples from 42 pedigrees segregating a congenital heart defect phenotype. CONCLUSIONS: CNGen is the first published software for the partitioning of copy number genotypes in pedigrees, making possible the use CNPs and CNVs for linkage analysis. It was implemented with the Python interpreter version 2.5.2. It was successfully tested on current Linux, Windows and Mac OS workstations.

Mutations in DCC cause congenital mirror movements.

Mirror movements are involuntary contralateral movements that mirror voluntary ones and are often associated with defects in midline crossing of the developing central nervous system. We studied two large families, one French Canadian and one Iranian, in which isolated congenital mirror movements were inherited as an autosomal dominant trait. We found that affected individuals carried protein-truncating mutations in DCC (deleted in colorectal carcinoma), a gene on chromosome 18q21.2 that encodes a receptor for netrin-1, a diffusible protein that helps guide axon growth across the midline. Functional analysis of the mutant DCC protein from the French Canadian family revealed a defect in netrin-1 binding. Thus, DCC has an important role in lateralization of the human nervous system.

Application of homozygosity haplotype analysis to genetic mapping with high-density SNP genotype data.

BACKGROUND: In families segregating a monogenic genetic disorder with a single disease gene introduction, patients share a mutation-carrying chromosomal interval with identity-by-descent (IBD). Such a shared chromosomal interval or haplotype, surrounding the actual pathogenic mutation, is typically detected and defined by multipoint linkage and phased haplotype analysis using microsatellite or SNP genotype data. High-density SNP genotype data presents a computational challenge for conventional genetic analyses. A novel non-parametric method termed Homozygosity Haplotype (HH) was recently proposed for the genome-wide search of the autosomal segments shared among patients using high density SNP genotype data. METHODOLOGY/PRINCIPAL FINDINGS: The applicability and the effectiveness of HH in identifying the potential linkage of disease causative gene with high-density SNP genotype data were studied with a series of monogenic disorders ascertained in eastern Canadian populations. The HH approach was validated using the genotypes of patients from a family affected with a rare autosomal dominant disease Schnyder crystalline corneal dystrophy. HH accurately detected the approximately 1 Mb genomic interval encompassing the causative gene UBIAD1 using the genotypes of only four affected subjects. The successful application of HH to identify the potential linkage for a family with pericentral retinal disorder indicates that HH can be applied to perform family-based association analysis by treating affected and unaffected family members as cases and controls respectively. A new strategy for the genome-wide screening of known causative genes or loci with HH was proposed, as shown the applications to a myoclonus dystonia and a renal failure cohort. CONCLUSIONS/SIGNIFICANCE: Our study of the HH approach demonstrates that HH is very efficient and effective in identifying potential disease linked region. HH has the potential to be used as an efficient alternative approach to sequencing or microsatellite-based fine mapping for screening the known causative genes in genetic disease study.

Autosomal-dominant locus for Restless Legs Syndrome in French-Canadians on chromosome 16p12.1.

We describe an autosomal-dominant locus for Restless Legs Syndrome (RLS) in a French-Canadian (FC) pedigree. Genome-wide microsatellite scan and linkage analysis were used in this study. The locus maps to chromosome 16p12.1 and spans 1.18 Mega bases. The maximum multipoint LOD scores are of 3.5 over the total of 10 markers. Evidence for the same locus was also found in a smaller FC pedigree sime095. The analysis of the sequence of 8 annotated genes within the region did not reveal any pathogenic mutations. Copy number variation and karyotype analyses did not reveal any chromosomal abnormality in the region. Further analyses of the region are necessary to find the genetic cause of RLS in this family.