Pharmacogenomics: Driving Personalized Medicine.

Personalized medicine tailors therapies, disease prevention, and health maintenance to the individual, with pharmacogenomics serving as a key tool to improve outcomes and prevent adverse effects. Advances in genomics have transformed pharmacogenetics, traditionally focused on single gene-drug pairs, into pharmacogenomics, encompassing all "-omics" fields (e.g., proteomics, transcriptomics, metabolomics, and metagenomics). This review summarizes basic genomics principles relevant to translation into therapies, assessing pharmacogenomics' central role in converging diverse elements of personalized medicine. We discuss genetic variations in pharmacogenes (drug-metabolizing enzymes, drug transporters, and receptors), their clinical relevance as biomarkers, and the legacy of decades of research in pharmacogenetics. All types of therapies, including proteins, nucleic acids, viruses, cells, genes, and irradiation, can benefit from genomics, expanding the role of pharmacogenomics across medicine. Food and Drug Administration approvals of personalized therapeutics involving biomarkers increase rapidly, demonstrating the growing impact of pharmacogenomics. A beacon for all therapeutic approaches, molecularly targeted cancer therapies highlight trends in drug discovery and clinical applications. To account for human complexity, multicomponent biomarker panels encompassing genetic, personal, and environmental factors can guide diagnosis and therapies, increasingly involving artificial intelligence to cope with extreme data complexities. However, clinical application encounters substantial hurdles, such as unknown validity across ethnic groups, underlying bias in health care, and real-world validation. This review address the underlying science and technologies germane to pharmacogenomics and personalized medicine, integrated with economic, ethical, and regulatory issues, providing insights into the current status and future direction of health care. SIGNIFICANCE STATEMENT: Personalized medicine aims to optimize health care for the individual patients with use of predictive biomarkers to improve outcomes and prevent adverse effects. Pharmacogenomics drives biomarker discovery and guides the development of targeted therapeutics. This review addresses basic principles and current trends in pharmacogenomics, with large-scale data repositories accelerating medical advances. The impact of pharmacogenomics is discussed, along with hurdles impeding broad clinical implementation, in the context of clinical care, ethics, economics, and regulatory affairs.

Ligand-Free Signaling of G-Protein-Coupled Receptors: Physiology, Pharmacology, and Genetics.

G-protein-coupled receptors (GPCRs) are ubiquitous sensors and regulators of cellular functions. Each GPCR exists in complex aggregates with multiple resting and active conformations. Designed to detect weak stimuli, GPCRs can also activate spontaneously, resulting in basal ligand-free signaling. Agonists trigger a cascade of events leading to an activated agonist-receptor G-protein complex with high agonist affinity. However, the ensuing signaling process can further remodel the receptor complex to reduce agonist affinity, causing rapid ligand dissociation. The acutely activated ligand-free receptor can continue signaling, as proposed for rhodopsin and µ opioid receptors, resulting in robust receptor activation at low agonist occupancy with enhanced agonist potency. Continued receptor stimulation can further modify the receptor complex, regulating sustained ligand-free signaling-proposed to play a role in opioid dependence. Basal, acutely agonist-triggered, and sustained elevated ligand-free signaling could each have distinct functions, reflecting multi-state conformations of GPCRs. This review addresses basal and stimulus-activated ligand-free signaling, its regulation, genetic factors, and pharmacological implications, focusing on opioid and serotonin receptors, and the growth hormone secretagogue receptor (GHSR). The hypothesis is proposed that ligand-free signaling of 5-HT2A receptors mediate therapeutic effects of psychedelic drugs. Research avenues are suggested to close the gaps in our knowledge of ligand-free GPCR signaling.

Ligand-Free Signaling of G-Protein-Coupled Receptors: Relevance to µ Opioid Receptors in Analgesia and Addiction.

