A program for real-time surveillance of SARS-CoV-2 genetics.

The COVID-19 pandemic brought forth an urgent need for widespread genomic surveillance for rapid detection and monitoring of emerging SARS-CoV-2 variants. It necessitated design, development, and deployment of a nationwide infrastructure designed for sequestration, consolidation, and characterization of patient samples that disseminates de-identified information to public authorities in tight turnaround times. Here, we describe our development of such an infrastructure, which sequenced 594,832 high coverage SARS-CoV-2 genomes from isolates we collected in the United States (U.S.) from March 13th 2020 to July 3rd 2023. Our sequencing protocol ('Virseq') utilizes wet and dry lab procedures to generate mutation-resistant sequencing of the entire SARS-CoV-2 genome, capturing all major lineages. We also characterize 379 clinically relevant SARS-CoV-2 multi-strain co-infections and ensure robust detection of emerging lineages via simulation. The modular infrastructure, sequencing, and analysis capabilities we describe support the U.S. Centers for Disease Control and Prevention national surveillance program and serve as a model for rapid response to emerging pandemics at a national scale.

Meta-analysis of genome-wide association studies of stable warfarin dose in patients of African ancestry.

Warfarin dose requirements are highly variable due to clinical and genetic factors. While genetic variants influencing warfarin dose have been identified in European and East Asian populations, more work is needed to identify African-specific genetic variants to help optimize warfarin dosing. We performed genome-wide association studies (GWAS) in four African cohorts from Uganda, South Africa, and Zimbabwe, totalling 989 warfarin-treated participants who reached stable dose and had international normalized ratios within therapeutic ranges. We also included two African American cohorts recruited by the International Warfarin Pharmacogenetics Consortium (n=316) and the University of Alabama at Birmingham (n=199). Following the GWAS, we performed standard error-weighted meta-analyses and then conducted stepwise conditional analyses to account for known loci (the CYP2C cluster SNP rs12777823 and CYP2C9 in chromosome 10; VKORC1 in chromosome 16). The genome-wide significance threshold was set at P<5×10-8. The meta-analysis, comprising 1,504 participants identified 242 significant SNPs across three genomic loci, with 99.6% of these located within known loci on chromosomes 10 (top SNP: rs58800757, P=4.27×10-13) and 16 (top SNP: rs9925964, P=9.97×10-16). Adjustment for the VKORC1 SNP -1639G>A revealed an additional locus on chromosome 2 (top SNPs rs116057875/rs115254730/rs115240773, P=3.64×10-8), implicating the MALL gene, that could indirectly influence warfarin response through interactions with caveolin-1. In conclusion, our meta-analysis of six cohorts of warfarin-treated patients of African ancestry reaffirmed the importance of CYP2C9 and VKORC1 in influencing warfarin dose requirements. We also identified a new locus (MALL), that still requires direct evidence of biological plausibility.

Jellyfish Venom Peptides Targeting Human Potassium Channels Identified through Ligand Screening: Morphometric and Molecular Identification of the Species and Antibiotic Potential.

The relative lack of marine venom could be attributed to the difficulty in dealing with venomous marine animals. Moreover, the venom of marine animals consists of various bioactive molecules, many of which are proteins with unique properties. In this study, we investigated the potential toxic proteins of jellyfish collected for ligand screening to understand the mechanism of the toxic effects of jellyfish. Since taxonomic identification is problematic due to the lack of proper keys, we conducted morphological and molecular mitochondrial DNA sequencing from COI and ITS regions. The venom extract from nematocysts found along the bell margins was used for protein characterization using SDS-gel electrophoresis and nano-liquid chromatography-tandem mass spectrometry. Ligand screening for the most potent toxin and antibacterial and cytotoxicity assays were carried out. The phylogenetic tree showed distinct clustering from other Catostylus sp. The proteomic analysis revealed venom with many bioactive proteins. Only 13 venom proteins were identified with molecular weights ranging from 4318 to 184,923 Da, exhibiting the venom's complexity. The overall toxin protein composition of Catostylus sp. venom was dominated by potassium channel toxin alpha-KTx. Molecular docking of toxin alpha-KTx 1.13 revealed high specificity towards the human voltage-gated potassium channel Kv3 with a high fitness score and a minimum energy barrier of -17.9 kcal/mol. Disc diffusion and cytotoxicity assays revealed potent antibacterial activity against Klebsiella pneumoniae with no cytotoxicity. Further studies on detailed characterization and therapeutic potentials are warranted.

