GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis.

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

GSEApy: a comprehensive package for performing gene set enrichment analysis in Python.

MOTIVATION: Gene set enrichment analysis (GSEA) is a commonly used algorithm for characterizing gene expression changes. However, the currently available tools used to perform GSEA have a limited ability to analyze large datasets, which is particularly problematic for the analysis of single-cell data. To overcome this limitation, we developed a GSEA package in Python (GSEApy), which could efficiently analyze large single-cell datasets. RESULTS: We present a package (GSEApy) that performs GSEA in either the command line or Python environment. GSEApy uses a Rust implementation to enable it to calculate the same enrichment statistic as GSEA for a collection of pathways. The Rust implementation of GSEApy is 3-fold faster than the Numpy version of GSEApy (v0.10.8) and uses >4-fold less memory. GSEApy also provides an interface between Python and Enrichr web services, as well as for BioMart. The Enrichr application programming interface enables GSEApy to perform over-representation analysis for an input gene list. Furthermore, GSEApy consists of several tools, each designed to facilitate a particular type of enrichment analysis. AVAILABILITY AND IMPLEMENTATION: The new GSEApy with Rust extension is deposited in PyPI: https://pypi.org/project/gseapy/. The GSEApy source code is freely available at https://github.com/zqfang/GSEApy. Also, the documentation website is available at https://gseapy.rtfd.io/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hepatic organoids move from adolescence to maturity.

Since organoids were developed 15 years ago, they are now in their adolescence as a research tool. The ability to generate 'tissue in a dish' has created enormous opportunities for biomedical research. We examine the contributions that hepatic organoids have made to three areas of liver research: as a source of cells and tissue for basic research, for drug discovery and drug safety testing, and for understanding disease pathobiology. We discuss the features that enable hepatic organoids to provide useful models for human liver diseases and identify four types of advances that will enable them to become a mature (i.e., adult) research tool over the next 5 years. During this period, advances in single-cell RNA sequencing and CRISPR technologies coupled with improved hepatic organoid methodology, which enables them to have a wider range of cell types that are present in liver and to be grown in microwells, will generate discoveries that will dramatically advance our understanding of liver development and the pathogenesis of liver diseases. It will generate also new approaches for treating liver fibrosis, which remains a major public health problem with few treatment options.

Analysis of structural variation among inbred mouse strains.

BACKGROUND: 'Long read' sequencing methods have been used to identify previously uncharacterized structural variants that cause human genetic diseases. Therefore, we investigated whether long read sequencing could facilitate genetic analysis of murine models for human diseases. RESULTS: The genomes of six inbred strains (BTBR T + Itpr3tf/J, 129Sv1/J, C57BL/6/J, Balb/c/J, A/J, SJL/J) were analyzed using long read sequencing. Our results revealed that (i) Structural variants are very abundant within the genome of inbred strains (4.8 per gene) and (ii) that we cannot accurately infer whether structural variants are present using conventional short read genomic sequence data, even when nearby SNP alleles are known. The advantage of having a more complete map was demonstrated by analyzing the genomic sequence of BTBR mice. Based upon this analysis, knockin mice were generated and used to characterize a BTBR-unique 8-bp deletion within Draxin that contributes to the BTBR neuroanatomic abnormalities, which resemble human autism spectrum disorder. CONCLUSION: A more complete map of the pattern of genetic variation among inbred strains, which is produced by long read genomic sequencing of the genomes of additional inbred strains, could facilitate genetic discovery when murine models of human diseases are analyzed.

An automated multi-modal graph-based pipeline for mouse genetic discovery.

