Early life exposure to vitamin D deficiency impairs molecular mechanisms that regulate liver cholesterol biosynthesis, energy metabolism, inflammation, and detoxification.

Introduction: Emerging data suggests liver disease may be initiated during development when there is high genome plasticity and the molecular pathways supporting liver function are being developed. Methods: Here, we leveraged our Collaborative Cross mouse model of developmental vitamin D deficiency (DVD) to investigate the role of DVD in dysregulating the molecular mechanisms underlying liver disease. We defined the effects on the adult liver transcriptome and metabolome and examined the role of epigenetic dysregulation. Given that the parental origin of the genome (POG) influences response to DVD, we used our established POG model [POG1-(CC011xCC001)F1 and POG2-(CC001xCC011)F1] to identify interindividual differences. Results: We found that DVD altered the adult liver transcriptome, primarily downregulating genes controlling liver development, response to injury/infection (detoxification & inflammation), cholesterol biosynthesis, and energy production. In concordance with these transcriptional changes, we found that DVD decreased liver cell membrane-associated lipids (including cholesterol) and pentose phosphate pathway metabolites. Each POG also exhibited distinct responses. POG1 exhibited almost 2X more differentially expressed genes (DEGs) with effects indicative of increased energy utilization. This included upregulation of lipid and amino acid metabolism genes and increased intermediate lipid and amino acid metabolites, increased energy cofactors, and decreased energy substrates. POG2 exhibited broader downregulation of cholesterol biosynthesis genes with a metabolomics profile indicative of decreased energy utilization. Although DVD primarily caused loss of liver DNA methylation for both POGs, only one epimutation was shared, and POG2 had 6.5X more differentially methylated genes. Differential methylation was detected at DEGs regulating developmental processes such as amino acid transport (POG1) and cell growth & differentiation (e.g., Wnt & cadherin signaling, POG2). Conclusions: These findings implicate a novel role for maternal vitamin D in programming essential offspring liver functions that are dysregulated in liver disease. Importantly, impairment of these processes was not rescued by vitamin D treatment at weaning, suggesting these effects require preventative measures. Substantial differences in POG response to DVD demonstrate that the parental genomic context of exposure determines offspring susceptibility.

The collaborative cross strains and their founders vary widely in cocaine-induced behavioral sensitization.

Cocaine use and overdose deaths attributed to cocaine have increased significantly in the United States in the last 10 years. Despite the prevalence of cocaine use disorder (CUD) and the personal and societal problems it presents, there are currently no approved pharmaceutical treatments. The absence of treatment options is due, in part, to our lack of knowledge about the etiology of CUDs. There is ample evidence that genetics plays a role in increasing CUD risk but thus far, very few risk genes have been identified in human studies. Genetic studies in mice have been extremely useful for identifying genetic loci and genes, but have been limited to very few genetic backgrounds, leaving substantial phenotypic, and genetic diversity unexplored. Herein we report the measurement of cocaine-induced behavioral sensitization using a 19-day protocol that captures baseline locomotor activity, initial locomotor response to an acute exposure to cocaine and locomotor sensitization across 5 exposures to the drug. These behaviors were measured in 51 genetically diverse Collaborative Cross (CC) strains along with their inbred founder strains. The CC was generated by crossing eight genetically diverse inbred strains such that each inbred CC strain has genetic contributions from each of the founder strains. Inbred CC mice are infinitely reproducible and provide a stable, yet diverse genetic platform on which to study the genetic architecture and genetic correlations among phenotypes. We have identified significant differences in cocaine locomotor sensitivity and behavioral sensitization across the panel of CC strains and their founders. We have established relationships among cocaine sensitization behaviors and identified extreme responding strains that can be used in future studies aimed at understanding the genetic, biological, and pharmacological mechanisms that drive addiction-related behaviors. Finally, we have determined that these behaviors exhibit relatively robust heritability making them amenable to future genetic mapping studies to identify addiction risk genes and genetic pathways that can be studied as potential targets for the development of novel therapeutics.

Cocaine-Induced Locomotor Activation Differs Across Inbred Mouse Substrains.