LA-GEM: imputation of gene expression with incorporation of Local Ancestry.

Gene imputation and TWAS have become a staple in the genomics medicine discovery space; helping to identify genes whose regulation effects may contribute to disease susceptibility. However, the cohorts on which these methods are built are overwhelmingly of European Ancestry. This means that the unique regulatory variation that exist in non-European populations, specifically African Ancestry populations, may not be included in the current models. Moreover, African Americans are an admixed population, with a mix of European and African segments within their genome. No gene imputation model thus far has incorporated the effect of local ancestry (LA) on gene expression imputation. As such, we created LA-GEM which was trained and tested on a cohort of 60 African American hepatocyte primary cultures. Uniquely, LA-GEM include local ancestry inference in its prediction of gene expression. We compared the performance of LA-GEM to PrediXcan trained the same dataset (with no inclusion of local ancestry) We were able to reliably predict the expression of 2559 genes (1326 in LA-GEM and 1236 in PrediXcan). Of these, 546 genes were unique to LA-GEM, including the CYP3A5 gene which is critical to drug metabolism. We conducted TWAS analysis on two African American clinical cohorts with pharmacogenomics phenotypic information to identity novel gene associations. In our IWPC warfarin cohort, we identified 17 transcriptome-wide significant hits. No gene reached are prespecified significance level in the clopidogrel cohort. We did see suggestive association with RAS3A to P2RY12 Reactivity Units (PRU), a clinical measure of response to anti-platelet therapy. This method demonstrated the need for the incorporation of LA into study in admixed populations.

Investigation of genomic and transcriptomic risk factors in clopidogrel response in African Americans.

Clopidogrel, an anti-platelet drug, used to prevent thrombosis after percutaneous coronary intervention. Clopidogrel resistance results in recurring ischemic episodes, with African Americans suffering disproportionately. The aim of this study was to identify biomarkers of clopidogrel resistance in African American patients. We conducted a genome-wide association study, including local ancestry adjustment, in 141 African Americans on clopidogrel to identify associations with high on-treatment platelet reactivity (HTPR). We validated genome-wide and suggestive hits in an independent cohort of African American clopidogrel patients (N = 823) from the Million Veteran's Program (MVP) along with in vitro functional follow up. We performed differential gene expression (DGE) analysis in whole blood with functional follow-up in MEG-01 cells. We identified rs7807369, within thrombospondin 7A (THSD7A), as significantly associated with increasing risk of HTPR (p = 4.56 × 10-9). Higher THSD7A expression was associated with HTPR in an independent gene expression cohort of clopidogrel treated patients (p = 0.004) and supported by increased gene expression on THSD7A in primary human endothelial cells carrying the risk haplotype. Two SNPs (rs1149515 and rs191786) were validated in the MVP cohort. DGE analysis identified an association with decreased LAIR1 expression to HTPR. LAIR1 knockdown in a MEG-01 cells resulted in increased expression of SYK and AKT1, suggesting an inhibitory role of LAIR1 in the Glycoprotein VI pathway. Notably, the CYP2C19 variants showed no association with clopidogrel response in the discovery or MVP cohorts. In summary, these finding suggest that other variants outside of CYP2C19 star alleles play an important role in clopidogrel response in African Americans.

Discovery of Ancestry-specific Variants Associated with Clopidogrel Response among Caribbean Hispanics.

