Twenty-Five Novel Loci for Carotid Intima-Media Thickness: A Genome-Wide Association Study in >45 000 Individuals and Meta-Analysis of >100 000 Individuals.

OBJECTIVE: Carotid artery intima-media thickness (cIMT) is a widely accepted marker of subclinical atherosclerosis. Twenty susceptibility loci for cIMT were previously identified and the identification of additional susceptibility loci furthers our knowledge on the genetic architecture underlying atherosclerosis. APPROACH AND RESULTS: We performed 3 genome-wide association studies in 45 185 participants from the UK Biobank study who underwent cIMT measurements and had data on minimum, mean, and maximum thickness. We replicated 15 known loci and identified 20 novel loci associated with cIMT at P<5×10-8. Seven novel loci (ZNF385D, ADAMTS9, EDNRA, HAND2, MYOCD, ITCH/EDEM2/MMP24, and MRTFA) were identified in all 3 phenotypes. An additional new locus (LOXL1) was identified in the meta-analysis of the 3 phenotypes. Sex interaction analysis revealed sex differences in 7 loci including a novel locus (SYNE3) in males. Meta-analysis of UK Biobank data with a previous meta-analysis led to identification of three novel loci (APOB, FIP1L1, and LOXL4). Transcriptome-wide association analyses implicated additional genes ARHGAP42, NDRG4, and KANK2. Gene set analysis showed an enrichment in extracellular organization and the PDGF (platelet-derived growth factor) signaling pathway. We found positive genetic correlations of cIMT with coronary artery disease rg=0.21 (P=1.4×10-7), peripheral artery disease rg=0.45 (P=5.3×10-5), and systolic blood pressure rg=0.30 (P=4.0×10-18). A negative genetic correlation between average of maximum cIMT and high-density lipoprotein was found rg=-0.12 (P=7.0×10-4). CONCLUSIONS: Genome-wide association meta-analyses in >100 000 individuals identified 25 novel loci associated with cIMT providing insights into genes and tissue-specific regulatory mechanisms of proatherosclerotic processes. We found evidence for shared biological mechanisms with cardiovascular diseases.

Genome-wide association analysis in dilated cardiomyopathy reveals two new players in systolic heart failure on chromosomes 3p25.1 and 22q11.23.

AIMS: Our objective was to better understand the genetic bases of dilated cardiomyopathy (DCM), a leading cause of systolic heart failure. METHODS AND RESULTS: We conducted the largest genome-wide association study performed so far in DCM, with 2719 cases and 4440 controls in the discovery population. We identified and replicated two new DCM-associated loci on chromosome 3p25.1 [lead single-nucleotide polymorphism (SNP) rs62232870, P = 8.7 × 10-11 and 7.7 × 10-4 in the discovery and replication steps, respectively] and chromosome 22q11.23 (lead SNP rs7284877, P = 3.3 × 10-8 and 1.4 × 10-3 in the discovery and replication steps, respectively), while confirming two previously identified DCM loci on chromosomes 10 and 1, BAG3 and HSPB7. A genetic risk score constructed from the number of risk alleles at these four DCM loci revealed a 3-fold increased risk of DCM for individuals with 8 risk alleles compared to individuals with 5 risk alleles (median of the referral population). In silico annotation and functional 4C-sequencing analyses on iPSC-derived cardiomyocytes identify SLC6A6 as the most likely DCM gene at the 3p25.1 locus. This gene encodes a taurine transporter whose involvement in myocardial dysfunction and DCM is supported by numerous observations in humans and animals. At the 22q11.23 locus, in silico and data mining annotations, and to a lesser extent functional analysis, strongly suggest SMARCB1 as the candidate culprit gene. CONCLUSION: This study provides a better understanding of the genetic architecture of DCM and sheds light on novel biological pathways underlying heart failure.

Association of Factor V Leiden With Subsequent Atherothrombotic Events: A GENIUS-CHD Study of Individual Participant Data.

