Allele-specific expression and high-throughput reporter assay reveal functional genetic variants associated with alcohol use disorders.

Genome-wide association studies (GWAS) of complex traits, such as alcohol use disorders (AUD), usually identify variants in non-coding regions and cannot by themselves distinguish whether the associated variants are functional or in linkage disequilibrium with the functional variants. Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with AUD may cause expression differences, we integrated data from deep RNA-seq and GWAS of four postmortem brain regions from 30 subjects with AUD and 30 controls to analyze allele-specific expression (ASE). We identified 88 genes with differential ASE in subjects with AUD compared to controls. Next, to test one potential mechanism contributing to the differential ASE, we analyzed single nucleotide polymorphisms (SNPs) in the 3' untranslated regions (3'UTR) of these genes. Of the 88 genes with differential ASE, 61 genes contained 437 SNPs in the 3'UTR with at least one heterozygote among the subjects studied. Using a modified PASSPORT-seq (parallel assessment of polymorphisms in miRNA target-sites by sequencing) assay, we identified 25 SNPs that affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these SNPs are in binding sites of miRNAs and RNA-binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE. In sum, we demonstrate that a combination of computational and experimental approaches provides a powerful strategy to uncover functionally relevant variants associated with the risk for AUD.

A new Suzuki synthesis of triphenylethylenes that inhibit aromatase and bind to estrogen receptors α and ß.

The design and synthesis of dual aromatase inhibitors/selective estrogen receptor modulators (AI/SERMs) is an attractive strategy for the discovery of new breast cancer therapeutic agents. Previous efforts led to the preparation of norendoxifen (4) derivatives with dual aromatase inhibitory activity and estrogen receptor binding activity. In the present study, some of the structural features of the potent AI letrozole were incorporated into the lead compound (norendoxifen) to afford a series of new dual AI/SERM agents based on a symmetrical diphenylmethylene substructure that eliminates the problem of E,Z isomerization encountered with norendoxifen-based AI/SERMs. Compound 12d had good aromatase inhibitory activity (IC50=62.2nM) while also exhibiting good binding activity to both ER-α (EC50=72.1nM) and ER-ß (EC50=70.8nM). In addition, a new synthesis was devised for the preparation of norendoxifen and its analogues through a bis-Suzuki coupling strategy.

Evidence toward a dual phosphatase mechanism that restricts Aurora A (Thr-295) phosphorylation during the early embryonic cell cycle.

The mitotic kinase Aurora A (AurA) is regulated by a complex network of factors that includes co-activator binding, autophosphorylation, and dephosphorylation. Dephosphorylation of AurA by PP2A (human, Ser-51; Xenopus, Ser-53) destabilizes the protein, whereas mitotic dephosphorylation of its T-loop (human, Thr-288; Xenopus, Thr-295) by PP6 represses AurA activity. However, AurA(Thr-295) phosphorylation is restricted throughout the early embryonic cell cycle, not just during M-phase, and how Thr-295 is kept dephosphorylated during interphase and whether or not this mechanism impacts the cell cycle oscillator were unknown. Titration of okadaic acid (OA) or fostriecin into Xenopus early embryonic extract revealed that phosphatase activity other than PP1 continuously suppresses AurA(Thr-295) phosphorylation during the early embryonic cell cycle. Unexpectedly, we observed that inhibiting a phosphatase activity highly sensitive to OA caused an abnormal increase in AurA(Thr-295) phosphorylation late during interphase that corresponded with delayed cyclin-dependent kinase 1 (CDK1) activation. AurA(Thr-295) phosphorylation indeed influenced this timing, because AurA isoforms retaining an intact Thr-295 residue further delayed M-phase entry. Using mathematical modeling, we determined that one phosphatase would be insufficient to restrict AurA phosphorylation and regulate CDK1 activation, whereas a dual phosphatase topology best recapitulated our experimental observations. We propose that two phosphatases target Thr-295 of AurA to prevent premature AurA activation during interphase and that phosphorylated AurA(Thr-295) acts as a competitor substrate with a CDK1-activating phosphatase in late interphase. These results suggest a novel relationship between AurA and protein phosphatases during progression throughout the early embryonic cell cycle and shed new light on potential defects caused by AurA overexpression.

