Epigenetics in drug disposition & drug therapy: symposium report of the 24th North American meeting of the International Society for the Study of Xenobiotics (ISSX).

The 24th North American International Society for the Study of Xenobiotics (ISSX) meeting, held virtually from September 13 to 17, 2021, embraced the theme of "Broadening Our Horizons." This reinforces a key mission of ISSX: striving to share innovative science related to drug discovery and development. Session speakers and the ISSX New Investigators Group, which supports the scientific and professional development of student and early career ISSX members, elected to highlight the scientific content presented during the captivating session titled, "Epigenetics in Drug Disposition & Drug Therapy." The impact genetic variation has on drug response is well established; however, this session underscored the importance of investigating the role of epigenetics in drug disposition and drug discovery. Session speakers, Drs. Ning, McClay, and Lazarus, detailed mechanisms by which epigenetic players including long non-coding RNA (lncRNAs), microRNA (miRNAs), DNA methylation, and histone acetylation can alter the expression of genes involved in pharmacokinetics, pharmacodynamics, and toxicity. Dr. Ning detailed current knowledge about miRNAs and lncRNAs and the mechanisms by which they can affect the expression of drug metabolizing enzymes (DMEs) and nuclear receptors. Dr. Lazarus discussed the potential role of miRNAs on UDP-glucuronosyltransferase (UGT) expression and activity. Dr. McClay provided evidence that aging alters methylation and acetylation of DMEs in the liver, affecting gene expression and activity. These topics, compiled by the symposium organizers, presenters, and the ISSX New Investigators Group, are herein discussed, along with exciting future perspectives for epigenetics in drug disposition and drug discovery research.

A comprehensive rat transcriptome built from large scale RNA-seq-based annotation.

The rat is an important model organism in biomedical research for studying human disease mechanisms and treatments, but its annotated transcriptome is far from complete. We constructed a Rat Transcriptome Re-annotation named RTR using RNA-seq data from 320 samples in 11 different organs generated by the SEQC consortium. Totally, there are 52 807 genes and 114 152 transcripts in RTR. Transcribed regions and exons in RTR account for ∼42% and ∼6.5% of the genome, respectively. Of all 73 074 newly annotated transcripts in RTR, 34 213 were annotated as high confident coding transcripts and 24 728 as high confident long noncoding transcripts. Different tissues rather than different stages have a significant influence on the expression patterns of transcripts. We also found that 11 715 genes and 15 852 transcripts were expressed in all 11 tissues and that 849 house-keeping genes expressed different isoforms among tissues. This comprehensive transcriptome is freely available at http://www.unimd.org/rtr/. Our new rat transcriptome provides essential reference for genetics and gene expression studies in rat disease and toxicity models.

Impact of Sequencing Depth and Library Preparation on Toxicological Interpretation of RNA-Seq Data in a "Three-Sample" Scenario.

While RNA-sequencing (RNA-seq) has emerged as a standard approach in toxicogenomics, its full potential in gaining underlying toxicological mechanisms is still not clear when only three biological replicates are used. This "three-sample" study design is common in toxicological research, particularly in animal studies during preclinical drug development. Sequencing depth (the total number of reads in an experiment) and library preparation are critical to the resolution and integrity of RNA-seq data and biological interpretation. We used aflatoxin b1 (AFB1), a model toxicant, to investigate the effect of sequencing depth and library preparation in RNA-seq on toxicological interpretation in the "three-sample" scenario. We also compared different gene profiling platforms (RNA-seq, TempO-seq, microarray, and qPCR) using identical liver samples. Well-established mechanisms of AFB1 toxicity served as ground truth for our comparative analyses. We found that a minimum of 20 million reads was sufficient to elicit key toxicity functions and pathways underlying AFB1-induced liver toxicity using three replicates and that identification of differentially expressed genes was positively associated with sequencing depth to a certain extent. Further, our results showed that RNA-seq revealed toxicological insights from pathway enrichment with overall higher statistical power and overlap ratio, compared with TempO-seq and microarray. Moreover, library preparation using the same methods was important to reproducing the toxicological interpretation.

MicroRNA hsa-miR-1301-3p Regulates Human ADH6, ALDH5A1 and ALDH8A1 in the Ethanol-Acetaldehyde-Acetate Metabolic Pathway.