BACKGROUND: Studies examining the role of factor V Leiden among patients at higher risk of atherothrombotic events, such as those with established coronary heart disease (CHD), are lacking. Given that coagulation is involved in the thrombus formation stage on atherosclerotic plaque rupture, we hypothesized that factor V Leiden may be a stronger risk factor for atherothrombotic events in patients with established CHD. METHODS: We performed an individual-level meta-analysis including 25 prospective studies (18 cohorts, 3 case-cohorts, 4 randomized trials) from the GENIUS-CHD (Genetics of Subsequent Coronary Heart Disease) consortium involving patients with established CHD at baseline. Participating studies genotyped factor V Leiden status and shared risk estimates for the outcomes of interest using a centrally developed statistical code with harmonized definitions across studies. Cox proportional hazards regression models were used to obtain age- and sex-adjusted estimates. The obtained estimates were pooled using fixed-effect meta-analysis. The primary outcome was composite of myocardial infarction and CHD death. Secondary outcomes included any stroke, ischemic stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality. RESULTS: The studies included 69 681 individuals of whom 3190 (4.6%) were either heterozygous or homozygous (n=47) carriers of factor V Leiden. Median follow-up per study ranged from 1.0 to 10.6 years. A total of 20 studies with 61 147 participants and 6849 events contributed to analyses of the primary outcome. Factor V Leiden was not associated with the combined outcome of myocardial infarction and CHD death (hazard ratio, 1.03 [95% CI, 0.92-1.16]; I2=28%; P-heterogeneity=0.12). Subgroup analysis according to baseline characteristics or strata of traditional cardiovascular risk factors did not show relevant differences. Similarly, risk estimates for the secondary outcomes including stroke, coronary revascularization, cardiovascular mortality, and all-cause mortality were also close to identity. CONCLUSIONS: Factor V Leiden was not associated with increased risk of subsequent atherothrombotic events and mortality in high-risk participants with established and treated CHD. Routine assessment of factor V Leiden status is unlikely to improve atherothrombotic events risk stratification in this population.

A comparison of two workflows for regulome and transcriptome-based prioritization of genetic variants associated with myocardial mass.

A typical task arising from main effect analyses in a Genome Wide Association Study (GWAS) is to identify single nucleotide polymorphisms (SNPs), in linkage disequilibrium with the observed signals, that are likely causal variants and the affected genes. The affected genes may not be those closest to associating SNPs. Functional genomics data from relevant tissues are believed to be helpful in selecting likely causal SNPs and interpreting implicated biological mechanisms, ultimately facilitating prevention and treatment in the case of a disease trait. These data are typically used post GWAS analyses to fine-map the statistically significant signals identified agnostically by testing all SNPs and applying a multiple testing correction. The number of tested SNPs is typically in the millions, so the multiple testing burden is high. Motivated by this, in this study we investigated an alternative workflow, which consists in utilizing the available functional genomics data as a first step to reduce the number of SNPs tested for association. We analyzed GWAS on electrocardiographic QRS duration using these two workflows. The alternative workflow identified more SNPs, including some residing in loci not discovered with the typical workflow. Moreover, the latter are corroborated by other reports on QRS duration. This indicates the potential value of incorporating functional genomics information at the onset in GWAS analyses.

Contribution of non-HLA incompatibility between donor and recipient to kidney allograft survival: genome-wide analysis in a prospective cohort.

