Identification of Translational microRNA Biomarker Candidates for Ketoconazole-Induced Liver Injury Using Next-Generation Sequencing.

Drug-induced liver injury (DILI) is a leading cause of acute liver failure. Reliable and translational biomarkers are needed for early detection of DILI. microRNAs (miRNAs) have received wide attention as a novel class of potential DILI biomarkers. However, it is unclear how DILI drugs other than acetaminophen may influence miRNA expression or which miRNAs could serve as useful biomarkers in humans. We selected ketoconazole (KCZ), a classic hepatotoxin, to study miRNA biomarkers for DILI as a proof of concept for a workflow that integrated in vivo, in vitro, and bioinformatics analyses. We examined hepatic miRNA expression in KCZ-treated rats at multiple doses and durations using miRNA-sequencing and correlated our results with conventional DILI biomarkers such as liver histology. Significant dysregulation of rno-miR-34a-5p, rno-miR-331-3p, rno-miR-15b-3p, and rno-miR-676 was associated with cytoplasmic vacuolization, a phenotype in rat livers with KCZ-induced injury, which preceded the elevation of serum liver transaminases (ALT and AST). Between rats and humans, miR-34a-5p, miR-331-3p, and miR-15b-3p were evolutionarily conserved with identical sequences, whereas miR-676 showed 73% sequence similarity. Using quantitative PCR, we found that the levels of hsa-miR-34a-5p, hsa-miR-331-3p, and hsa-miR-15b-3p were significantly elevated in the culture media of HepaRG cells treated with 100 µM KCZ (a concentration that induced cytotoxicity). Additionally, we computationally characterized the miRNA candidates for their gene targeting, target functions, and miRNA/target evolutionary conservation. In conclusion, we identified miR-34a-5p, miR-331-3p, and miR-15b-3p as translational biomarker candidates for early detection of KCZ-induced liver injury with a workflow applicable to computational toxicology studies.

Long noncoding RNA LINC00844-mediated molecular network regulates expression of drug metabolizing enzymes and nuclear receptors in human liver cells.

Noncoding RNAs, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), regulate gene expression in many physiological and pathological processes, including drug metabolism. Drug metabolizing enzymes (DMEs) are critical components in drug-induced liver toxicity. In this study, we used human hepatic HepaRG cells treated with 5 or 10 mM acetaminophen (APAP) as a model system and identified LINC00844 as a toxicity-responsive lncRNA. We analyzed the expression profiles of LINC00844 in different human tissues. In addition, we examined the correlations between the levels of LINC00844 and those of key DMEs and nuclear receptors (NRs) for APAP metabolism in humans. Our results showed that lncRNA LINC00844 is enriched in the liver and its expression correlates positively with mRNA levels of CYP3A4, CYP2E1, SULT2A1, pregnane X receptor (PXR), and hepatocyte nuclear factor (HNF) 4α. We demonstrated that LINC00844 regulates the expression of these five genes in HepaRG cells using gain- and loss-of-function assays. Further, we discovered that LINC00844 is localized predominantly in the cytoplasm and acts as an hsa-miR-486-5p sponge, via direct binding, to protect SULT2A1 from miRNA-mediated gene silencing. Our data also demonstrated a functional interaction between LINC00844 and hsa-miR-486-5p in regulating DME and NR expression in HepaRG cells and primary human hepatocytes. We depicted a LINC00844-mediated regulatory network that involves miRNA and NRs and influences DME expression in response to APAP toxicity.

Coordinated Regulation of UGT2B15 Expression by Long Noncoding RNA LINC00574 and hsa-miR-129-5p in HepaRG Cells.