MOTIVATION: Our ability to identify causative genetic factors for mouse genetic models of human diseases and biomedical traits has been limited by the difficulties associated with identifying true causative factors, which are often obscured by the many false positive genetic associations produced by a GWAS. RESULTS: To accelerate the pace of genetic discovery, we developed a graph neural network (GNN)-based automated pipeline (GNNHap) that could rapidly analyze mouse genetic model data and identify high probability causal genetic factors for analyzed traits. After assessing the strength of allelic associations with the strain response pattern; this pipeline analyzes 29M published papers to assess candidate gene-phenotype relationships; and incorporates the information obtained from a protein-protein interaction network and protein sequence features into the analysis. The GNN model produces markedly improved results relative to that of a simple linear neural network. We demonstrate that GNNHap can identify novel causative genetic factors for murine models of diabetes/obesity and for cataract formation, which were validated by the phenotypes appearing in previously analyzed gene knockout mice. The diabetes/obesity results indicate how characterization of the underlying genetic architecture enables new therapies to be discovered and tested by applying 'precision medicine' principles to murine models. AVAILABILITY AND IMPLEMENTATION: The GNNHap source code is freely available at https://github.com/zqfang/gnnhap, and the new version of the HBCGM program is available at https://github.com/zqfang/haplomap. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Dynamic landscape and regulation of RNA editing in mammals.

Adenosine-to-inosine (A-to-I) RNA editing is a conserved post-transcriptional mechanism mediated by ADAR enzymes that diversifies the transcriptome by altering selected nucleotides in RNA molecules. Although many editing sites have recently been discovered, the extent to which most sites are edited and how the editing is regulated in different biological contexts are not fully understood. Here we report dynamic spatiotemporal patterns and new regulators of RNA editing, discovered through an extensive profiling of A-to-I RNA editing in 8,551 human samples (representing 53 body sites from 552 individuals) from the Genotype-Tissue Expression (GTEx) project and in hundreds of other primate and mouse samples. We show that editing levels in non-repetitive coding regions vary more between tissues than editing levels in repetitive regions. Globally, ADAR1 is the primary editor of repetitive sites and ADAR2 is the primary editor of non-repetitive coding sites, whereas the catalytically inactive ADAR3 predominantly acts as an inhibitor of editing. Cross-species analysis of RNA editing in several tissues revealed that species, rather than tissue type, is the primary determinant of editing levels, suggesting stronger cis-directed regulation of RNA editing for most sites, although the small set of conserved coding sites is under stronger trans-regulation. In addition, we curated an extensive set of ADAR1 and ADAR2 targets and showed that many editing sites display distinct tissue-specific regulation by the ADAR enzymes in vivo. Further analysis of the GTEx data revealed several potential regulators of editing, such as AIMP2, which reduces editing in muscles by enhancing the degradation of the ADAR proteins. Collectively, our work provides insights into the complex cis- and trans-regulation of A-to-I editing.

Genetic basis of thermal nociceptive sensitivity and brain weight in a BALB/c reduced complexity cross.

Thermal nociception involves the transmission of temperature-related noxious information from the periphery to the CNS and is a heritable trait that could predict transition to persistent pain. Rodent forward genetics complement human studies by controlling genetic complexity and environmental factors, analysis of end point tissue, and validation of variants on appropriate genetic backgrounds. Reduced complexity crosses between nearly identical inbred substrains with robust trait differences can greatly facilitate unbiased discovery of novel genes and variants. We found BALB/cByJ mice showed enhanced sensitivity on the 53.5°C hot plate and mechanical stimulation in the von Frey test compared to BALB/cJ mice and replicated decreased gross brain weight in BALB/cByJ versus BALB/cJ. We then identified a quantitative trait locus (QTL) on chromosome 13 for hot plate sensitivity (LOD = 10.7; p < 0.001; peak = 56 Mb) and a QTL for brain weight on chromosome 5 (LOD = 8.7; p < 0.001). Expression QTL mapping of brain tissues identified H2afy (56.07 Mb) as the top transcript with the strongest association at the hot plate locus (FDR = 0.0002) and spliceome analysis identified differential exon usage within H2afy associated with the same locus. Whole brain proteomics further supported decreased H2AFY expression could underlie enhanced hot plate sensitivity, and identified ACADS as a candidate for reduced brain weight. To summarize, a BALB/c reduced complexity cross combined with multiple-omics approaches facilitated identification of candidate genes underlying thermal nociception and brain weight. These substrains provide a powerful, reciprocal platform for future validation of candidate variants.

Zhx2 Is a Candidate Gene Underlying Oxymorphone Metabolite Brain Concentration Associated with State-Dependent Oxycodone Reward.