BACKGROUND: The introduction of HLA matching of donors and recipients was a breakthrough in kidney transplantation. However, half of all transplanted kidneys still fail within 15 years after transplantation. Epidemiological data suggest a fundamental role of non-HLA alloimmunity. METHODS: We genotyped 477 pairs of deceased donors and first kidney transplant recipients with stable graft function at three months that were transplanted between Dec 1, 2005, and April 30, 2015. Genome-wide genetic mismatches in non-synonymous single nucleotide polymorphisms (nsSNPs) were calculated to identify incompatibilities in transmembrane and secreted proteins. We estimated the association between nsSNP mismatch and graft loss in a Cox proportional hazard model, adjusting for HLA mismatch and clinical covariates. Customised peptide arrays were generated to screen for antibodies against genotype-derived mismatched epitopes in 25 patients with biopsy-confirmed chronic antibody-mediated rejection. FINDINGS: 59 268 nsSNPs affecting a transmembrane or secreted protein were analysed. The median number of nsSNP mismatches in immune-accessible transmembrane and secreted proteins between donors and recipients was 1892 (IQR 1850-1936). The degree of nsSNP mismatch was independently associated with graft loss in a multivariable model adjusted for HLA eplet mismatch (HLA-A, HLA-B, HLA-C, HLA-DP, HLA-DQ, and HLA-DR). Each increase by a unit of one IQR had an HR of 1·68 (95% CI 1·17-2·41, p=0·005). 5-year death censored graft survival was 98% in the quartile with the lowest mismatch, 91% in the second quartile, 89% in the third quartile, and 82% in the highest quartile (p=0·003, log-rank test). Customised peptide arrays verified a donor-specific alloimmune response to genetically predicted mismatched epitopes. INTERPRETATION: Genetic mismatch of non-HLA haplotypes coding for transmembrane or secreted proteins is associated with an increased risk of functional graft loss independently of HLA incompatibility. As in HLA alloimmunity, donor-specific alloantibodies can be identified against genotype derived non-HLA epitopes. FUNDING: Austrian Science Fund, WWTF (Vienna Science and Technology Fund), and Ministry of Health of the Czech Republic.

Pharmacogenomics in kidney transplant recipients and potential for integration into practice.

WHAT IS KNOWN AND OBJECTIVE: Pharmacogenomic biomarkers are now used in many clinical care settings and represent one of the successes of precision medicine. Genetic variants are associated with pharmacokinetic and pharmacodynamic changes leading to medication adverse effects and changes in clinical response. Actionable pharmacogenomic variants are common in transplant recipients and have implications for medications used in transplant, but yet are not broadly incorporated into practice. METHODS: From the Clinical Pharmacogenetics Implementation Consortium and Dutch Pharmacogenetics Working Group guidelines, and PharmGKB databases, 12 pharmacogenomic genes with 30 variants were selected and used to create diplotypes and actionable pharmacogenomic phenotypes. A total of 853 kidney allograft recipients who had genomic information available from a genome-wide association study were included. RESULTS: Each recipient had at least one actionable pharmacogenomic diplotype/phenotype, whereas the majority (58%) had three or four actionable diplotypes/phenotypes and 17.4% had five or more among the 12 genes. The participants carried actionable diplotypes/phenotypes for multiple medications, including tacrolimus, azathioprine, clopidogrel, warfarin, simvastatin, voriconazole, antidepressants and proton-pump inhibitors. WHAT IS NEW AND CONCLUSION: Pharmacogenomic variants are common in transplant recipients, and transplant recipients receive medications that have actionable variants. CLINICAL TRIAL: Genomics of Transplantation, clinicaltrials.gov (NCT01714440).

Discovery of novel heart rate-associated loci using the Exome Chip.

Resting heart rate is a heritable trait, and an increase in heart rate is associated with increased mortality risk. Genome-wide association study analyses have found loci associated with resting heart rate, at the time of our study these loci explained 0.9% of the variation. This study aims to discover new genetic loci associated with heart rate from Exome Chip meta-analyses.Heart rate was measured from either elecrtrocardiograms or pulse recordings. We meta-analysed heart rate association results from 104 452 European-ancestry individuals from 30 cohorts, genotyped using the Exome Chip. Twenty-four variants were selected for follow-up in an independent dataset (UK Biobank, N = 134 251). Conditional and gene-based testing was undertaken, and variants were investigated with bioinformatics methods.We discovered five novel heart rate loci, and one new independent low-frequency non-synonymous variant in an established heart rate locus (KIAA1755). Lead variants in four of the novel loci are non-synonymous variants in the genes C10orf71, DALDR3, TESK2 and SEC31B. The variant at SEC31B is significantly associated with SEC31B expression in heart and tibial nerve tissue. Further candidate genes were detected from long-range regulatory chromatin interactions in heart tissue (SCD, SLF2 and MAPK8). We observed significant enrichment in DNase I hypersensitive sites in fetal heart and lung. Moreover, enrichment was seen for the first time in human neuronal progenitor cells (derived from embryonic stem cells) and fetal muscle samples by including our novel variants.Our findings advance the knowledge of the genetic architecture of heart rate, and indicate new candidate genes for follow-up functional studies.