Background: High on-treatment platelet reactivity (HTPR) with clopidogrel is predictive of ischemic events in adults with coronary artery disease. Despite strong data suggesting HTPR varies with ethnicity, including clinical and genetic variables, no genome-wide association study (GWAS) of clopidogrel response has been performed among Caribbean Hispanics. This study aimed to identify genetic predictors of HTPR in a cohort of Caribbean Hispanic cardiovascular patients from Puerto Rico. Methods: Local Ancestry inference (LAI) and traditional GWASs were performed on a cohort of 511 clopidogrel-treated patients, stratified based on their P2Y12 reaction units (PRU) into responders and non-responders (HTPR). Results: The LAI GWAS identified variants within the CYP2C19 region associated with HTPR, predominantly driven by individuals of European ancestry and absent in those with native ancestry. Incorporating local ancestry adjustment notably enhanced our ability to detect associations. While no loci reached traditional GWAS significance, three variants showed suggestive significance at chromosomes 3, 14 and 22 (OSBPL10 rs1376606, DERL3 rs5030613, and RGS6 rs9323567). In addition, a variant in the UNC5C gene on chromosome 4 was associated with an increased risk of HTPR. These findings were not identified in other cohorts, highlighting the unique genetic landscape of Caribbean Hispanics. Conclusion: This is the first GWAS of clopidogrel response in Hispanics, confirming the relevance of the CYP2C19 cluster, particularly among those with European ancestry, and also identifying novel markers in a diverse patient population. Further studies are warranted to replicate our findings in other diverse cohorts and meta-analyses.

Molecular Identification and Evolutionary Divergence of the Sri Lankan Sambar Deer, Rusa unicolor (Kerr 1792).

The Sambar is one of the largest deer species distributed mainly in Asia, and it has been listed as a vulnerable species. Taxonomy based on morphological characterization has been the gold standard method used to identify the Sambar deer species. Yet, morphological identification is challenging and requires expertise. To conduct species identification and taxonomic decisions, we performed the molecular identification of R. unicolor found in Sri Lanka using DNA barcodes, COI, and Cyt b to compare the Sri Lankan R. unicolor with the Indian R. unicolor and other R. unicolor subspecies. We obtained mitochondrial DNA sequences from COI and Cyt b from blood samples collected from the wet zone in Sri Lanka. A phylogenetic tree was constructed based on the Bayesian analyses using MrBayes 3.2.7. Molecular dating was implemented in Bayesian Evolutionary Analysis Sampling Trees (BEAST v1.8.2) on the concatenated sequence using a log-normal relaxed clock and Yule species tree prior, with four categories. The results showed that the Sri Lankan R. unicolor is genetically different from the Indian R. unicolor and other R. unicolor subspecies. The divergence occurred approximately 1.1 MYA (million years ago) in the Pleistocene era. The results are essential for designing new conservation platforms for these Sambar deer species.

Prostate cancer cell-platelet bidirectional signaling promotes calcium mobilization, invasion and apoptotic resistance via distinct receptor-ligand pairs.

