RNA sequencing reveals dynamic changes of mRNA abundance of cytochromes P450 and their alternative transcripts during mouse liver development

Epigenetic Mechanisms Contribute to Intraindividual Variations of Drug Metabolism Mediated by Cytochrome P450 Enzymes

Significant interindividual and intraindividual variations on cytochrome P450 (CYP)-mediated drug metabolism exist in the general population globally. Genetic polymorphisms are one of the major contribution factors for interindividual variations, but epigenetic mechanisms mainly contribute to intraindividual variations, including DNA methylation, histone modifications, microRNAs, and long non-coding RNAs. The current review provides analysis of advanced knowledge in the last decade on contributions of epigenetic mechanisms to intraindividual variations on CYP-mediated drug metabolism in several situations, including (1) ontogeny, the developmental changes of CYP expression in individuals from neonates to adults; (2) increased activities of CYP enzymes induced by drug treatment; (3) increased activities of CYP enzymes in adult ages induced by drug treatment at neonate ages; and (4) decreased activities of CYP enzymes in individuals with drug-induced liver injury (DILI). Furthermore, current challenges, knowledge gaps, and future perspective of the epigenetic mechanisms in development of CYP pharmacoepigenetics are discussed. In conclusion, epigenetic mechanisms have been proven to contribute to intraindividual variations of drug metabolism mediated by CYP enzymes in age development, drug induction, and DILI conditions. The knowledge has helped understanding how intraindividual variation are generated. Future studies are needed to develop CYP-based pharmacoepigenetics to guide clinical applications for precision medicine with improved therapeutic efficacy and reduced risk of adverse drug reactions and toxicity. SIGNIFICANCE STATEMENT: Understanding epigenetic mechanisms in contribution to intraindividual variations of CYP-mediated drug metabolism may help to develop CYP-based pharmacoepigenetics for precision medicine to improve therapeutic efficacy and reduce adverse drug reactions and toxicity for drugs metabolized by CYP enzymes.

Prenatal Lipopolysaccharide Exposure Alters Hepatic Drug-Metabolizing Enzyme Expression in Mouse Offspring via Histone Modifications

Inflammation is a major regulator of drug-metabolizing enzymes (DMEs), therefore contributing to the interindividual variability of drug effects. However, whether prenatal inflammation affects DMEs expression in offspring remains obscure. This study investigated the effects of prenatal lipopolysaccharide (LPS) exposure on hepatic expression of inflammatory-related genes, nuclear receptors, and DMEs in offspring mice. Prenatal LPS exposure on gestational day (GD) 10 led to higher expression of NF-κB, Pxr, and Cyp2b10, while lower expression of Car, Ahr, Cyp3a11, and Ugt1a1 in postnatal day (PD) 30 offspring. However, multiple doses of LPS exposure on GD10-14 resulted in higher levels of inflammatory-related genes, Cyp1a2, and Cyp2b10, and lower levels of Pxr and Cyp3a11 in PD30 offspring liver. For PD60 offspring, decreased hepatic expression of NF-κB and IL-6, and increased expression of Pxr and Cyp3a11 were seen in single-dose LPS groups, whereas opposite results were observed in the multiple-dose LPS groups. Notably, enhanced H3K4me3 levels in the PXR response elements of the Cyp3a11 promoter were observed in the liver of PD60 offspring mice from dams treated with multiple doses of LPS during pregnancy. Overall, this study suggests that parental LPS exposure could persistently alter the hepatic expression of DMEs, and histone modifications may contribute to the long-term effects.

Dynamic Transcriptomic Profiling During Liver Development in Schizothorax Prenanti

Liver is an important organ for glucose and lipid metabolism, immunity, and detoxification in fish. However, the gene regulatory network of postnatal liver development still remains unknown in teleost fish. In this study, we performed transcriptome analysis on the liver of S. prenanti at three stages. A total of 1692 differentially expressed genes (DGEs) were identified across three liver developmental stages. The oil red O staining and PAS staining revealed that the lipid content of liver was increased and the glycogen content of liver was decreased during liver development. The fatty acids biosynthesis related genes were upregulated in adult and young stages compared with juvenile stage, while lipid degradation related genes were downregulated. The genes related to glycolysis, gluconeogenesis and glycogenolysis were upregulated in juvenile or young stages compared with adult stage. Further pathway analysis indicated that the CYP450 pathway, cell cycle and amino acid metabolic pathway were induced in the process of liver maturation. Our study presents the gene expression pattern in different liver development stages of S. prenanti and may guide future studies on metabolism of S. prenanti liver.

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

Aristolochic acids exposure was not the main cause of liver tumorigenesis in adulthood

Aristolochic acids (AAs) have long been considered as a potent carcinogen due to its nephrotoxicity. Aristolochic acid I (AAI) reacts with DNA to form covalent aristolactam (AL)–DNA adducts, leading to subsequent A to T transversion mutation, commonly referred as AA mutational signature. Previous research inferred that AAs were widely implicated in liver cancer throughout Asia. In this study, we explored whether AAs exposure was the main cause of liver cancer in the context of HBV infection in mainland China. Totally 1256 liver cancer samples were randomly retrieved from 3 medical centers and a refined bioanalytical method was used to detect AAI–DNA adducts. 5.10% of these samples could be identified as AAI positive exposure. Whole genome sequencing suggested 8.41% of 107 liver cancer patients exhibited the dominant AA mutational signature, indicating a relatively low overall AAI exposure rate. In animal models, long-term administration of AAI barely increased liver tumorigenesis in adult mice, opposite from its tumor-inducing role when subjected to infant mice. Furthermore, AAI induced dose-dependent accumulation of AA–DNA adduct in target organs in adult mice, with the most detected in kidney instead of liver. Taken together, our data indicate that AA exposure was not the major threat of liver cancer in adulthood.

Effects of Overexpression of Fibroblast Growth Factor 15/19 on Hepatic Drug Metabolizing Enzymes

Fibroblast growth factors 15 (FGF15) and 19 (FGF19) are endocrine growth factors that play an important role in maintaining bile acid homeostasis. FGF15/19-based therapies are currently being tested in clinical trials for the treatment of nonalcoholic steatohepatitis and cholestatic liver diseases. To determine the physiologic impact of long-term elevations of FGF15/19, a transgenic mouse model with overexpression of Fgf15 (Fgf15 Tg) was used in the current study. The RNA sequencing (RNA-seq) analysis revealed elevations of the expression of several genes encoding phase I drug metabolizing enzymes (DMEs), including Cyp2b10 and Cyp3a11, in Fgf15 Tg mice. We found that the induction of several Cyp2b isoforms resulted in increased function of CYP2B in microsomal metabolism and pharmacokinetics studies. Because the CYP2B family is known to be induced by constitutive androstane receptor (CAR), to determine the role of CAR in the observed inductions, we crossed Fgf15 Tg mice with CAR knockout mice and found that CAR played a minor role in the observed alterations in DME expression. Interestingly, we found that the overexpression of Fgf15 in male mice resulted in a phenotypical switch from the male hepatic expression pattern of DMEs to that of female mice. Differences in secretion of growth hormone (GH) between male and female mice are known to drive sexually dimorphic, STAT5b-dependent expression patterns of hepatic genes. We found that male Fgf15 Tg mice presented with many features similar to GH deficiency, including lowered body length and weight, Igf-1 and Igfals expression, and STAT5 signaling. SIGNIFICANCE STATEMENT: The overexpression of Fgf15 in mice causes an alteration in DMEs at the mRNA, protein, and functional levels, which is not entirely due to CAR activation but associated with lower GH signaling.

Targeting Peroxisome Proliferator-Activated Receptor-β/δ (PPARβ/δ) for the Treatment or Prevention of Alcoholic Liver Disease

Excessive, chronic alcohol consumption can lead to alcoholic liver disease. The etiology of alcoholic liver disease is multifactorial and is influenced by alterations in gene expression and changes in fatty acid metabolism, oxidative stress, and insulin resistance. These events can lead to steatosis, fibrosis, and eventually to cirrhosis and liver cancer. Many of these functions are regulated by peroxisome proliferator-activated receptors (PPARs). Thus, it is not surprising that PPARs can modulate the mechanisms that cause alcoholic liver disease. While the roles of PPARα and PPARγ are clearer, the role of PPARβ/δ in alcoholic liver disease requires further clarification. This review summarizes the current understanding based on recent studies that indicate that PPARβ/δ can likely be targeted for the treatment and/or the prevention of alcoholic liver disease.

Ontogeny of Drug-Metabolizing Enzymes

Almost 50% of prescription drugs lack age-appropriate dosing guidelines and therefore are used "off-label." Only ~10% drugs prescribed to neonates and infants have been studied for safety or efficacy. Immaturity of drug metabolism in children is often associated with drug toxicity. This chapter summarizes data on the ontogeny of major human metabolizing enzymes involved in oxidation, reduction, hydrolysis, and conjugation of drugs. The ontogeny data of individual drug-metabolizing enzymes are important for accurate prediction of drug pharmacokinetics and toxicity in children. This information is critical for designing clinical studies to appropriately test pharmacological hypotheses and develop safer pediatric drugs, and to replace the long-standing practice of body weight- or surface area-normalized drug dosing. The application of ontogeny data in physiologically based pharmacokinetic model and regulatory submission are discussed.

Alternative RNA Splicing in Fatty Liver Disease