Alcohol dehydrogenases (ADHs) and aldehyde dehydrogenases (ALDHs) are vital enzymes involved in the metabolism of a variety of alcohols. Differences in the expression and enzymatic activity of human ADHs and ALDHs correlate with individual variability in metabolizing alcohols and drugs and in the susceptibility to alcoholic liver disease. MicroRNAs (miRNAs) function as epigenetic modulators to regulate the expression of drug-metabolizing enzymes. To characterize miRNAs that target ADHs and ALDHs in human liver cells, we carried out a systematic bioinformatics analysis to analyze free energies of the interaction between miRNAs and their cognate sequences in ADH and ALDH transcripts and then calculated expression correlations between miRNAs and their targeting ADH and ALDH genes using a public data base. Candidate miRNAs were selected to evaluate bioinformatic predictions using a series of biochemical assays. Our results showed that 11 miRNAs have the potential to modulate the expression of two ADH and seven ALDH genes in the human liver. We found that hsa-miR-1301-3p suppressed the expression of ADH6, ALDH5A1, and ALDH8A1 in liver cells and blocked their induction by ethanol. In summary, our results revealed that hsa-miR-1301-3p plays an important role in ethanol metabolism by regulating ADH and ALDH gene expression. SIGNIFICANCE STATEMENT: Systematic bioinformatics analysis showed that 11 microRNAs might play regulatory roles in the expression of two alcohol dehydrogenase (ADH) and seven aldehyde dehydrogenase (ALDH) genes in the human liver. Experimental evidences proved that hsa-miR-1301-3p suppressed the expression of ADH6, ALDH5A1, and ALDH8A1 in liver cells and decreased their inducibility by ethanol.

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.

The role of hepatic cytochrome P450s in the cytotoxicity of sertraline.

Sertraline, an antidepressant, is commonly used to manage mental health symptoms related to depression, anxiety disorders, and obsessive-compulsive disorder. The use of sertraline has been associated with rare but severe hepatotoxicity. Previous research demonstrated that mitochondrial dysfunction, apoptosis, and endoplasmic reticulum stress were involved in sertraline-associated cytotoxicity. In this study, we reported that after a 24-h treatment in HepG2 cells, sertraline caused cytotoxicity, suppressed topoisomerase I and IIα, and damaged DNA in a concentration-dependent manner. We also investigated the role of cytochrome P450 (CYP)-mediated metabolism in sertraline-induced toxicity using our previously established HepG2 cell lines individually expressing 14 CYPs (1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, 3A5, and 3A7). We demonstrated that CYP2D6, 2C19, 2B6, and 2C9 metabolize sertraline, and sertraline-induced cytotoxicity was significantly decreased in the cells expressing these CYPs. Western blot analysis demonstrated that the induction of É£H2A.X (a hallmark of DNA damage) and topoisomerase inhibition were partially reversed in CYP2D6-, 2C19-, 2B6-, and 2C9-overexpressing HepG2 cells. These data indicate that DNA damage and topoisomerase inhibition are involved in sertraline-induced cytotoxicity and that CYPs-mediated metabolism plays a role in decreasing the toxicity of sertraline.

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.

Regulation of cytochrome P450 expression by microRNAs and long noncoding RNAs: Epigenetic mechanisms in environmental toxicology and carcinogenesis.

Environmental exposures to hazardous chemicals are associated with a variety of human diseases and disorders, including cancers. Phase I metabolic activation and detoxification reactions catalyzed by cytochrome P450 enzymes (CYPs) affect the toxicities of many xenobiotic compounds. Proper regulation of CYP expression influences their biological effects. Noncoding RNAs (ncRNAs) are involved in regulating CYP expression, and ncRNA expression is regulated in response to environmental chemicals. The mechanistic interactions between ncRNAs and CYPs associated with the toxicity and carcinogenicity of environmental chemicals are described in this review, focusing on microRNA-dependent CYP regulation. The role of long noncoding RNAs in regulating CYP expression is also presented and new avenues of research concerning this regulatory mechanism are described.

The role of hepatic cytochrome P450s in the cytotoxicity of dronedarone.