Cocaine use disorders (CUD) are devastating for affected individuals and impose a significant societal burden, but there are currently no FDA-approved therapies. The development of novel and effective treatments has been hindered by substantial gaps in our knowledge about the etiology of these disorders. The risk for developing a CUD is influenced by genetics, the environment and complex interactions between the two. Identifying specific genes and environmental risk factors that increase CUD risk would provide an avenue for the development of novel treatments. Rodent models of addiction-relevant behaviors have been a valuable tool for studying the genetics of behavioral responses to drugs of abuse. Traditional genetic mapping using genetically and phenotypically divergent inbred mice has been successful in identifying numerous chromosomal regions that influence addiction-relevant behaviors, but these strategies rarely result in identification of the causal gene or genetic variant. To overcome this challenge, reduced complexity crosses (RCC) between closely related inbred mouse strains have been proposed as a method for rapidly identifying and validating functional variants. The RCC approach is dependent on identifying phenotypic differences between substrains. To date, however, the study of addiction-relevant behaviors has been limited to very few sets of substrains, mostly comprising the C57BL/6 lineage. The present study expands upon the current literature to assess cocaine-induced locomotor activation in 20 inbred mouse substrains representing six inbred strain lineages (A/J, BALB/c, FVB/N, C3H/He, DBA/2 and NOD) that were either bred in-house or supplied directly by a commercial vendor. To our knowledge, we are the first to identify significant differences in cocaine-induced locomotor response in several of these inbred substrains. The identification of substrain differences allows for the initiation of RCC populations to more rapidly identify specific genetic variants associated with acute cocaine response. The observation of behavioral profiles that differ between mice generated in-house and those that are vendor-supplied also presents an opportunity to investigate the influence of environmental factors on cocaine-induced locomotor activity.

Bayesian modeling of skewed X inactivation in genetically diverse mice identifies a novel Xce allele associated with copy number changes.

Female mammals are functional mosaics of their parental X-linked gene expression due to X chromosome inactivation (XCI). This process inactivates one copy of the X chromosome in each cell during embryogenesis and that state is maintained clonally through mitosis. In mice, the choice of which parental X chromosome remains active is determined by the X chromosome controlling element (Xce), which has been mapped to a 176-kb candidate interval. A series of functional Xce alleles has been characterized or inferred for classical inbred strains based on biased, or skewed, inactivation of the parental X chromosomes in crosses between strains. To further explore the function structure basis and location of the Xce, we measured allele-specific expression of X-linked genes in a large population of F1 females generated from Collaborative Cross (CC) strains. Using published sequence data and applying a Bayesian "Pólya urn" model of XCI skew, we report two major findings. First, inter-individual variability in XCI suggests mouse epiblasts contain on average 20-30 cells contributing to brain. Second, CC founder strain NOD/ShiLtJ has a novel and unique functional allele, Xceg, that is the weakest in the Xce allelic series. Despite phylogenetic analysis confirming that NOD/ShiLtJ carries a haplotype almost identical to the well-characterized C57BL/6J (Xceb), we observed unexpected patterns of XCI skewing in females carrying the NOD/ShiLtJ haplotype within the Xce. Copy number variation is common at the Xce locus and we conclude that the observed allelic series is a product of independent and recurring duplications shared between weak Xce alleles.

Behavioral characterization of a novel Cisd2 mutant mouse.

Wolfram syndrome (WFS) is a rare autosomal recessive disorder characterized by diabetes mellitus and insipidus, progressive optic atrophy and sensorineural deafness. An increased incidence of psychiatric disorders has also been reported in WFS patients. There are two subtypes of WFS. Type 1 (WFS1) is caused by mutations in the WFS1 gene and type 2 (WFS2) results from mutations in the CISD2 gene. Existing Wfs1 knockout mice exhibit many WFS1 cardinal symptoms including diabetic nephropathy, metabolic disruptions and optic atrophy. Far fewer studies have examined loss of Cisd2 function in mice. We identified B6.DDY-Cisd2m1Lmt, a mouse model with a spontaneous mutation in the Cisd2 gene. B6.DDY-Cisd2m1Lmt mice were initially identified based on the presence of audible sonic vocalizations as well as decreased body size and weight compared to unaffected wildtype littermates. Although Wfs1 knockout mice have been characterized for numerous behavioral phenotypes, similar studies have been lacking for Cisd2 mutant mice. We tested B6.DDY-Cisd2m1Lmt mice in a battery of behavioral assays that model phenotypes related to neurological and psychiatric disorders including anxiety, sensorimotor gating, stress response, social interaction and learning and memory. B6.DDY-Cisd2m1Lmt mice displayed hypoactivity across several behavioral tests, exhibited increased stress response and had deficits in spatial learning and memory and sensorimotor gating compared to wildtype littermates. Our data indicate that the B6.DDY-Cisd2m1Lmt mouse strain is a useful model to investigate potential mechanisms underlying the neurological and psychiatric symptoms observed in WFS.