Tacrolimus troughs and genetic determinants of metabolism in kidney transplant recipients: A comparison of four ancestry groups.

Tacrolimus trough and dose requirements vary dramatically between individuals of European and African American ancestry. These differences are less well described in other populations. We conducted an observational, prospective, multicenter study from which 2595 kidney transplant recipients of European, African, Native American, and Asian ancestry were studied for tacrolimus trough, doses, and genetic determinants of metabolism. We studied the well-known variants and conducted a CYP3A4/5 gene-wide analysis to identify new variants. Daily doses, and dose-normalized troughs were significantly different between the four groups (P < .001). CYP3A5*3 (rs776746) was associated with higher dose-normalized tacrolimus troughs in all groups but occurred at different allele frequencies and had differing effect sizes. The CYP3A5*6 (rs10264272) and *7 (rs413003343) variants were only present in African Americans. CYP3A4*22 (rs35599367) was not found in any of the Asian ancestry samples. We identified seven suggestive variants in the CYP3A4/5 genes associated with dose-normalized troughs in Native Americans (P = 1.1 × 10-5 -8.8 × 10-6 ) and one suggestive variant in Asian Americans (P = 5.6 × 10-6 ). Tacrolimus daily doses and dose-normalized troughs vary significantly among different ancestry groups. We identified potential new variants important in Asians and Native Americans. Studies with larger populations should be conducted to assess the importance of the identified suggestive variants.

Polygenic risk score as a determinant of risk of non-melanoma skin cancer in a European-descent renal transplant cohort.

Renal transplant recipients have an increased risk of non-melanoma skin cancer (NMSC) compared to in the general population. Here, we show polygenic risk scores (PRS) calculated from genome-wide association studies (GWAS) of NMSC in a general, nontransplant setting, can predict risk of, and time to posttransplant skin cancer. Genetic variants, reaching predefined P-value thresholds were chosen from published squamous cell carcinoma (SCC) and basal cell carcinoma (BCC) nontransplant GWAS. Using these GWAS, BCC and SCC PRS were calculated for each sample across three European ancestry renal transplant cohorts (n = 889) and tested as predictors of case:control status and time to NMSC posttransplant. BCC PRS calculated at P-value threshold 1 × 10-5 was the most significant predictor of case:control status of NMSC posttransplant (OR = 1.61; adjusted P = .0022; AUC [full model adjusted for clinical predictors and PRS] = 0.81). SCC PRS at P-value threshold 1 × 10-5 was the most significant predictor of time to posttransplant NMSC (adjusted P = 9.39 × 10-7 ; HR = 1.41, concordance [full model] = 0.74). PRS of nontransplant NMSC is predictive of case:control status and time to NMSC posttransplant. These results are relevant to how genomics can risk stratify patients to help develop personalized treatment regimens.

The impact of donor and recipient common clinical and genetic variation on estimated glomerular filtration rate in a European renal transplant population.