Numerous G-protein-coupled receptors (GPCRs) display ligand-free basal signaling with potential physiological functions, a target in drug development. As an example, the µ opioid receptor (MOR) signals in ligand-free form (MOR-µ*), influencing opioid responses. In addition, agonists bind to MOR but can dissociate upon MOR activation, with ligand-free MOR-µ* carrying out signaling. Opioid pain therapy is effective but incurs adverse effects (ADRs) and risk of opioid use disorder (OUD). Sustained opioid agonist exposure increases persistent basal MOR-µ* activity, which could be a driving force for OUD and ADRs. Antagonists competitively prevent resting MOR (MOR-µ) activation to MOR-µ*, while common antagonists, such as naloxone and naltrexone, also bind to and block ligand-free MOR-µ*, acting as potent inverse agonists. A neutral antagonist, 6ß-naltrexol (6BN), binds to but does not block MOR-µ*, preventing MOR-µ activation only competitively with reduced potency. We hypothesize that 6BN gradually accelerates MOR-µ* reversal to resting-state MOR-µ. Thus, 6BN potently prevents opioid dependence in rodents, at doses well below those blocking antinociception or causing withdrawal. Acting as a 'retrograde addiction modulator', 6BN could represent a novel class of therapeutics for OUD. Further studies need to address regulation of MOR-µ* and, more broadly, the physiological and pharmacological significance of ligand-free signaling in GPCRs.

Phagosome-regulated mTOR signalling during sarcoidosis granuloma biogenesis.

INTRODUCTION: Sarcoidosis and tuberculosis are granulomatous pulmonary diseases characterised by heightened immune reactivity to Mycobacterium tuberculosis antigens. We hypothesised that an unsupervised analysis comparing the molecular characteristics of granulomas formed in response to M. tuberculosis antigens in patients with sarcoidosis or latent tuberculosis infection (LTBI) would provide novel insights into the pathogenesis of sarcoidosis. METHODS: A genomic analysis identified differentially expressed genes in granuloma-like cell aggregates formed by sarcoidosis (n=12) or LTBI patients (n=5) in an established in vitro human granuloma model wherein peripheral blood mononuclear cells were exposed to M. tuberculosis antigens (beads coated with purified protein derivative) and cultured for 7 days. Pathway analysis of differentially expressed genes identified canonical pathways, most notably antigen processing and presentation via phagolysosomes, as a prominent pathway in sarcoidosis granuloma formation. The phagolysosomal pathway promoted mechanistic target of rapamycin complex 1 (mTORc1)/STAT3 signal transduction. Thus, granuloma formation and related immune mediators were evaluated in the absence or presence of various pre-treatments known to prevent phagolysosome formation (chloroquine) or phagosome acidification (bafilomycin A1) or directly inhibit mTORc1 activation (rapamycin). RESULTS: In keeping with genomic analyses indicating enhanced phagolysosomal activation and predicted mTORc1 signalling, it was determined that sarcoidosis granuloma formation and related inflammatory mediator release was dependent upon phagolysosome assembly and acidification and mTORc1/S6/STAT3 signal transduction. CONCLUSIONS: Sarcoidosis granulomas exhibit enhanced and sustained intracellular antigen processing and presentation capacities, and related phagolysosome assembly and acidification are required to support mTORc1 signalling to promote sarcoidosis granuloma formation.

Functional CYP3A variants affecting tacrolimus trough blood concentrations in Chinese renal transplant recipients.

The aim of this study was to identify novel genetic variants affecting tacrolimus trough blood concentrations. We analyzed the association between 58 single nucleotide polymorphisms (SNPs) across the CYP3A gene cluster and the log-transformed tacrolimus concentration/dose ratio (log (C0/D)) in 819 renal transplant recipients (Discovery cohort). Multivariate linear regression was used to test for associations between tacrolimus log (C0/D) and clinical factors. Luciferase reporter gene assays were used to evaluate the functions of select SNPs. Associations of putative functional SNPs with log (C0/D) were further tested in 631 renal transplant recipients (Replication cohort). Nine SNPs were significantly associated with tacrolimus log (C0/D) after adjustment for CYP3A5*3 and clinical factors. Dual luciferase reporter assays indicated that the rs4646450 G allele and rs3823812 T allele were significantly associated with increased normalized luciferase activity ratios (p < 0.01). Moreover, CYP3A7*2 was associated with higher TAC log(C0/D) in the group of CYP3A5 expressers. Age, serum creatinine and hematocrit were significantly associated with tacrolimus log (C0/D). CYP3A7*2, rs4646450, and rs3823812 are proposed as functional SNPs affecting tacrolimus trough blood concentrations in Chinese renal transplant recipients. Clinical factors also significantly affect tacrolimus metabolism.