Recent studies have shown that microRNAs and long noncoding RNAs (lncRNAs) regulate the expression of drug metabolizing enzymes (DMEs) in human hepatic cells and that a set of DMEs, including UDP glucuronosyltransferase (UGT) 2B15, is down-regulated dramatically in liver cells by toxic acetaminophen (APAP) concentrations. In this study we analyzed mRNA, microRNA, and lncRNA expression profiles in APAP-treated HepaRG cells to explore noncoding RNA-dependent regulation of DME expression. The expression of UGT2B15 and lncRNA LINC00574 was decreased in APAP-treated HepaRG cells. UGT2B15 levels were diminished by LINC00574 suppression using antisense oligonucleotides or small interfering RNA. Furthermore, we found that hsa-miR-129-5p suppressed LINC00574 and decreased UGT2B15 expression via LINC00574 in HepaRG cells. In conclusion, our results indicate that LINC00574 acts as an important regulator of UGT2B15 expression in human hepatic cells, providing experimental evidence and new clues to understand the role of cross-talk between noncoding RNAs. SIGNIFICANCE STATEMENT: We showed a molecular network that displays the cross-talk and consequences among mRNA, micro RNA, long noncoding RNA, and proteins in acetaminophen (APAP)-treated HepaRG cells. APAP treatment increased the level of hsa-miR-129-5p and decreased that of LINC00574, ultimately decreasing the production of UDP glucuronosyltransferase (UGT) 2B15. The proposed regulatory network suppresses UGT2B15 expression through interaction of hsa-miR-129-5p and LINC00574, which may be achieved potentially by recruiting RNA binding proteins.

FREMSA: A Method That Provides Direct Evidence of the Interaction between microRNA and mRNA.

microRNAs (miRNAs) modulate the expression of enzymes responsible for activation or detoxification of xenobiotics and toxicants. miRNAs are dysregulated in response to environmental exposure and have been implicated in toxicological events. Many in vivo and in vitro experimental approaches have been employed to delineate the mechanisms by which miRNAs regulate target genes; however, all these methods provide only indirect evidence for the interaction between miRNAs and their counterpart mRNA molecules. In this chapter, we describe a novel approach-a fluorescent-based RNA electrophoretic mobility shift assay (FREMSA) that is a sensitive and time-saving method, with a high specificity, to visualize the interactions among miRNAs, mRNAs, and proteins, as direct evidence of mRNA/miRNA complex formation.

MicroRNAs hsa-miR-495-3p and hsa-miR-486-5p suppress basal and rifampicin-induced expression of human sulfotransferase 2A1 (SULT2A1) by facilitating mRNA degradation.

Drug metabolizing enzymes mediate biotransformation of drugs and play an essential role in drug efficacy and toxicity. Human sulfotransferases are a superfamily of Phase II detoxification enzymes that metabolize a wide spectrum of endogenous compounds and xenobiotics. SULT2A1 is one of the most abundant hepatic sulfotransferases and it catalyzes the sulfate conjugation of many endogenous substrates, such as bile acids and steroids. In the current study, we utilized a systematic approach by combining a series of computational analyses and in vitro methods to identify miRNAs that repress SULT2A1 expression post-transcriptionally. Our in silico analyses predicted miRNA response elements for hsa-miR-495-3p and hsa-miR-486-5p within the 3'-UTR of SULT2A1 mRNA and the levels of these miRNAs were inversely correlated with that of SULT2A1 mRNA in human liver. Using fluorescence-based RNA electrophoretic mobility shift assays, we found that hsa-miR-495-3p and hsa-miR-486-5p interacted directly with the SULT2A1 3'-UTR. The activity of a luciferase reporter gene construct containing sequences from the SULT2A1 3-UTR was suppressed by hsa-miR-486-5p and hsa-miR-495-3p. Furthermore, gain- and loss-of-function assays demonstrated that hsa-miR-486-5p and hsa-miR-495-3p negatively modulate basal and rifampicin-induced expression of SULT2A1 in HepG2 cells by decreasing mRNA stability.

Epigenome-wide association study of peripheral blood mononuclear cells in systemic lupus erythematosus: Identifying DNA methylation signatures associated with interferon-related genes based on ethnicity and SLEDAI.