Understanding the pharmacogenomics of opioid metabolism and behavior is vital to therapeutic success, as mutations can dramatically alter therapeutic efficacy and addiction liability. We found robust, sex-dependent BALB/c substrain differences in oxycodone behaviors and whole brain concentration of oxycodone metabolites. BALB/cJ females showed robust state-dependent oxycodone reward learning as measured via conditioned place preference when compared with the closely related BALB/cByJ substrain. Accordingly, BALB/cJ females also showed a robust increase in brain concentration of the inactive metabolite noroxycodone and the active metabolite oxymorphone compared with BALB/cByJ mice. Oxymorphone is a highly potent, full agonist at the mu opioid receptor that could enhance drug-induced interoception and state-dependent oxycodone reward learning. Quantitative trait locus (QTL) mapping in a BALB/c F2 reduced complexity cross revealed one major QTL on chromosome 15 underlying brain oxymorphone concentration that explained 32% of the female variance. BALB/cJ and BALB/cByJ differ by fewer than 10,000 variants, which can greatly facilitate candidate gene/variant identification. Hippocampal and striatal cis-expression QTL (eQTL) and exon-level eQTL analysis identified Zhx2, a candidate gene coding for a transcriptional repressor with a private BALB/cJ retroviral insertion that reduces Zhx2 expression and sex-dependent dysregulation of cytochrome P450 enzymes. Whole brain proteomics corroborated the Zhx2 eQTL and identified upregulated CYP2D11 that could increase brain oxymorphone in BALB/cJ females. To summarize, Zhx2 is a highly promising candidate gene underlying brain oxycodone metabolite levels. Future studies will validate Zhx2 and its site of action using reciprocal gene editing and tissue-specific viral manipulations in BALB/c substrains. SIGNIFICANCE STATEMENT: Our findings show that genetic variation can result in sex-specific alterations in whole brain concentration of a bioactive opioid metabolite after oxycodone administration, reinforcing the need for sex as a biological factor in pharmacogenomic studies. The cooccurrence of female-specific increased oxymorphone and state-dependent reward learning suggests that this minor yet potent and efficacious metabolite of oxycodone could increase opioid interoception and drug-cue associative learning of opioid reward, which has implications for cue-induced relapse of drug-seeking behavior and for precision pharmacogenetics.

The phosphatidylethanolamine biosynthesis pathway provides a new target for cancer chemotherapy.

BACKGROUND & AIMS: Since human induced pluripotent stem cells (iPSCs) develop into hepatic organoids through stages that resemble human embryonic liver development, they can be used to study developmental processes and disease pathology. Therefore, we examined the early stages of hepatic organoid formation to identify key pathways affecting early liver development. METHODS: Single-cell RNA-sequencing and metabolomic analysis was performed on developing organoid cultures at the iPSC, hepatoblast (day 9) and mature organoid stage. The importance of the phosphatidylethanolamine biosynthesis pathway to early liver development was examined in developing organoid cultures using iPSC with a CRISPR-mediated gene knockout and an over the counter medication (meclizine) that inhibits the rate-limiting enzyme in this pathway. Meclizine's effect on the growth of a human hepatocarcinoma cell line in a xenotransplantation model and on the growth of acute myeloid leukemia cells in vitro was also examined. RESULTS: Transcriptomic and metabolomic analysis of organoid development indicated that the phosphatidylethanolamine biosynthesis pathway is essential for early liver development. Unexpectedly, early hepatoblasts were selectively sensitive to the cytotoxic effect of meclizine. We demonstrate that meclizine could be repurposed for use in a new synergistic combination therapy for primary liver cancer: a glycolysis inhibitor reprograms cancer cell metabolism to make it susceptible to the cytotoxic effect of meclizine. This combination inhibited the growth of a human liver carcinoma cell line in vitro and in a xenotransplantation model, without causing significant side effects. This drug combination was also highly active against acute myeloid leukemia cells. CONCLUSION: Our data indicate that phosphatidylethanolamine biosynthesis is a targetable pathway for cancer; meclizine may have clinical efficacy as a repurposed anti-cancer drug when used as part of a new combination therapy. LAY SUMMARY: The early stages of human liver development were modeled using human hepatic organoids. We identified a pathway that was essential for early liver development. Based upon this finding, a novel combination drug therapy was identified that could be used to treat primary liver cancer and possibly other types of cancer.

