Progressive multifocal leukoencephalopathy genetic risk variants for pharmacovigilance of immunosuppressant therapies

Background

Progressive multifocal leukoencephalopathy (PML) is a rare and often lethal brain disorder caused by the common, typically benign polyomavirus 2, also known as JC virus (JCV). In a small percentage of immunosuppressed individuals, JCV is reactivated and infects the brain, causing devastating neurological defects. A wide range of immunosuppressed groups can develop PML, such as patients with: HIV/AIDS, hematological malignancies (e.g., leukemias, lymphomas, and multiple myeloma), autoimmune disorders (e.g., psoriasis, rheumatoid arthritis, and systemic lupus erythematosus), and organ transplants. In some patients, iatrogenic (i.e., drug-induced) PML occurs as a serious adverse event from exposure to immunosuppressant therapies used to treat their disease (e.g., hematological malignancies and multiple sclerosis). While JCV infection and immunosuppression are necessary, they are not sufficient to cause PML.

Methods

We hypothesized that patients may also have a genetic susceptibility from the presence of rare deleterious genetic variants in immune-relevant genes (e.g., those that cause inborn errors of immunity). In our prior genetic study of 184 PML cases, we discovered 19 candidate PML risk variants. In the current study of another 152 cases, we validated 4 of 19 variants in both population controls (gnomAD 3.1) and matched controls (JCV+ multiple sclerosis patients on a PML-linked drug ≥ 2 years).

Results

The four variants, found in immune system genes with strong biological links, are: C8B, 1-57409459-C-A, rs139498867; LY9 (alias SLAMF3), 1-160769595-AG-A, rs763811636; FCN2, 9-137779251-G-A, rs76267164; STXBP2, 19-7712287-G-C, rs35490401. Carriers of any one of these variants are shown to be at high risk of PML when drug-exposed PML cases are compared to drug-exposed matched controls: P value = 3.50E-06, OR = 8.7 [3.7-20.6]. Measures of clinical validity and utility compare favorably to other genetic risk tests, such as BRCA1 and BRCA2 screening for breast cancer risk and HLA-B^*15:02 pharmacogenetic screening for pharmacovigilance of carbamazepine to prevent Stevens-Johnson Syndrome and Toxic Epidermal Necrolysis.

Conclusion

For the first time, a PML genetic risk test can be implemented for screening patients taking or considering treatment with a PML-linked drug in order to decrease the incidence of PML and enable safer use of highly effective therapies used to treat their underlying disease.

A principal component meta-analysis on multiple anthropometric traits identifies novel loci for body shape

Large consortia have revealed hundreds of genetic loci associated with anthropometric traits, one trait at a time. We examined whether genetic variants affect body shape as a composite phenotype that is represented by a combination of anthropometric traits. We developed an approach that calculates averaged PCs (AvPCs) representing body shape derived from six anthropometric traits (body mass index, height, weight, waist and hip circumference, waist-to-hip ratio). The first four AvPCs explain >99% of the variability, are heritable, and associate with cardiometabolic outcomes. We performed genome-wide association analyses for each body shape composite phenotype across 65 studies and meta-analysed summary statistics. We identify six novel loci: LEMD2 and CD47 for AvPC1, RPS6KA5/C14orf159 and GANAB for AvPC3, and ARL15 and ANP32 for AvPC4. Our findings highlight the value of using multiple traits to define complex phenotypes for discovery, which are not captured by single-trait analyses, and may shed light onto new pathways.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study