Systemic lupus erythematosus (SLE or lupus) is a heterogeneous autoimmune disease characterized by the involvement of multiple organs and the production of antinuclear antibodies. DNA methylation plays an important role in the pathogenesis of lupus. We have performed an epigenome-wide DNA methylation study in lupus and healthy control (non-lupus) subjects to identify epigenetic patterns in lupus characterized ethnicity and SLE disease activity index (SLEDAI). A total of fifty-seven lupus patients (39 African American (AA) and 18 European American (EA)) and 33 healthy controls (17 AA and 16 EA) were studied. Differential DNA methylation between lupus patients and controls was assessed for approximately 485,000 CpG sites across the genome. We identified 41 differentially methylated sites (associated with 30 genes) between lupus and control s subjects, 85% of which were hypomethylated. Significant hypomethylation of differentially methylated sites was associated with several interferon-related genes, including MX1, IFI44L, PARP9, DT3XL, IFIT1, IFI44, RSAD2, PLSCR1, and IRF7. Several of these associated genes were also hypomethylated in comparisons between AA lupus and AA non-lupus subjects and between lupus patients with SLEDAI>6 and non-lupus subjects. Our analysis of gene expression data through RT-PCR confirmed these findings. Thus, the results indicate epigenetics susceptibility in lupus, which may be associated with SLEDAI score and ethnicity. In addition, our findings support the importance of the Type 1 interferon pathway in lupus pathogenesis.

A functional SNP in the 3'-UTR of TAP2 gene interacts with microRNA hsa-miR-1270 to suppress the gene expression.

The transporter associated with antigen processing 2 (TAP2) is involved in the development of multidrug resistance and the etiology of immunological diseases. In this study, we investigated whether the expression of TAP2 can be perturbed by single nucleotide polymorphisms (SNPs) located in 3'-untranslated region (3'-UTR) of the gene via interactions with microRNAs. Using a series of in silico assays, we selected the candidate microRNAs (miRNAs) with the potential to interact with functional SNPs of TAP2. The SNP rs241456-located in the 3'-UTR of TAP2-resides in a potential binding site for hsa-miR-1270 and hsa-miR-620. HEK 293 cells, from a human kidney cell line, were used to characterize the extent of binding of miRNAs to each polymorphic allele of the SNP by a luciferase reporter gene assay. RNA electrophoretic mobility shift assays were used to evaluate the interaction between the miRNAs and each allele sequence of the SNP. We found that hsa-miR-1270 inhibited luciferase activity by binding to the T allele of the SNP in an allele-specific manner. A negative correlation was also found between the expression of hsa-miR-1270 and the T allele of the SNP in kidney tissues. Our findings support the hypothesis that hsa-miR-1270 suppresses the production of TAP2 by binding to this SNP in the 3'-UTR of this gene. Environ. Mol. Mutagen. 59:134-143, 2018. © 2017 Wiley Periodicals, Inc.

Multiple microRNAs function as self-protective modules in acetaminophen-induced hepatotoxicity in humans.

Acetaminophen (APAP) overdose is the leading cause of acute liver failure. Yet the mechanisms underlying adaptive tolerance toward APAP-induced liver injury are not fully understood. To better understand molecular mechanisms contributing to adaptive tolerance to APAP is an underpinning foundation for APAP-related precision medicine. In the current study, the mRNA and microRNA (miRNA) expression profiles derived from next generation sequencing data for APAP-treated (5 and 10 mM) HepaRG cells and controls were analyzed systematically. Putative miRNAs targeting key dysregulated genes involved in APAP hepatotoxicity were selected using in silico prediction algorithms, un-biased gene ontology, and network analyses. Luciferase reporter assays, RNA electrophoresis mobility shift assays, and miRNA pull-down assays were performed to investigate the role of miRNAs affecting the expression of dysregulated genes. Levels of selected miRNAs were measured in serum samples obtained from children with APAP overdose (58.6-559.4 mg/kg) and from healthy controls. As results, 2758 differentially expressed genes and 47 miRNAs were identified. Four of these miRNAs (hsa-miR-224-5p, hsa-miR-320a, hsa-miR-449a, and hsa-miR-877-5p) suppressed drug metabolizing enzyme (DME) levels involved in APAP-induced liver injury by downregulating HNF1A, HNF4A and NR1I2 expression. Exogenous transfection of these miRNAs into HepaRG cells effectively rescued them from APAP toxicity, as indicated by decreased alanine aminotransferase levels. Importantly, hsa-miR-320a and hsa-miR-877-5p levels were significantly elevated in serum samples obtained from children with APAP overdose compared to health controls. Collectively, these data indicate that hsa-miR-224-5p, hsa-miR-320a, hsa-miR-449a, and hsa-miR-877-5p suppress DME expression involved in APAP-induced hepatotoxicity and they contribute to an adaptive response in hepatocytes.