Biased Opioid Antagonists as Modulators of Opioid Dependence: Opportunities to Improve Pain Therapy and Opioid Use Management.

Opioid analgesics are effective pain therapeutics but they cause various adverse effects and addiction. For safer pain therapy, biased opioid agonists selectively target distinct µ opioid receptor (MOR) conformations, while the potential of biased opioid antagonists has been neglected. Agonists convert a dormant receptor form (MOR-µ) to a ligand-free active form (MOR-µ*), which mediates MOR signaling. Moreover, MOR-µ converts spontaneously to MOR-µ* (basal signaling). Persistent upregulation of MOR-µ* has been invoked as a hallmark of opioid dependence. Contrasting interactions with both MOR-µ and MOR-µ* can account for distinct pharmacological characteristics of inverse agonists (naltrexone), neutral antagonists (6ß-naltrexol), and mixed opioid agonist-antagonists (buprenorphine). Upon binding to MOR-µ*, naltrexone but not 6ß-naltrexol suppresses MOR-µ*signaling. Naltrexone blocks opioid analgesia non-competitively at MOR-µ*with high potency, whereas 6ß-naltrexol must compete with agonists at MOR-µ, accounting for ~100-fold lower in vivo potency. Buprenorphine's bell-shaped dose-response curve may also result from opposing effects on MOR-µ and MOR-µ*. In contrast, we find that 6ß-naltrexol potently prevents dependence, below doses affecting analgesia or causing withdrawal, possibly binding to MOR conformations relevant to opioid dependence. We propose that 6ß-naltrexol is a biased opioid antagonist modulating opioid dependence at low doses, opening novel avenues for opioid pain therapy and use management.

IL-13-regulated Macrophage Polarization during Granuloma Formation in an In Vitro Human Sarcoidosis Model.

The mechanisms underlying abnormal granuloma formation in patients with sarcoidosis are complex and remain poorly understood. A novel in vitro human granuloma model was used to determine the molecular mechanisms of granuloma genesis in patients with sarcoidosis in response to putative disease-causing mycobacterial antigens. Peripheral blood mononuclear cells (PBMCs) from patients with active sarcoidosis and from normal, disease-free control subjects were incubated for 7 days with purified protein derivative-coated polystyrene beads. Molecular responses, as reflected by differential expression of genes, extracellular cytokine patterns, and cell surface receptor expression, were analyzed. Unbiased systems biology approaches were used to identify signaling pathways engaged during granuloma formation. Model findings were compared with human lung and mediastinal lymph node gene expression profiles. Compared with identically treated PBMCs of control subjects (n = 5), purified protein derivative-treated sarcoidosis PBMCs (n = 6) were distinguished by the formation of cellular aggregates resembling granulomas. Ingenuity Pathway Analysis of differential expression gene patterns identified molecular pathways that are primarily regulated by IL-13, which promotes alternatively activated (M2) macrophage polarization. M2 polarization was further demonstrated by immunohistochemistry performed on the in vitro sarcoidosis granuloma-like structures. IL-13-regulated gene pathways were confirmed in human sarcoidosis lung and mediastinal lymph node tissues. The in vitro human sarcoidosis granuloma model provides novel insights into early granuloma formation, particularly IL-13 regulation of molecular networks that regulate M2 macrophage polarization. M2 macrophages are predisposed to aggregation and multinucleated giant cell formation, which are characteristic features of sarcoidosis granulomas. Clinical trial registered with www.clinicaltrials.gov (NCT01857401).

Ligand-Free Estrogen Receptor α (ESR1) as Master Regulator for the Expression of CYP3A4 and Other Cytochrome P450 Enzymes in the Human Liver.