Genome-wide association studies (GWAS) have identified more than 100 genetic variants contributing to BMI, a measure of body size, or waist-to-hip ratio (adjusted for BMI, WHR_adjBMI), a measure of body shape. Body size and shape change as people grow older and these changes differ substantially between men and women. To systematically screen for age- and/or sex-specific effects of genetic variants on BMI and WHR_adjBMI, we performed meta-analyses of 114 studies (up to 320,485 individuals of European descent) with genome-wide chip and/or Metabochip data by the Genetic Investigation of Anthropometric Traits (GIANT) Consortium. Each study tested the association of up to ~2.8M SNPs with BMI and WHR_adjBMI in four strata (men ≤50y, men >50y, women ≤50y, women >50y) and summary statistics were combined in stratum-specific meta-analyses. We then screened for variants that showed age-specific effects (G x AGE), sex-specific effects (G x SEX) or age-specific effects that differed between men and women (G x AGE x SEX). For BMI, we identified 15 loci (11 previously established for main effects, four novel) that showed significant (FDR<5%) age-specific effects, of which 11 had larger effects in younger (<50y) than in older adults (≥50y). No sex-dependent effects were identified for BMI. For WHR_adjBMI, we identified 44 loci (27 previously established for main effects, 17 novel) with sex-specific effects, of which 28 showed larger effects in women than in men, five showed larger effects in men than in women, and 11 showed opposite effects between sexes. No age-dependent effects were identified for WHR_adjBMI. This is the first genome-wide interaction meta-analysis to report convincing evidence of age-dependent genetic effects on BMI. In addition, we confirm the sex-specificity of genetic effects on WHR_adjBMI. These results may provide further insights into the biology that underlies weight change with age or the sexually dimorphism of body shape.

Adult body size and body shape differ substantially between men and women and change over time. More than 100 genetic variants that influence body mass index (measure of body size) or waist-to-hip ratio (measure of body shape) have been identified. While there is evidence that some genetic loci affect body shape differently in men than in women, little is known about whether genetic effects differ in older compared to younger adults, and whether such changes differ between men and women. Therefore, we conducted a systematic genome-wide search, including 114 studies (>320,000 individuals), to specifically identify genetic loci with age- and or sex-dependent effects on body size and shape. We identified 15 loci of which the effect on BMI was different in older compared to younger adults, whereas we found no evidence for loci with different effects in men compared to women. The opposite was seen for body shape as we identified 44 loci of which the effect on waist-to-hip ratio differed between men and women, but no difference between younger and older adults were observed. Our observations may provide new insights into the biology that underlies weight change with age or the sexual dimorphism of body shape.

Interactions between vascular wall and perivascular adipose tissue reveal novel roles for adiponectin in the regulation of endothelial nitric oxide synthase function in human vessels

BACKGROUND

Adiponectin is an adipokine with potentially important roles in human cardiovascular disease states. We studied the role of adiponectin in the cross-talk between adipose tissue and vascular redox state in patients with atherosclerosis.

METHODS AND RESULTS

The study included 677 patients undergoing coronary artery bypass graft surgery. Endothelial function was evaluated by flow-mediated dilation of the brachial artery in vivo and by vasomotor studies in saphenous vein segments ex vivo. Vascular superoxide (O2(-)) and endothelial nitric oxide synthase (eNOS) uncoupling were quantified in saphenous vein and internal mammary artery segments. Local adiponectin gene expression and ex vivo release were quantified in perivascular (saphenous vein and internal mammary artery) subcutaneous and mesothoracic adipose tissue from 248 patients. Circulating adiponectin was independently associated with nitric oxide bioavailability and O2(-) production/eNOS uncoupling in both arteries and veins. These findings were supported by a similar association between functional polymorphisms in the adiponectin gene and vascular redox state. In contrast, local adiponectin gene expression/release in perivascular adipose tissue was positively correlated with O2(-) and eNOS uncoupling in the underlying vessels. In ex vivo experiments with human saphenous veins and internal mammary arteries, adiponectin induced Akt-mediated eNOS phosphorylation and increased tetrahydrobiopterin bioavailability, improving eNOS coupling. In ex vivo experiments with human saphenous veins/internal mammary arteries and adipose tissue, we demonstrated that peroxidation products produced in the vascular wall (ie, 4-hydroxynonenal) upregulate adiponectin gene expression in perivascular adipose tissue via a peroxisome proliferator-activated receptor-γ-dependent mechanism.

CONCLUSIONS

