A dataset of proteomic changes during human heat stress and heat acclimation.

Hotter climates have important impacts on human health and performance. Yet, the cellular and molecular responses involved in human heat stress and acclimation remain understudied. This dataset includes physiological measurements and the plasma concentration of 2,938 proteins collected from 10 healthy adults, before and during passive heat stress that was performed both prior to and after a 7-day heat acclimation protocol. Physiological measurements included body temperatures, sweat rate, cutaneous vascular conductance, blood pressure, and skin sympathetic nerve activity. The proteomic dataset was generated using the Olink Explore 3072 assay, enabling a high-multiplex antibody-based assessment of protein changes based on proximity extension assay technology. The data need to be interpreted in the context of the moderate level of body hyperthermia attained and the specific demographic of young, healthy adults. We have made this dataset publicly available to facilitate research into the cellular and molecular mechanisms involved in human heat stress and acclimation, crucial for addressing the health and performance challenges posed by rising temperatures.

Impact of amiodarone use on metoprolol concentrations, α-OH-metoprolol concentrations, metoprolol dosing and heart rate: A cross-sectional study.

Small studies suggest that amiodarone is a weak inhibitor of cytochrome P450 (CYP) 2D6. Inhibition of CYP2D6 leads to increases in concentrations of drugs metabolized by the enzyme, such as metoprolol. Considering that both metoprolol and amiodarone have ß-adrenergic blocking properties and that the modest interaction between the two drugs would result in increased metoprolol concentrations, this could lead to a higher risk of bradycardia and atrioventricular block. The primary objective of this study was to evaluate whether metoprolol plasma concentrations collected at random timepoints from patients enrolled in the Montreal Heart Institute Hospital Cohort could be useful in identifying the modest pharmacokinetic interaction between amiodarone and metoprolol. We performed an analysis of a cross-sectional study, conducted as part of the Montreal Heart Institute Hospital Cohort. All participants were self-described "White" adults with metoprolol being a part of their daily pharmacotherapy regimen. Of the 999 patients being treated with metoprolol, 36 were also taking amiodarone. Amiodarone use was associated with higher metoprolol concentrations following adjustment for different covariates (p = .0132). Consistently, the association between amiodarone use and lower heart rate was apparent and significant after adjustment for all covariates under study (p = .0001). Our results highlight that single randomly collected blood samples can be leveraged to detect modest pharmacokinetic interactions.

Pharmacogenomic markers of metoprolol and α-OH-metoprolol concentrations: a genome-wide association study.

Aim: 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. Patients & methods: The authors performed a GWAS of a cross-sectional study of 993 patients from the Montreal Heart Institute Biobank taking metoprolol. Results: A total of 391 and 444 SNPs reached the significance threshold of 5 × 10-8 for metoprolol and α-OH-metoprolol concentrations, respectively. All were located on chromosome 22 at or near the CYP2D6 gene, encoding CYP450 2D6, metoprolol's main metabolizing enzyme. Conclusion: The results reinforce previous findings of the importance of the CYP2D6 locus for metoprolol concentrations and confirm that large biobanks can be used to identify genetic determinants of drug pharmacokinetics at a GWAS significance level.

Females present higher dose-adjusted drug concentrations of metoprolol and allopurinol/oxypurinol than males.

Females present a higher risk of adverse drug reactions. Sex-related differences in drug concentrations may contribute to these observations but they remain understudied given the underrepresentation of females in clinical trials. The aim of this study was to investigate whether anthropometric and socioeconomic factors and comorbidities could explain sex-related differences in concentrations and dosing for metoprolol and oxypurinol, the active metabolite of allopurinol. We conducted an analysis of two cross-sectional studies. Participants were self-described "White" adults taking metoprolol or allopurinol selected from the Montreal Heart Institute Hospital Cohort. A total of 1007 participants were included in the metoprolol subpopulation and 459 participants in the allopurinol subpopulation; 73% and 86% of the participants from the metoprolol and allopurinol subpopulations were males, respectively. Females presented higher age- and dose-adjusted concentrations of both metoprolol and oxypurinol (both p < 0.03). Accordingly, females presented higher unadjusted and age-adjusted concentration:dose ratio of both metoprolol and allopurinol/oxypurinol compared to males (all p < 3.0 × 10-4 ). Sex remained an independent predictor of metoprolol concentrations (p < 0.01), but not of oxypurinol concentrations, after adjusting for other predictors. In addition to sex, age, daily dose, use of moderate to strong CYP2D6 inhibitors, weight, and CYP2D6 genotype-inferred phenotype were associated with concentrations of metoprolol (all p < 0.01). Daily dose, weight, estimated glomerular filtration rate (eGFR), and employment status were associated with oxypurinol concentrations (all p < 0.01). Females present higher dose-adjusted concentrations of metoprolol and oxypurinol than males. This suggests the need for sex-specific dosing requirements for these drugs, although this hypothesis should be validated in prospective studies.