Dronedarone is used to treat patients with cardiac arrhythmias and has been reported to be associated with liver injury. Our previous mechanistic work demonstrated that DNA damage-induced apoptosis contributes to the cytotoxicity of dronedarone. In this study, we examined further the underlying mechanisms and found that after a 24-h treatment of HepG2 cells, dronedarone caused cytotoxicity, G1-phase cell cycle arrest, suppression of topoisomerase II, and DNA damage in a concentration-dependent manner. We also investigated the role of cytochrome P450s (CYPs)-mediated metabolism in the dronedarone-induced toxicity using our previously established HepG2 cell lines expressing individually 14 human CYPs (1A1, 1A2, 1B1, 2A6, 2B6, 2C8, 2C9, 2C18, 2C19, 2D6, 2E1, 3A4, 3A5, and 3A7). We demonstrated that CYP3A4, 3A5, and 2D6 were the major enzymes that metabolize dronedarone, and that CYP3A7, 2E1, 2C19, 2C18, 1A1, and 2B6 also metabolize dronedarone, but to a lesser extent. Our data showed that the cytotoxicity of dronedarone was decreased in CYP3A4-, 3A5-, or 2D6-overexpressing cells compared to the control HepG2 cells, indicating that the parent dronedarone has higher potency than the metabolites to induce cytotoxicity in these cells. In contrast, cytotoxicity was increased in CYP1A1-overexpressing cells, demonstrating that CYP1A1 exerts an opposite effect in dronedarone's toxicity, comparing to CYP3A4, 3A5, or 2D6. We also studied the involvement of topoisomerase II in dronedarone-induced toxicity, and demonstrated that the overexpression of topoisomerase II caused an increase in cell viability and a decrease in γ-H2A.X induction, suggesting that suppression of topoisomerase II may be one of the mechanisms involved in dronedarone-induced liver toxicity.

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 Development of a Database for Herbal and Dietary Supplement Induced Liver Toxicity.

The growing use of herbal dietary supplements (HDS) in the United States provides compelling evidence for risk of herbal-induced liver injury (HILI). Information on HDS products was retrieved from MedlinePlus of the U.S. National Library of Medicine and the herbal monograph of the European Medicines Agency. The hepatotoxic potential of HDS was ascertained by considering published case reports. Other relevant data were collected from governmental documents, public databases, web sources, and the literature. We collected information for 296 unique HDS products. Evidence of hepatotoxicity was reported for 67, that is 1 in 5, of these HDS products. The database revealed an apparent gender preponderance with women representing 61% of HILI cases. Culprit hepatotoxic HDS were mostly used for weight control, followed by pain and inflammation, mental stress, and mood disorders. Commonly discussed mechanistic events associated with HILI are reactive metabolites and oxidative stress, mitochondrial injury, as well as inhibition of transporters. HDS⁻drug interactions, causing both synergistic and antagonizing effects of drugs, were also reported for certain HDS. The database contains information for nearly 300 commonly used HDS products to provide a single-entry point for better comprehension of their impact on public health.

Activation of the Nrf2 signaling pathway in usnic acid-induced toxicity in HepG2 cells.

Many usnic acid-containing dietary supplements have been marketed as weight loss agents, although severe hepatotoxicity and acute liver failure have been associated with their overuse. Our previous mechanistic studies revealed that autophagy, disturbance of calcium homeostasis, and ER stress are involved in usnic acid-induced toxicity. In this study, we investigated the role of oxidative stress and the Nrf2 signaling pathway in usnic acid-induced toxicity in HepG2 cells. We found that a 24-h treatment with usnic acid caused DNA damage and S-phase cell cycle arrest in a concentration-dependent manner. Usnic acid also triggered oxidative stress as demonstrated by increased reactive oxygen species generation and glutathione depletion. Short-term treatment (6 h) with usnic acid significantly increased the protein level for Nrf2 (nuclear factor erythroid 2-related factor 2), promoted Nrf2 translocation to the nucleus, up-regulated antioxidant response element (ARE)-luciferase reporter activity, and induced the expression of Nrf2-regulated targets, including glutathione reductase, glutathione S-transferase, and NAD(P)H quinone oxidoreductase-1 (NQO1). Furthermore, knockdown of Nrf2 with shRNA potentiated usnic acid-induced DNA damage and cytotoxicity. Taken together, our results show that usnic acid causes cell cycle dysregulation, DNA damage, and oxidative stress and that the Nrf2 signaling pathway is activated in usnic acid-induced cytotoxicity.

Single nucleotide polymorphism rs11669203 in TGFBR3L is associated with the risk of neuroblastoma in a Chinese population.

With a primary mortality, neuroblastoma (NB) is the most common extracranial solid tumor in childhood. Amplification of the MYCN (v-myc avian myelocytomatosis viral oncogene neuroblastoma derived homolog) oncogene is observed in 20-30 % of NB cases, a feature which also characterizes a highly aggressive subtype of the disease. However, the systematic study of association between single nucleotide polymorphisms (SNPs) in MYCN-regulated genes and the risk of NB has not been investigated. In the current study, we scanned a set of 16 SNPs located within known or predicted MYCN binding sites in a cohort of 247 patients of Chinese origin with neuroblastic family tumors, including neuroblastoma (NB), ganglioneuroma (GN), and ganglioneuroblastoma (GNB), and in 290 cancer-free controls to determine whether any of the tested SNPs are associated with neuroblastic family tumors. We found that the rs11669203 G>C polymorphism, located in TGFBR3L promoter, is significantly associated with the risk of NB. Further, we found that this association is site specific to adrenal NB compared to non-adrenal NB. In addition, transcriptome analysis indicated that increased expression of TGFBR3L is strongly correlated with poor survival. The SNP rs11669203 located at the MYCN binding site of TGFBR3L is significantly associated with elevated risk of NB, and abnormal MYCN-regulated TGFBR3L expression may contribute to NB oncogenesis.

