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

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.2 nM) while also exhibiting good binding activity to both ER-α (EC50 = 72.1 nM) and ER-β (EC50 = 70.8 nM). In addition, a new synthesis was devised for the preparation of norendoxifen and its analogues through a bis-Suzuki coupling strategy.

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

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

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

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

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

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 (C(max) ) 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% (R(2) ) 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

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

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