Alternative RNA splicing is a process by which introns are removed and exons are assembled to construct different RNA transcript isoforms from a single pre-mRNA. Previous studies have demonstrated an association between dysregulation of RNA splicing and a number of clinical syndromes, but the generality to common disease has not been established. Non-alcoholic fatty liver disease (NAFLD) is the most common liver disease affecting one-third of adults worldwide, increasing the risk of cirrhosis and hepatocellular carcinoma (HCC). In this review we focus on the change in alternative RNA splicing in fatty liver disease and the role for splicing regulation in disease progression.

Ontogeny of Hepatic Transporters and Drug-Metabolizing Enzymes in Humans and in Nonclinical Species

The liver represents a major eliminating and detoxifying organ, determining exposure to endogenous compounds, drugs, and other xenobiotics. Drug transporters (DTs) and drug-metabolizing enzymes (DMEs) are key determinants of disposition, efficacy, and toxicity of drugs. Changes in their mRNA and protein expression levels and associated functional activity between the perinatal period until adulthood impact drug disposition. However, high-resolution ontogeny profiles for hepatic DTs and DMEs in nonclinical species and humans are lacking. Meanwhile, increasing use of physiologically based pharmacokinetic (PBPK) models necessitates availability of underlying ontogeny profiles to reliably predict drug exposure in children. In addition, understanding of species similarities and differences in DT/DME ontogeny is crucial for selecting the most appropriate animal species when studying the impact of development on pharmacokinetics. Cross-species ontogeny mapping is also required for adequate translation of drug disposition data in developing nonclinical species to humans. This review presents a quantitative cross-species compilation of the ontogeny of DTs and DMEs relevant to hepatic drug disposition. A comprehensive literature search was conducted on PubMed Central: Tables and graphs (often after digitization) in original manuscripts were used to extract ontogeny data. Data from independent studies were standardized and normalized before being compiled in graphs and tables for further interpretation. New insights gained from these high-resolution ontogeny profiles will be indispensable to understand cross-species differences in maturation of hepatic DTs and DMEs. Integration of these ontogeny data into PBPK models will support improved predictions of pediatric hepatic drug disposition processes. SIGNIFICANCE STATEMENT: Hepatic drug transporters (DTs) and drug-metabolizing enzymes (DMEs) play pivotal roles in hepatic drug disposition. Developmental changes in expression levels and activities of these proteins drive age-dependent pharmacokinetics. This review compiles the currently available ontogeny profiles of DTs and DMEs expressed in livers of humans and nonclinical species, enabling robust interpretation of age-related changes in drug disposition and ultimately optimization of pediatric drug therapy.

Placental BCRP/ABCG2 Transporter Prevents Fetal Exposure to the Estrogenic Mycotoxin Zearalenone

In the placenta, the breast cancer resistance protein (BCRP)/ABCG2 efflux transporter limits the maternal-to-fetal transfer of drugs and chemicals. Previous research has pointed to the estrogenic mycotoxin zearalenone as a potential substrate for BCRP. Here, we sought to assess the role of the BCRP transporter in the transplacental disposition of zearalenone during pregnancy. In vitro transwell transport assays employing BCRP/Bcrp-transfected Madine-Darby canine kidney cells and BeWo trophoblasts with reduced BCRP expression were used to characterize the impact of BCRP on the bidirectional transport of zearalenone. In both models, the presence of BCRP protein increased the basolateral-to-apical transport and reduced the apical-to-basolateral transport of zearalenone over a 2-h period. In vivo pharmacokinetic analyses were then performed using pregnant wild-type and Bcrp-/- mice after a single tail vein injection of zearalenone. Zearalenone and its metabolite α-zearalenol were detectable in serum, placentas, and fetuses from all animals, and β-zearalenol was detected in serum and fetuses, but not placentas. There were no significant differences in the maternal serum concentrations of any analytes between the two genotypes. In Bcrp-/- mice, the free fetal concentrations of zearalenone, α-zearalenol, and β-zearalenol were increased by 115%, 84%, and 150%, respectively, when compared with wild-type mice. Concentrations of free zearalenone and α-zearalenol were elevated 145% and 78% in Bcrp-/- placentas, respectively, when compared with wild-type placentas. Taken together, these data indicate that the placental BCRP transporter functions to reduce the fetal accumulation of zearalenone, which may impact susceptibility to developmental toxicities associated with in utero zearalenone exposure.

Acetaminophen-Induced Liver Injury Alters Expression and Activities of Cytochrome P450 Enzymes in an Age-Dependent Manner in Mouse Liver

Drug-induced liver injury (DILI) is a global medical problem. The risk of DILI is often related to expression and activities of drug-metabolizing enzymes, especially cytochrome P450s (P450s). However, changes on expression and activities of P450s after DILI have not been determined. The aim of this study is to fill this knowledge gap. Acetaminophen (APAP) was used as a model drug to induce DILI in C57BL/6J mice at different ages of days 10 (infant), 22 (child), and 60 (adult). DILI was assessed by levels of alanine aminotransferase and aspartate aminotransferase in plasma with a confirmation by H&E staining on liver tissue sections. The expression of selected P450s at mRNA and protein levels was measured by real-time polymerase chain reaction and liquid chromatography-tandem mass spectrometry, respectively. The activities of these P450s were determined by the formation of metabolites from probe drugs for each P450 using ultraperformance liquid chromatography-quadrupole time of flight mass spectrometry. DILI was induced at mild to severe levels in a dose-dependent manner in 200, 300, and 400 mg/kg APAP-treated groups at child and adult ages, but not at the infant age. Significantly decreased expression at mRNA and protein levels as well as enzymatic activities of CYP2E1, 3A11, 1A2, and 2C29 were found at child and adult ages. Adult male mice were more susceptible to APAP-induced liver injury than female mice with more decreased expression of P450s. These results suggest that altered levels of P450s in livers severely injured by drugs may affect the therapeutic efficacy of drugs, which are metabolized by P450s, more particularly for males. SIGNIFICANCE STATEMENT: The current study in an animal model demonstrates that acetaminophen-induced liver injury results in decreased expression and enzyme activities of several examined drug-metabolizing cytochrome P450s (P450s). The extent of such decreases is correlated to the degree of liver injury severity. The generated data may be translated to human health for patients who have drug-induced liver injury with decreased capability to metabolize drugs by certain P450s.

A time-resolved multi-omic atlas of the developing mouse liver

Liver organogenesis and development are composed of a series of complex, well-orchestrated events. Identifying key factors and pathways governing liver development will help elucidate the physiological and pathological processes including those of cancer. We conducted multidimensional omics measurements including protein, mRNA, and transcription factor (TF) DNA-binding activity for mouse liver tissues collected from embryonic day 12.5 (E12.5) to postnatal week 8 (W8), encompassing major developmental stages. These data sets reveal dynamic changes of core liver functions and canonical signaling pathways governing development at both mRNA and protein levels. The TF DNA-binding activity data set highlights the importance of TF activity in early embryonic development. A comparison between mouse liver development and human hepatocellular carcinoma (HCC) proteomic profiles reveal that more aggressive tumors are characterized with the activation of early embryonic development pathways, whereas less aggressive ones maintain liver function-related pathways that are elevated in the mature liver. This work offers a panoramic view of mouse liver development and provides a rich resource to explore in-depth functional characterization.

Impact of obese levels on the hepatic expression of nuclear receptors and drug-metabolizing enzymes in adult and offspring mice

The prevalence of obesity-associated conditions raises new challenges in clinical medication. Although altered expression of drug-metabolizing enzymes (DMEs) has been shown in obesity, the impacts of obese levels (overweight, obesity, and severe obesity) on the expression of DMEs have not been elucidated. Especially, limited information is available on whether parental obese levels affect ontogenic expression of DMEs in children. Here, a high-fat diet (HFD) and three feeding durations were used to mimic different obese levels in C57BL/6 mice. The hepatic expression of five nuclear receptors (NRs) and nine DMEs was examined. In general, a trend of induced expression of NRs and DMEs (except for Cyp2c29 and 3a11) was observed in HFD groups compared to low-fat diet (LFD) groups. Differential effects of HFD on the hepatic expression of DMEs were found in adult mice at different obese levels. Family-based dietary style of an HFD altered the ontogenic expression of DMEs in the offspring older than 15 days. Furthermore, obese levels of parental mice affected the hepatic expression of DMEs in offspring. Overall, the results indicate that obese levels affected expression of the DMEs in adult individuals and that of their children. Drug dosage might need to be optimized based on the obese levels.

Chemical Activation of the Constitutive Androstane Receptor Leads to Activation of Oxidant-Induced Nrf2