Platelets play a crucial role in cancer and thrombosis. However, the receptor-ligand repertoire mediating prostate cancer (PCa) cell-platelet interactions and ensuing consequences have not been fully elucidated. Microvilli emanating from the plasma membrane of PCa cell lines (RC77 T/E, MDA PCa 2b) directly contacted individual platelets and platelet aggregates. PCa cell-platelet interactions were associated with calcium mobilization in platelets, and translocation of P-selectin and integrin αIIbß3 onto the platelet surface. PCa cell-platelet interactions reciprocally promoted PCa cell invasion and apoptotic resistance, and these events were insensitive to androgen receptor blockade by bicalutamide. PCa cells were exceedingly sensitive to activation by platelets in vitro, occurring at a PCa cell:platelet coculture ratio as low as 1:10 (whereas PCa patient blood contains 1:2,000,000 per ml). Conditioned medium from cocultures stimulated PCa cell invasion but not apoptotic resistance nor platelet aggregation. Candidate transmembrane signaling proteins responsible for PCa cell-platelet oncogenic events were identified by RNA-Seq and broadly divided into 4 major categories: (1) integrin-ligand, (2) EPH receptor-ephrin, (3) immune checkpoint receptor-ligand, and (4) miscellaneous receptor-ligand interactions. Based on antibody neutralization and small molecule inhibitor assays, PCa cell-stimulated calcium mobilization in platelets was found to be mediated by a fibronectin1 (FN1)-αIIbß3 signaling axis. Platelet-stimulated PCa cell invasion was facilitated by a CD55-adhesion G protein coupled receptor E5 (ADGRE5) axis, with contribution from platelet cytokines CCL3L1 and IL32. Platelet-stimulated PCa cell apoptotic resistance relied on ephrin-EPH receptor and lysophosphatidic acid (LPA)-LPA receptor (LPAR) signaling. Of participating signaling partners, FN1 and LPAR3 overexpression was observed in PCa specimens compared to normal prostate, while high expression of CCR1 (CCL3L1 receptor), EPHA1 and LPAR5 in PCa was associated with poor patient survival. These findings emphasize that non-overlapping receptor-ligand pairs participate in oncogenesis and thrombosis, highlighting the complexity of any contemplated clinical intervention strategy.

Implementation of pharmacogenomics into inpatient general medicine.

Pharmacogenomics is a crucial piece of personalized medicine. Preemptive pharmacogenomic testing is only used sparsely in the inpatient setting and there are few models to date for fostering the adoption of pharmacogenomic treatment in the inpatient setting. We created a multi-institutional project in Chicago to enable the translation of pharmacogenomics into inpatient practice. We are reporting our implementation process and barriers we encountered with solutions. This study, 'Implementation of Point-of-Care Pharmacogenomic Decision Support Accounting for Minority Disparities', sought to implement pharmacogenomics into inpatient practice at three sites: The University of Chicago, Northwestern Memorial Hospital, and the University of Illinois at Chicago. This study involved enrolling African American adult patients for preemptive genotyping across a panel of actionable germline variants predicting drug response or toxicity risk. We report our approach to implementation and the barriers we encountered engaging hospitalists and general medical providers in the inpatient pharmacogenomic intervention. Our strategies included: a streamlined delivery system for pharmacogenomic information, attendance at hospital medicine section meetings, use of physician and pharmacist champions, focus on hospitalists' care and optimizing system function to fit their workflow, hand-offs, and dealing with hospitalists turnover. Our work provides insights into strategies for the initial engagement of inpatient general medicine providers that we hope will benefit other institutions seeking to implement pharmacogenomics in the inpatient setting.

Leveraging drug perturbation to reveal genetic regulators of hepatic gene expression in African Americans.

Expression quantitative locus (eQTL) studies have paved the way in identifying genetic variation impacting gene expression levels. African Americans (AAs) are disproportionately underrepresented in eQTL studies, resulting in a lack of power to identify population-specific regulatory variants especially related to drug response. Specific drugs are known to affect the biosynthesis of drug metabolism enzymes as well as other genes. We used drug perturbation in cultured primary hepatocytes derived from AAs to determine the effect of drug treatment on eQTL mapping and to identify the drug response eQTLs (reQTLs) that show altered effect size following drug treatment. Whole-genome genotyping (Illumina MEGA array) and RNA sequencing were performed on 60 primary hepatocyte cultures after treatment with six drugs (Rifampin, Phenytoin, Carbamazepine, Dexamethasone, Phenobarbital, and Omeprazole) and at baseline (no treatment). eQTLs were mapped by treatment and jointly with Meta-Tissue. We found varying transcriptional changes across different drug treatments and identified Nrf2 as a potential general transcriptional regulator. We jointly mapped eQTLs with gene expression data across all drug treatments and baseline, which increased our power to detect eQTLs by 2.7-fold. We also identified 2,988 reQTLs (eQTLs with altered effect size after drug treatment). reQTLs were more likely to overlap transcription factor binding sites, and we uncovered reQTLs for drug metabolizing genes such as CYP3A5. Our results provide insights into the genetic regulation of gene expression in hepatocytes through drug perturbation and provide insight into SNPs that effect the liver's ability to respond to transcription upregulation.