Quantification of spatial pharmacogene expression heterogeneity in breast tumors.

BACKGROUND: Chemotherapeutic drug concentrations vary across different regions of tumors and this is thought to be involved in development of chemotherapy resistance. Insufficient drug delivery to some regions of the tumor may be due to spatial differences in expression of genes involved in the disposition, transport, and detoxification of drugs (pharmacogenes). Therefore, in this study, we analyzed the spatial expression of 286 pharmacogenes in six breast cancer tissues using the recently developed Visium spatial transcriptomics platform to (1) determine if these pharmacogenes are expressed heterogeneously across tumor tissue and (2) to determine which pharmacogenes have the most spatial expression heterogeneity. METHODS AND RESULTS: The spatial transcriptomics technology sequences the transcriptome of 55 um diameter barcoded sections (spots) across a tissue sample. We analyzed spatial gene expression profiles of four biobank-sourced breast tumor samples in addition to two breast tumor sample datasets from 10× Genomics. We define heterogeneity as the interquartile range of read counts. Collectively, we identified 8887 spots in tumor regions, 3814 in stroma, 44 in lymphocytes, and 116 in normal regions based on pathologist annotation of the tissues. We showed statistically significant differences in expression of pharmacogenes in tumor regions compared to surrounding non-tumor regions. We also observed that the most heterogeneously expressed genes within tumor regions were involved in reactive oxygen species (ROS) handling and detoxification mechanisms. GPX4, GSTP1, MGST3, SOD1, CYP4Z1, CYB5R3, GSTK1, and NAT1 showed the most heterogeneous expression within tumor regions. CONCLUSIONS: The heterogeneous expression of these pharmacogenes may have important implications for cancer therapy due to their ability to impact drug distribution and efficacy throughout the tumor. Our results suggest that chemoresistance caused by expression of GPX4, GSTP1, MGST3, and SOD1 may be intrinsic, not acquired, since the heterogeneity is not specific to chemotherapy-treated samples or cell type. Additionally, we identified candidate chemoresistance pharmacogenes that can be further tested through focused follow-up studies.

Clinically important alterations in pharmacogene expression in histologically severe nonalcoholic fatty liver disease.

Polypharmacy is common in patients with nonalcoholic fatty liver disease (NAFLD) and previous reports suggest that NAFLD is associated with altered drug disposition. This study aims to determine if patients with NAFLD are at risk for altered drug response by characterizing changes in hepatic mRNA expression of genes mediating drug disposition (pharmacogenes) across the histological NAFLD severity spectrum. We utilize RNA-seq for 93 liver biopsies with histologically staged NAFLD Activity Score (NAS), fibrosis stage, and steatohepatitis (NASH). We identify 37 significant pharmacogene-NAFLD severity associations including CYP2C19 downregulation. We chose to validate CYP2C19 due to its actionability in drug prescribing. Meta-analysis of 16 independent studies demonstrate that CYP2C19 is significantly downregulated to 46% in NASH, to 58% in high NAS, and to 43% in severe fibrosis. Our data demonstrate the downregulation of CYP2C19 in NAFLD which supports developing personalized medicine approaches for drugs sensitive to metabolism by the CYP2C19 enzyme.

MicroRNA sequencing in patients with coronary artery disease - considerations for use as biomarker for thrombotic risk.