Exposure to environmentally relevant chemicals that activate the xenobiotic receptors aryl hydrocarbon receptor (AhR), constitutive androstane receptor (CAR), and peroxisome proliferator-activated receptor alpha (PPARα) in rodent test systems often leads to increases in oxidative stress (OS) that contributes to liver cancer induction. We hypothesized that activation of the oxidant-induced transcription factor Nrf2 could be used as a surrogate endpoint for increases in OS. We examined the relationships between activation of xenobiotic receptors and Nrf2 using previously characterized gene expression biomarkers that accurately predict modulation. Using a correlation approach (Running Fisher Test), the biomarkers were compared with microarray profiles in a mouse liver gene expression compendium. Out of the 163 chemicals examined, 47% from 53 studies activated Nrf2. We found consistent coupling between CAR and Nrf2 activation. Out of the 41 chemicals from 32 studies that activated CAR, 90% also activated Nrf2. CAR was activated earlier and at lower doses than Nrf2, indicating CAR activation preceded Nrf2 activation. Nrf2 activation by 2 CAR activators was abolished in CAR-null mice. We hypothesized that Nrf2 is activated by reactive oxygen species from the increased activity of enzymes encoded by Cyp2b family members. However, Nrf2 was similarly activated in the livers of both TCPOBOP-treated wild-type and Cyp2b9/10/13-null mice. This study provides evidence that Nrf2 activation (1) often occurs after exposure to xenobiotic chemicals, (2) is tightly linked to activation of CAR, and (3) does not require induction of 3 Cyp2b genes secondary to CAR activation.

Epigenetic Memory Is Involved in the Persistent Alterations of Drug-Processing Genes in Adult Mice due to PCN-Activated PXR during Early Life

Pregnane X receptor (PXR), which can be activated by xenobiotic chemicals (including pediatric drugs), plays a key role in the regulation of drug-processing genes (DPGs). The induction of DPGs due to PXR activation may reduce therapeutic efficacy or cause toxicity. This work aims to demonstrate the impact of pregnenolone 16α-carbonitrile (PCN)-mediated PXR activation during early life on DPGs expression and drug sensitivity in adulthood, as well as the underlying mechanism. In this study, mice were sacrificed at postnatal day 60 to detect the hepatic expression of selected DPGs and histone modifications in the Cyp3a11 promoter. We found that all doses of PCN treatment (50-200 mg/kg/day) at postnatal days 5-8 resulted in persistently increased CYP2B10 expression, whereas only high doses of PCN treatment (150 and 200 mg/kg/day) persistently induced the expression of CYP3A11, 1A2, and UGTA1A1. We also demonstrated that PCN treatment before postnatal day 15 had a long-term impact on the expression of CYP3A11, 2B10, ABCC4, and PAPSS2. Additionally, elevated expression of CYP3A11, SULT2A1, UGT1A1, and PAPSS2 was observed in PCN-treated groups at days 25-28. Attenuated inducibility of CYP3A11 by PCN was seen in the primary hepatocytes derived from PCN-pretreated mice. Moreover, enhanced H3K4me3 level and reduced H3K27me3 level in the PXR response elements (PXREs) of the Cyp3a11 promoter may contribute to the persistent up-regulation of CYP3A11 by neonatal PCN treatment. Overall, our study suggests that PXR activation during early life could persistently alter the hepatic expression of DPGs and epigenetic memory may be an underlying mechanism in mice.

Cyp2b-null male mice are susceptible to diet-induced obesity and perturbations in lipid homeostasis

Obesity is an endemic problem in the United States and elsewhere, and data indicate that in addition to overconsumption, exposure to specific chemicals enhances obesity. CYP2B metabolizes multiple endo- and xenobiotics, and recent data suggests that repression of Cyp2b activity increases dyslipidemia and age-onset obesity, especially in males. To investigate the role played by Cyp2b in lipid homeostasis and obesity, we treated wildtype and Cyp2b-null mice with a normal (ND) or 60% high-fat diet (HFD) for 10 weeks and determined metabolic and molecular changes. Male HFD-fed Cyp2b-null mice weigh 15% more than HFD-fed wildtype mice, primarily due to an increase in white adipose tissue (WAT); however, Cyp2b-null female mice did not demonstrate greater body mass or WAT. Serum parameters indicate increased ketosis, leptin and cholesterol in HFD-fed Cyp2b-null male mice compared to HFD-fed wildtype mice. Liver triglycerides and liver:serum triglyceride ratios were higher than their similarly treated wildtype counterparts in Cyp2b-null male mice, indicating a role for Cyp2b in fatty acid metabolism regardless of diet. Furthermore, RNAseq demonstrates that hepatic gene expression in ND-fed Cyp2b-null male mice is similar to HFD-fed WT male mice, suggestive of fatty liver disease progression and a role for Cyp2b in lipid homeostasis. Females did not show as demonstrative changes in liver health, and significantly fewer changes in gene expression, as well as gene expression associated with liver disease. Overall our data indicates that the repression or inhibition of CYP2B may exacerbate metabolic disorders and cause obesity by perturbing fatty acid metabolism, especially in males.

RNA-Seq provides new insights on the relative mRNA abundance of antioxidant components during mouse liver development

Mammals have developed a variety of antioxidant systems to protect them from the oxygen environment and toxic stimuli. Little is known about the mRNA abundance of antioxidant components during postnatal development of the liver. Therefore, the purpose of this study was to compare the mRNA abundance of antioxidant components during liver development. Livers from male C57BL/6J mice were collected at 12 ages from prenatal to adulthood. The transcriptome was determined by RNA-Seq with transcript abundance estimated by Cufflinks. RNA-Seq provided a complete, more accurate, and unbiased quantification of the transcriptome. Among 33 known antioxidant components examined, three ontogeny patterns of liver antioxidant components were observed: (1) Prenatal-enriched, in which the mRNAs decreased from fetal livers to adulthood, such as metallothionein and heme oxygenase-1; (2) adolescent-rich and relatively stable expression, such as peroxiredoxins; and (3) adult-rich, in which the mRNA increased with age, such as catalase and superoxide dismutase. Alternative splicing of several antioxidant genes, such as Keap1, Glrx2, Gpx3, and Txnrd1, were also detected by RNA-Seq. In summary, RNA-Seq revealed the relative abundance of hepatic antioxidant enzymes, which are important in protecting against the deleterious effects of oxidative stress.

Molecular regulation of mammalian hepatic architecture

The essential liver exocrine and endocrine functions require a precise spatial arrangement of the hepatic lobule consisting of the central vein, portal vein, hepatic artery, intrahepatic bile duct system, and hepatocyte zonation. This allows blood to be carried through the liver parenchyma sampled by all hepatocytes and bile produced by the hepatocytes to be carried out of the liver through the intrahepatic bile duct system composed of cholangiocytes. The molecular orchestration of multiple signaling pathways and epigenetic factors is required to set up lineage restriction of the bipotential hepatoblast progenitor into the hepatocyte and cholangiocyte cell lineages, and to further refine cell fate heterogeneity within each cell lineage reflected in the functional heterogeneity of hepatocytes and cholangiocytes. In addition to the complex molecular regulation, there is a complicated morphogenetic choreography observed in building the refined hepatic epithelial architecture. Given the multifaceted molecular and cellular regulation, it is not surprising that impairment of any of these processes can result in acute and chronic hepatobiliary diseases. To enlighten the development of potential molecular and cellular targets for therapeutic options, an understanding of how the intricate hepatic molecular and cellular interactions are regulated is imperative. Here, we review the signaling pathways and epigenetic factors regulating hepatic cell lineages, fates, and epithelial architecture.

Cyp2b-Knockdown Mice Poorly Metabolize Corn Oil and Are Age-Dependent Obese

We previously made a RNAi-based cytochrome P450 2b (Cyp2b)-knockdown (Cyp2b-KD) mouse to determine the in vivo role of the Cyp2b subfamily in xenobiotic detoxification. Further studies reported here indicate a role for Cyp2b in unsaturated fatty-acid (UFA) metabolism and in turn obesity. Mice were treated intraperitoneally (i.p.) with 100 μL corn oil as a carrier or the potent Cyp2b-inducer 3,3',5,5'-Tetrachloro-1,4-bis(pyridyloxy)benzene (TCPOBOP (TC)) dissolved in corn oil. Surprisingly, female Cyp2b-KD mice but not male mice showed increased liver lipid accumulation. Male Cyp2b-KD mice had higher serum triacylglycerols, cholesterol, very low-density lipoprotein (VLDL), low-density lipoprotein (LDL), and high-density lipoprotein (HDL) than wildtype (WT) mice; females had higher cholesterol, LDL, and HDL. Thus, Cyp2b-KD mice are unable to clear a high bolus dose of corn oil, potentially because the Cyp2b-KD mice were unable to metabolize the UFA in the corn oil. Therefore, WT and Cyp2b-KD mice were housed for 35 weeks and necropsies performed to test whether Cyp2b-KD mice develop age onset obesity. Cyp2b-KD mice exhibited a significant increase in body weight caused by an increase in white adipose tissue deposition relative to WT mice. Serum cholesterol, triacylglycerol, LDL, and VLDL were significantly greater in 35-week-old Cyp2b-KD males compared to WT males; only serum triacylglycerol and LDL were higher in females. In conclusion, changes in Cyp2b expression led to perturbation in lipid metabolism and depuration in Cyp2b-KD mice. This suggests that Cyp2b is more than a detoxification enzyme, but also involved in the metabolism of UFA, as Cyp2b-KD mice have increased the body weight, fat deposition, and serum lipids.

Identification and characterization of long noncoding RNAs and mRNAs expression profiles related to postnatal liver maturation of breeder roosters using Ribo-zero RNA sequencing

BACKGROUND