Can Humanized Mice Predict Drug "Behavior" in Humans?

Most of what we know about a drug prior to human clinical studies is derived from animal testing. Because animals and humans have substantial differences in their physiology and in their drug metabolism pathways, we do not know very much about the pharmacokinetic and pharmacodynamic behavior of a drug in humans until after it is administered to many people. Hence, drug-induced liver injury has become a significant public health problem, and we have a very inefficient drug development process with a high failure rate. Because the human liver is at the heart of these problems, chimeric mice with humanized livers could be used to address these issues. We examine recent evidence indicating that drug testing in chimeric mice could provide better information about a drug's metabolism, disposition, and toxicity (i.e., its "behavior") in humans and could aid in developing personalized medicine strategies, which would improve drug efficacy and safety.

Profiling of ARDS pulmonary edema fluid identifies a metabolically distinct subset.

There is considerable biological and physiological heterogeneity among patients who meet standard clinical criteria for acute respiratory distress syndrome (ARDS). In this study, we tested the hypothesis that there exists a subgroup of ARDS patients who exhibit a metabolically distinct profile. We examined undiluted pulmonary edema fluid obtained at the time of endotracheal intubation from 16 clinically phenotyped ARDS patients and 13 control patients with hydrostatic pulmonary edema. Nontargeted metabolic profiling was carried out on the undiluted edema fluid. Univariate and multivariate statistical analyses including principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) were conducted to find discriminant metabolites. Seven-hundred and sixty unique metabolites were identified in the pulmonary edema fluid of these 29 patients. We found that a subset of ARDS patients (6/16, 38%) presented a distinct metabolic profile with the overrepresentation of 235 metabolites compared with edema fluid from the other 10 ARDS patients, whose edema fluid metabolic profile was indistinguishable from those of the 13 control patients with hydrostatic edema. This "high metabolite" endotype was characterized by higher concentrations of metabolites belonging to all of the main metabolic classes including lipids, amino acids, and carbohydrates. This distinct group with high metabolite levels in the edema fluid was also associated with a higher mortality rate. Thus metabolic profiling of the edema fluid of ARDS patients supports the hypothesis that there is considerable biological heterogeneity among ARDS patients who meet standard clinical and physiological criteria for ARDS.

Palonosetron and hydroxyzine pre-treatment reduces the objective signs of experimentally-induced acute opioid withdrawal in humans: a double-blinded, randomized, placebo-controlled crossover study.

BACKGROUND: Treatments for reducing opioid withdrawal are limited and prone to problematic side effects. Laboratory studies, clinical observations, and limited human trial data suggest 5-HT3-receptor antagonists and antihistamines may be effective. OBJECTIVES: This double-blind, crossover, placebo-controlled study employing an acute physical dependence model evaluated whether (i) treatment with a 5-HT3-receptor antagonist (palonosetron) would reduce opioid withdrawal symptoms, and (ii) co-administration of an antihistamine (hydroxyzine) would enhance any treatment effect. METHODS: At timepoint T = 0, healthy (non-opioid dependent, non-substance abuser) male volunteers (N = 10) were pre-treated with either a) placebo, b) palonosetron IV (0.75 mg), or c) palonosetron IV (0.75 mg) and hydroxyzine PO (100 mg) in a crossover study design. This was followed at T = 30 by intravenous morphine (10 mg/70kg). At T = 165, 10 mg/70kg naloxone IV was given to precipitate opioid withdrawal. The objective opioid withdrawal score (OOWS) and subjective opioid withdrawal score (SOWS) were determined 5 and 15 minutes after naloxone administration (T = 170, 180, respectively). Baseline measurements were recorded at T = -30 and T = -15. RESULTS: Comparison of average baseline OOWS scores with OOWS scores obtained 15 minutes after naloxone was significant (p = 0.0001). Scores from 15 minutes post-naloxone infusion showed significant differences in OOWS scores between treatment groups: placebo, 3.7 ± 2.4; palonosetron, 1.5 ± 0.97; and palonosetron with hydroxyzine, 0.2 ± 0.1333. CONCLUSIONS: Pretreatment with palonosetron significantly reduced many signs of experimentally-induced opioid withdrawal. Co-administration with hydroxyzine further reduced opioid withdrawal severity. These results suggest that 5-HT3 receptor antagonists, alone or in combination with an antihistamine, may be useful in the treatment of opioid withdrawal.