We demonstrate for the first time that adiponectin improves the redox state in human vessels by restoring eNOS coupling, and we identify a novel role of vascular oxidative stress in the regulation of adiponectin expression in human perivascular adipose tissue.

Normal passive viscoelasticity but abnormal myofibrillar force generation in human hypertrophic cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is characterized by left ventricular hypertrophy, increased ventricular stiffness and impaired diastolic filling. We investigated to what extent myocardial functional defects can be explained by alterations in the passive and active properties of human cardiac myofibrils. Skinned ventricular myocytes were prepared from patients with obstructive HCM (two patients with MYBPC3 mutations, one with a MYH7 mutation, and three with no mutation in either gene) and from four donors. Passive stiffness, viscous properties, and titin isoform expression were similar in HCM myocytes and donor myocytes. Maximal Ca²⁺-activated force was much lower in HCM myocytes (14 ± 1 kN/m²) than in donor myocytes (23 ± 3 kN/m²; P < 0.01), though cross-bridge kinetics (k_tr) during maximal Ca²⁺ activation were 10% faster in HCM myocytes. Myofibrillar Ca²⁺ sensitivity in HCM myocytes (pCa₅₀ = 6.40 ± 0.05) was higher than for donor myocytes (pCa₅₀ = 6.09 ± 0.02; P < 0.001) and was associated with reduced phosphorylation of troponin-I (ser-23/24) and MyBP-C (ser-282) in HCM myocytes. These characteristics were common to all six HCM patients and may therefore represent a secondary consequence of the known and unknown underlying genetic variants. Some HCM patients did however exhibit an altered relationship between force and cross-bridge kinetics at submaximal Ca²⁺ concentrations, which may reflect the primary mutation. We conclude that the passive viscoelastic properties of the myocytes are unlikely to account for the increased stiffness of the HCM ventricle. However, the low maximum Ca²⁺-activated force and high Ca²⁺ sensitivity of the myofilaments are likely to contribute substantially to any systolic and diastolic dysfunction, respectively, in hearts of HCM patients.

Evidence from human myectomy samples that MYBPC3 mutations cause hypertrophic cardiomyopathy through haploinsufficiency

RATIONALE

Most sarcomere gene mutations that cause hypertrophic cardiomyopathy are missense alleles that encode dominant negative proteins. The potential exceptions are mutations in the MYBPC3 gene (encoding cardiac myosin-binding protein-C [MyBP-C]), which frequently encode truncated proteins.

OBJECTIVE

We sought to determine whether there was evidence of haploinsufficiency in hypertrophic cardiomyopathy caused by MYBPC3 mutations by comparing left ventricular muscle from patients undergoing surgical myectomy with samples from donor hearts.

METHODS AND RESULTS

MyBP-C protein and mRNA levels were quantitated using immunoblotting and RT-PCR. Nine of 37 myectomy samples had mutations in MYBPC3: 2 missense alleles (Glu258Lys, Arg502Trp) and 7 premature terminations. No specific truncated MyBP-C peptides were detected in whole muscle homogenates of hypertrophic cardiomyopathy tissue. However, the overall level of MyBP-C in myofibrils was significantly reduced (P<0.0005) in tissue containing either a truncation or missense MYBPC3 mutation: 0.76+/-0.03 compared with 1.00+/-0.05 in donor and 1.01+/-0.06 in non-MYBPC3 mutant myectomies.

CONCLUSIONS

The absence of any detectable truncated MyBP-C argues against its incorporation in the myofiber and any dominant negative effect. In contrast, the lowered relative level of full length protein in both truncation and missense MYBPC3 mutations argues strongly that haploinsufficiency is sufficient to cause the disease.

The molecular phenotype of human cardiac myosin associated with hypertrophic obstructive cardiomyopathy

AIM

The aim of the study was to compare the functional and structural properties of the motor protein, myosin, and isolated myocyte contractility in heart muscle excised from hypertrophic cardiomyopathy patients by surgical myectomy with explanted failing heart and non-failing donor heart muscle.