The liver is mainly hematopoietic in the embryo, and converts into a major metabolic organ in the adult. Therefore, it is intensively remodeled after birth to adapt and perform adult functions. Long non-coding RNAs (lncRNAs) are involved in organ development and cell differentiation, likely they have potential roles in regulating postnatal liver development. Herein, in order to understand the roles of lncRNAs in postnatal liver maturation, we analyzed the lncRNAs and mRNAs expression profiles in immature and mature livers from one-day-old and adult (40 weeks of age) breeder roosters by Ribo-Zero RNA-Sequencing.

RESULTS

Around 21,939 protein-coding genes and 2220 predicted lncRNAs were expressed in livers of breeder roosters. Compared to protein-coding genes, the identified chicken lncRNAs shared fewer exons, shorter transcript length, and significantly lower expression levels. Notably, in comparison between the livers of newborn and adult breeder roosters, a total of 1570 mRNAs and 214 lncRNAs were differentially expressed with the criteria of log2fold change > 1 or < - 1 and P values < 0.05, which were validated by qPCR using randomly selected five mRNAs and five lncRNAs. Further GO and KEGG analyses have revealed that the differentially expressed mRNAs were involved in the hepatic metabolic and immune functional changes, as well as some biological processes and pathways including cell proliferation, apoptotic and cell cycle that are implicated in the development of liver. We also investigated the cis- and trans- regulatory effects of differentially expressed lncRNAs on its target genes. GO and KEGG analyses indicated that these lncRNAs had their neighbor protein coding genes and trans-regulated genes associated with adapting of adult hepatic functions, as well as some pathways involved in liver development, such as cell cycle pathway, Notch signaling pathway, Hedgehog signaling pathway, and Wnt signaling pathway.

CONCLUSIONS

This study provides a catalog of mRNAs and lncRNAs related to postnatal liver maturation of chicken, and will contribute to a fuller understanding of biological processes or signaling pathways involved in significant functional transition during postnatal liver development that differentially expressed genes and lncRNAs could take part in.

Consequences of Phenytoin Exposure on Hepatic Cytochrome P450 Expression during Postnatal Liver Maturation in Mice

The induction of cytochrome P450 (P450) enzymes in response to drug treatment is a significant contributing factor to drug-drug interactions, which may reduce therapeutic efficacy and/or cause toxicity. Since most studies on P450 induction are performed in adults, enzyme induction at neonatal, infant, and adolescent ages is not well understood. Previous work defined the postnatal ontogeny of drug-metabolizing P450s in human and mouse livers; however, there are limited data on the ontogeny of the induction potential of each enzyme in response to drug treatment. Induction of P450s at the neonatal age may also cause permanent alterations in P450 expression in adults. The goal of this study was to investigate the short- and long-term effects of phenytoin treatment on mRNA and protein expressions and enzyme activities of CYP2B10, 2C29, 3A11, and 3A16 at different ages during postnatal liver maturation in mice. Induction of mRNA immediately following phenytoin treatment appeared to depend on basal expression of the enzyme at a specific age. While neonatal mice showed the greatest fold changes in CYP2B10, 2C29, and 3A11 mRNA expression following treatment, the levels of induced protein expression and enzymatic activity were much lower than that of induced levels in adults. The expression of fetal CYP3A16 was repressed by phenytoin treatment. Neonatal treatment with phenytoin did not permanently induce enzyme expression in adulthood. Taken together, our data suggest that inducibility of drug-metabolizing P450s is much lower in neonatal mice than it is in adults and neonatal induction by phenytoin is not permanent.

Relationship between low midazolam metabolism by cytochrome P450 3A in mice and the high incidence of birth defects

The use of midazolam in early stages of pregnancy has resulted in a high incidence of birth defects; however, the underlying reason is unknown. We investigated expression changes of the CYP3A molecular species and focused on its midazolam metabolizing activity from the foetal period to adulthood. CYP3A16 was the only CYP3A species found to be expressed in the liver during the foetal period. However, CYP3A11 is upregulated in adult mice, but has been found to be downregulated during the foetal period and to gradually increase after birth. When CYP3A16 expression was induced in a microsomal fraction of cells used to study midazolam metabolism by CYP3A16, its activity was suppressed. These results showed that the capacity to metabolize midazolam in the liver during the foetal period is very low, which could hence result in a high incidence of birth defects associated with the use of midazolam during early stages of pregnancy.

De novo formation of the biliary system by TGFβ-mediated hepatocyte transdifferentiation

Transdifferentiation is a complete and stable change in cell identity that serves as an alternative to stem-cell-mediated organ regeneration. In adult mammals, findings of transdifferentiation have been limited to the replenishment of cells lost from preexisting structures, in the presence of a fully developed scaffold and niche¹. Here we show that transdifferentiation of hepatocytes in the mouse liver can build a structure that failed to form in development-the biliary system in a mouse model that mimics the hepatic phenotype of human Alagille syndrome (ALGS)². In these mice, hepatocytes convert into mature cholangiocytes and form bile ducts that are effective in draining bile and persist after the cholestatic liver injury is reversed, consistent with transdifferentiation. These findings redefine hepatocyte plasticity, which appeared to be limited to metaplasia, that is, incomplete and transient biliary differentiation as an adaptation to cell injury, based on previous studies in mice with a fully developed biliary system^3-6. In contrast to bile duct development^7-9, we show that de novo bile duct formation by hepatocyte transdifferentiation is independent of NOTCH signalling. We identify TGFβ signalling as the driver of this compensatory mechanism and show that it is active in some patients with ALGS. Furthermore, we show that TGFβ signalling can be targeted to enhance the formation of the biliary system from hepatocytes, and that the transdifferentiation-inducing signals and remodelling capacity of the bile-duct-deficient liver can be harnessed with transplanted hepatocytes. Our results define the regenerative potential of mammalian transdifferentiation and reveal opportunities for the treatment of ALGS and other cholestatic liver diseases.

Effect of the active ingredient of Kaempferia parviflora, 5,7-dimethoxyflavone, on the pharmacokinetics of midazolam

5,7-Dimethoxyflavone (5,7-DMF), one of the major components of Kaempferia parviflora, has anti-obesity, anti-inflammatory, and antineoplastic effects. On the other hand, in vitro studies have reported that it directly inhibits the drug metabolizing enzyme family cytochrome P450 (CYP) 3As. In this study, its safety was evaluated from a pharmacokinetic point of view, based on daily ingestion of 5,7-DMF. Midazolam, a substrate of CYP3As, was orally administered to mice treated with 5,7-DMF for 10 days, and its pharmacokinetic properties were investigated. In the group administered 5,7-DMF, the area under the curve (AUC) of midazolam increased by 130% and its biological half-life was extended by approximately 100 min compared to the control group. Compared to the control group, 5,7-DMF markedly decreased the expression of CYP3A11 and CYP3A25 in the liver. These results suggest that continued ingestion of 5,7-DMF decreases the expression of CYP3As in the liver, consequently increasing the blood concentrations of drugs metabolized by CYP3As.

A CYP2B6-humanized mouse model and its potential applications

CYP2B6 is a human microsomal cytochrome P450 enzyme with broad substrate selectivity. CYP2B6 is the only functional member of the human CYP2B gene subfamily, which differs from the situation in rodents, such as mouse, where multiple functional Cyp2b genes are expressed. Recent studies with Cyp2b knockout or knockdown mouse models have yielded insights into the in vivo roles of mouse CYP2B enzymes in drug disposition and xenobiotic toxicity. A CYP2B6-humanized mouse model (CYP2A13/2B6/2F1-transgenic/Cyp2abfgs-null), which expresses human CYP2B6 in the liver, and human CYP2A13 and CYP2F1 in the respiratory tract, but not any of the mouse Cyp2b genes, has also been established. In the CYP2B6-humanized mouse, the CYP2B6 transgene is expressed primarily in the liver, where it was found to be active toward prototype CYP2B6 substrate drugs. The regulatory elements of the CYP2B6 transgene appear to be compatible with mouse nuclear receptors that mediate CYP2B induction. Therefore, the CYP2B6-humanized mouse is a valuable animal model for studying the impact of CYP2B6 expression or induction on drug metabolism, drug efficacy, drug-drug interaction, and drug/xenobiotic toxicity. In this mini-review, we provide a brief background on CYP2B6 and the Cyp2b-knockout and CYP2B6-humanized mice, and discuss the potential applications and limitations of the current models.

The role of H19, a long non-coding RNA, in mouse liver postnatal maturation

H19 RNA is highly expressed at early postnatal ages and precipitously decreases at a specific time corresponding with increases in expression of genes important for mature liver function, such as drug metabolizing enzymes. H19’s role in the regulation of liver maturation is currently unknown. Using an H19 knockout mouse model to determine the role of H19 in liver development, we quantified gene expression for insulin growth factor signaling, Wnt signaling, key cytochrome P450 (P450) enzymes known to change as the liver develops, and fetal and adult plasma protein produced in liver. In mice lacking H19 expression, liver weights were significantly increased immediately after birth and significant increases were found in the number of actively proliferating cells. Increases in cell proliferation may be due to increases in β-catenin protein affecting Wnt signaling, increases in insulin-like growth factor 2 (IGF2) expression, and/or increases in insulin-like growth factor 1 receptor (IGF1R) expression at the protein level. Loss of targeted inhibition of IGF1R by microRNA 675 (miR-675) may be the cause of IGF1R increases, as miR-675 expression is also abrogated with loss of H19 expression in our model. P450 expression patterns were largely unchanged. No change in the production of plasma proteins was found, indicating H19 may not be important for liver maturation despite its role in controlling cell proliferation during liver growth. H19 may be important for normal liver development, and understanding how the liver matures will assist in predicting drug efficacy and toxicity in pediatric populations.