The Role of Global and Local Ancestry on Clopidogrel Response in African Americans.

Pharmacogenomics has long lacked dedicated studies in African Americans, resulting in a lack of indepth data in this populations. The ACCOuNT consortium has collected a cohort of 167 African American patients on steady state clopidogrel with the goal of discovering population specific variation that may contribute to the response of this anti-platelet agent. Here we analyze the role of both global and local ancestry on the clinical phenotypes of P2Y12 reaction units (PRU) and high on-treatment platelet reactivity (HTPR) in this cohort. We found that local ancestry at the TSS of three genes, IRS-1, ABCB1 and KDR were nominally associated with PRU, and local ancestry-adjusted SNP association identified variants in ITGA2 associated to increased PRU. These finding help to explain the variability in drug response seen in African Americans, especially as few studies on genes outside of CYP2C19 has been conducted in this population.

Multi-Omics Studies in Historically Excluded Populations: The Road to Equity.

Over the past few decades, genomewide association studies (GWASs) have identified the specific genetics variants contributing to many complex diseases by testing millions of genetic variations across the human genome against a variety of phenotypes. However, GWASs are limited in their ability to uncover mechanistic insight given that most significant associations are found in non-coding region of the genome. Furthermore, the lack of diversity in studies has stymied the advance of precision medicine for many historically excluded populations. In this review, we summarize most popular multi-omics approaches (genomics, transcriptomics, proteomics, and metabolomics) related to precision medicine and highlight if diverse populations have been included and how their findings have advance biological understanding of disease and drug response. New methods that incorporate local ancestry have been to improve the power of GWASs for admixed populations (such as African Americans and Latinx). Because most signals from GWAS are in the non-coding region, other machine learning and omics approaches have been developed to identify the potential causative single-nucleotide polymorphisms and genes that explain these phenotypes. These include polygenic risk scores, expression quantitative trait locus mapping, and transcriptome-wide association studies. Analogous protein methods, such as proteins quantitative trait locus mapping, proteome-wide association studies, and metabolomic approaches provide insight into the consequences of genetic variation on protein abundance. Whereas, integrated multi-omics studies have improved our understanding of the mechanisms for genetic association, we still lack the datasets and cohorts for historically excluded populations to provide equity in precision medicine and pharmacogenomics.

Genomewide Association Study Identifies Copy Number Variants Associated With Warfarin Dose Response and Risk of Venous Thromboembolism in African Americans.

The anticoagulant warfarin is commonly used to control and prevent thrombotic disorders, such as venous thromboembolism (VTE), which disproportionately afflicts African Americans. Despite the importance of copy number variants (CNVs), few studies have focused on characterizing and understanding their role in drug response and disease risk among African Americans. In this study, we conduct the first genome-wide analysis of CNVs to more comprehensively account for the contribution of genetic variation in warfarin dose requirement and VTE risk among African Americans. We used hidden Markov models to detect CNVs from high-throughput single-nucleotide polymorphism arrays for 340 African American participants in the International Warfarin Pharmacogenetics Consortium. We identified 11,570 CNVs resulting in 2,038 copy number variable regions (CNVRs) and found 3 CNVRs associated with warfarin dose requirement and 3 CNVRs associated with VTE risk in African Americans. CNVRs 1q31.2del and 6q14.1del were associated with increased warfarin dose requirement (ß = 11.18 and 4.94, respectively; Pemp  = < 0.002); CNVR 19p13.31del was associated with decreased warfarin dose requirement (ß = -1.41, Pemp  = 0.0004); CNVRs (2p22.1del and 5q35.1-q35.2del) were found to be associated with increased risk of VTE (odds ratios (ORs) = 1.88 and 14.9, respectively; Pemp ≤0.02); and CNVR 10q26.12del was associated with a decreased risk of VTE (OR = 0.6; Pemp  = 0.05). Modeling of the 10q26.12del in HepG2 cells revealed that this deletion results in decreased fibrinogen gene expression, decreased fibrinogen and WDR11 protein levels, and decreased secretion of fibrinogen into the extracellular matrix. We found robust evidence that CNVRs could contribute to warfarin dose requirement and risk of VTE in African Americans and for 10q26.3del describe a possible pathogenic mechanism.