The multiple PDZ domain protein Mpdz/MUPP1 regulates opioid tolerance and opioid-induced hyperalgesia.

BACKGROUND: Opioids are a mainstay for the treatment of chronic pain. Unfortunately, therapy-limiting maladaptations such as loss of treatment effect (tolerance), and paradoxical opioid-induced hyperalgesia (OIH) can occur. The objective of this study was to identify genes responsible for opioid tolerance and OIH. RESULTS: These studies used a well-established model of ascending morphine administration to induce tolerance, OIH and other opioid maladaptations in 23 strains of inbred mice. Genome-wide computational genetic mapping was then applied to the data in combination with a false discovery rate filter. Transgenic mice, gene expression experiments and immunoprecipitation assays were used to confirm the functional roles of the most strongly linked gene. The behavioral data processed using computational genetic mapping and false discovery rate filtering provided several strongly linked biologically plausible gene associations. The strongest of these was the highly polymorphic Mpdz gene coding for the post-synaptic scaffolding protein Mpdz/MUPP1. Heterozygous Mpdz +/- mice displayed reduced opioid tolerance and OIH. Mpdz gene expression and Mpdz/MUPP1 protein levels were lower in the spinal cords of low-adapting 129S1/Svlm mice than in high-adapting C57BL/6 mice. Morphine did not alter Mpdz expression levels. In addition, association of Mpdz/MUPP1 with its known binding partner CaMKII did not differ between these high- and low-adapting strains. CONCLUSIONS: The degrees of maladaptive changes in response to repeated administration of morphine vary greatly across inbred strains of mice. Variants of the multiple PDZ domain gene Mpdz may contribute to the observed inter-strain variability in tolerance and OIH by virtue of changes in the level of their expression.

The role of Abcb5 alleles in susceptibility to haloperidol-induced toxicity in mice and humans.

BACKGROUND: We know very little about the genetic factors affecting susceptibility to drug-induced central nervous system (CNS) toxicities, and this has limited our ability to optimally utilize existing drugs or to develop new drugs for CNS disorders. For example, haloperidol is a potent dopamine antagonist that is used to treat psychotic disorders, but 50% of treated patients develop characteristic extrapyramidal symptoms caused by haloperidol-induced toxicity (HIT), which limits its clinical utility. We do not have any information about the genetic factors affecting this drug-induced toxicity. HIT in humans is directly mirrored in a murine genetic model, where inbred mouse strains are differentially susceptible to HIT. Therefore, we genetically analyzed this murine model and performed a translational human genetic association study. METHODS AND FINDINGS: A whole genome SNP database and computational genetic mapping were used to analyze the murine genetic model of HIT. Guided by the mouse genetic analysis, we demonstrate that genetic variation within an ABC-drug efflux transporter (Abcb5) affected susceptibility to HIT. In situ hybridization results reveal that Abcb5 is expressed in brain capillaries, and by cerebellar Purkinje cells. We also analyzed chromosome substitution strains, imaged haloperidol abundance in brain tissue sections and directly measured haloperidol (and its metabolite) levels in brain, and characterized Abcb5 knockout mice. Our results demonstrate that Abcb5 is part of the blood-brain barrier; it affects susceptibility to HIT by altering the brain concentration of haloperidol. Moreover, a genetic association study in a haloperidol-treated human cohort indicates that human ABCB5 alleles had a time-dependent effect on susceptibility to individual and combined measures of HIT. Abcb5 alleles are pharmacogenetic factors that affect susceptibility to HIT, but it is likely that additional pharmacogenetic susceptibility factors will be discovered. CONCLUSIONS: ABCB5 alleles alter susceptibility to HIT in mouse and humans. This discovery leads to a new model that (at least in part) explains inter-individual differences in susceptibility to a drug-induced CNS toxicity.

A better prognosis for genetic association studies in mice.