Role of Farnesoid X Receptor in the Determination of Liver Transcriptome during Postnatal Maturation in Mice

The liver is a vital organ with critical functions in metabolism of various biologically useful materials, synthesis of several vital proteins, detoxification of toxic substances, and immune defense. Most liver functions are not mature at birth and many changes happen during postnatal liver development, which lead to differential vulnerabilities of the liver at different developmental stages. However, the details of what changes occur in liver after birth, at what developmental stages they occur, and molecular mechanisms in the regulation of the developmental process are not clearly known. The nuclear receptor Farnesoid X receptor (FXR) is an important transcriptional regulator in liver. Here, we used RNA-Sequencing to analyze the transcriptome of mouse liver from perinatal to adult ages in both C57BL/6 and Fxr-/- mice. We have defined a clear timeline of functional transition from prenatal through neonatal and adolescent to adult in C57BL/6 mice. Without FXR, activation of neonatal-specific pathways was prolonged and maturation of multiple metabolic pathways was delayed. The loss of FXR also led to increased expression of 27 other transcription regulators. Our data support a conclusion that developmental transcriptome revealed significant functional transition during postnatal liver development and FXR plays an important role in control of postnatal liver maturation.

In Vivo Hepatic Reprogramming of Myofibroblasts with AAV Vectors as a Therapeutic Strategy for Liver Fibrosis

Liver fibrosis, a form of scarring, develops in chronic liver diseases when hepatocyte regeneration cannot compensate for hepatocyte death. Initially, collagen produced by myofibroblasts (MFs) functions to maintain the integrity of the liver, but excessive collagen accumulation suppresses residual hepatocyte function, leading to liver failure. As a strategy to generate new hepatocytes and limit collagen deposition in the chronically injured liver, we developed in vivo reprogramming of MFs into hepatocytes using adeno-associated virus (AAV) vectors expressing hepatic transcription factors. We first identified the AAV6 capsid as effective in transducing MFs in a mouse model of liver fibrosis. We then showed in lineage-tracing mice that AAV6 vector-mediated in vivo hepatic reprogramming of MFs generates hepatocytes that replicate function and proliferation of primary hepatocytes, and reduces liver fibrosis. Because AAV vectors are already used for liver-directed human gene therapy, our strategy has potential for clinical translation into a therapy for liver fibrosis.

Effectiveness of the liver micronucleus assay using juvenile mice

This study investigated the effectiveness of the liver micronucleus (MN) assay using juvenile mice. Therefore, we analyzed various hepatic cytochrome P450 (CYP)- mediated activities of ethoxyresorufin O-deethylation, pentoxyresorufin O-dealkylation, tolbutamide hydroxylation, bufuralol 1’-hydroxylation, aniline hydroxylation and midazolam 4-hydroxylation by CYP1A, CYP2B, CYP2C, CYP2D, CYP2E and CYP3A, respectively, in non-treated male ICR mice aged between 3 and 8 weeks. The enzyme efficiency levels in 3- and 4-week-old mice were approximately similar to or higher than those in 8-week-old mice, except for CYP1A and CYP2E in 3- and 4-week-old mice, respectively. Since these results suggest that juvenile mice have sufficient activities for most CYP enzymes, we also conducted a liver MN assay using diethylnitrosamine (DEN), a rodent hepatocarcinogen, on male ICR mice aged between 3 and 6 weeks. A peripheral blood (PB) MN assay was performed simultaneously in 4-week-old mice. Assays incorporating DEN produced positive results in 3- and 4-week-old mice and showed a dose-dependent increase in the micronucleated hepatocyte frequencies at 4 weeks. Both the liver MN assay in 5- and 6-week-old mice and the PB MN assay had negative results when using DEN. These results suggest that 3- and 4-week-old mice have micronuclei-inducing potential in the liver to detect genotoxic compounds using the liver MN assay.

Inhibition of the all-trans Retinoic Acid (atRA) Hydroxylases CYP26A1 and CYP26B1 Results in Dynamic, Tissue-Specific Changes in Endogenous atRA Signaling

All-trans retinoic acid (atRA), the active metabolite of vitamin A, is a ligand for several nuclear receptors and acts as a critical regulator of many physiologic processes. The cytochrome P450 family 26 (CYP26) enzymes are responsible for atRA clearance, and are potential drug targets to increase concentrations of endogenous atRA in a tissue-specific manner. Talarozole is a potent inhibitor of CYP26A1 and CYP26B1, and has shown some success in clinical trials. However, it is not known what magnitude of change is needed in tissue atRA concentrations to promote atRA signaling changes. The aim of this study was to quantify the increase in endogenous atRA concentrations necessary to alter atRA signaling in target organs, and to establish the relationship between CYP26 inhibition and altered atRA concentrations in tissues. Following a single 2.5-mg/kg dose of talarozole to mice, atRA concentrations increased up to 5.7-, 2.7-, and 2.5-fold in serum, liver, and testis, respectively, resulting in induction of Cyp26a1 in the liver and testis and Rar β and Pgc 1β in liver. The increase in atRA concentrations was well predicted from talarozole pharmacokinetics and in vitro data of CYP26 inhibition. After multiple doses of talarozole, a significant increase in atRA concentrations was observed in serum but not in liver or testis. This lack of increase in atRA concentrations correlated with an increase in CYP26A1 expression in the liver. The increased atRA concentrations in serum without a change in liver suggest that CYP26B1 in extrahepatic sites plays a key role in regulating systemic atRA exposure.

Compensatory changes in CYP expression in three different toxicology mouse models: CAR-null, Cyp3a-null, and Cyp2b9/10/13-null mice

Targeted mutant models are common in mechanistic toxicology experiments investigating the absorption, metabolism, distribution, or elimination (ADME) of chemicals from individuals. Key models include those for xenosensing transcription factors and cytochrome P450s (CYP). Here we investigated changes in transcript levels, protein expression, and steroid hydroxylation of several xenobiotic detoxifying CYPs in constitutive androstane receptor (CAR)-null and two CYP-null mouse models that have subfamily members regulated by CAR; the Cyp3a-null and a newly described Cyp2b9/10/13-null mouse model. Compensatory changes in CYP expression that occur in these models may also occur in polymorphic humans, or may complicate interpretation of ADME studies performed using these models. The loss of CAR causes significant changes in several CYPs probably due to loss of CAR-mediated constitutive regulation of these CYPs. Expression and activity changes include significant repression of Cyp2a and Cyp2b members with corresponding drops in 6α- and 16β-testosterone hydroxylase activity. Further, the ratio of 6α-/15α-hydroxylase activity, a biomarker of sexual dimorphism in the liver, indicates masculinization of female CAR-null mice, suggesting a role for CAR in the regulation of sexually dimorphic liver CYP profiles. The loss of Cyp3a causes fewer changes than CAR. Nevertheless, there are compensatory changes including gender-specific increases in Cyp2a and Cyp2b. Cyp2a and Cyp2b were down-regulated in CAR-null mice, suggesting activation of CAR and potentially PXR following loss of the Cyp3a members. However, the loss of Cyp2b causes few changes in hepatic CYP transcript levels and almost no significant compensatory changes in protein expression or activity with the possible exception of 6α-hydroxylase activity. This lack of a compensatory response in the Cyp2b9/10/13-null mice is probably due to low CYP2B hepatic expression, especially in male mice. Overall, compensatory and regulatory CYP changes followed the order CAR-null > Cyp3a-null > Cyp2b-null mice.

RNA Sequencing Reveals Age and Species Differences of Constitutive Androstane Receptor-Targeted Drug-Processing Genes in the Liver

The constitutive androstane receptor (CAR/Nr1i3) is an important xenobiotic-sensing nuclear receptor that is highly expressed in the liver and is well known to have species differences. During development, age-specific activation of CAR may lead to modified pharmacokinetics and toxicokinetics of drugs and environmental chemicals, leading to higher risks for adverse drug reactions in newborns and children. The goal of this study was to systematically investigate the age- and species-specific regulation of various drug-processing genes (DPGs) after neonatal or adult CAR activation in the livers of wild-type, CAR-null, and humanized CAR transgenic mice. At either 5 or 60 days of age, the three genotypes of mice were administered a species-appropriate CAR ligand or vehicle once daily for 4 days (i.p.). The majority of DPGs were differentially regulated by age and/or CAR activation. Thirty-six DPGs were commonly upregulated by CAR activation regardless of age or species of CAR. Although the cumulative mRNAs of uptake transporters were not readily altered by CAR, the cumulative phase I and phase II enzymes as well as efflux transporters were all increased after CAR activation in both species. In general, mouse CAR activation produced comparable or even greater fold increases of many DPGs in newborns than in adults; conversely, humanized CAR activation produced weaker induction in newborns than in adults. Western blotting and enzyme activity assays confirmed the age and species specificities of selected CAR-targeted DPGs. In conclusion, this study systematically compared the effect of age and species of CAR proteins on the regulation of DPGs in the liver and demonstrated that the regulation of xenobiotic biotransformation by CAR is profoundly modified by age and species.

Alternative RNA Splicing in the Pathogenesis of Liver Disease