The expression, induction and pharmacological activity of CYP1A2 are post-transcriptionally regulated by microRNA hsa-miR-132-5p.

Cytochrome P450 1A2 (CYP1A2) is one of the most abundant and important drug metabolizing enzymes in human liver. However, little is known about the post-transcriptional regulation of CYP1A2, especially the mechanisms involving microRNAs (miRNAs). This study applied a systematic approach to investigate the post-transcriptional regulation of CYP1A2 by miRNAs. Candidate miRNAs targeting the 3'-untranslated region (3'-UTR) of CYP1A2 were screened in silico, resulting in the selection of sixty-two potential miRNAs for further analysis. The levels of two miRNAs, hsa-miR-132-5p and hsa-miR-221-5p, were inversely correlated with the expression of CYP1A2 mRNA transcripts in normal human liver tissue samples represented in The Cancer Genome Atlas (TCGA) dataset. The interactions between these miRNAs and cognate CYP1A2 mRNA sequences were evaluated using luciferase reporter gene studies and electrophoretic mobility shift assays, by which a direct interaction was confirmed involving hsa-miR-132-5p and a cognate binding site present in the CYP1A2 3'-UTR. Experiments by which hsa-miR-132-5p or random miRNA controls were introduced into HepG2, Huh-7 and HepaRG hepatic cell lines showed that only hsa-miR-132-5p suppressed the endogenous and lansoprazole-induced expression of CYP1A2, at biological activity, protein production, and mRNA transcript levels. Furthermore, 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), and lactate dehydrogenase (LDH) assays showed that hsa-miR-132-5p attenuates CYP1A2-mediated, lansoprazole-enhanced, flutamide-induced hepatic cell toxicity. Results from multilayer experiments demonstrate that hsa-miR-132-5p suppresses the expression of CYP1A2 and that this suppression is able to decrease the extent of an adverse drug-drug interaction involving lansoprazole and flutamide.

MicroRNA hsa-miR-370-3p suppresses the expression and induction of CYP2D6 by facilitating mRNA degradation.

Cytochrome P450 2D6 (CYP2D6) participates in the metabolism of approximately 20-25% of prescribed drugs. Genetic polymorphisms influence the expression and/or activity of CYP2D6, and inter-individual differences in drug activation and elimination caused by CYP2D6 genetic variants were reported. However, little is known about the potential modulation of CYP2D6 expression by microRNAs (miRNAs). In the current study, by using in silico prediction of the stabilities of miRNA/mRNA complexes, we screened 38 miRNA candidates that may interact with the transcript of CYP2D6. An inverse correlation between the expression of miRNA hsa-miR-370-3p and the expression of CYP2D6 was observed in human liver tissue samples. Electrophoretic mobility shift assays confirmed that hsa-miR-370-3p was able to directly bind to its cognate target within the coding region of the CYP2D6 transcript. The transfection of hsa-miR-370-3p mimics into the HepG2CYP2D6 cell line, a genetically modified cell line that overexpresses exogenous CYP2D6, was able to suppress the expression of CYP2D6 significantly at both mRNA and protein levels. The transfection of hsa-miR-370-3p mimics was also able to inhibit endogenous mRNA expression and/or protein production of CYP2D6 in HepaRG cells. Furthermore, in HepaRG, HepG2, and Huh7 cells, dexamethasone-induced expression of CYP2D6 was inhibited by hsa-miR-370-3p mimics. To investigate whether the miRNA mediated suppression is caused by inhibiting protein translation or promoting mRNA degradation, an actinomycin D assay was used to measure the stability of CYP2D6 transcripts. The results indicated that hsa-miR-370-3p mimics facilitated significantly the degradation of CYP2D6 mRNA. In addition, proteomics analyses of proteins isolated from the miRNA/mRNA/protein complex suggested that a group of multifunctional proteins facilitated the interaction between hsa-miR-370-3p and CYP2D6, thereby promoting mRNA degradation.