Genetic variation across the human leukocyte antigen loci is known to influence renal-transplant outcome. However, the impact of genetic variation beyond the human leukocyte antigen loci is less clear. We tested the association of common genetic variation and clinical characteristics, from both the donor and recipient, with posttransplant eGFR at different time-points, out to 5 years posttransplantation. We conducted GWAS meta-analyses across 10 844 donors and recipients from five European ancestry cohorts. We also analyzed the impact of polygenic risk scores (PRS), calculated using genetic variants associated with nontransplant eGFR, on posttransplant eGFR. PRS calculated using the recipient genotype alone, as well as combined donor and recipient genotypes were significantly associated with eGFR at 1-year posttransplant. Thirty-two percent of the variability in eGFR at 1-year posttransplant was explained by our model containing clinical covariates (including weights for death/graft-failure), principal components and combined donor-recipient PRS, with 0.3% contributed by the PRS. No individual genetic variant was significantly associated with eGFR posttransplant in the GWAS. This is the first study to examine PRS, composed of variants that impact kidney function in the general population, in a posttransplant context. Despite PRS being a significant predictor of eGFR posttransplant, the effect size of common genetic factors is limited compared to clinical variables.

NPHP1 (Nephrocystin-1) Gene Deletions Cause Adult-Onset ESRD.

Background Nephronophthisis (NPH) is the most prevalent genetic cause for ESRD in children. However, little is known about the prevalence of NPH in adult-onset ESRD. Homozygous full gene deletions of the NPHP1 gene encoding nephrocystin-1 are a prominent cause of NPH. We determined the prevalence of NPH in adults by assessing homozygous NPHP1 full gene deletions in adult-onset ESRD.Methods Adult renal transplant recipients from five cohorts of the International Genetics and Translational Research in Transplantation Network (iGeneTRAiN) underwent single-nucleotide polymorphism genotyping. After quality control, we determined autosomal copy number variants (such as deletions) on the basis of median log2 ratios and B-allele frequency patterns. The findings were independently validated in one cohort. Patients were included in the analysis if they had adult-onset ESRD, defined as start of RRT at ≥18 years old.Results We included 5606 patients with adult-onset ESRD; 26 (0.5%) showed homozygous NPHP1 deletions. No donor controls showed homozygosity for this deletion. Median age at ESRD onset was 30 (range, 18-61) years old for patients with NPH, with 54% of patients age ≥30 years old. Notably, only three (12%) patients were phenotypically classified as having NPH, whereas most patients were defined as having CKD with unknown etiology (n=11; 42%).Conclusions Considering that other mutation types in NPHP1 or mutations in other NPH-causing genes were not analyzed, NPH is a relatively frequent monogenic cause of adult-onset ESRD. Because 88% of patients had not been clinically diagnosed with NPH, wider application of genetic testing in adult-onset ESRD may be warranted.

Human Validation of Genes Associated With a Murine Atherosclerotic Phenotype.

OBJECTIVE: The genetically modified mouse is the most commonly used animal model for studying the pathogenesis of atherosclerotic disease. We aimed to assess if mice atherosclerosis-related genes could be validated in human disease through examination of results from genome-wide association studies. APPROACH AND RESULTS: We performed a systematic review to identify atherosclerosis-causing genes in mice and carried out gene-based association tests of their human orthologs for an association with human coronary artery disease and human large artery ischemic stroke. Moreover, we investigated the association of these genes with human atherosclerotic plaque characteristics. In addition, we assessed the presence of tissue-specific cis-acting expression quantitative trait loci for these genes in humans. Finally, using pathway analyses we show that the putative atherosclerosis-causing genes revealed few associations with human coronary artery disease, large artery ischemic stroke, or atherosclerotic plaque characteristics, despite the fact that the majority of these genes have cis-acting expression quantitative trait loci. CONCLUSIONS: A role for genes that has been observed in mice for atherosclerotic lesion development could scarcely be confirmed by studying associations of disease development with common human genetic variants. The value of murine atherosclerotic models for selection of therapeutic targets in human disease remains unclear.

No evidence for shared genetic basis of common variants in multiple sclerosis and amyotrophic lateral sclerosis.