An association study of ABCG2 rs2231142 on the concentrations of allopurinol and its metabolites.

ABCG2 is a gene that codes for the human breast cancer resistance protein (BCRP). It is established that rs2231142 G>T, a single nucleotide polymorphism of the ABCG2 gene, is associated with gout and poor response to allopurinol, a uric acid-lowering agent used to treat this condition. It has also been suggested that oxypurinol, the primary active metabolite of allopurinol, is a substrate of the BCRP. We thus hypothesized that carrying the rs2231142 variant would be associated with decreased oxypurinol concentrations, which would explain the lower reduction in uric acid. We performed a cross-sectional study to investigate the association between the ABCG2 rs2231142 variant and oxypurinol, allopurinol, and allopurinol riboside concentrations in 459 participants from the Montreal Heart Institute Hospital Cohort. Age, sex, weight, use of diuretics, and estimated glomerular filtration rate were all significantly associated with oxypurinol plasma concentration. No association was found between rs2231142 and oxypurinol, allopurinol and allopurinol riboside plasma concentrations. Rs2231142 was not significantly associated with daily allopurinol dose in the overall population, but an association was observed in men, with T carriers receiving higher doses. Our results do not support a major role of ABCG2 in the pharmacokinetics of allopurinol or its metabolites. The underlying mechanism of the association between rs2231142 and allopurinol efficacy requires further investigation.

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.

Leveraging large observational studies to discover genetic determinants of drug concentrations: A proof-of-concept study.

Large, observational genetic studies are commonly used to identify genetic factors associated with diseases and disease-related traits. Such cohorts have not been commonly used to identify genetic predictors of drug dosing or concentrations, perhaps because of the heterogeneity in drug dosing and formulation, and the random timing of blood sampling. We hypothesized that large sample sizes relative to traditional pharmacokinetic studies would compensate for this variability and enable the identification of pharmacogenetic predictors of drug concentrations. We performed a cross-sectional, proof-of-concept association study to replicate the well-established association between metoprolol concentrations and CYP2D6 genotype-inferred metabolizer phenotypes in participants from the Montreal Heart Institute Hospital Cohort undergoing metoprolol therapy. Plasma concentrations of metoprolol and α-hydroxymetoprolol (α-OH-metoprolol) were measured in samples collected randomly regarding the previous metoprolol dose. A total of 999 individuals were included. The metoprolol daily dose ranged from 6.25 to 400 mg (mean 84.3 ± 57.1 mg). CYP2D6-inferred phenotype was significantly associated with both metoprolol and α-OH-metoprolol in unadjusted and adjusted models (all p < 10-14 ). Models for metoprolol daily dose showed consistent results. Our study suggests that randomly drawn blood samples from biobanks can serve as a new approach to discover genetic associations related to drug concentrations and dosing, with potentially broader implications for genomewide association studies on the pharmacogenomics of drug metabolism.

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.

The Human Immunopeptidome Project: A Roadmap to Predict and Treat Immune Diseases.