Local Ancestry-Informed Candidate Pathway Analysis of Warfarin Stable Dose in Latino Populations.

Accuracy of warfarin dose prediction algorithms may be improved by including data from diverse populations in genetic studies of dose variability. Here, we surveyed single nucleotide polymorphisms in vitamin K-related genetic pathways for association with warfarin dose requirements in two admixed Latino populations in standard-principal component adjusted and contemporary-local ancestry adjusted regression models. A total of five variants from vitamin K-related genes/pathways were associated with warfarin dose in both cohorts (P < 0.0125) in standard models. Local ancestry-adjusted analysis unveiled 35 associated variants with absolute effects ranging from ß = 9.04 ( ±2.23) to 39.18 ( ±10.89) per ancestral allele in the discovery cohort and ß = 6.47 (± 2.02) to 17.82 (± 6.83) in the replication cohort. Importantly, we demonstrate the technical validity of the Tractor model in cohorts with admixed ancestry from three founder populations and bring attention to the technical hurdles obstructing the inclusion of diverse, especially admixed, populations in pharmacogenomic research.

Warfarin Dosing in Patients With CYP2C9*5 Variant Alleles.

Pharmacogenetic dosing improves the accuracy of warfarin dosing, but current pharmacogenetic dosing algorithms are less accurate in populations of African ancestry. The cytochrome P450 2C9*5 (CYP2C9*5) allele is found almost exclusively in populations of African ancestry, and in vitro studies suggest CYP2C9*5 is associated with reduced clearance of warfarin. The clinical relevance of this single-nucleotide variation (SNV) (formerly SNP) is uncertain. In this multicentered study of 2,298 patients (49% female, 35% Black) taking warfarin, we quantified the association between the CYP2C9*5 allele and warfarin requirements. The CYP2C9*5 SNV was present in 2.3% of Black and 0.07% of White patients. Without taking CYP2C9*5 into account, pharmacogenetic algorithms that include other SNVs overestimated the warfarin dose by 30% (95% confidence interval (19-40%), P < 0.001), an average of 1.87 mg/day (SD 1.64) in heterozygotes (P < 0.001). Noncarriers required a slightly (0.23 mg/day, SD 2.09) higher than predicted dose. Genotyping for CYP2C9*5 corrected the potential overdose and halved overall dosing error in heterozygotes. Patients carrying CYP2C9*5 require a clinically relevant reduction in warfarin dose. Given the potential to improve the accuracy and safety of warfarin dosing in populations of African ancestry, we have incorporated this SNV into a nonprofit website to assist warfarin initiation (www.WarfarinDosing.org).

Apixaban Concentrations in Routine Clinical Care of Older Adults With Nonvalvular Atrial Fibrillation.