MicroRNAs (miRNAs) are small RNAs integral in the regulation of gene expression. Analysis of circulating miRNA levels may identify patients with coronary artery disease (CAD) at risk for recurrent myocardial infarction (MI) after percutaneous coronary interventions (PCIs). Subjects with CAD were selected from the GENCATH cardiac catheterization biobank. Subjects with recurrent MI after PCI were compared with those without recurrent MI during follow-up in the initial (n = 48) and replication cohort (n = 67). Next generation MiRNA sequencing was performed on plasma samples and whole blood samples fixed with PAXGENE tubes upon collection. Overall, 164 miRNAs derived from whole blood were differentially expressed in the replication cohort between subjects with and without recurrent MI events (p < 0.05), with 69 remaining significant after false-discovery rate (FDR) correction. None of the miRNAs in plasma was significantly different by FDR among subjects with and without MI. Overall, correlation between direction of effects between plasma and whole blood assays was variable, and only two miRNAs were concordant and significant in both. Associations of miRNA with vascular disease, MI, and thrombosis were further explored. MiRNA profiling has potential as the future biomarker for disease prognosis and treatment response marker in secondary treatment of patients with CAD after PCI. Whole blood may be the preferred sample source as compared to plasma.

Mapping the miRNA-mRNA Interactome in Human Hepatocytes and Identification of Functional mirSNPs in Pharmacogenes.

MiRNAs regulate the expression of hepatic genes involved in pharmacokinetics and pharmacodynamics. Genetic variants affecting miRNA binding (mirSNPs) have been associated with altered drug response, but previously used methods to identify miRNA binding sites and functional mirSNPs in pharmacogenes are indirect and limited by low throughput. We utilized the high-throughput chimeric-eCLIP assay to directly map thousands of miRNA-mRNA interactions and define the miRNA binding sites in primary hepatocytes. We then used the high-throughput PASSPORT-seq assay to functionally test 262 potential mirSNPs with coordinates overlapping the identified miRNA binding sites. Using chimeric-eCLIP, we identified a network of 448 miRNAs that collectively target 11,263 unique genes in primary hepatocytes pooled from 100 donors. Our data provide an extensive map of miRNA binding of each gene, including pharmacogenes, expressed in primary hepatocytes. For example, we identified the hsa-mir-27b-DPYD interaction at a previously validated binding site. A second example is our identification of 19 unique miRNAs that bind to CYP2B6 across 20 putative binding sites on the transcript. Using PASSPORT-seq, we then identified 24 mirSNPs that functionally impacted reporter mRNA levels. To our knowledge, this is the most comprehensive identification of miRNA binding sites in pharmacogenes. Combining chimeric-eCLIP with PASSPORT-seq successfully identified functional mirSNPs in pharmacogenes that may affect transcript levels through altered miRNA binding. These results provide additional insights into potential mechanisms contributing to interindividual variability in drug response.

Variability of Dosing and Number of Medications Needed to Achieve Adequate Sedation in Mechanically Ventilated Pediatric Intensive Care Patients.

Children admitted to the pediatric intensive care unit (PICU) often require multiple medications to achieve comfort and sedation. Although starting doses are available, these medications are typically titrated to the desired effect. Both oversedation and undersedation are associated with adverse events. The aim of this retrospective study was to evaluate cumulative medication burden necessary to achieve comfort in patients in the PICU and determine relevant predictors of medication needs. In order to account for all of the sedative medications, z-scores were used to assess the population average dose of each medication and compare each patient day to this population average. Sedation regimens for 130 patients in the PICU were evaluated. Mean overall infusion rates of fentanyl, morphine, and hydromorphone were 1.67 ± 0.81 µg/kg/hour, 0.12 ± 0.08 mg/kg/hour, and 17.84 ± 13.4 µg/kg/hour, respectively. The mean infusion rate of dexmedetomidine was 0.59 ± 0.28 µg/kg/hour, and midazolam was 0.14 ± 0.1 mg/kg/hour. Summation z-sores were used to rank the amount of sedation medication needed to achieve comfort for each individual patient for his/her PICU stay in relation to the entire sample. Patient age, weight, and length of mechanical ventilation were all significant predictors of sedation requirement. This study will provide data necessary to develop a model of cumulative medication burden needed to achieve appropriate sedation in this population. This descriptive model in appropriately ranking patients based on sedative needs is the first step in exploring potential genetic factors that may provide an insight into homing in on the appropriate sedation regimen.