Systems view of adipogenesis via novel omics-driven and tissue-specific activity scoring of network functional modules.

The investigation of the complex processes involved in cellular differentiation must be based on unbiased, high throughput data processing methods to identify relevant biological pathways. A number of bioinformatics tools are available that can generate lists of pathways ranked by statistical significance (i.e. by p-value), while ideally it would be desirable to functionally score the pathways relative to each other or to other interacting parts of the system or process. We describe a new computational method (Network Activity Score Finder - NASFinder) to identify tissue-specific, omics-determined sub-networks and the connections with their upstream regulator receptors to obtain a systems view of the differentiation of human adipocytes. Adipogenesis of human SBGS pre-adipocyte cells in vitro was monitored with a transcriptomic data set comprising six time points (0, 6, 48, 96, 192, 384 hours). To elucidate the mechanisms of adipogenesis, NASFinder was used to perform time-point analysis by comparing each time point against the control (0 h) and time-lapse analysis by comparing each time point with the previous one. NASFinder identified the coordinated activity of seemingly unrelated processes between each comparison, providing the first systems view of adipogenesis in culture. NASFinder has been implemented into a web-based, freely available resource associated with novel, easy to read visualization of omics data sets and network modules.

MicroRNA hsa-miR-25-3p suppresses the expression and drug induction of CYP2B6 in human hepatocytes.

Cytochrome P450 2B6 (CYP2B6), mainly expressed in the liver and brain, is important for processing a number of widely used drugs. Variations in CYP2B6 expression are associated with decreased drug efficacy or adverse effects in some patients. Although CYP2B6 genetic variants are associated with its differential expression, epigenetic mechanisms affecting CYP2B6 gene regulation have not been established. Sequence analysis identified 29 domains in the CYP2B6 mRNA transcript that could be subject to regulation by microRNAs. Inverse correlations were found in human hepatocytes for the levels of the microRNAs hsa-miR-504-5p and hsa-miR-25-3p compared with CYP2B6 mRNA. Reporter gene assays showed that hsa-miR-25-3p suppresses CYP2B6 expression by targeting a specific sequence in the 3'-untranslated region of the mRNA transcript. Electrophoretic mobility shift assays confirmed that hsa-miR-25-3p forms stable complexes with its cognate mRNA sequence and that it recruits cellular factors, including Ago-4. Transfection of HepaRG cells with hsa-miR-25-3p mimics inhibited expression of the endogenous CYP2B6 gene and it also decreased rifampicin-dependent induction of CYP2B6 at the mRNA and protein levels. In summary, in silico and in vitro analyses show that hsa-miR-25-3p suppresses CYP2B6 expression in human liver cells via an epigenetic mechanism.

Whole genome and normalized mRNA sequencing reveal genetic status of TK6, WTK1, and NH32 human B-lymphoblastoid cell lines.

Closely related TK6, WTK1, and NH32 human B-lymphoblastoid cell lines differ in their p53 functional status. These lines are used frequently in genotoxicity studies and in studies aimed at understanding the role of p53 in DNA repair. Despite their routine use, little is known about the genetic status of these cells. To provide insight into their genetic composition, we sequenced and analyzed the entire genome of TK6 cells, as well as the normalized transcriptomes of TK6, WTK1, and NH32 cells. Whole genome sequencing (WGS) identified 21,561 genes and 5.17×10(6) small variants. Within the small variants, 50.54% were naturally occurring single nucleotide polymorphisms (SNPs) and 49.46% were mutations. The mutations were comprised of 92.97% single base-pair substitutions and 7.03% insertions or deletions (indels). The number of predicted genes, SNPs, and small mutations are similar to frequencies observed in the human population in general. Normalized mRNA-seq analysis identified the expression of transcripts bearing SNPs or mutations for TK6, WTK1, and NH32 as 2.88%, 2.04%, and 1.71%, respectively, and several of the variant transcripts identified appear to have important implications in genetic toxicology. These include a single base deletion mutation in the ferritin heavy chain gene (FTH1) resulting in a frame shift and protein truncation in TK6 that impairs iron metabolism. SNPs in the thiopurine S-methyltransferase (TPMT) gene (TPMT*3A SNP), and in the xenobiotic metabolizing enzyme, NADPH quinine oxidoreductase 1 (NQO1) gene (NQO1*2 SNP), are both associated with decreased enzyme activity. The clinically relevant TPMT*3A and NQO1*2 SNPs can make these cell lines useful in pharmacogenetic studies aimed at improving or tailoring drug treatment regimens that minimize toxicity and enhance efficacy.

