Search

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 twenty-four 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.

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.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.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.

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.

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.

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.

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.

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.

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.

Background Thrombelastography (TEG), an ex-vivo clotting assay can identify subjects at high risk of subsequent coronary thrombotic events. Synthesis of clotting factors is subject to post-translational regulation, which is modulated at least in part by miRNA. Hypothesis We hypothesized that miRNA sequencing may identify specific miRNA linked with measures of hypercoaguability by TEG. Methods Kaolin activated thrombelastography was performed in platelet poor citrate plasma from 61 subjects referred for cardiac catheterization. Time to clot formation (R), clot stabilization time (K), and maximal fibrin clot strength (MA) was measured. Next generation miRNA sequencing was done from RNA isolated from whole blood samples, which includes miRNA derived from leukocytes and platelets. Prediction of miRNA gene targets was performed with targetscan. Results Sequencing resulted in quantification of 371 distinct miRNA from whole blood samples. We found 13 miRNA correlating with alteration in TEG-R, 33 miRNA correlating with TEG-K, and 21 miRNA correlating with TEG-MA. Coagulation factors or genes associated with coagulation were found to be among predicted targets in 49 out of these 67 miRNA. Most common predicted targets included factors II, V, VII, X, XIII, fibrinogen, plasminogen-activator inhibitor, and tissue factor. Factor XIIIA1 was highly conserved gene target by miR-96-5p (one of only 3miRNA predicted for this gene). MiR-96-5p correlated with clot stabilization time (ρ=-0.26, p=0.047) which has been shown to be dependent on FXIIIa activity. MiR-22-3p was significantly correlated with TEG-K (ρ=0.28, p=0.034) and was only miRNA with highly conserved target site for prothrombin (Factor II). Conclusions In patients with cardiovascular disease miRNA sequencing combined with coagulation phenotype in silico analysis may identify novel links to coagulation that are associated with increased thrombotic risk. Acknowledgement/Funding Charles Fisch Cardiovascular Research Award endowed by S. Knoebel

A large number of single nucleotide variants (SNVs) in the human genome are known to be responsible for inherited disease. An even larger number of SNVs, particularly those located in introns, have yet to be investigated for their pathogenic potential. Using known pathogenic and neutral intronic SNVs (iSNVs), we developed the regSNPs-intron algorithm based on a random forest classifier that integrates RNA splicing, protein structure and evolutionary conservation features. regSNPs-intron showed high accuracy in computing disease-causing probabilities of iSNVs. Using a high-throughput functional reporter assay called ASSET-seq (ASsay for Splicing using ExonTrap and sequencing), we validated regSNPs-intron predictions by measuring the impact of iSNVs on splicing outcome. Together, regSNPs-intron and ASSET-seq enable effective prioritization of iSNVs for disease pathogenesis. regSNPs-intron is available at https://regsnps-intron.ccbb.iupui.edu.

Transcriptome studies can identify genes whose expression differs between alcoholics and controls. To test which variants associated with alcohol use disorder (AUDs) may cause expression differences, we integrated deep RNA-seq and genome-wide association studies (GWAS) data from four postmortem brain regions of 30 AUDs subjects and 30 controls (social/non-drinkers) and analyzed allele-specific expression (ASE). We identified 90 genes with differential ASE in subjects with AUDs compared to controls. Of these, 61 genes contained 437 single nucleotide polymorphisms (SNPs) in the 3’ untranslated regions (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 showed affected RNA levels in a consistent manner in two neuroblastoma cell lines, SH-SY5Y and SK-N-BE(2). Many of these are in binding sites of miRNAs and RNA binding proteins, indicating that these SNPs are likely causal variants of AUD-associated differential ASE.

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 33% decreased CYP2B6 activity among normal metabolizers (AG/AG vs. CA/CA [p

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.

145 Background: Abiraterone is a CYP17 inhibitor approved for treatment of men with mCRPC. The NCI 9012 trial evaluated abiraterone alone with or without the PARP inhibitor veliparib in mCRPC patients. We hypothesized that germline genetic variation in the androgen axis and other metabolic enzymes would predict response to veliparib + abiraterone vs. abiraterone alone. Methods: A randomized trial cohort of (148) men with advanced mCRPC treated with abiraterone with or without veliparib was genotyped for 120 DNA polymorphisms in genes involved in androgen metabolism using Lifetech Open array chips. Blood for pharmacogenomic SNP analysis were collected at pre-treatment from each subject into 10-mL EDTA tube. Polymorphisms were tested using Cox models without treatment for prognostic testing and with treatment arm for predictive testing. Results: Genotyping was completed in 143 of 148 men; all were treated with abiraterone; 72 without veliparib (Median PFS: 10.3m) and 71 with veliparib (Median PFS: 11.3m). Polymorphisms in separate genes (SLCO2B1, KIF3C CYP19A, ESR1) were significantly (P ≤ .025) associated with progression-free survival (PFS) during abiraterone (q-value < 0.69). Polymorphisms in (CYP11A1, HSD17B4, ABHD13;LIG4, CYP19A1, HSD17B4, TRMT11) were predictive for PFS in patients treated with combination of abiraterone/veliparib compared to abiraterone alone (p-value < 0.025; q-value < 0.28). Conclusions: This analysis examines the influence of inherited variations on the efficacy of abiraterone, establishing the importance of pharmacogenomics on individual’s response to this therapy. Genotyping patients at these loci could be predictive of improved PFS to valiparib in combination with abiraterone. Further analysis of the association of more than one polymorphisms compared to zero or one with PFS associated with improved TTP demonstrated a better response to therapy than individuals carrying zero or one is ongoing. Clinical trial information: NCT01576172.

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.

We evaluated in this prospective study the effectiveness of continuous thoracic epidural anesthesia (TEA) and postoperative analgesia with ropivacaine and compared it with general anesthesia (GA) and opioids for pain relief, side effects, postanesthesia recovery, and hospital discharge after modified radical mastectomy. Sixty ASA physical status II and III patients undergoing mastectomy were randomly assigned to two study groups of 30 patients each. In the TEA group, an epidural catheter was inserted at T6-7, and 5--10 mL of 0.2% ropivacaine was injected to maintain anesthesia and to continuously administer adequate analgesia for 48 h. GA was induced with IV 1--2 mg of midazolam or 50--100 microg/mL of fentanyl followed by 50--150 mg of propofol and was maintained with sevoflurane and 50% N(2)O in oxygen. The Aldrete score system was used to evaluate postanesthesia recovery, a verbal rating scale was used for assessment of pain intensity, and a postanesthesia discharge scoring system was used for discharge home. The demographic data and side effects (except for nausea and vomiting) (GA 43%, TEA 10%, P = 0.0074) and discharge home were similar in both groups. However, the number of patients ready for discharge from the recovery room during the first postanesthesia hour (Aldrete score of 10) was significantly larger after TEA (80%) than after GA (33%) (P = 0.0006). GA patients experienced significantly more (P < 0.001) substantial pain than TEA patients on Day 0 (70%), Day 1 (53%), and Day 2 (27%) after the surgery. Patient satisfaction was greater with TEA (70%) than with GA (30%) (P < 0.001). We conclude that TEA with ropivacaine provides better postoperative pain relief and less nausea and vomiting, facilitates postanesthesia recovery, and gives greater patient satisfaction than GA.

IMPLICATIONS During a right lobectomy operation, a patient with carcinoma of the lung developed postoperative headache caused by a leak of cerebrospinal fluid from an area of dura injured during the procedure. Conservative treatment was unsuccessful. Injection of 10 mL of the patient's own blood into the epidural space relieved the headache.