Circulating miRNAs as Biomarkers for CYP2B6 Enzyme Activity.

The CYP2B6 gene is highly polymorphic and its activity shows wide interindividual variability. However, substantial variability in CYP2B6 activity remains unexplained by the known CYP2B6 genetic variations. Circulating, cell-free micro RNAs (miRNAs) may serve as biomarkers of hepatic enzyme activity. CYP2B6 activity in 72 healthy volunteers was determined using the disposition of efavirenz as a probe drug. Circulating miRNA expression was quantified from baseline plasma samples. A linear model consisting of the effects of miRNA expression, genotype-determined metabolizer status, and demographic information was developed to predict CYP2B6 activity. Expression of 2,510 miRNAs were quantified out of which 7 miRNAs, together with the CYP2B6-genotypic metabolizer status and demographics, was shown to be predictive markers for CYP2B6 activity. The reproducibility of the model was evaluated by cross-validation. The average Pearson's correlation (R) between the predicted and observed maximum plasma concentration (Cmax ) ratios of efavirenz and its metabolite-8-OH efavirenz using the linear model with all features (7 miRNA + metabolizer status + age + sex + race) was 0.6702. Similar results were also observed using area under the curve (AUC) ratios (Pearson correlation's R = 0.6035). Thus, at least 36% (R2 ) of the variability of in vivo CYP2B6 activity was explained using this model. This is a significant improvement over the models using only the genotype-based metabolizer status or the demographic information, which explained only 6% or less of the variability of in vivo CYP2B6 activity. Our results, therefore, demonstrate that circulating plasma miRNAs can be valuable biomarkers for in vivo CYP2B6 activity.

RegSNPs-intron: a computational framework for predicting pathogenic impact of intronic single nucleotide variants.

Single nucleotide variants (SNVs) in intronic regions have yet to be systematically investigated for their disease-causing potential. Using known pathogenic and neutral intronic SNVs (iSNVs) as training data, we develop the RegSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure, and evolutionary conservation features. RegSNPs-intron showed excellent performance in evaluating the pathogenic impacts of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we evaluate the impact of RegSNPs-intron predictions on splicing outcome. Together, RegSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis.

PASSPORT-seq: A Novel High-Throughput Bioassay to Functionally Test Polymorphisms in Micro-RNA Target Sites.

Next-generation sequencing (NGS) studies have identified large numbers of genetic variants that are predicted to alter miRNA-mRNA interactions. We developed a novel high-throughput bioassay, PASSPORT-seq, that can functionally test in parallel 100s of these variants in miRNA binding sites (mirSNPs). The results are highly reproducible across both technical and biological replicates. The utility of the bioassay was demonstrated by testing 100 mirSNPs in HEK293, HepG2, and HeLa cells. The results of several of the variants were validated in all three cell lines using traditional individual luciferase assays. Fifty-five mirSNPs were functional in at least one of three cell lines (FDR ≤ 0.05); 11, 36, and 27 of them were functional in HEK293, HepG2, and HeLa cells, respectively. Only four of the variants were functional in all three cell lines, which demonstrates the cell-type specific effects of mirSNPs and the importance of testing the mirSNPs in multiple cell lines. Using PASSPORT-seq, we functionally tested 111 variants in the 3' UTR of 17 pharmacogenes that are predicted to alter miRNA regulation. Thirty-three of the variants tested were functional in at least one cell line.

Next generation MicroRNA sequencing to identify coronary artery disease patients at risk of recurrent myocardial infarction.