Cytochrome P450 3A4 isoform (CYP3A4) transcription is controlled by hepatic transcription factors (TFs), but how TFs dynamically interact remains uncertain. We hypothesize that several TFs form a regulatory network with nonlinear, dynamic, and hierarchical interactions. To resolve complex interactions, we have applied a computational approach for estimating Sobol's sensitivity indices (SSI) under generalized linear models to existing liver RNA expression microarray data (GSE9588) and RNA-seq data from genotype-tissue expression (GTEx), generating robust importance ranking of TF effects and interactions. The SSI-based analysis identified TFs and interacting TF pairs, triplets, and quadruplets involved in CYP3A4 expression. In addition to known CYP3A4 TFs, estrogen receptor α (ESR1) emerges as key TF with the strongest main effect and as the most frequently included TF interacting partner. Model predictions were validated using small interfering RNA (siRNA)/short hairpin RNA (shRNA) gene knockdown and clustered regularly interspaced short palindromic repeats (CRISPR)-mediated transcriptional activation of ESR1 in biliary epithelial Huh7 cells and human hepatocytes in the absence of estrogen. Moreover, ESR1 and known CYP3A4 TFs mutually regulate each other. Detectable in both male and female hepatocytes without added estrogen, the results demonstrate a role for unliganded ESR1 in CYP3A4 expression consistent with unliganded ESR1 signaling reported in other cell types. Added estrogen further enhances ESR1 effects. We propose a hierarchical regulatory network for CYP3A4 expression directed by ESR1 through self-regulation, cross regulation, and TF-TF interactions. We also demonstrate that ESR1 regulates the expression of other P450 enzymes, suggesting broad influence of ESR1 on xenobiotics metabolism in human liver. Further studies are required to understand the mechanisms underlying role of ESR1 in P450 regulation. SIGNIFICANCE STATEMENT: This study focuses on identifying key transcription factors and regulatory networks for CYP3A4, the main drug metabolizing enzymes in liver. We applied a new computational approach (Sobol's sensitivity analysis) to existing hepatic gene expression data to determine the role of transcription factors in regulating CYP3A4 expression, and used molecular genetics methods (siRNA/shRNA gene knockdown and CRISPR-mediated transcriptional activation) to test these interactions in life cells. This approach reveals a robust network of TFs, including their putative interactions and the relative impact of each interaction. We find that ESR1 serves as a key transcription factor function in regulating CYP3A4, and it appears to be acting at least in part in a ligand-free fashion.

Significant association of DRD2 enhancer variant rs12364283 with heroin addiction in a Pakistani population.

The dopamine D2 receptor encoded by DRD2 has been implicated in multiple psychiatric disorders, mediated at least in part by two intronic variants affecting mRNA splicing, rs1076560 and rs2283265, and a less frequent enhancer variant, rs12364283, which increases DRD2 mRNA expression. This study tests whether these functionally validated variants confer susceptibility toward heroin addiction in a Pakistani population. A total of 540 heroin addicts and 467 healthy controls were genotyped, basic allele and genotype tests were performed. Neither rs1076560 nor rs2283265 significantly associated with heroin addiction. The enhancer rs12364283 occurs more frequently in heroin-dependent cases than controls (MAF 13% vs. 7%, respectively), revealing significant association with heroin addiction (p = 3.0E-06, OR 2.1). This study identifies rs12364283 of DRD2 as a potential risk factor for heroin addiction in the Pakistani study population. This enhancer variant had been shown to increase DRD2 mRNA expression, a possible factor in increased vulnerability to heroin addiction. Further studies are needed to validate this association of rs12364283.

CYP2D6 haplotypes with enhancer single-nucleotide polymorphism rs5758550 and rs16947 (*2 allele): implications for CYP2D6 genotyping panels.

INTRODUCTION: CYP2D6 metabolizes ∼25% of all clinically used drugs, with numerous genetic polymorphisms affecting enzyme activity and drug response. Clinical utility of current CYP2D6 genotyping is partially compromised the unresolved complex haplotype structure of the CYP2D6 locus. We have identified a distal enhancer single-nucleotide polymorphism rs5758550 that robustly increases CYP2D6 expression, whereas rs16947 (CYP2D6*2), previously considered inert, reduces correct mRNA splicing and expression, thereby affecting presumed activity of other alleles on the *2 haplotype. OBJECTIVE: This study aims to determine the structure and frequency of haplotypes containing either rs5758550 or rs16947, or both, together with other relevant CYP2D6 alleles, assigning predictive enzyme activity scores to each, and addressing ambiguities in estimating diplotypes in different populations. METHODS: The structure and frequency of haplotypes containing rs5758550 and/or rs16947 in different populations were determined by using phased genotype data from 'The 1000 Genomes Project'. The assigned haplotype-phenotype relationship was tested by associating assigned CYP2D6 activity score with CYP2D6 enzyme activity in a cohort of 122 human liver microsomes. RESULTS: Addition of enhancer single-nucleotide polymorphism rs5758550 and *2 to a CYP2D6 panel improves prediction of CYP2D6 activity. Moreover, the haplotype containing rs5758550 and rs16947 predict extensive CYP2D6 activity more accurately than CYP2D6*2A, a surrogate marker for extensive activity. CONCLUSION: With further studies, the results support possible incorporation of rs5758550 and rs16947 into CYP2D6 biomarker panels for more accurate prediction of CYP2D6 metabolizer status.