Non-alcoholic fatty liver disease (NAFLD) is becoming increasingly prevalent due to the worldwide obesity epidemic and currently affects one-third of adults or about one billion people worldwide. NAFLD is predicted to affect over 50% of the world's population by the end of the next decade. It is the most common form of liver disease and is associated with increased risk for progression to a more severe form non-alcoholic steatohepatitis, as well as insulin resistance, type 2 diabetes mellitus, cirrhosis, and eventually hepatocellular carcinoma. This review article will focus on the role of alternative splicing in normal liver physiology and dysregulation in liver disease.

Role of farnesoid X receptor in establishment of ontogeny of phase-I drug metabolizing enzyme genes in mouse liver

The expression of phase-I drug metabolizing enzymes in liver changes dramatically during postnatal liver maturation. Farnesoid X receptor (FXR) is critical for bile acid and lipid homeostasis in liver. However, the role of FXR in regulating ontogeny of phase-I drug metabolizing genes is not clear. Hence, we applied RNA-sequencing to quantify the developmental expression of phase-I genes in both Fxr-null and control (C57BL/6) mouse livers during development. Liver samples of male C57BL/6 and Fxr-null mice at 6 different ages from prenatal to adult were used. The Fxr-null showed an overall effect to diminish the “day-1 surge” of phase-I gene expression, including cytochrome P450s at neonatal ages. Among the 185 phase-I genes from 12 different families, 136 were expressed, and differential expression during development occurred in genes from all 12 phase-I families, including hydrolysis: carboxylesterase (Ces), paraoxonase (Pon), and epoxide hydrolase (Ephx); reduction: aldoketo reductase (Akr), quinone oxidoreductase (Nqo), and dihydropyrimidine dehydrogenase (Dpyd); and oxidation: alcohol dehydrogenase (Adh), aldehyde dehydrogenase (Aldh), flavin monooxygenases (Fmo), molybdenum hydroxylase (Aox and Xdh), cytochrome P450 (P450), and cytochrome P450 oxidoreductase (Por). The data also suggested new phase-I genes potentially targeted by FXR. These results revealed an important role of FXR in regulation of ontogeny of phase-I genes.

Editor's Highlight: Neonatal Activation of the Xenobiotic-Sensors PXR and CAR Results in Acute and Persistent Down-regulation of PPARα-Signaling in Mouse Liver

Safety concerns have emerged regarding the potential long-lasting effects due to developmental exposure to xenobiotics. The pregnane X receptor (PXR) and constitutive androstane receptor (CAR) are critical xenobiotic-sensing nuclear receptors that are highly expressed in liver. The goal of this study was to test our hypothesis that neonatal exposure to PXR- or CAR-activators not only acutely but also persistently regulates the expression of drug-processing genes (DPGs). A single dose of the PXR-ligand PCN (75 mg/kg), CAR-ligand TCPOBOP (3 mg/kg), or vehicle (corn oil) was administered intraperitoneally to 3-day-old neonatal wild-type mice. Livers were collected 24 h post-dose or from adult mice at 60 days of age, and global gene expression of these mice was determined using Affymetrix Mouse Transcriptome Assay 1.0. In neonatal liver, PCN up-regulated 464 and down-regulated 449 genes, whereas TCPOBOP up-regulated 308 and down-regulated 112 genes. In adult liver, there were 15 persistently up-regulated and 22 persistently down-regulated genes following neonatal exposure to PCN, as well as 130 persistently up-regulated and 18 persistently down-regulated genes following neonatal exposure to TCPOBOP. Neonatal exposure to both PCN and TCPOBOP persistently down-regulated multiple Cyp4a members, which are prototypical-target genes of the lipid-sensor PPARα, and this correlated with decreased PPARα-binding to the Cyp4a gene loci. RT-qPCR, western blotting, and enzyme activity assays in livers of wild-type, PXR-null, and CAR-null mice confirmed that the persistent down-regulation of Cyp4a was PXR and CAR dependent. In conclusion, neonatal exposure to PXR- and CAR-activators both acutely and persistently regulates critical genes involved in xenobiotic and lipid metabolism in liver.

RNA-Seq reveals common and unique PXR- and CAR-target gene signatures in the mouse liver transcriptome

The pregnane X receptor (PXR) and constitutive androstane receptor (CAR) are well-known xenobiotic-sensing nuclear receptors with overlapping functions. However, there lacks a quantitative characterization to distinguish between the PXR and CAR target genes and signaling pathways in the liver. The present study performed a transcriptomic comparison of the PXR- and CAR-targets using RNA-Seq in livers of adult wild-type mice that were treated with the prototypical PXR ligand PCN (200mg/kg, i.p. once daily for 4days in corn oil) or the prototypical CAR ligand TCPOBOP (3mg/kg, i.p., once daily for 4days in corn oil). At the given doses, TCPOBOP differentially regulated many more genes (2125) than PCN (212), and 147 of the same genes were differentially regulated by both chemicals. As expected, the top pathways differentially regulated by both PCN and TCPOBOP were involved in xenobiotic metabolism, and they also up-regulated genes involved in retinoid metabolism, but down-regulated genes involved in inflammation and iron homeostasis. Regarding unique pathways, PXR activation appeared to overlap with the aryl hydrocarbon receptor signaling, whereas CAR activation appeared to overlap with the farnesoid X receptor signaling, acute-phase response, and mitochondrial dysfunction. The mRNAs of differentially regulated drug-processing genes (DPGs) partitioned into three patterns, namely TCPOBOP-induced, PCN-induced, as well as TCPOBOP-suppressed gene clusters. The cumulative mRNAs of the differentially regulated DPGs, phase-I and -II enzymes, as well as efflux transporters were all up-regulated by both PCN and TCPOBOPOP, whereas the cumulative mRNAs of the uptake transporters were down-regulated only by TCPOBOP. The absolute mRNA abundance in control and receptor-activated conditions was examined in each DPG category to predict the contribution of specific DPG genes in the PXR/CAR-mediated pharmacokinetic responses. The preferable differential regulation by TCPOBOP in the entire hepatic transcriptome correlated with a marked change in the expression of many DNA and histone epigenetic modifiers. In conclusion, the present study has revealed known and novel, as well as common and unique targets of PXR and CAR in mouse liver following pharmacological activation using their prototypical ligands. Results from this study will further support the role of these receptors in regulating the homeostasis of xenobiotic and intermediary metabolism in the liver, and aid in distinguishing between PXR and CAR signaling at various physiological and pathophysiological conditions. This article is part of a Special Issue entitled: Xenobiotic nuclear receptors: New Tricks for An Old Dog, edited by Dr. Wen Xie.

RNA Sequencing Quantification of Xenobiotic-Processing Genes in Various Sections of the Intestine in Comparison to the Liver of Male Mice

Previous reports on tissue distribution of xenobiotic-processing genes (XPGs) have limitations, because many non-cytochrome P450 phase I enzymes have not been investigated, and one cannot compare the real mRNA abundance of multiple XPGs using conventional quantification methods. Therefore, this study aimed to quantify and compare the mRNA abundance of all major XPGs in the liver and intestine using RNA sequencing. The mRNA profiles of 304 XPGs, including phase I, phase II enzymes, phase II cosubstrate synthetic enzymes, xenobiotic transporters, as well as xenobiotic-related transcription factors, were systematically examined in the liver and various sections of the intestine in adult male C57BL/6J mice. By two-way hierarchical clustering, over 80% of the XPGs had tissue-divergent expression, which partitioned into liver-predominant, small intestine-predominant, and large intestine-predominant patterns. Among the genes, 54% were expressed highest in the liver, 21% in the duodenum, 4% in the jejunum, 6% in the ileum, and 15% in the large intestine. The highest-expressed XPG in the liver was Mgst1; in the duodenum, Cyp3a11; in the jejunum and ileum, Ces2e; and in the large intestine, Cyp2c55. Interestingly, XPGs in the same family usually exhibited highly different tissue distribution patterns, and many XPGs were almost exclusively expressed in one tissue and minimally expressed in others. In conclusion, the present study is among the first and the most comprehensive investigations of the real mRNA abundance and tissue-divergent expression of all major XPGs in mouse liver and intestine, which aids in understanding the tissue-specific biotransformation and toxicity of drugs and other xenobiotics.

Interindividual Variability in Cytochrome P450-Mediated Drug Metabolism

The cytochrome P450 (P450) enzymes are the predominant enzyme system involved in human drug metabolism. Alterations in the expression and/or activity of these enzymes result in changes in pharmacokinetics (and consequently the pharmacodynamics) of drugs that are metabolized by this set of enzymes. Apart from changes in activity as a result of drug-drug interactions (by P450 induction or inhibition), the P450 enzymes can exhibit substantial interindividual variation in basal expression and/or activity, leading to differences in the rates of drug elimination and response. This interindividual variation can result from a myriad of factors, including genetic variation in the promoter or coding regions, variation in transcriptional regulators, alterations in microRNA that affect P450 expression, and ontogenic changes due to exposure to xenobiotics during the developmental and early postnatal periods. Other than administering a probe drug or cocktail of drugs to obtain the phenotype or conducting a genetic analysis to determine genotype, methods to determine interindividual variation are limited. Phenotyping via a probe drug requires exposure to a xenobiotic, and genotyping is not always well correlated with phenotype, making both methodologies less than ideal. This article describes recent work evaluating the effect of some of these factors on interindividual variation in human P450-mediated metabolism and the potential utility of endogenous probe compounds to assess rates of drug metabolism among individuals.