Genome-wide association studies have been successful in identifying common variants that influence the susceptibility to complex diseases. From these studies, it has emerged that there is substantial overlap in susceptibility loci between diseases. In line with those findings, we hypothesized that shared genetic pathways may exist between multiple sclerosis (MS) and amyotrophic lateral sclerosis (ALS). While both diseases may have inflammatory and neurodegenerative features, epidemiological studies have indicated an increased co-occurrence within individuals and families. To this purpose, we combined genome-wide data from 4088 MS patients, 3762 ALS patients and 12 030 healthy control individuals in whom 5 440 446 single-nucleotide polymorphisms (SNPs) were successfully genotyped or imputed. We tested these SNPs for the excess association shared between MS and ALS and also explored whether polygenic models of SNPs below genome-wide significance could explain some of the observed trait variance between diseases. Genome-wide association meta-analysis of SNPs as well as polygenic analyses fails to provide evidence in favor of an overlap in genetic susceptibility between MS and ALS. Hence, our findings do not support a shared genetic background of common risk variants in MS and ALS.

A genome-wide association study identifies a functional ERAP2 haplotype associated with birdshot chorioretinopathy.

Birdshot chorioretinopathy (BSCR) is a rare form of autoimmune uveitis that can lead to severe visual impairment. Intriguingly, >95% of cases carry the HLA-A29 allele, which defines the strongest documented HLA association for a human disease. We have conducted a genome-wide association study in 96 Dutch and 27 Spanish cases, and 398 unrelated Dutch and 380 Spanish controls. Fine-mapping the primary MHC association through high-resolution imputation at classical HLA loci, identified HLA-A*29:02 as the principal MHC association (odds ratio (OR) = 157.5, 95% CI 91.6-272.6, P = 6.6 × 10(-74)). We also identified two novel susceptibility loci at 5q15 near ERAP2 (rs7705093; OR = 2.3, 95% CI 1.7-3.1, for the T allele, P = 8.6 × 10(-8)) and at 14q32.31 in the TECPR2 gene (rs150571175; OR = 6.1, 95% CI 3.2-11.7, for the A allele, P = 3.2 × 10(-8)). The association near ERAP2 was confirmed in an independent British case-control samples (combined meta-analysis P = 1.7 × 10(-9)). Functional analyses revealed that the risk allele of the polymorphism near ERAP2 is strongly associated with high mRNA and protein expression of ERAP2 in B cells. This study further defined an extremely strong MHC risk component in BSCR, and detected evidence for a novel disease mechanism that affects peptide processing in the endoplasmic reticulum.

Loci influencing blood pressure identified using a cardiovascular gene-centric array.

Blood pressure (BP) is a heritable determinant of risk for cardiovascular disease (CVD). To investigate genetic associations with systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP) and pulse pressure (PP), we genotyped ∼50 000 single-nucleotide polymorphisms (SNPs) that capture variation in ∼2100 candidate genes for cardiovascular phenotypes in 61 619 individuals of European ancestry from cohort studies in the USA and Europe. We identified novel associations between rs347591 and SBP (chromosome 3p25.3, in an intron of HRH1) and between rs2169137 and DBP (chromosome1q32.1 in an intron of MDM4) and between rs2014408 and SBP (chromosome 11p15 in an intron of SOX6), previously reported to be associated with MAP. We also confirmed 10 previously known loci associated with SBP, DBP, MAP or PP (ADRB1, ATP2B1, SH2B3/ATXN2, CSK, CYP17A1, FURIN, HFE, LSP1, MTHFR, SOX6) at array-wide significance (P < 2.4 × 10(-6)). We then replicated these associations in an independent set of 65 886 individuals of European ancestry. The findings from expression QTL (eQTL) analysis showed associations of SNPs in the MDM4 region with MDM4 expression. We did not find any evidence of association of the two novel SNPs in MDM4 and HRH1 with sequelae of high BP including coronary artery disease (CAD), left ventricular hypertrophy (LVH) or stroke. In summary, we identified two novel loci associated with BP and confirmed multiple previously reported associations. Our findings extend our understanding of genes involved in BP regulation, some of which may eventually provide new targets for therapeutic intervention.