METHODS

Myosin was isolated and studied using an in vitro motility assay. The distribution of myosin light chain-1 isoforms was measured by two-dimensional electrophoresis. Myosin light chain-2 phosphorylation was measured by sodium dodecyl sulphate-polyacrylamide gel electrophoresis using Pro-Q Diamond phosphoprotein stain.

RESULTS

The fraction of actin filaments moving when powered by myectomy myosin was 21% less than with donor myosin (P = 0.006), whereas the sliding speed was not different (0.310 +/- 0.034 for myectomy myosin vs. 0.305 +/- 0.019 microm/s for donor myosin in six paired experiments). Failing heart myosin showed 18% reduced motility. One myectomy myosin sample produced a consistently higher sliding speed than donor heart myosin and was identified with a disease-causing heavy chain mutation (V606M). In myectomy myosin, the level of atrial light chain-1 relative to ventricular light chain-1 was 20 +/- 5% compared with 11 +/- 5% in donor heart myosin and the level of myosin light chain-2 phosphorylation was decreased by 30-45%. Isolated cardiomyocytes showed reduced contraction amplitude (1.61 +/- 0.25 vs. 3.58 +/- 0.40%) and reduced relaxation rates compared with donor myocytes (TT(50%) = 0.32 +/- 0.09 vs. 0.17 +/- 0.02 s).

CONCLUSION

Contractility in myectomy samples resembles the hypocontractile phenotype found in end-stage failing heart muscle irrespective of the primary stimulus, and this phenotype is not a direct effect of the hypertrophy-inducing mutation. The presence of a myosin heavy chain mutation causing hypertrophic cardiomyopathy can be predicted from a simple functional assay.

Polymorphisms in the cardiac sodium channel promoter displaying variant in vitro expression activity

Variable transcription of the cardiac sodium channel gene is a candidate mechanism determining arrhythmia susceptibility. We have previously cloned and characterized the core promoter and flanking region of SCN5A, encoding the cardiac sodium channel. Loss-of-function mutations in this gene have been reported in approximately 20% of patients with Brugada syndrome, an inherited cardiac electrical disorder associated with a high incidence of life-threatening arrhythmias. In this study, we identified DNA variants in the proximal 2.8 kb promoter region of SCN5A and determined their frequency in 1,121 subjects. This population consisted of 88 Brugada syndrome patients with no SCN5A coding region mutation, and 1,033 anonymized subjects from various ethnicities. Variant promoter activity was assayed in CHO cells and neonatal cardiomyocytes by transient transfection of promoter-reporter constructs. Single-nucleotide polymorphisms (SNPs) were identified at approximately 1/200 base pairs which are: 11 in the 5'-flanking region, 1 in exon 1, and 5 in intron 1. In addition, a haplotype consisting of two SNPs in complete linkage disequilibrium was identified. Minor allele frequencies were >5% in at least one ethnic panel at 5/19 polymorphic sites. In vitro functional analysis in cardiomyocytes identified four variants with significantly (P<0.05) reduced reporter activity (up to 63% reduction). The largest changes were seen with c.-225-1790 G>A, which reduced reporter activity by 62.8% in CHO cells and 55% in cardiomyocytes. From these results, we can conclude that the SCN5A core promoter includes multiple DNA polymorphisms with altered in vitro activity, further supporting the concept of interindividual variability in transcription of this cardiac ion channel gene.

Genome-wide association study of restless legs syndrome identifies common variants in three genomic regions

Restless legs syndrome (RLS) is a frequent neurological disorder characterized by an imperative urge to move the legs during night, unpleasant sensation in the lower limbs, disturbed sleep and increased cardiovascular morbidity. In a genome-wide association study we found highly significant associations between RLS and intronic variants in the homeobox gene MEIS1, the BTBD9 gene encoding a BTB(POZ) domain as well as variants in a third locus containing the genes encoding mitogen-activated protein kinase MAP2K5 and the transcription factor LBXCOR1 on chromosomes 2p, 6p and 15q, respectively. Two independent replications confirmed these association signals. Each genetic variant was associated with a more than 50% increase in risk for RLS, with the combined allelic variants conferring more than half of the risk. MEIS1 has been implicated in limb development, raising the possibility that RLS has components of a developmental disorder.