Age-Specific Regulation of Drug-Processing Genes in Mouse Liver by Ligands of Xenobiotic-Sensing Transcription Factors

The xenobiotic-sensing transcription factors (xeno-sensors) AhR, CAR, and PXR upregulate the expression of many drug-processing genes (DPGs) in liver. Previous studies have unveiled profound changes in the basal expression of DPGs during development; however, knowledge on the ontogeny of the inducibility of DPGs in response to pharmacological activation of xeno-sensors is still limited. The goal of this study was to investigate the age-specific regulation of DPGs by prototypical xeno-sensor ligands: 2,3,7,8-tetrachlorodibenzodioxin (TCDD) for AhR; 1,4-bis [2-(3,5-dichloropyridyloxy)] benzene (TCPOBOP) for CAR; and pregnane-16α-carbonitrile (PCN) for PXR during mouse liver development. The basal mRNAs of most DPGs were low during neonatal age, but gradually increased to adult levels, whereas some DPGs (Cyp1a2, Cyp2b10, Cyp3a11, Gstm2, Gstm3, Papss2, and Oatp1a4) exhibited an adolescent-predominant expression pattern. The inducibility of DPGs was age-specific: 1) during neonatal age, the highest fold increase in the mRNA expression was observed for Cyp1a2, Sult5a1, and Ugt1a9 by TCDD; Cyp3a11 and Mrp2 by TCPOBOP; as well as Gstm2 and Gstm3 by PCN; 2) during adolescent age, the highest fold increase in the mRNA expression was observed for Ugt1a6 and Mrp4 by TCDD, Cyp2b10, Ugt2b34, and Ugt2b35 by TCPOBOP, as well as Gsta1, Gsta4, Sult1e1, Ugt1a1, Mrp3, and Mrp4 by PCN; 3) in adults, the highest fold increase in the mRNA expression was observed for Aldh1a1, Aldh1a7, and Ugt2b36 by TCPOBOP, as well as Papss2 and Oatp1a4 by PCN. In conclusion, the inducibility of hepatic DPGs following the pharmacological activation of xeno-sensors is age specific.

Deciphering the Developmental Dynamics of the Mouse Liver Transcriptome

During development, liver undergoes a rapid transition from a hematopoietic organ to a major organ for drug metabolism and nutrient homeostasis. However, little is known on a transcriptome level of the genes and RNA-splicing variants that are differentially regulated with age, and which up-stream regulators orchestrate age-specific biological functions in liver. We used RNA-Seq to interrogate the developmental dynamics of the liver transcriptome in mice at 12 ages from late embryonic stage (2-days before birth) to maturity (60-days after birth). Among 21,889 unique NCBI RefSeq-annotated genes, 9,641 were significantly expressed in at least one age, 7,289 were differently regulated with age, and 859 had multiple (> = 2) RNA splicing-variants. Factor analysis showed that the dynamics of hepatic genes fall into six distinct groups based on their temporal expression. The average expression of cytokines, ion channels, kinases, phosphatases, transcription regulators and translation regulators decreased with age, whereas the average expression of peptidases, enzymes and transmembrane receptors increased with age. The average expression of growth factors peak between Day-3 and Day-10, and decrease thereafter. We identified critical biological functions, upstream regulators, and putative transcription modules that seem to govern age-specific gene expression. We also observed differential ontogenic expression of known splicing variants of certain genes, and 1,455 novel splicing isoform candidates. In conclusion, the hepatic ontogeny of the transcriptome ontogeny has unveiled critical networks and up-stream regulators that orchestrate age-specific biological functions in liver, and suggest that age contributes to the complexity of the alternative splicing landscape of the hepatic transcriptome.

Dose of Phenobarbital and Age of Treatment at Early Life are Two Key Factors for the Persistent Induction of Cytochrome P450 Enzymes in Adult Mouse Liver

Drug treatment of neonates and infants and its long-term consequences on drug responses have emerged in recent years as a major challenge for health care professionals. In the current study, we use phenobarbital as a model drug and mouse as an in vivo model to demonstrate that the dose of phenobarbital and age of treatment are two key factors for the persistent induction of gene expression and consequential increases of enzyme activities of Cyp2b, Cyp2c, and Cyp3a in adult livers. We show that phenobarbital treatment at early life of day 5 after birth with a low dose (<100 mg/kg) does not change expression and enzyme activities of Cyp2b, Cyp2c, and Cyp3a in adult mouse liver, whereas phenobarbital treatment with a high dose (>200 mg/kg) significantly increases expression and enzyme activities of these P450s in adult liver. We also demonstrate that phenobarbital treatment before day 10 after birth, but not at later ages, significantly increases mRNAs, proteins, and enzyme activities of the tested P450s. Such persistent induction of P450 gene expression and enzyme activities in adult livers by phenobarbital treatment only occurs within a sensitive age window early in life. The persistent induction in gene expression and enzyme activities is higher in female mice than in male mice for Cyp2b10 but not for Cyp2c29 and Cyp3a11. These results will stimulate studies to evaluate the long-term impacts of drug treatment with different doses at neonatal and infant ages on drug metabolism, therapeutic efficacy, and drug-induced toxicity throughout the rest of life.

Moderate activation of IKK2-NF-kB in unstressed adult mouse liver induces cytoprotective genes and lipogenesis without apparent signs of inflammation or fibrosis

Background

The NF-kB signaling, regulated by IKK1-p52/RelB and IKK2-p65, is activated by various stresses to protect or damage the liver, in context-specific manners. Two previous studies of liver-specific expression of constitutive active IKK2 (IKK2ca) showed that strong activation of IKK2-NF-kB in mouse livers caused inflammation, insulin resistance, and/or fibrosis. The purpose of this study was to understand how moderate activation of IKK2-NF-kB in adult mouse livers alters hepatic gene expression and pathophysiology.

Method

We generated mice with adult hepatocyte-specific activation of Ikk2 (Liv-Ikk2ca) using Alb-cre mice and Ikk2ca Rosa26 knockin mice in which a moderate expression of Ikk2ca transgene was driven by the endogenous Rosa26 promoter.

Results

Surprisingly, compared to wild-type mice, adult male Liv-Ikk2ca mice had higher hepatic mRNA expression of Ikk2 and classical NF-kB targets (e.g. Lcn2 and A20), as well as IKK1, NIK, and RelB, but no changes in markers of inflammation or fibrosis. Blood levels of IL-6 and MCP-1 remained unchanged, and histology analysis showed a lack of injury or infiltration of inflammatory cells in livers of Liv-Ikk2ca mice. Moreover, Liv-Ikk2ca mice had lower mRNA expression of prooxidative enzymes Cyp2e1 and Cyp4a14, higher expression of antioxidative enzymes Sod2, Gpx1, and Nqo1, without changes in key enzymes for fatty acid oxidation, glucose utilization, or gluconeogenesis. In parallel, Liv-Ikk2ca mice and wild-type mice had similar levels of hepatic reduced glutathione, endogenous reactive oxygen species, and lipid peroxidation. Additionally, Liv-Ikk2ca mice had higher Cyp3a11 without down-regulation of most drug processing genes. Regarding nuclear proteins of NF-kB subunits, Liv-Ikk2ca mice had moderately higher p65 and p50 but much higher RelB. Results of ChIP-qPCR showed that the binding of p50 to multiple NF-kB-target genes was markedly increased in Liv-Ikk2ca mice. Additionally, Liv-Ikk2ca mice had moderate increase in triglycerides in liver, which was associated with higher lipogenic factors Pparγ, Lxr, Fasn, Scd1, and CD36.

Conclusion

In summary, moderate activation of IKK2-NF-kB in unstressed adult mouse hepatocytes produces a cytoprotective gene expression profile and induces lipogenesis without apparent signs of inflammation or fibrosis, likely due to strong activation of the anti-inflammatory IKK1-RelB alternative NF-kB pathway as well as the Lxr.

Electronic supplementary material

The online version of this article (doi:10.1186/s12876-015-0325-z) contains supplementary material, which is available to authorized users.

Developmental Regulation of Drug-Processing Genes in Livers of Germ-Free Mice

Very little is known about the effect of gut microbiota on the ontogeny of drug-processing genes (DPGs) in liver. In this study, livers were harvested from conventional (CV) and germ-free (GF) male and female mice from 1 to 90 days of age. RNA-Seq in livers of 90-day-old male mice showed that xenobiotic metabolism was the most downregulated pathway within the mRNA transcriptome in absence of intestinal bacteria. In male livers, the mRNAs of 67 critical DPGs partitioned into 4 developmental patterns (real-time-quantitative polymerase chain reaction): Pattern-1 gradually increased to adult levels in livers of CV mice and were downregulated in livers of GF mice, as exemplified by the major drug-metabolizing enzymes cytochrome 3a (Cyp3a) family, which are prototypical pregnane X receptor (PXR)-target genes. Genes in Pattern-2 include Cyp1a2 (aryl hydrocarbon receptor-target gene), Cyp2c family, and Cyp2e1, which were all upregulated mainly at 90 days of age; as well as the peroxisome proliferator-activated receptor α (PPARα)-target genes Cyp4a family and Aldh3a2, which were upregulated not only in 90-days adult age, but also between neonatal and adolescent ages (from 1 to 30 days of age). Genes in Pattern-3 were enriched predominantly in livers of 15-day-old mice, among which the sterol-efflux transporter dimers Abcg5/Abcg8 were downregulated in GF mice. Genes in Pattern-4 were neonatal-enriched, among which the transporter Octn1 mRNA tended to be lower in GF mice at younger ages but higher in adult GF mice as compared with age-matched CV mice. Protein assays confirmed the downregulation of the PXR-target gene Cyp3a protein (Western-blot and liquid chromatography tandem mass spectroscopy), and decreased Cyp3a enzyme activities in male GF livers. Increased microsomal-Cyp4a proteins and nuclear-PPARα were also observed in male GF livers. Interestingly, in contrast to male livers, the mRNAs of Cyp2c or Cyp4a were not readily upregulated in female GF livers approaching adult age, suggesting the maturation of female-specific hormones interferes with the interactions between intestinal microbiota and DPG ontogeny. In conclusion, intestinal microbiota markedly impacts the ontogeny of many hepatic DPGs in a gender-specific manner.