BACKGROUND: Direct-acting oral anticoagulants are first-line agents for prevention of stroke in patients with nonvalvular atrial fibrillation (NVAF), but data are limited for the oldest patients, and with reduced dosing. OBJECTIVES: To determine steady-state apixaban peak and trough concentrations during routine care of older adults with NVAF, compare concentrations to clinical trial concentrations, and explore factors associated with concentrations. METHODS: A cross-sectional study of medically stable older adults with NVAF (≥75 years or ≥70 years if Black) receiving apixaban. Peak (2-4.4 hours post-dose) and trough (before next dose) concentrations were determined by anti-Xa activity calibrated chromogenic assay. Patient characteristics associated with concentrations were determined by multivariate modeling. RESULTS: The median age of patients (n = 115) was 80 (interquartile range: 77-84) years. The cohort comprised 46 women and 69 men; of which 98 are White, 11 Black, and 6 Asian. With 5 mg twice daily per labelling (n = 88), peak concentrations were higher in women: 248 ± 105 vs 174 ± 67 ng/mL in men (P < 0.001) and exceeded expected 95% range in 6 of 30 vs 0 of 55 men (P = 0.002). With 2.5 mg twice daily per label (n = 11), concentrations were <5 mg twice daily (peak: 136 ± 87 vs 201 ± 90 ng/mL, P = 0.026; trough: 65 ± 28 vs 109 ± 56 ng/mL, P < 0.001), but not different than 2.5 mg twice daily without reduction criteria (n = 13; peak: 132 ± 88; trough: 65 ± 31 ng/mL). Covariates associated with concentrations included sex, number of daily medications, and creatinine clearance. CONCLUSIONS: Older women had higher than expected peak apixaban concentrations, and 2.5 mg twice daily produced lower concentrations than standard dosing. Factors not currently included in dosing recommendations affected concentrations. The impact of apixaban concentrations on outcomes needs evaluation.

Genetic association of primary nonresponse to anti-TNFα therapy in patients with inflammatory bowel disease.

OBJECTIVES: Primary nonresponse (PNR) to antitumor necrosis factor-α (TNFα) biologics is a serious concern in patients with inflammatory bowel disease (IBD). We aimed to identify the genetic variants associated with PNR. PATIENTS AND METHODS: Patients were recruited from outpatient GI clinics and PNR was determined using both clinical and endoscopic findings. A case-control genome-wide association study was performed in 589 IBD patients and associations were replicated in an independent cohort of 293 patients. Effect of the associated variant on gene expression and TNFα secretion was assessed by cell-based assays. Pleiotropic effects were investigated by Phenome-wide association study (PheWAS). RESULTS: We identified rs34767465 as associated with PNR to anti-TNFα therapy (odds ratio: 2.07, 95% CI, 1.46-2.94, P = 2.43 × 10-7, [replication odds ratio: 1.8, 95% CI, 1.04-3.16, P = 0.03]). rs34767465 is a multiple-tissue expression quantitative trait loci for FAM114A2. Using RNA-sequencing and protein quantification from HapMap lymphoblastoid cell lines (LCLs), we found a significant decrease in FAM114A2 mRNA and protein expression in both heterozygous and homozygous genotypes when compared to wild type LCLs. TNFα secretion was significantly higher in THP-1 cells [differentiated into macrophages] with FAM114A2 knockdown versus controls. Immunoblotting experiments showed that depletion of FAM114A2 impaired autophagy-related pathway genes suggesting autophagy-mediated TNFα secretion as a potential mechanism. PheWAS showed rs34767465 was associated with comorbid conditions found in IBD patients (derangement of joints [P = 3.7 × 10-4], pigmentary iris degeneration [P = 5.9 × 10-4], diverticulum of esophagus [P = 7 × 10-4]). CONCLUSIONS: We identified a variant rs34767465 associated with PNR to anti-TNFα biologics, which increases TNFα secretion through mechanism related to autophagy. rs34767465 may also explain the comorbidities associated with IBD.

Applicability of Pharmacogenomically Guided Medication Treatment during Hospitalization of At-Risk Minority Patients.