The science that investigates the ensembles of all peptides associated to human leukocyte antigen (HLA) molecules is termed "immunopeptidomics" and is typically driven by mass spectrometry (MS) technologies. Recent advances in MS technologies, neoantigen discovery and cancer immunotherapy have catalyzed the launch of the Human Immunopeptidome Project (HIPP) with the goal of providing a complete map of the human immunopeptidome and making the technology so robust that it will be available in every clinic. Here, we provide a long-term perspective of the field and we use this framework to explore how we think the completion of the HIPP will truly impact the society in the future. In this context, we introduce the concept of immunopeptidome-wide association studies (IWAS). We highlight the importance of large cohort studies for the future and how applying quantitative immunopeptidomics at population scale may provide a new look at individual predisposition to common immune diseases as well as responsiveness to vaccines and immunotherapies. Through this vision, we aim to provide a fresh view of the field to stimulate new discussions within the community, and present what we see as the key challenges for the future for unlocking the full potential of immunopeptidomics in this era of precision medicine.

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.

Pharmacogenetic content of commercial genome-wide genotyping arrays.

AIM: We have evaluated the pharmacogenetic content of commercial human genome-wide genotyping arrays, as it is a critical determinant to enabling pharmacogenomic discoveries. METHODS: Using bioinformatics approaches, we assessed 27,811 genetic variants in 3146 genes for their presence in 18 Illumina and 15 Affymetrix genome-wide arrays. RESULTS: 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. CONCLUSION: 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.

A prospective study of the impact of AGTR1 A1166C on the effects of candesartan in patients with heart failure.

AIM: To evaluate the impact of AGTR1 A1166C (rs5186) on the response to candesartan in patients with heart failure. MATERIALS & METHODS: Prospective, multicentre, open-label study. We studied 299 symptomatic patients with heart failure presenting a left ventricular ejection fraction ≤40%. RESULTS: Reductions in the primary end points of natriuretic peptides were not significantly associated with AGTR1 A1166C. Nevertheless, carrying the 1166C allele was associated with a greater compensatory increase in renin activity (p = 0.037) after 16 weeks of treatment with candesartan and a more modest effect on aldosterone concentrations (p = 0.022). CONCLUSION: AGTR1 1166C carriers may experience a greater long-term compensatory renin-angiotensin-aldosterone system activation following treatment with candesartan. Whether these associations ultimately influence clinical outcomes requires investigation. Clinicaltrials.gov : NCT00400582.

Rationale, design, and preliminary results of the Quebec Warfarin Cohort Study.

Over- and undercoagulation with warfarin are associated with hemorrhagic and thromboembolic events, respectively. Genetic and clinical factors affect warfarin response, and the causes of this variability remain unclear. We present descriptive statistics and test for predictors of poor anticoagulation control. The Quebec Warfarin Cohort (QWC) comprises 1059 new warfarin users, with prospective follow-up using telephone questionnaires every 3 months for 1 year, and using healthcare administrative databases (RAMQ and Med-Echo) for 5 years prior to cohort entry and up to 10 years following active patient participation. Genetic material was collected, and genotyping of CYP2C9 and VKORC1 genes was conducted. Measured outcomes included the percentage of time patients spent within therapeutic range, anticoagulation control, warfarin dose, bleeding, and thromboembolic events. We report baseline characteristics and outcomes after 1 year of follow-up. Poor anticoagulation control was defined as time in therapeutic range <60% in the 3- to 12-month interval. Participants had a mean age of 71 years, and 62% were men. The most common indication for warfarin was atrial fibrillation (87%). Mean time in therapeutic range was 56% (±25%) in the 3 months following warfarin initiation, and 70% (±21%) in the 3- to 12-month interval. During follow-up, the rate of stroke or systemic embolism was 1.8 events per 100 person-years; for major bleeding events, 3.3 events per 100 person-years. Independent predictors of poor anticoagulation control were chronic kidney disease, heart failure, dyslipidemia, and age. The QWC represents a good research cohort to investigate clinical and genetic factors in a warfarin-anticoagulated population.

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.

Genotype-Dependent Effects of Dalcetrapib on Cholesterol Efflux and Inflammation: Concordance With Clinical Outcomes.