The association between a genetic variant in the SULF2 gene, metabolic parameters and vascular disease in patients at high cardiovascular risk.

Clearance of triglyceride-rich lipoproteins (TRLs) is mediated by several receptors, including heparan sulfate proteoglycans (HSPGs). Sulfate glucosamine-6-O-endosulfatase-2 is a gene related to the regulation of HSPG. A variant in this gene, rs2281279, has been shown to be associated with triglycerides and insulin resistance. Objective: To determine the relationship between rs2281279, metabolic parameters and vascular events, and type 2 diabetes mellitus (T2DM) in patients at high cardiovascular risk and whether APOE genotype modifies this relationship. Methods: Patients (n = 4386) at high cardiovascular risk from the Utrecht Cardiovascular Cohort-Second Manifestations of Arterial Disease study were stratified according to their imputed rs2281279 genotype: AA (n = 2438), AG (n = 1642) and GG (n = 306). Effects of rs2281279 on metabolic parameters, vascular events and T2DM were analyzed with linear regression and Cox models. Results: There was no relationship between imputed rs2281279 genotype and triglycerides, non-high-density lipoprotein (HDL)-cholesterol, insulin and quantitative insulin sensitivity check index. During a median follow-up of 11.8 (IQR, 9.3-15.5) years, 1026 cardiovascular events and 320 limb events occurred. The presence of the G allele in rs2281279 did not affect the risk of vascular events [hazard ratio (HR), 1.03; 95% confidence interval (CI), 0.94-1.14] or limb events (HR, 0.92; 95% CI, 0.77-1.10). The presence of the G allele in rs2281279 did not affect the risk of T2DM (HR, 1.09; 95% CI, 0.94-1.27). The presence of the minor G allele of rs2281279 was associated with a beneficial risk profile in ε2ε2 patients, but not in ε3ε3 patients. Conclusions: Imputed rs2281279 genotype is not associated with metabolic parameters and does not increase the risk of vascular events or T2DM in patients at high risk for cardiovascular disease.

Cardiac MRI to guide heart failure and atrial fibrillation drug discovery: a Mendelian randomization analysis.