Although inbred mouse strains have been the premier model organism used in biomedical research, multiple studies and analyses have indicated that genome-wide association studies (GWAS) cannot be productively performed using inbred mouse strains. However, there is one type of GWAS in mice that has successfully identified the genetic basis for many biomedical traits of interest: haplotype-based computational genetic mapping (HBCGM). Here, we describe how the methodological basis for a HBCGM study significantly differs from that of a conventional murine GWAS, and how an integrative analysis of its output within the context of other 'omic' information can enable genetic discovery. Consideration of these factors will substantially improve the prognosis for the utility of murine genetic association studies for biomedical discovery.

Humanized thymidine kinase-NOG mice can be used to identify drugs that cause animal-specific hepatotoxicity: a case study with furosemide.

Interspecies differences have limited the predictive utility of toxicology studies performed using animal species. A drug that could be a safe and effective treatment in humans could cause toxicity in animals, preventing it from being used in humans. We investigated whether the use of thymidine kinase (TK)-NOG mice with humanized livers could prevent this unfortunate outcome (i.e., "rescue" a drug for use in humans). A high dose of furosemide is known to cause severe liver toxicity in mice, but it is a safe and effective treatment in humans. We demonstrate that administration of a high dose of furosemide (200 mg/kg i.p.) causes extensive hepatotoxicity in control mice but not in humanized TK-NOG mice. This interspecies difference results from a higher rate of production of the toxicity-causing metabolite by mouse liver. Comparison of their survival curves indicated that the humanized mice were more resistant than control mice to the hepatotoxicity caused by high doses of furosemide. In this test case, humanized TK-NOG mouse studies indicate that humans could be safely treated with a high dose of furosemide.

Chimeric TK-NOG mice: a predictive model for cholestatic human liver toxicity.

Due to the substantial interspecies differences in drug metabolism and disposition, drug-induced liver injury (DILI) in humans is often not predicted by studies performed in animal species. For example, a drug (bosentan) used to treat pulmonary artery hypertension caused unexpected cholestatic liver toxicity in humans, which was not predicted by preclinical toxicology studies in multiple animal species. In this study, we demonstrate that NOG mice expressing a thymidine kinase transgene (TK-NOG) with humanized livers have a humanized profile of biliary excretion of a test (cefmetazole) drug, which was shown by an in situ perfusion study to result from interspecies differences in the rate of biliary transport and in liver retention of this drug. We also found that readily detectable cholestatic liver injury develops in TK-NOG mice with humanized livers after 1 week of treatment with bosentan (160, 32, or 6 mg/kg per day by mouth), whereas liver toxicity did not develop in control mice after 1 month of treatment. The laboratory and histologic features of bosentan-induced liver toxicity in humanized mice mirrored that of human subjects. Because DILI has become a significant public health problem, drug safety could be improved if preclinical toxicology studies were performed using humanized TK-NOG.

The Netrin-1 receptor DCC is a regulator of maladaptive responses to chronic morphine administration.

BACKGROUND: Opioids are the cornerstone of treatment for moderate to severe pain, but chronic use leads to maladaptations that include: tolerance, dependence and opioid-induced hyperalgesia (OIH). These responses limit the utility of opioids, as well as our ability to control chronic pain. Despite decades of research, we have no therapies or proven strategies to overcome this problem. However, murine haplotype based computational genetic mapping and a SNP data base generated from analysis of whole-genome sequence data (whole-genome HBCGM), provides a hypothesis-free method for discovering novel genes affecting opioid maladaptive responses. RESULTS: Whole genome-HBCGM was used to analyze phenotypic data on morphine-induced tolerance, dependence and hyperalgesia obtained from 23 inbred strains. The robustness of the genetic mapping results was analyzed using strain subsets. In addition, the results of analyzing all of the opioid-related traits together were examined. To characterize the functional role of the leading candidate gene, we analyzed transgenic animals, mRNA and protein expression in behaviorally divergent mouse strains, and immunohistochemistry in spinal cord tissue. Our mapping procedure identified the allelic pattern within the netrin-1 receptor gene (Dcc) as most robustly associated with OIH, and it was also strongly associated with the combination of the other maladaptive opioid traits analyzed. Adult mice heterozygous for the Dcc gene had significantly less tendency to develop OIH, become tolerant or show evidence of dependence after chronic exposure to morphine. The difference in opiate responses was shown not to be due to basal or morphine-stimulated differences in the level of Dcc expression in spinal cord tissue, and was not associated with nociceptive neurochemical or anatomical alterations in the spinal cord or dorsal root ganglia in adult animals. CONCLUSIONS: Whole-genome HBCGM is a powerful tool for identifying genes affecting biomedical traits such as opioid maladaptations. We demonstrate that Dcc affects tolerance, dependence and OIH after chronic opioid exposure, though not through simple differences in expression in the adult spinal cord.