BACKGROUND AND AIMS: Variation in micro-RNA (miRNA) levels in blood has been associated with alterations of physiological functions of the cardiovascular system. Circulating miRNA have the potential to become reliable biomarkers for risk stratification and early detection of cardiovascular events. Recurrent thrombotic events in patients with established coronary artery disease (CAD) demonstrate the need for personalized approaches to secondary prevention, especially in light of recent novel treatment approaches. METHODS: In a single center cohort study, whole blood samples were collected from 437 subjects undergoing cardiac catheterization, who were followed for recurrent cardiovascular events during a mean follow up of 1.5 years. We selected a case cohort (n = 22) with recurrent thrombotic events on standard medical therapy (stent thrombosis (n = 6) or spontaneous myocardial infarction (MI) (n = 16)) and a matched cohort with CAD, but uneventful clinical follow up (n = 26), as well as a control group with cardiovascular risk factors, but without angiographic CAD (n = 24). We performed complete miRNA next generation sequencing of RNA extracted from whole blood samples (including leukocytes and platelets). RESULTS: A differential pattern of miRNA expression was found among controls, CAD patients with no events, and CAD patients with recurrent events. MiRNA previously associated with MI, CAD, endothelial function, vascular smooth muscle cells, platelets, angiogenesis, heart failure, cardiac hypertrophy, arrhythmia, and stroke were found variably expressed in our case-control cohorts. Seventy miRNA (FDR <0.05) were linked to the risk of recurrent myocardial infarction and future stent thrombosis, as compared to CAD patients with subsequently uneventful follow up. CONCLUSIONS: MiRNA next generation sequencing demonstrates altered fingerprint profile of whole blood miRNA expression among subjects with subsequent recurrent thrombotic events on standard medical therapy ('non-responders'), as compared to subjects with no recurrent cardiovascular events. MiRNA profiling may be useful to identify high risk subjects and provide additional insights into disease mechanisms not currently attenuated with standard medical therapy used in CAD treatment.

Variants in the CYP2B6 3'UTR Alter In Vitro and In Vivo CYP2B6 Activity: Potential Role of MicroRNAs.

CYP2B6*6 and CYP2B6*18 are the most clinically important variants causing reduced CYP2B6 protein expression and activity. However, these variants do not account for all variability in CYP2B6 activity. Emerging evidence has shown that genetic variants in the 3'UTR may explain variable drug response by altering microRNA regulation. Five 3'UTR variants were associated with significantly altered efavirenz AUC0-48 (8-OH-EFV/EFV) ratios in healthy human volunteers. The rs70950385 (AG>CA) variant, predicted to create a microRNA binding site for miR-1275, was associated with a 33% decreased CYP2B6 activity among normal metabolizers (AG/AG vs. CA/CA (P < 0.05)). In vitro luciferase assays were used to confirm that the CA on the variant allele created a microRNA binding site causing an 11.3% decrease in activity compared to the AG allele when treated with miR-1275 (P = 0.0035). Our results show that a 3'UTR variant contributes to variability in CYP2B6 activity.

In Vivo siRNA Delivery and Rebound of Renal LRP2 in Mice.

siRNA stabilized for in vivo applications is filtered and reabsorbed in the renal proximal tubule (PT), reducing mRNA expression transiently. Prior siRNA efforts have successfully prevented upregulation of mRNA in response to injury. We proposed reducing constitutive gene and protein expression of LRP2 (megalin) in order to understand its molecular regulation in mice. Using siRNA targeting mouse LRP2 (siLRP2), reduction of LRP2 mRNA expression was compared to scrambled siRNA (siSCR) in mouse PT cells. Mice received siLRP2 administration optimized for dose, administration site, carrier solution, administration frequency, and administration duration. Kidney cortex was collected upon sacrifice. Renal gene and protein expression were compared by qRT-PCR, immunoblot, and immunohistochemistry (IHC). Compared to siSCR, siLRP2 reduced mRNA expression in PT cells to 16.6% ± 0.6%. In mouse kidney cortex, siLRP2 reduced mRNA expression to 74.8 ± 6.3% 3 h and 70.1 ± 6.3% 6 h after administration. mRNA expression rebounded at 12 h (160.6 ± 11.2%). No megalin renal protein expression reduction was observed by immunoblot or IHC, even after serial twice daily dosing for 3.5 days. Megalin is a constitutively expressed protein. Although LRP2 renal mRNA expression reduction was achieved, siRNA remains a costly and inefficient intervention to reduce in vivo megalin protein expression.