Known disparities exist in the availability of pharmacogenomic information for minority populations, amplifying uncertainty around clinical utility for these groups. We conducted a multi-site inpatient pharmacogenomic implementation program among self-identified African-Americans (AA; n = 135) with numerous rehospitalizations (n = 341) from 2017 to 2020 (NIH-funded ACCOuNT project/clinicaltrials.gov#NCT03225820). We evaluated the point-of-care availability of patient pharmacogenomic results to healthcare providers via an electronic clinical decision support tool. Among newly added medications during hospitalizations and at discharge, we examined the most frequently utilized medications with associated pharmacogenomic results. The population was predominantly female (61%) with a mean age of 53 years (range 19-86). On average, six medications were newly prescribed during each individual hospital admission. For 48% of all hospitalizations, clinical pharmacogenomic information was applicable to at least one newly prescribed medication. Most results indicated genomic favorability, although nearly 29% of newly prescribed medications indicated increased genomic caution (increase in toxicity risk/suboptimal response). More than one of every five medications prescribed to AA patients at hospital discharge were associated with cautionary pharmacogenomic results (most commonly pantoprazole/suboptimal antacid effect). Notably, high-risk pharmacogenomic results (genomic contraindication) were exceedingly rare. We conclude that the applicability of pharmacogenomic information during hospitalizations for vulnerable populations at-risk for experiencing health disparities is substantial and warrants continued prospective investigation.

Differences in the Platelet mRNA Landscape Portend Racial Disparities in Platelet Function and Suggest Novel Therapeutic Targets.

The African American (AA) population displays a 1.6 to 3-fold higher incidence of thrombosis and stroke mortality compared with European Americans (EAs). Current antiplatelet therapies target the ADP-mediated signaling pathway, which displays significant pharmacogenetic variation for platelet reactivity. The focus of this study was to define underlying population differences in platelet function in an effort to identify novel molecular targets for future antiplatelet therapy. We performed deep coverage RNA-Seq to compare gene expression levels in platelets derived from a cohort of healthy volunteers defined by ancestry determination. We identified > 13,000 expressed platelet genes of which 480 were significantly differentially expressed genes (DEGs) between AAs and EAs. DEGs encoding proteins known or predicted to modulate platelet aggregation, morphology, or platelet count were upregulated in AA platelets. Numerous G-protein coupled receptors, ion channels, and pro-inflammatory cytokines not previously associated with platelet function were likewise differentially expressed. Many of the signaling proteins represent potential pharmacologic targets of intervention. Notably, we confirmed the differential expression of cytokines IL32 and PROK2 in an independent cohort by quantitative real-time polymerase chain reaction, and provide functional validation of the opposing actions of these two cytokines on collagen-induced AA platelet aggregation. Using Genotype-Tissue Expression whole blood data, we identified 516 expression quantitative trait locuses with Fst values > 0.25, suggesting that population-differentiated alleles may contribute to differences in gene expression. This study identifies gene expression differences at the population level that may affect platelet function and serve as potential biomarkers to identify cardiovascular disease risk. Additionally, our analysis uncovers candidate novel druggable targets for future antiplatelet therapies.

Incorporation of DNA methylation into eQTL mapping in African Americans.

Epigenetics is a reversible molecular mechanism that plays a critical role in many developmental, adaptive, and disease processes. DNA methylation has been shown to regulate gene expression and the advent of high throughput technologies has made genome-wide DNA methylation analysis possible. We investigated the effect of DNA methylation on eQTL mapping (methylation-adjusted eQTLs), by incorporating DNA methylation as a SNP-based covariate in eQTL mapping in African American derived hepatocytes. We found that the addition of DNA methylation uncovered new eQTLs and eGenes. Previously discovered eQTLs were significantly altered by the addition of DNA methylation data suggesting that methylation may modulate the association of SNPs to gene expression. We found that methylation-adjusted eQTLs that were less significant compared to PC-adjusted eQTLs were enriched in lipoprotein measurements (FDR=0.0040), immune system disorders (FDR = 0.0042), and liver enzyme measurements (FDR=0.047), suggesting that DNA methylation modulates the genetic regulation of these phenotypes. Our methylation-adjusted eQTL analysis also uncovered novel SNP-gene pairs. For example, we found that the SNP, rs1332018, was associated to GSTM3. GSTM3 expression has been linked to Hepatitis B which African Americans suffer from disproportionately. Our methylation-adjusted method adds new understanding to the genetic basis of complex diseases that disproportionally affect African Americans.