Characterization of genetically complex Collaborative Cross mouse strains that model divergent locomotor activating and reinforcing properties of cocaine.

RATIONALE: Few effective treatments exist for cocaine use disorders due to gaps in knowledge about its complex etiology. Genetically defined animal models provide a useful tool for advancing our understanding of the biological and genetic underpinnings of addiction-related behavior and evaluating potential treatments. However, many attempts at developing mouse models of behavioral disorders were based on overly simplified single gene perturbations, often leading to inconsistent and misleading results in pre-clinical pharmacology studies. A genetically complex mouse model may better reflect disease-related behaviors. OBJECTIVES: Screening defined, yet genetically complex, intercrosses of the Collaborative Cross (CC) mice revealed two lines, RIX04/17 and RIX41/51, with extreme high and low behavioral responses to cocaine. We characterized these lines as well as their CC parents, CC004/TauUnc and CC041/TauUnc, to evaluate their utility as novel model systems for studying the biological and genetic mechanisms underlying behavioral responses to cocaine. METHODS: Behavioral responses to acute (initial locomotor sensitivity) and repeated (behavioral sensitization, conditioned place preference, intravenous self-administration) exposures to cocaine were assessed. We also examined the monoaminergic system (striatal tissue content and in vivo fast scan cyclic voltammetry), HPA axis reactivity, and circadian rhythms as potential mechanisms for the divergent phenotypic behaviors observed in the two strains, as these systems have a previously known role in mediating addiction-related behaviors. RESULTS: RIX04/17 and 41/51 show strikingly divergent initial locomotor sensitivity to cocaine with RIX04/17 exhibiting very high and RIX41/51 almost no response. The lines also differ in the emergence of behavioral sensitization with RIX41/51 requiring more exposures to exhibit a sensitized response. Both lines show conditioned place preference for cocaine. We determined that the cocaine sensitivity phenotype in each RIX line was largely driven by the genetic influence of one CC parental strain, CC004/TauUnc and CC041/TauUnc. CC004 demonstrates active operant cocaine self-administration and CC041 is unable to acquire under the same testing conditions, a deficit which is specific to cocaine as both strains show operant response for a natural food reward. Examination of potential mechanisms driving differential responses to cocaine show strain differences in molecular and behavioral circadian rhythms. Additionally, while there is no difference in striatal dopamine tissue content or dynamics, there are selective differences in striatal norepinephrine and serotonergic tissue content. CONCLUSIONS: These CC strains offer a complex polygenic model system to study underlying mechanisms of cocaine response. We propose that CC041/TauUnc and CC004/TauUnc will be useful for studying genetic and biological mechanisms underlying resistance or vulnerability to the stimulatory and reinforcing effects of cocaine.

National Institute on Drug Abuse genomics consortium white paper: Coordinating efforts between human and animal addiction studies.

The National Institute on Drug Abuse Genetics and Epigenetics Cross-Cutting Research Team convened a diverse group of researchers, clinicians, and healthcare providers on the campus of the University of California, San Diego, in June 2018. The goal was to develop strategies to integrate genetics and phenotypes across species to achieve a better understanding of substance use disorders through associations between genotypes and addictive behaviors. This conference (a) discussed progress in harmonizing large opioid genetics cohorts, (b) discussed phenotypes that are used for genetics studies in humans, (c) examined phenotypes that are used for genetics studies in animal models, (d) identified synergies and gaps in phenotypic analyses of human and animal models and (e) identified strategies to integrate genetics and genomics data with phenotypes across species. The meeting consisted of panels that focused on phenotype harmonization (Dr. Laura Bierut, Dr. Olivier George, Dr. Dan Larach and Dr. Sesh Mudumbai), translating genetic findings between species (Dr. Elissa Chesler, Dr. Gary Peltz and Dr. Abraham Palmer), interpreting and understanding allelic variations (Dr. Vanessa Troiani and Dr. Tamara Richards) and pathway conservation in animal models and human studies (Dr. Robert Hitzemann, Dr. Huda Akil and Dr. Laura Saba). There were also updates that were provided by large consortia (Dr. Susan Tapert, Dr. Danielle Dick, Dr. Howard Edenberg and Dr. Eric Johnson). Collectively, the conference was convened to discuss progress and changes in genome-wide association studies.