The common non-synonymous variant G38S of the KCNE1-(minK)-gene is not associated to QT interval in Central European Caucasians: results from the KORA study

AIMS

The QT interval in the general population is a complex trait with 30-50% heritability. QT prolongation is associated with an increased risk of sudden death. A recent family-based study found an association between QT interval and the common non-synonymous Glycin 38 Serine variant (G38S, rs1805127) of the KCNE1 gene coding for the minK-potassium channel subunit. We intended to replicate this finding in a large population sample of central European Caucasian ancestry as part of our ongoing search for genetic variants predisposing to arrhythmias.

METHODS AND RESULTS

We studied 3966 unrelated individuals from the KORA S4 population-based study without atrial fibrillation, pacemaker implant, or pregnancy. Individuals were genotyped by MALDI-TOF mass spectrometry. We did not detect any significant association between the genotypes of the G38S variant and the QT interval in the entire population or in any gender.

CONCLUSION

Unlike the common Lysine 897 Threonine variant of KCNH2 (K897T, rs1805123) the G38S variant of KCNE1 does not appear to have a strong modifying effect on QT interval. However, we cannot rule out an effect of G38S on QT in other ethnic groups, under exercise or medications or on the risk for arrhythmias and sudden death.

A common genetic variant in the NOS1 regulator NOS1AP modulates cardiac repolarization

Extremes of the electrocardiographic QT interval, a measure of cardiac repolarization, are associated with increased cardiovascular mortality. We identified a common genetic variant influencing this quantitative trait through a genome-wide association study on 200 subjects at the extremes of a population-based QT interval distribution of 3,966 subjects from the KORA cohort in Germany, with follow-up screening of selected markers in the remainder of the cohort. We validated statistically significant findings in two independent samples of 2,646 subjects from Germany and 1,805 subjects from the US Framingham Heart Study. This genome-wide study identified NOS1AP (CAPON), a regulator of neuronal nitric oxide synthase, as a new target that modulates cardiac repolarization. Approximately 60% of subjects of European ancestry carry at least one minor allele of the NOS1AP genetic variant, which explains up to 1.5% of QT interval variation.

Common Variants in Myocardial Ion Channel Genes Modify the QT Interval in the General Population

Altered myocardial repolarization is one of the important substrates of ventricular tachycardia and fibrillation. The influence of rare gene variants on repolarization is evident in familial long QT syndrome. To investigate the influence of common gene variants on the QT interval we performed a linkage disequilibrium based SNP association study of four candidate genes. Using a two-step design we analyzed 174 SNPs from the KCNQ1, KCNH2, KCNE1, and KCNE2 genes in 689 individuals from the population-based KORA study and 14 SNPs with results suggestive of association in a confirmatory sample of 3277 individuals from the same survey. We detected association to a gene variant in intron 1 of the KCNQ1 gene (rs757092, +1.7 ms/allele, P =0.0002) and observed weaker association to a variant upstream of the KCNE1 gene (rs727957, +1.2 ms/allele, P =0.0051). In addition we detected association to two SNPs in the KCNH2 gene, the previously described K897T variant (rs1805123, −1.9 ms/allele, P =0.0006) and a gene variant that tags a different haplotype in the same block (rs3815459, +1.7 ms/allele, P =0.0004). The analysis of additive effects by an allelic score explained a 10.5 ms difference in corrected QT interval length between extreme score groups and 0.95% of trait variance ( P <0.00005). These results confirm previous heritability studies indicating that repolarization is a complex trait with a significant heritable component and demonstrate that high-resolution SNP-mapping in large population samples can detect and fine map quantitative trait loci even if locus specific heritabilities are small.