Association of Regulatory Genetic Variants for Protein Kinase Cα with Mortality and Drug Efficacy in Patients with Heart Failure.

PURPOSE: Protein kinase C alpha (gene: PRKCA) is a key regulator of cardiac contractility. Two genetic variants have recently been discovered to regulate PRKCA expression in failing human heart tissue (rs9909004 [T → C] and rs9303504 [C → G]). The association of those variants with clinical outcomes in patients with heart failure (HF), and their interaction with HF drug efficacy, is unknown. METHODS: Patients with HF in a prospective registry starting in 2007 were genotyped by whole genome array (n = 951). The primary outcome was all-cause mortality. Cox proportional hazards models adjusted for established clinical risk factors and genomic ancestry tested the independent association of rs9909004 or rs9303504 and the variant interactions with cornerstone HF pharmacotherapies (beta-blockers or angiotensin-converting enzyme inhibitors/angiotensin receptor blockers) in additive genetic models. RESULTS: The minor allele of rs9909004, but not of rs9303504, was independently associated with a decreased risk for all-cause mortality: adjusted HR = 0.81 (95% CI = 0.67-0.98), p = 0.032. The variants did not significantly interact with mortality benefit associated with cornerstone HF pharmacotherapies (p > 0.1 for all). CONCLUSIONS: A recently discovered cardiac-specific regulatory variant for PRKCA (rs9909004) was independently associated with a decreased risk for all-cause mortality in patients with HF. The variant did not interact with mortality benefit associated with cornerstone HF pharmacotherapies.

Combined genetic influence of the nicotinic receptor gene cluster CHRNA5/A3/B4 on nicotine dependence.

BACKGROUND: The CHRNA5/A3/B4 gene locus is associated with nicotine dependence and other smoking related disorders. While the non-synonymous CHRNA5 variant rs16969968 appears to be the main risk factor, linkage disequilibrium (LD) bins in the gene cluster carry frequent variants that regulate expression. Pairwise LD and haplotype analyses had identified at least three haplotype tagging SNPs including rs16969968 as main genetic risk factors. Searching for variants with evidence of regulatory functions, we have reported interactions between CHRNA5 and CHRNA3 enhancer variants (tagged by rs880395 and rs1948, respectively) and rs16969968, forming 3-SNP haplotypes and diplotypes that may more accurately reflect the cluster's combined effects on nicotine dependence (Barrie et al., Hum Mutat 38:112-9, 2017). Here we address further contributions by variants affecting CHRNB4, a possibly limiting component of nicotinic receptors. RESULTS: We identify an LD bin (tagged by rs4887074) associated with expression of CHRNB4. Additive logistic regression models indicate that rs4887074 is associated with nicotine dependence and modulates the effect of rs16969968 in GWAS datasets (COGEND, UW-TTURC, SAGE). 4-SNP haplotype and diplotype analyses (rs880395-rs16969968-rs1948 -rs4887074) yield nicotine dependence risk values that further differentiate those obtained with the 3-SNP model. Moreover, both the main G allele of rs16969968 and the minor G allele of rs4887074 (associated with reduced expression of CHRNB4), residing predominantly on common haplotypes that are protective, represent significant allele-specific variance QTLs, indicating that they interact with each other. CONCLUSIONS: These results indicate rs4887074 is associated with CHRNB4 expression, and along with two regulatory variants of CHRNA3 and CHRNA5, modulates the effect of rs16969968 on nicotine dependence risk. Assignable to individuals because of strong LD structures, 4-SNP haplotypes and diplotypes serve to assess the combined genetic influence of this multi-gene cluster on complex traits, accounting for complex LD relationships and tissue-specific genetic effects (CHRNA5/3) relevant to the traits analyzed. The 4-SNP haplotypes account at least in part for previous tagging SNPs, including the highly GWAS-significant rs6495308, located in a distinct pair-wise LD bin but included in protective 4-SNP haplotypes. Our approach refines and integrates the cluster's overall genetic influence, an important variable when integrating the genetics of multiple genomic loci.