Fialuridine induces acute liver failure in chimeric TK-NOG mice: a model for detecting hepatic drug toxicity prior to human testing.

BACKGROUND: Seven of 15 clinical trial participants treated with a nucleoside analogue (fialuridine [FIAU]) developed acute liver failure. Five treated participants died, and two required a liver transplant. Preclinical toxicology studies in mice, rats, dogs, and primates did not provide any indication that FIAU would be hepatotoxic in humans. Therefore, we investigated whether FIAU-induced liver toxicity could be detected in chimeric TK-NOG mice with humanized livers. METHODS AND FINDINGS: Control and chimeric TK-NOG mice with humanized livers were treated orally with FIAU 400, 100, 25, or 2.5 mg/kg/d. The response to drug treatment was evaluated by measuring plasma lactate and liver enzymes, by assessing liver histology, and by electron microscopy. After treatment with FIAU 400 mg/kg/d for 4 d, chimeric mice developed clinical and serologic evidence of liver failure and lactic acidosis. Analysis of liver tissue revealed steatosis in regions with human, but not mouse, hepatocytes. Electron micrographs revealed lipid and mitochondrial abnormalities in the human hepatocytes in FIAU-treated chimeric mice. Dose-dependent liver toxicity was detected in chimeric mice treated with FIAU 100, 25, or 2.5 mg/kg/d for 14 d. Liver toxicity did not develop in control mice that were treated with the same FIAU doses for 14 d. In contrast, treatment with another nucleotide analogue (sofosbuvir 440 or 44 mg/kg/d po) for 14 d, which did not cause liver toxicity in human trial participants, did not cause liver toxicity in mice with humanized livers. CONCLUSIONS: FIAU-induced liver toxicity could be readily detected using chimeric TK-NOG mice with humanized livers, even when the mice were treated with a FIAU dose that was only 10-fold above the dose used in human participants. The clinical features, laboratory abnormalities, liver histology, and ultra-structural changes observed in FIAU-treated chimeric mice mirrored those of FIAU-treated human participants. The use of chimeric mice in preclinical toxicology studies could improve the safety of candidate medications selected for testing in human participants. Please see later in the article for the Editors' Summary.

Preconception opioids interact with mouse strain to alter morphine withdrawal in the next generation.

RATIONALE: Transgenerational effects of preconception morphine exposure in female rats have been reported which suggest that epigenetic modifications triggered by female opioid exposure, even when that exposure ends several weeks prior to pregnancy, has significant ramifications for their future offspring. OBJECTIVE: The current study compares two mouse strains with well-established genetic variation in their response to mu opioid receptor agonists, C57BL/6J (BL6) and 129S1/svlmJ (129) to determine whether genetic background modifies the impact of preconception opioid exposure. METHODS: Adolescent females from both strains were injected daily with morphine for a total of 10 days using an increasing dosing regimen with controls receiving saline. Several weeks after their final injection, aged-matched BL6 and 129 morphine (Mor-F0) or saline (Sal-F0) females were mated with drug naïve males to generate Mor-F1 and Sal-F1 offspring, respectively. As adults, F1 mice were made morphine dependent using thrice daily morphine injections for 4 days. On day 5, mice were administered either saline or morphine followed 3 h later by naloxone. Behavioral and physiological signs of withdrawal were then measured. RESULTS: Regardless of strain or sex, morphine-dependent Mor-F1 mice had significantly lower levels of withdrawal-induced corticosterone but significantly higher glucose levels when compared to Sal-F1 controls. In contrast, both strain- and preconception opioid exposure effects on physical signs of morphine dependence were observed.