Incorporating the human gene annotations in different databases significantly improved transcriptomic and genetic analyses.

Human gene annotation is crucial for conducting transcriptomic and genetic studies; however, the impacts of human gene annotations in diverse databases on related studies have been less evaluated. To enable full use of various human annotation resources and better understand the human transcriptome, here we systematically compare the human annotations present in RefSeq, Ensembl (GENCODE), and AceView on diverse transcriptomic and genetic analyses. We found that the human gene annotations in the three databases are far from complete. Although Ensembl and AceView annotated more genes than RefSeq, more than 15,800 genes from Ensembl (or AceView) are within the intergenic and intronic regions of AceView (or Ensembl) annotation. The human transcriptome annotations in RefSeq, Ensembl, and AceView had distinct effects on short-read mapping, gene and isoform expression profiling, and differential expression calling. Furthermore, our findings indicate that the integrated annotation of these databases can obtain a more complete gene set and significantly enhance those transcriptomic analyses. We also observed that many more known SNPs were located within genes annotated in Ensembl and AceView than in RefSeq. In particular, 1033 of 3041 trait/disease-associated SNPs involved in about 200 human traits/diseases that were previously reported to be in RefSeq intergenic regions could be relocated within Ensembl and AceView genes. Our findings illustrate that a more complete transcriptome generated by incorporating human gene annotations in diverse databases can strikingly improve the overall results of transcriptomic and genetic studies.

Toxicogenomics and cancer susceptibility: advances with next-generation sequencing.

The aim of this review is to comprehensively summarize the recent achievements in the field of toxicogenomics and cancer research regarding genetic-environmental interactions in carcinogenesis and detection of genetic aberrations in cancer genomes by next-generation sequencing technology. Cancer is primarily a genetic disease in which genetic factors and environmental stimuli interact to cause genetic and epigenetic aberrations in human cells. Mutations in the germline act as either high-penetrance alleles that strongly increase the risk of cancer development, or as low-penetrance alleles that mildly change an individual's susceptibility to cancer. Somatic mutations, resulting from either DNA damage induced by exposure to environmental mutagens or from spontaneous errors in DNA replication or repair are involved in the development or progression of the cancer. Induced or spontaneous changes in the epigenome may also drive carcinogenesis. Advances in next-generation sequencing technology provide us opportunities to accurately, economically, and rapidly identify genetic variants, somatic mutations, gene expression profiles, and epigenetic alterations with single-base resolution. Whole genome sequencing, whole exome sequencing, and RNA sequencing of paired cancer and adjacent normal tissue present a comprehensive picture of the cancer genome. These new findings should benefit public health by providing insights in understanding cancer biology, and in improving cancer diagnosis and therapy.

Cytotoxicity and inhibitory effects of low-concentration triclosan on adipogenic differentiation of human mesenchymal stem cells.

Humans at all ages are continually exposed to triclosan (TCS), a widely used antimicrobial agent that can be found in many daily hygiene products, such as toothpastes and shampoos; however, the toxicological and biological effects of TCS in the human body after long-term and low-concentration exposure are far from being well understood. In the current study, we investigated the effects of TCS on the differentiation of human mesenchymal stem cells (hMSCs) by measuring the cytotoxicity, morphological changes, lipid accumulation, and the expression of adipocyte differentiation biomarkers during 21-day adipogenesis. Significant cytotoxicity was observed in un-induced hMSCs treated with high-concentration TCS (≥ 5.0 µM TCS), but not with low-concentration treatments (≤ 2.5 µM TCS). TCS inhibited adipocyte differentiation of hMSCs in a concentration-dependent manner in the 0.156 to 2.5 µM range as indicated by morphological changes with Oil Red O staining, which is an index of lipid accumulation. The inhibitory effect was confirmed by a decrease in gene expression of specific adipocyte differentiation biomarkers including adipocyte protein 2, lipoprotein lipase, and adiponectin. Our study demonstrates that TCS inhibits adipocyte differentiation of hMSCs under concentrations that are not cytotoxic and in the range observed in human blood.

Differential promoter activities of functional haplotypes in the 5'-flanking region of human sulfotransferase 1A1.