The CHRNA5/CHRNA3/CHRNB4 Nicotinic Receptor Regulome: Genomic Architecture, Regulatory Variants, and Clinical Associations.

Functionally related genes often cluster into a genome region under coordinated regulation, forming a local regulome. To understand regulation of the CHRNA5/CHRNA3/CHRNB4 nicotinic receptor gene cluster, we integrate large-scale RNA expression data (brain and peripheral) from GTEx (Genotype Tissue Expression), clinical associations (GRASP), and linkage disequilibrium data (1000 Genomes) to find candidate SNPs representing independent regulatory variants. CHRNA3, CHRNA5, CHRNB4 mRNAs, and a well-expressed CHRNA5 antisense RNA (RP11-650L12.2) are co-expressed in many human tissues, suggesting common regulatory elements. The CHRNA5 enhancer haplotype tagged by rs880395 not only increases CHRNA5 mRNA expression in all tissues, but also enhances RP11-650L12.2 and CHRNA3 expression, suggesting DNA looping to multiple promoters. However, in nucleus accumbens and putamen, but not other brain regions, CHRNA3 expression associates uniquely with a haplotype tagged by rs1948 (located in the CHRNB4 3'UTR). Haplotype/diplotype analysis of rs880395 and rs1948 plus rs16969968 (a nonsynonymous CHRNA5 risk variant) in GWAS (COGEND, UW-TTURC, SAGE) yields a nicotine dependence risk profile only partially captured by rs16969968 alone. An example of local gene clusters, this nicotinic regulome is controlled by complex genetic variation, with broad implications for interpreting GWAS.

A Novel In Vitro Human Granuloma Model of Sarcoidosis and Latent Tuberculosis Infection.

Many aspects of pathogenic granuloma formation are poorly understood, requiring new relevant laboratory models that represent the complexity (genetics and diversity) of human disease. To address this need, we developed an in vitro model of granuloma formation using human peripheral blood mononuclear cells (PBMCs) derived from patients with active sarcoidosis, latent tuberculosis (TB) infection (LTBI), or normal healthy control subjects. PBMCs were incubated for 7 days with uncoated polystyrene beads or beads coated with purified protein derivative (PPD) or human serum albumin. In response to PPD-coated beads, PBMCs from donors with sarcoidosis and LTBI formed robust multicellular aggregates resembling granulomas, displaying a typical T-helper cell type 1 immune response, as assessed by cytokine analyses. In contrast, minimal PBMC aggregation occurred when control PBMCs were incubated with PPD-coated beads, whereas the response to uncoated beads was negligible in all groups. Sarcoidosis PBMCs responded to human serum albumin-coated beads with modest cellular aggregation and inflammatory cytokine release. Whereas the granuloma-like aggregates formed in response to PPD-coated beads were similar for sarcoidosis and LTBI, molecular profiles differed significantly. mRNA expression patterns revealed distinct pathways engaged in early granuloma formation in sarcoidosis and LTBI, and they resemble molecular patterns reported in diseased human tissues. This novel in vitro human granuloma model is proposed as a tool to investigate mechanisms of early granuloma formation and for preclinical drug discovery research of human granulomatous disorders. Clinical trial registered with www.clinicaltrials.gov (NCT01857401).

Regulatory Variants Modulate Protein Kinase C α (PRKCA) Gene Expression in Human Heart.

PURPOSE: Protein kinase C α (PRKCA) is involved in multiple functions and has been implicated in heart failure risks and treatment outcomes. This study aims to identify regulatory variants affecting PRKCA expression in human heart, and evaluate attributable risk of heart disease. METHODS: mRNA expression quantitative trait loci (eQTLs) were extracted from the Genotype and Tissue Expression Project (GTEx). Allelic mRNA ratios were measured in 51 human heart tissues to identify cis-acting regulatory variants. Potential regulatory regions were tested with luciferase reporter gene assays and further evaluated in GTEx and genome-wide association studies. RESULTS: Located in a region with robust enhancer activity in luciferase reporter assays, rs9909004 (T > C, minor allele frequency =0.47) resides in a haplotype displaying strong eQTLs for PRKCA in heart (p = 1.2 × 10-23). The minor C allele is associated with both decreased PRKCA mRNA expression and decreased risk of phenotypes characteristic of heart failure in GWAS analyses (QT interval p = 3.0 × 10-14). While rs9909004 is the likely regulatory variant, other variants in high linkage disequilibrium cannot be excluded. Distinct regulatory variants appear to affect expression in other tissues. CONCLUSIONS: The haplotype carrying rs9909004 influences PRKCA expression in the heart and is associated with traits linked to heart failure, potentially affecting therapy of heart failure.

CYP2D6 Genetic Variation and Beta-Blocker Maintenance Dose in Patients with Heart Failure.

PURPOSE: This study examined whether a CYP2D6 polymorphism (CYP2D6*4) was related to beta-blocker maintenance dose in patients with heart failure. METHODS: Logistic regression modeling was utilized in a retrospective chart-review analysis of heart-failure patients (60% Male, 90% of European descent) to assess whether CYP2D6*4 (non-functional CYP2D6 allele present in 1 of 5 individuals of European descent) is associated with maintenance dose of carvedilol (n = 65) or metoprolol (n = 33). RESULTS: CYP2D6*4 was associated with lower maintenance dose of metoprolol (OR 0.13 [95% CI 0.02-0.75] p = 0.023), and a trend was observed between CYP2D6*4 and higher maintenance dose of carvedilol (OR 2.94 [95% CI 0.84-10.30] p = 0.093). None of the patients that carried CYP2D6*4 achieved the recommended target dose of metoprolol (200 mg/day). CONCLUSION: Consistent with the role of CYP2D6 in the metabolism of metoprolol, the tolerated maintenance dose of metoprolol was lower in CYP2D6*4 carriers compared to non-carriers. Consistent with the role of CYP2D6 in activation of carvedilol, tolerated maintenance dose of carvedilol was higher in CYP2D6*4 carriers compared to non-carriers. Further investigation is warranted to ascertain the potential of CYP2D6 as a potential predictive biomarker of beta-blocker maintenance dose in heart failure patients.

Elemental Ingredients in the Macrophage Cocktail: Role of ZIP8 in Host Response to Mycobacterium tuberculosis.

Tuberculosis (TB) is a global epidemic caused by the infection of human macrophages with the world's most deadly single bacterial pathogen, Mycobacterium tuberculosis (M.tb). M.tb resides in a phagosomal niche within macrophages, where trace element concentrations impact the immune response, bacterial metal metabolism, and bacterial survival. The manipulation of micronutrients is a critical mechanism of host defense against infection. In particular, the human zinc transporter Zrt-/Irt-like protein 8 (ZIP8), one of 14 ZIP family members, is important in the flux of divalent cations, including zinc, into the cytoplasm of macrophages. It also has been observed to exist on the membrane of cellular organelles, where it can serve as an efflux pump that transports zinc into the cytosol. ZIP8 is highly inducible in response to M.tb infection of macrophages, and we have observed its localization to the M.tb phagosome. The expression, localization, and function of ZIP8 and other divalent cation transporters within macrophages have important implications for TB prevention and dissemination and warrant further study. In particular, given the importance of zinc as an essential nutrient required for humans and M.tb, it is not yet clear whether ZIP-guided zinc transport serves as a host protective factor or, rather, is targeted by M.tb to enable its phagosomal survival.

Non-linear interactions between candidate genes of myocardial infarction revealed in mRNA expression profiles.

BACKGROUND: Alterations in gene expression are key events in disease etiology and risk. Poor reproducibility in detecting differentially expressed genes across studies suggests individual genes may not be sufficiently informative for complex diseases, such as myocardial infarction (MI). Rather, dysregulation of the 'molecular network' may be critical for pathogenic processes. Such a dynamic network can be built from pairwise non-linear interactions. RESULTS: We investigate non-linear interactions represented in mRNA expression profiles that integrate genetic background and environmental factors. Using logistic regression, we test the association of individual GWAS-based candidate genes and non-linear interaction terms (between these mRNA expression levels) with MI. Based on microarray data in CATHGEN (CATHeterization in GENetics) and FHS (Framingham Heart Study), we find individual genes and pairs of mRNAs, encoded by 41 MI candidate genes, with significant interaction terms in the logistic regression model. Two pairs replicate between CATHGEN and FHS (CNNM2|GUCY1A3 and CNNM2|ZEB2). Analysis of RNAseq data from GTEx (Genotype-Tissue Expression) shows that 20 % of these disease-associated RNA pairs are co-expressed, further prioritizing significant interactions. Because edges in sparse co-expression networks formed solely by the 41 candidate genes are unlikely to represent direct physical interactions, we identify additional RNAs as links between network pairs of candidate genes. This approach reveals additional mRNAs and interaction terms significant in the context of MI, for example, the path CNNM2|ACSL5|SCARF1|GUCY1A3, characterized by the common themes of magnesium and lipid processing. CONCLUSIONS: The results of this study support a role for non-linear interactions between genes in MI and provide a basis for further study of MI systems biology. mRNA expression profiles encoded by a limited number of candidate genes yield sparse networks of MI-relevant interactions that can be expanded to include additional candidates by co-expression analysis. The non-linear interactions observed here inform our understanding of the clinical relevance of gene-gene interactions in the pathophysiology of MI, while providing a new strategy in developing clinical biomarker panels.

Common CYP2D6 polymorphisms affecting alternative splicing and transcription: long-range haplotypes with two regulatory variants modulate CYP2D6 activity.

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 25% of clinically used drugs. Genetic polymorphisms cause substantial variation in CYP2D6 activity and serve as biomarkers guiding drug therapy. However, genotype-phenotype relationships remain ambiguous except for poor metabolizers carrying null alleles, suggesting the presence of yet unknown genetic variants. Searching for regulatory CYP2D6 polymorphisms, we find that a SNP defining the CYP2D6*2 allele, rs16947 [R296C, 17-60% minor allele frequency (MAF)], previously thought to convey normal activity, alters exon 6 splicing, thereby reducing CYP2D6 expression at least 2-fold. In addition, two completely linked SNPs (rs5758550/rs133333, MAF 13-42%) increase CYP2D6 transcription more than 2-fold, located in a distant downstream enhancer region (>100 kb) that interacts with the CYP2D6 promoter. In high linkage disequilibrium (LD) with each other, rs16947 and the enhancer SNPs form haplotypes that affect CYP2D6 enzyme activity in vivo. In a pediatric cohort of 164 individuals, rs16947 alone (minor haplotype frequency 28%) was associated with reduced CYP2D6 metabolic activity (measured as dextromethorphan/metabolite ratios), whereas rs5758550/rs133333 alone (frequency 3%) resulted in increased CYP2D6 activity, while haplotypes containing both rs16947 and rs5758550/rs133333 were similar to the wild-type. Other alleles used in biomarker panels carrying these variants such as CYP2D6*41 require re-evaluation of independent effects on CYP2D6 activity. The occurrence of two regulatory variants of high frequency and in high LD, residing on a long haplotype, highlights the importance of gene architecture, likely shaped by evolutionary selection pressures, in determining activity of encoded proteins.

CYP3A4 intronic SNP rs35599367 (CYP3A4*22) alters RNA splicing.

Cytochrome P450 3A4 (CYP3A4) metabolizes 30-50% of clinically used drugs. Large interperson variability in CYP3A4 activity affects response to CYP3A4 substrate drugs. We had demonstrated that an intronic single nucleotide polymorphism rs35599367 (CYP3A4*22, located in intron 6) reduces mRNA/protein expression; however, the underlying mechanism remained unknown. Here we show that CYP3A4*22 is associated with a two-fold or greater increase in formation of a nonfunctional CYP3A4 alternative splice variant with partial intron 6 retention in human liver (P=0.006), but not in small intestines. Consistent with this observation, in-vitro transfection experiments with a CYP3A4 minigene (spanning from intron 5 to intron 7) demonstrated that plasmids carrying the rs35599367 minor T allele caused significantly greater intron 6 retention than the C allele in liver derived HepG2 cells, but not in intestine-derived LS-174T cells. These results indicate that tissue-specific increased formation of nonfunctional alternative splice variant causes reduced CYP3A4 mRNA/protein expression in CYP3A4*22 carriers.