BACKGROUND: Dalcetrapib effects on cardiovascular outcomes are determined by adenylate cyclase 9 gene polymorphisms. Our aim was to determine whether these clinical end point results are also associated with changes in reverse cholesterol transport and inflammation. METHODS AND RESULTS: Participants of the dal-OUTCOMES and dal-PLAQUE-2 trials were randomly assigned to receive dalcetrapib or placebo in addition to standard care. High-sensitivity C-reactive protein was measured at baseline and at end of study in 5243 patients from dal-OUTCOMES also genotyped for the rs1967309 polymorphism in adenylate cyclase 9. Cholesterol efflux capacity of high-density lipoproteins from J774 macrophages after cAMP stimulation was determined at baseline and 12 months in 171 genotyped patients from dal-PLAQUE-2. Treatment with dalcetrapib resulted in placebo-adjusted geometric mean percent increases in high-sensitivity C-reactive protein from baseline to end of trial of 18.1% (P=0.0009) and 18.7% (P=0.00001) in participants with the GG and AG genotypes, respectively, but the change was -1.0% (P=0.89) in those with the protective AA genotype. There was an interaction between the treatment arm and the genotype groups (P=0.02). Although the mean change in cholesterol efflux was similar among study arms in patients with GG genotype (mean: 7.8% and 7.4%), increases were 22.3% and 3.5% with dalcetrapib and placebo for those with AA genotype (P=0.005). There was a significant genetic effect for change in efflux for dalcetrapib (P=0.02), but not with placebo. CONCLUSIONS: Genotype-dependent effects on C-reactive protein and cholesterol efflux are supportive of dalcetrapib benefits on atherosclerotic cardiovascular outcomes in patients with the AA genotype at polymorphism rs1967309. CLINICAL TRIALS REGISTRATION: ClinicalTrials.gov; Unique Identifiers: NCT00658515 and NCT01059682.

Whole-genome sequencing in French Canadians from Quebec.

Genome-wide association studies (GWAS) have had a tremendous success in the identification of common DNA sequence variants associated with complex human diseases and traits. However, because of their design, GWAS are largely inappropriate to characterize the role of rare and low-frequency DNA variants on human phenotypic variation. Rarer genetic variation is geographically more restricted, supporting the need for local whole-genome sequencing (WGS) efforts to study these variants in specific populations. Here, we present the first large-scale low-pass WGS of the French-Canadian population. Specifically, we sequenced at ~5.6× coverage the whole genome of 1970 French Canadians recruited by the Montreal Heart Institute Biobank and identified 29 million bi-allelic variants (31 % novel), including 19 million variants with a minor allele frequency (MAF) <0.5 %. Genotypes from the WGS data are highly concordant with genotypes obtained by exome array on the same individuals (99.8 %), even when restricting this analysis to rare variants (MAF <0.5, 99.9 %) or heterozygous sites (98.9 %). To further validate our data set, we showed that we can effectively use it to replicate several genetic associations with myocardial infarction risk and blood lipid levels. Furthermore, we analyze the utility of our WGS data set to generate a French-Canadian-specific imputation reference panel and to infer population structure in the Province of Quebec. Our results illustrate the value of low-pass WGS to study the genetics of human diseases in the founder French-Canadian population.

Testing the role of predicted gene knockouts in human anthropometric trait variation.

Although the role of complete gene inactivation by two loss-of-function mutations inherited in trans is well-established in recessive Mendelian diseases, we have not yet explored how such gene knockouts (KOs) could influence complex human phenotypes. Here, we developed a statistical framework to test the association between gene KOs and quantitative human traits. Our method is flexible, publicly available, and compatible with common genotype format files (e.g. PLINK and vcf). We characterized gene KOs in 4498 participants from the NHLBI Exome Sequence Project (ESP) sequenced at high coverage (>100×), 1976 French Canadians from the Montreal Heart Institute Biobank sequenced at low coverage (5.7×), and >100 000 participants from the Genetic Investigation of ANthropometric Traits (GIANT) Consortium genotyped on an exome array. We tested associations between gene KOs and three anthropometric traits: body mass index (BMI), height and BMI-adjusted waist-to-hip ratio (WHR). Despite our large sample size and multiple datasets available, we could not detect robust associations between specific gene KOs and quantitative anthropometric traits. Our results highlight several limitations and challenges for future gene KO studies in humans, in particular when there is no prior knowledge on the phenotypes that might be affected by the tested gene KOs. They also suggest that gene KOs identified with current DNA sequencing methodologies probably do not strongly influence normal variation in BMI, height, and WHR in the general human population.