Background: drug development and disease prevention of heart failure (HF) and atrial fibrillation (AF) are impeded by a lack of robust early-stage surrogates. We determined to what extent cardiac magnetic resonance (CMR) measurements act as surrogates for the development of HF or AF in healthy individuals. Methods: Genetic data was sourced on the association with 22 atrial and ventricular CMR measurements. Mendelian randomization was used to determine CMR associations with atrial fibrillation (AF), heart failure (HF), non-ischemic cardiomyopathy (CMP), and dilated cardiomyopathy (DCM). Additionally, for the CMR surrogates of AF and HF, we explored their association with non-cardiac traits. Results: In total we found that 9 CMR measures were associated with the development of HF, 7 with development of non-ischemic CMR, 6 with DCM, and 12 with AF. biventricular ejection fraction (EF), biventricular or end-systolic volumes (ESV) and left-ventricular (LV) end diastolic volume (EDV) were associated with all 4 cardiac outcomes. Increased LV-MVR (mass to volume ratio) affected HF (odds ratio (OR) 0.83, 95%CI 0.79; 0.88), and DCM (OR 0.26, 95%CI 0.20; 0.34. We were able to identify 9 CMR surrogates for HF and/or AF (including LV-MVR, biventricular EDV, ESV, and right-ventricular EF) which associated with non-cardiac traits such as blood pressure, lung function traits, BMI, cardioembolic stroke, and late-onset Alzheimer's disease. Conclusion: CMR measurements may act as surrogate endpoints for the development of HF (including non-ischemic CMP and DCM) or AF. Additionally, we show that changes in cardiac function and structure measured through CMR, may affect diseases of other organs leading to lung disease or late-onset Alzheimer's disease.

Effect of tissue-grouped regulatory variants associated to type 2 diabetes in related secondary outcomes.

Genome-wide association studies have identified over five hundred loci that contribute to variation in type 2 diabetes (T2D), an established risk factor for many diseases. However, the mechanisms and extent through which these loci contribute to subsequent outcomes remain elusive. We hypothesized that combinations of T2D-associated variants acting on tissue-specific regulatory elements might account for greater risk for tissue-specific outcomes, leading to diversity in T2D disease progression. We searched for T2D-associated variants acting on regulatory elements and expression quantitative trait loci (eQTLs) in nine tissues. We used T2D tissue-grouped variant sets as genetic instruments to conduct 2-Sample Mendelian Randomization (MR) in ten related outcomes whose risk is increased by T2D using the FinnGen cohort. We performed PheWAS analysis to investigate whether the T2D tissue-grouped variant sets had specific predicted disease signatures. We identified an average of 176 variants acting in nine tissues implicated in T2D, and an average of 30 variants acting on regulatory elements that are unique to the nine tissues of interest. In 2-Sample MR analyses, all subsets of regulatory variants acting in different tissues were associated with increased risk of the ten secondary outcomes studied on similar levels. No tissue-grouped variant set was associated with an outcome significantly more than other tissue-grouped variant sets. We did not identify different disease progression profiles based on tissue-specific regulatory and transcriptome information. Bigger sample sizes and other layers of regulatory information in critical tissues may help identify subsets of T2D variants that are implicated in certain secondary outcomes, uncovering system-specific disease progression.

LoFTK: a framework for fully automated calculation of predicted Loss-of-Function variants and genes.

BACKGROUND: Loss-of-Function (LoF) variants in human genes are important due to their impact on clinical phenotypes and frequent occurrence in the genomes of healthy individuals. The association of LoF variants with complex diseases and traits may lead to the discovery and validation of novel therapeutic targets. Current approaches predict high-confidence LoF variants without identifying the specific genes or the number of copies they affect. Moreover, there is a lack of methods for detecting knockout genes caused by compound heterozygous (CH) LoF variants. RESULTS: We have developed the Loss-of-Function ToolKit (LoFTK), which allows efficient and automated prediction of LoF variants from genotyped, imputed and sequenced genomes. LoFTK enables the identification of genes that are inactive in one or two copies and provides summary statistics for downstream analyses. LoFTK can identify CH LoF variants, which result in LoF genes with two copies lost. Using data from parents and offspring we show that 96% of CH LoF genes predicted by LoFTK in the offspring have the respective alleles donated by each parent. CONCLUSIONS: LoFTK is a command-line based tool that provides a reliable computational workflow for predicting LoF variants from genotyped and sequenced genomes, identifying genes that are inactive in 1 or 2 copies. LoFTK is an open software and is freely available to non-commercial users at https://github.com/CirculatoryHealth/LoFTK .

Druggable proteins influencing cardiac structure and function: Implications for heart failure therapies and cancer cardiotoxicity.

Dysfunction of either the right or left ventricle can lead to heart failure (HF) and subsequent morbidity and mortality. We performed a genome-wide association study (GWAS) of 16 cardiac magnetic resonance (CMR) imaging measurements of biventricular function and structure. Cis-Mendelian randomization (MR) was used to identify plasma proteins associating with CMR traits as well as with any of the following cardiac outcomes: HF, non-ischemic cardiomyopathy, dilated cardiomyopathy (DCM), atrial fibrillation, or coronary heart disease. In total, 33 plasma proteins were prioritized, including repurposing candidates for DCM and/or HF: IL18R (providing indirect evidence for IL18), I17RA, GPC5, LAMC2, PA2GA, CD33, and SLAF7. In addition, 13 of the 25 druggable proteins (52%; 95% confidence interval, 0.31 to 0.72) could be mapped to compounds with known oncological indications or side effects. These findings provide leads to facilitate drug development for cardiac disease and suggest that cardiotoxicities of several cancer treatments might represent mechanism-based adverse effects.