Establishment of metabolism and transport pathways in the rodent and human fetal liver

The ultimate fate of drugs and chemicals in the body is largely regulated by hepatic uptake, metabolism, and excretion. The liver acquires the functional ability to metabolize and transport chemicals during the perinatal period of development. Research using livers from fetal and juvenile rodents and humans has begun to reveal the timing, key enzymes and transporters, and regulatory factors that are responsible for the establishment of hepatic phase I and II metabolism as well as transport. The majority of this research has been limited to relative mRNA and protein quantification. However, the recent utilization of novel technology, such as RNA-Sequencing, and the improved availability and refinement of functional activity assays, has begun to provide more definitive information regarding the extent of hepatic drug disposition in the developing fetus. The goals of this review are to provide an overview of the early regulation of the major phase I and II enzymes and transporters in rodent and human livers and to highlight potential mechanisms that control the ontogeny of chemical metabolism and excretion pathways.

RNA-sequencing quantification of hepatic ontogeny of phase-I enzymes in mice

Phase-I drug metabolizing enzymes catalyze reactions of hydrolysis, reduction, and oxidation of drugs and play a critical role in drug metabolism. However, the functions of most phase-I enzymes are not mature at birth, which markedly affects drug metabolism in newborns. Therefore, characterization of the expression profiles of phase-I enzymes and the underlying regulatory mechanisms during liver maturation is needed for better estimation of using drugs in pediatric patients. The mouse is an animal model widely used for studying the mechanisms in the regulation of developmental expression of phase-I genes. Therefore, we applied RNA sequencing to provide a "true quantification" of the mRNA expression of phase-I genes in the mouse liver during development. Liver samples of male C57BL/6 mice at 12 different ages from prenatal to adulthood were used for defining the ontogenic mRNA profiles of phase-I families, including hydrolysis: carboxylesterase (Ces), paraoxonase (Pon), and epoxide hydrolase (Ephx); reduction: aldo-keto reductase (Akr), quinone oxidoreductase (Nqo), and dihydropyrimidine dehydrogenase (Dpyd); and oxidation: alcohol dehydrogenase (Adh), aldehyde dehydrogenase (Aldh), flavin monooxygenases (Fmo), molybdenum hydroxylase (Aox and Xdh), cytochrome P450 (P450), and cytochrome P450 oxidoreductase (Por). Two rapidly increasing stages of total phase-I gene expression after birth reflect functional transition of the liver during development. Diverse expression patterns were identified, and some large gene families contained the mRNA of genes that are enriched at different stages of development. Our study reveals the mRNA abundance of phase-I genes in the mouse liver during development and provides a valuable foundation for mechanistic studies in the future.

Altered subcellular localization of heat shock protein 90 is associated with impaired expression of the aryl hydrocarbon receptor pathway in dogs

The aryl hydrocarbon receptor (AHR) mediates biological responses to toxic chemicals. An unexpected role for AHR in vascularization was suggested when mice lacking AHR displayed impaired closure of the ductus venosus after birth, as did knockout mice for aryl hydrocarbon receptor interacting protein (AIP) and aryl hydrocarbon receptor nuclear translocator (ARNT). The resulting intrahepatic portosystemic shunts (IHPSS) are frequently diagnosed in specific dog breeds, such as the Irish wolfhound. We compared the expression of components of the AHR pathway in healthy Irish wolfhounds and dogs with IHPSS. To this end, we analyzed the mRNA expression in the liver of AHR,AIP, ARNT, and other genes involved in this pathway, namely, those for aryl hydrocarbon receptor nuclear translocator 2 (ARNT2), hypoxia inducible factor 1alpha (HIF1A), heat shock protein 90AA1 (HSP90AA1), cytochromes P450 (CYP1A1, CYP1A2, and CYP1B1), vascular endothelial growth factor A (VEGFA), nitric oxide synthesase 3 (NOS3), and endothelin (EDN1). The observed low expression of AHR mRNA in the Irish wolfhounds is in associated with a LINE-1 insertion in intron 2, for which these dogs were homozygous. Down regulation in Irish wolfhounds was observed for AIP, ARNT2, CYP1A2, CYP1B1 and HSP90AA1 expression, whereas the expression of HIF1A was increased. Immunohistochemistry revealed lower levels of AHR, HIF1A, and VEGFA protein in the nucleus and lower levels of ARNT and HSP90AA1 protein in the cytoplasm of the liver cells of Irish wolfhounds. The impaired expression of HSP90AA1 could trigger the observed differences in mRNA and protein levels and therefore explain the link between two very different functions of AHR: regulation of the closure of the ductus venosus and the response to toxins.

RNA-sequencing quantification of hepatic ontogeny and tissue distribution of mRNAs of phase II enzymes in mice

Phase II conjugating enzymes play key roles in the metabolism of xenobiotics. In the present study, RNA sequencing was used to elucidate hepatic ontogeny and tissue distribution of mRNA expression of all major known Phase II enzymes, including enzymes involved in glucuronidation, sulfation, glutathione conjugation, acetylation, methylation, and amino acid conjugation, as well as enzymes for the synthesis of Phase II cosubstrates, in male C57BL/6J mice. Livers from male C57BL/6J mice were collected at 12 ages from prenatal to adulthood. Many of these Phase II enzymes were expressed at much higher levels in adult livers than in perinatal livers, such as Ugt1a6b, -2a3, -2b1, -2b5, -2b36, -3a1, and -3a2; Gsta1, -m1, -p1, -p2, and -z1; mGst1; Nat8; Comt; Nnmt; Baat; Ugdh; and Gclc. In contrast, hepatic mRNA expression of a few Phase II enzymes decreased during postnatal liver development, such as mGst2, mGst3, Gclm, and Mat2a. Hepatic expression of certain Phase II enzymes peaked during the adolescent stage, such as Ugt1a1, Sult1a1, Sult1c2, Sult1d1, Sult2as, Sult5a1, Tpmt, Glyat, Ugp2, and Mat1a. In adult mice, the total transcripts for Phase II enzymes were comparable in liver, kidney, and small intestine; however, individual Phase II enzymes displayed marked tissue specificity among the three organs. In conclusion, this study unveils for the first time developmental changes in mRNA abundance of all major known Phase II enzymes in mouse liver, as well as their tissue-specific expression in key drug-metabolizing organs. The age- and tissue-specific expression of Phase II enzymes indicate that the detoxification of xenobiotics is highly regulated by age and cell type.

Hepatic ontogeny and tissue distribution of mRNAs of epigenetic modifiers in mice using RNA-sequencing

Developmental regulation of gene expression is controlled by distinct epigenetic signatures catalyzed by various epigenetic modifiers. Little is known about the ontogeny and tissue distribution of these epigenetic modifiers. In the present study, we used a novel approach of RNA-sequencing to elucidate hepatic ontogeny and tissue distribution of mRNA expression of 142 epigenetic modifiers, including enzymes involved in DNA methylation/demethylation, histone acetylation/deacetylation, histone methylation/demethylation, histone phosphorylation and chromosome remodeling factors in male C57BL/6 mice. Livers from male C57BL/6 mice were collected at 12 ages from prenatal to adulthood. Many of these epigenetic modifiers were expressed at much higher levels in perinatal livers than adult livers, such as Dnmt1, Dnmt3a, Dnmt3b, Apobec3, Kat1, Ncoa4, Setd8, Ash2l, Dot1l, Cbx1, Cbx3, Cbx5, Cbx6, Ezh2, Suz12, Eed, Suv39h1, Suv420h2, Dek, Hdac1, Hdac2, Hdac7, Kdm2b, Kdm5c, Kdm7, Prmt1-5, Prmt7, Smarca4, Smarcb1, Chd4 and Ino80e. In contrast, hepatic mRNA expression of a few epigenetic modifiers increased during postnatal liver development, such as Smarca2, Kdm1b, Cbx7 and Chd3. In adult mice (60 d of age), most epigenetic modifiers were expressed at moderately (1-3-fold) higher levels in kidney and/or small intestine than liver. In conclusion, this study, for the first time, unveils developmental changes in mRNA abundance of all major known epigenetic modifiers in mouse liver. These data suggest that ontogenic changes in mRNA expression of epigenetic modifiers may play important roles in determining the addition and/or removal of corresponding epigenetic signatures during liver development.