Previously, we reported five common single nucleotide polymorphisms (SNPs), -624G>C, -396G>A, -358A>C, -341C>G, and -294T>C, and six common haplotypes (CGACT, GAACT, GGAGC, GGACC, CAACT, and GAACC) in the 5'-flanking region of the SULT1A1 gene that were associated with altered enzymatic activity. In the present study, we performed in vitro assays to determine the functional impact of these genetic variations on the promoter activity. Dual luciferase reporter assays revealed that these SNPs are located in a negative regulatory fragment of the SULT1A1 gene. Further experiments demonstrated that these SNPs and haplotypes affected promoter activities of SULT1A1. Electrophoretic mobility shift assays showed distinctive binding patterns for the SNPs -396G>A and -294T>C, due to differential binding affinities of the G/A alleles and the T/C alleles to nuclear proteins extracted from the liver carcinoma cell lines, HepG2 and Huh7.

Integrative approaches for studying the role of noncoding RNAs in influencing drug efficacy and toxicity.

INTRODUCTION: Drug efficacy and toxicity are important factors for evaluation in drug development. Drug metabolizing enzymes and transporters (DMETs) play an essential role in drug efficacy and toxicity. Noncoding RNAs (ncRNAs) have been implicated to influence inter-individual variations in drug efficacy and safety by regulating DMETs. An efficient strategy is urgently needed to identify and functionally characterize ncRNAs that mediate drug efficacy and toxicity through regulating DMETs. AREAS COVERED: We outline an integrative strategy to identify ncRNAs that modulate DMETs. We include reliable tools and databases for computational prediction of ncRNA targets with regard to their advantages and limitations. Various biochemical, molecular, and cellular assays are discussed for in vitro experimental verification of the regulatory function of ncRNAs. In vivo approaches for association of ncRNAs with drug treatment and toxicity are also reviewed. EXPERT OPINION: A streamlined integration of computational prediction and wet-lab validation is important to elucidate mechanisms of ncRNAs in the regulation of DMETs related to drug efficacy and safety. Bioinformatic analyses using open-access tools and databases serve as a powerful booster for ncRNA Research in toxicology. Further refinement of computational algorithms and experimental technologies is needed to improve accuracy and efficiency in ncRNA target identification and characterization.

Assessing reproducibility of inherited variants detected with short-read whole genome sequencing.

BACKGROUND: Reproducible detection of inherited variants with whole genome sequencing (WGS) is vital for the implementation of precision medicine and is a complicated process in which each step affects variant call quality. Systematically assessing reproducibility of inherited variants with WGS and impact of each step in the process is needed for understanding and improving quality of inherited variants from WGS. RESULTS: To dissect the impact of factors involved in detection of inherited variants with WGS, we sequence triplicates of eight DNA samples representing two populations on three short-read sequencing platforms using three library kits in six labs and call variants with 56 combinations of aligners and callers. We find that bioinformatics pipelines (callers and aligners) have a larger impact on variant reproducibility than WGS platform or library preparation. Single-nucleotide variants (SNVs), particularly outside difficult-to-map regions, are more reproducible than small insertions and deletions (indels), which are least reproducible when > 5 bp. Increasing sequencing coverage improves indel reproducibility but has limited impact on SNVs above 30×. CONCLUSIONS: Our findings highlight sources of variability in variant detection and the need for improvement of bioinformatics pipelines in the era of precision medicine with WGS.

Similarities and differences in the expression of drug-metabolizing enzymes between human hepatic cell lines and primary human hepatocytes.

In addition to primary human hepatocytes, hepatoma cell lines, and transfected nonhepatoma, hepatic cell lines have been used for pharmacological and toxicological studies. However, a systematic evaluation and a general report of the gene expression spectra of drug-metabolizing enzymes and transporters (DMETs) in these in vitro systems are not currently available. To fill this information gap and to provide references for future studies, we systematically characterized the basal gene expression profiles of 251 drug-metabolizing enzymes in untreated primary human hepatocytes from six donors, four commonly used hepatoma cell lines (HepG2, Huh7, SK-Hep-1, and Hep3B), and one transfected human liver epithelial cell line. A large variation in DMET expression spectra was observed between hepatic cell lines and primary hepatocytes, with the complete absence or much lower abundance of certain DMETs in hepatic cell lines. Furthermore, the basal DMET expression spectra of five hepatic cell lines are summarized, providing references for researchers to choose carefully appropriate in vitro models for their studies of drug metabolism and toxicity, especially for studies with drugs in which toxicities are mediated through the formation of reactive metabolites.