Reciprocal F1 Hybrids of Two Inbred Mouse Strains Reveal Parent-of-Origin and Perinatal Diet Effects on Behavior and Expression.

Parent-of-origin effects (POE) in mammals typically arise from maternal effects or imprinting. In some instances, such POE have been associated with psychiatric disorders, as well as with changes in a handful of animal behaviors. However, POE on complex traits such as behavior remain largely uncharacterized. Moreover, although both behavior and epigenetic effects are known to be modified by perinatal environmental exposures such as nutrient deficiency, the architecture of such environment-by-POE is mostly unexplored. To study POE and environment-by-POE, we employ a relatively neglected but especially powerful experimental system for POE-detection: reciprocal F1 hybrids (RF1s). We exposed female NOD/ShiLtJ×C57Bl/6J and C57Bl/6J×NOD/ShiLtJ mice, perinatally, to one of four different diets, then after weaning recorded a set of behaviors that model psychiatric disease. Whole-brain microarray expression data revealed an imprinting-enriched set of 15 genes subject to POE. The most-significant expression POE, on the non-imprinted gene Carmil1 (a.k.a. Lrrc16a), was validated using qPCR in the same and in a new set of mice. Several behaviors, especially locomotor behaviors, also showed POE. Bayesian mediation analysis suggested Carmil1 expression suppresses behavioral POE, and that the imprinted gene Airn suppresses POE on Carmil1 expression. A suggestive diet-by-POE was observed on percent center time in the open field test, and a significant diet-by-POE was observed on one imprinted gene, Mir341, and on 16 non-imprinted genes. The relatively small, tractable set of POE and diet-by-POE detected on behavior and expression here motivates further studies examining such effects across RF1s on multiple genetic backgrounds.

Maternal vitamin D depletion alters DNA methylation at imprinted loci in multiple generations.

BACKGROUND: Environmental perturbation of epigenetic mechanisms is linked to a growing number of diseases. Characterizing the role environmental factors play in modifying the epigenome is important for disease etiology. Vitamin D is an essential nutrient affecting brain, bone, heart, immune and reproductive health. Vitamin D insufficiency is a global issue, and the role in maternal and child health remains under investigation. METHODS: We used Collaborative Cross (CC) inbred mice to characterize the effect of maternal vitamin D depletion on offspring phenotypic and epigenetic outcomes at imprinted domains (H19/Igf2, Snrpn, Dlk1/Gtl2, and Grb10) in the soma (liver) and germline (sperm). We assessed outcomes in two generations of offspring to determine heritability. We used reciprocal crosses between lines CC001/Unc and CC011/Unc to investigate parent of origin effects. RESULTS: Maternal vitamin D deficiency led to altered body weight and DNA methylation in two generations of offspring. Loci assayed in adult liver and sperm were mostly hypomethylated, but changes were few and small in effect size (<7 % difference on average). There was no change in total expression of genes adjacent to methylation changes in neonatal liver. Methylation changes were cell type specific such that changes at IG-DMR were present in sperm but not in liver. Some methylation changes were distinct between generations such that methylation changes at the H19ICR in second-generation liver were not present in first-generation sperm or liver. Interestingly, some diet-dependent changes in body weight and methylation were seemingly influenced by parent of origin such that reciprocal crosses exhibited inverse effects. CONCLUSIONS: These findings demonstrate that maternal vitamin D status plays a role in determining DNA methylation state in the germline and soma. Detection of methylation changes in the unexposed second-generation demonstrates that maternal vitamin D depletion can have long-term effects on the epigenome of subsequent generations. Differences in vitamin D-dependent epigenetic state between cell types and generations indicate perturbation of the epigenetic landscape rather than a targeted, locus-specific effect. While the biological importance of these subtle changes remains unclear, they warrant an investigation of epigenome-wide effects of maternal vitamin D depletion.