We have evaluated the pharmacogenetic content of commercial human genome-wide genotyping arrays, as it is a critical determinant to enabling pharmacogenomic discoveries.
Using bioinformatics ...approaches, we assessed 27,811 genetic variants in 3146 genes for their presence in 18 Illumina and 15 Affymetrix genome-wide arrays.
The pharmacogenetic content of the arrays varied greatly. The combination of the Affymetrix precision medicine array and PharmacoScan arrays (Affymetrix) had the highest coverage for a set of clinically actionable absorption, distribution, metabolism and excretion (ADME) variants, single nucleotide ADME variants and ADME insertions/deletions, with a physical coverage of 125/130 (96.2%), 9924/24,138 (41.1%) and 2252/3994 (56.4%), respectively.
The combination of the Affymetrix precision medicine array and PharmacoScan arrays provided both genome-wide and pharmacogene coverage, which is crucial in the discovering of new variants responsible for drug adverse effects. These results will help in the design of pharmacogenomic studies and will enable a critical review of results from past studies.
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.
Naturally occurring human genetic variants provide a valuable tool to identify drug targets and guide drug prioritization and clinical trial design. Ivabradine is a heart rate lowering drug with ...protective effects on heart failure despite increasing the risk of atrial fibrillation. In patients with coronary artery disease without heart failure, the drug does not protect against major cardiovascular adverse events prompting questions about the ability of genetics to have predicted those effects. This study evaluates the effect of a variant in HCN4, ivabradine's drug target, on safety and efficacy endpoints.
We used genetic association testing and Mendelian randomization to predict the effect of ivabradine and heart rate lowering on cardiovascular outcomes.
Using data from the UK Biobank and large GWAS consortia, we evaluated the effect of a heart rate-reducing genetic variant at the HCN4 locus encoding ivabradine's drug target. These genetic association analyses showed increases in risk for atrial fibrillation (OR 1.09, 95% CI: 1.06-1.13, P = 9.3 ×10-9) in the UK Biobank. In a cause-specific competing risk model to account for the increased risk of atrial fibrillation, the HCN4 variant reduced incident heart failure in participants that did not develop atrial fibrillation (HR 0.90, 95% CI: 0.83-0.98, P = 0.013). In contrast, the same heart rate reducing HCN4 variant did not prevent a composite endpoint of myocardial infarction or cardiovascular death (OR 0.99, 95% CI: 0.93-1.04, P = 0.61).
Genetic modelling of ivabradine recapitulates its benefits in heart failure, promotion of atrial fibrillation, and neutral effect on myocardial infarction.
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).
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.
Hematological traits are important clinical parameters. To test the effects of rare and low-frequency coding variants on hematological traits, we analyzed hemoglobin concentration, hematocrit levels, ...white blood cell (WBC) counts and platelet counts in 31,340 individuals genotyped on an exome array. We identified several missense variants in CXCR2 associated with reduced WBC count (gene-based P = 2.6 × 10(-13)). In a separate family-based resequencing study, we identified a CXCR2 frameshift mutation in a pedigree with congenital neutropenia that abolished ligand-induced CXCR2 signal transduction and chemotaxis. We also identified missense or splice-site variants in key hematopoiesis regulators (EPO, TFR2, HBB, TUBB1 and SH2B3) associated with blood cell traits. Finally, we were able to detect associations between a rare somatic JAK2 mutation (encoding p.Val617Phe) and platelet count (P = 3.9 × 10(-22)) as well as hemoglobin concentration (P = 0.002), hematocrit levels (P = 9.5 × 10(-7)) and WBC count (P = 3.1 × 10(-5)). In conclusion, exome arrays complement genome-wide association studies in identifying new variants that contribute to complex human traits.
The advent of high throughput sequencing methods breeds an important amount of technical challenges. Among those is the one raised by the discovery of copy-number variations (CNVs) using whole-genome ...sequencing data. CNVs are genomic structural variations defined as a variation in the number of copies of a large genomic fragment, usually more than one kilobase. Here, we aim to compare different CNV calling methods in order to assess their ability to consistently identify CNVs by comparison of the calls in 9 quartets of identical twin pairs. The use of monozygotic twins provides a means of estimating the error rate of each algorithm by observing CNVs that are inconsistently called when considering the rules of Mendelian inheritance and the assumption of an identical genome between twins. The similarity between the calls from the different tools and the advantage of combining call sets were also considered.
ERDS and CNVnator obtained the best performance when considering the inherited CNV rate with a mean of 0.74 and 0.70, respectively. Venn diagrams were generated to show the agreement between the different algorithms, before and after filtering out familial inconsistencies. This filtering revealed a high number of false positives for CNVer and Breakdancer. A low overall agreement between the methods suggested a high complementarity of the different tools when calling CNVs. The breakpoint sensitivity analysis indicated that CNVnator and ERDS achieved better resolution of CNV borders than the other tools. The highest inherited CNV rate was achieved through the intersection of these two tools (81%).
This study showed that ERDS and CNVnator provide good performance on whole genome sequencing data with respect to CNV consistency across families, CNV breakpoint resolution and CNV call specificity. The intersection of the calls from the two tools would be valuable for CNV genotyping pipelines.
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
, P = 5.81 × 10
, and P = 5.94 × 10
, 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.
Genetic variants in drug targets can be used to predict the long‐term, on‐target effect of drugs. Here, we extend this principle to assess how sex and body mass index may modify the effect of ...genetically predicted lower CETP levels on biomarkers and cardiovascular outcomes. We found sex and body mass index (BMI) to be modifiers of the association between genetically predicted lower CETP and lipid biomarkers in UK Biobank participants. Female sex and lower BMI were associated with higher high‐density lipoprotein cholesterol and lower low‐density lipoprotein cholesterol for the same genetically predicted reduction in CETP concentration. We found that sex also modulated the effect of genetically lower CETP on cholesterol efflux capacity in samples from the Montreal Heart Institute Biobank. However, these modifying effects did not extend to sex differences in cardiovascular outcomes in our data. Our results provide insight into the clinical effects of CETP inhibitors in the presence of effect modification based on genetic data. The approach can support precision medicine applications and help assess the external validity of clinical trials.
Few genome-wide association studies (GWASs) have been conducted to identify predictors of drug concentrations. The authors therefore sought to discover the pharmacogenomic markers involved in ...metoprolol pharmacokinetics.
The authors performed a GWAS of a cross-sectional study of 993 patients from the Montreal Heart Institute Biobank taking metoprolol.
A total of 391 and 444 SNPs reached the significance threshold of 5 × 10
for metoprolol and α-OH-metoprolol concentrations, respectively. All were located on chromosome 22 at or near the
gene, encoding CYP450 2D6, metoprolol's main metabolizing enzyme.
The results reinforce previous findings of the importance of the
locus for metoprolol concentrations and confirm that large biobanks can be used to identify genetic determinants of drug pharmacokinetics at a GWAS significance level.
Since the beginning of the COVID-19 pandemic, many countries, including Canada, have adopted unprecedented physical distancing measures such as closure of schools and non-essential businesses, and ...restrictions on gatherings and household visits. We described time trends in social contacts for the pre-pandemic and pandemic periods in Quebec, Canada.
CONNECT is a population-based study of social contacts conducted shortly before (2018/2019) and during the COVID-19 pandemic (April 2020 - February 2021), using the same methodology for both periods. We recruited participants by random digit dialing and collected data by self-administered web-based questionnaires. Questionnaires documented socio-demographic characteristics and social contacts for two assigned days. A contact was defined as a two-way conversation at a distance ≤ 2 m or as a physical contact, irrespective of masking. We used weighted generalized linear models with a Poisson distribution and robust variance (taking possible overdispersion into account) to compare the mean number of social contacts over time and by socio-demographic characteristics.
A total of 1291 and 5516 Quebecers completed the study before and during the pandemic, respectively. Contacts significantly decreased from a mean of 8 contacts/day prior to the pandemic to 3 contacts/day during the spring 2020 lockdown. Contacts remained lower than the pre-COVID period thereafter (lowest = 3 contacts/day during the Christmas 2020/2021 holidays, highest = 5 in September 2020). Contacts at work, during leisure activities/in other locations, and at home with visitors showed the greatest decreases since the beginning of the pandemic. All sociodemographic subgroups showed significant decreases of contacts since the beginning of the pandemic. The mixing matrices illustrated the impact of public health measures (e.g. school closure, gathering restrictions) with fewer contacts between children/teenagers and fewer contacts outside of the three main diagonals of contacts between same-age partners/siblings and between children and their parents.
Physical distancing measures in Quebec significantly decreased social contacts, which most likely mitigated the spread of COVID-19.