Modulation of ALDH5A1 and SLC22A7 by microRNA hsa-miR-29a-3p in human liver cells.

Observed variations in drug responses among patients may result from differences in heritable genetic traits or from alterations in the epigenetic regulation of drug metabolizing enzymes and transporters (DMETs). MicroRNAs (miRNAs), a group of small non-coding RNAs, provide an epigenetic mechanism for fine-tuning the expression of targeted DMET genes by regulating the efficiency of protein translation and by decreasing mRNA stability via enhanced degradation. In the current study we systematically screened 374 important genes encoding DMETs for potential response elements to hsa-miR-29a-3p, a highly abundant miRNA in human liver. RNA electrophoresis mobility shift assays displayed direct interactions between hsa-miR-29a-3p and its cognate targets within the mRNA transcripts for the ABCC6, SLC22A7 and ALDH5A1 genes. The expression of luciferase reporter genes containing the 3'-UTRs of SLC22A7 or ALDH5A1 and the expression of endogenous SLC22A7 and ALDH5A1 were each suppressed by transfection with hsa-miR-29a-3p mimics. Importantly, chemically-induced up-regulation of hsa-miR-29a-3p correlated inversely with the expression of SLC22A7 and ALDH5A1. However, our studies failed to detect suppressive effects of hsa-miR-29a-3p on ABCC6 expression, which might be explained by the notion that the interaction of hsa-miR-29a-3p and ABCC6 mRNA was unable to recruit ribonucleoproteins to form a RNA-induced silencing complex.

Expression of drug transporters in human kidney: impact of sex, age, and ethnicity.

BACKGROUND: Differences in expression of drug transporters in human kidney contribute to changes in pharmacokinetics and toxicokinetics of a variety of drug compounds. The basal expression levels of genes involved in drug transport processes in the kidney introduces differences in bioavailability, distribution, and clearance of drugs, possibly influencing drug efficacy and adverse reactions. Sex differences in gene expression of transporters are a key cause of differences in sex-dependent pharmacokinetics, which may characterize many drugs and contribute to individual differences in drug efficacy and toxicity. Therefore, evaluating the expression of drug transporters in normal human kidneys is important to better understand differences in drug bioavailability, distribution, and clearance of drugs in humans. Other factors such as age and ethnicity may also contribute to individual differences in gene expression of drug transporters in the human kidney. METHODS: Quantitative real-time PCR (QRT-PCR) was performed to determine the gene expression of 30 drug transporters in 95 age-matched normal human kidney tissues. Multiple Student's t-tests (Sidak-Bonferroni correction) and two-way ANOVA (Bonferroni correction) analyses were used to determine statistically significant differences. RESULTS: In the 30 transporter genes examined, sex, ethnicity, and age differences in gene expression were exhibited in normal human kidney tissue. These changes in expression were not found to be differentially significant. However, sex-age and sex-ethnicity interactions were found to be statistically significant. For sex-age interactions, SCL22A12 was found to be significantly higher expressed in females <50 years compared to males <50 years. Expression levels of SLC22A2, SLC22A12, SLC6A16, and ABCB6 were significantly higher in females <50 years compared to females ≥50 years. In sex-ethnicity interactions, expression levels of ATP7B and KCNJ8 were found to be significantly higher in African American females compared to European American females. Also, the expression of SLC31A2 was significantly higher in European American males compared to European American females. CONCLUSIONS: Sex, age, and ethnic differences impacted the expression of drug transporters in normal human kidneys, which suggests that the analysis of gene expression of drug transporters will aid in improving the usage/dosage of drug therapies influencing personalized medicine and susceptibility to adverse drug reactions.