L-NAC and L-NAC methyl ester prevent and overcome physical dependence to fentanyl in male rats

N-acetyl-L-cysteine (L-NAC) is a proposed therapeutic for opioid use disorder. This study determined whether co-injections of L-NAC (500 μmol/kg, IV) or its highly cell-penetrant analogue, L-NAC methyl ester (L-NACme, 500 μmol/kg, IV), prevent acquisition of acute physical dependence induced by twice-daily injections of fentanyl (125 μg/kg, IV), and overcome acquired dependence to these injections in freely-moving male Sprague Dawley rats. The injection of the opioid receptor antagonist, naloxone HCl (NLX; 1.5 mg/kg, IV), elicited a series of withdrawal phenomena (i.e. behavioral and cardiorespiratory responses, hypothermia and body weight loss) in rats that received 5 or 10 injections of fentanyl and similar numbers of vehicle co-injections. With respect to the development of dependence, the NLX-precipitated withdrawal phenomena were reduced in rats that received had co-injections of L-NAC, and more greatly reduced in rats that received co-injections of L-NACme. In regard to overcoming established dependence, the NLX-precipitated withdrawal phenomena in rats that had received 10 injections of fentanyl (125 μg/kg, IV) were reduced in rats that had received co-injections of L-NAC, and more greatly reduced in rats that received co-injections of L-NACme beginning with injection 6 of fentanyl. This study provides compelling evidence that co-injections of L-NAC and L-NACme prevent the acquisition of physical dependence and overcome acquired dependence to fentanyl in male rats. The higher efficacy of L-NACme is likely due to its greater cell penetrability in brain regions mediating dependence to fentanyl and interaction with intracellular signaling cascades, including redox-dependent processes, responsible for the acquisition of physical dependence to fentanyl.

Lipophilic analogues of D-cysteine prevent and reverse physical dependence to fentanyl in male rats

We examined whether co-injections of the cell-permeant D-cysteine analogues, D-cysteine ethyl ester (D-CYSee) and D-cysteine ethyl amide (D-CYSea), prevent acquisition of physical dependence induced by twice-daily injections of fentanyl, and reverse acquired dependence to these injections in freely-moving male Sprague Dawley rats. Injection of the opioid receptor antagonist, naloxone HCl (NLX, 1.5 mg/kg, IV), elicited a series of withdrawal phenomena that included cardiorespiratory and behavioral responses, and falls in body weight and body temperature, in rats that received 5 or 10 injections of fentanyl (125 μg/kg, IV), and the same number of vehicle co-injections. Regarding the development of physical dependence, the NLX-precipitated withdrawal phenomena were markedly reduced in fentanyl-injected rats that had received co-injections of D-CYSee (250 μmol/kg, IV) or D-CYSea (100 μmol/kg, IV), but not D-cysteine (250 μmol/kg, IV). Regarding reversal of established dependence to fentanyl, the NLX-precipitated withdrawal phenomena in rats that had received 10 injections of fentanyl (125 μg/kg, IV) was markedly reduced in rats that received co-injections of D-CYSee (250 μmol/kg, IV) or D-CYSea (100 μmol/kg, IV), but not D-cysteine (250 μmol/kg, IV), starting with injection 6 of fentanyl. This study provides evidence that co-injections of D-CYSee and D-CYSea prevent the acquisition of physical dependence, and reverse acquired dependence to fentanyl in male rats. The lack of effect of D-cysteine suggests that the enhanced cell-penetrability of D-CYSee and D-CYSea into cells, particularly within the brain, is key to their ability to interact with intracellular signaling events involved in acquisition to physical dependence to fentanyl.

Discovery and validation of genes driving drug-intake and related behavioral traits in mice

Substance use disorders are heritable disorders characterized by compulsive drug use, the biological mechanisms for which remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference and sensation seeking, are predictive of drug-use phenotypes, thereby implicating shared genetic mechanisms. High-throughput behavioral screening in knockout (KO) mice allows efficient discovery of the function of genes. We used this strategy in two rounds of candidate prioritization in which we identified 33 drug-use candidate genes based upon predisposing drug-naïve phenotypes and ultimately validated the perturbation of 22 genes as causal drivers of substance intake. We selected 19/221 KO strains (8.5%) that had a difference from control on at least one drug-naïve predictive behavioral phenotype and determined that 15/19 (~80%) affected the consumption or preference for alcohol, methamphetamine or both. No mutant exhibited a difference in nicotine consumption or preference which was possibly confounded with saccharin. In the second round of prioritization, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15/401 KO strains (3.7%, which included one gene from the first cohort) that differed most from controls for the predisposing phenotypes. 8 of 15 gene deletions (53%) affected intake or preference for alcohol, methamphetamine or both. Using multivariate and bioinformatic analyses, we observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

L-cysteine ethyl ester prevents and reverses acquired physical dependence on morphine in male Sprague Dawley rats

The molecular mechanisms underlying the acquisition of addiction/dependence on morphine may result from the ability of the opioid to diminish the transport of L-cysteine into neurons via inhibition of excitatory amino acid transporter 3 (EAA3). The objective of this study was to determine whether the co-administration of the cell-penetrant L-thiol ester, L-cysteine ethyl ester (L-CYSee), would reduce physical dependence on morphine in male Sprague Dawley rats. Injection of the opioid-receptor antagonist, naloxone HCl (NLX; 1.5 mg/kg, IP), elicited pronounced withdrawal phenomena in rats which received a subcutaneous depot of morphine (150 mg/kg) for 36 h and were receiving a continuous infusion of saline (20 μL/h, IV) via osmotic minipumps for the same 36 h period. The withdrawal phenomena included wet-dog shakes, jumping, rearing, fore-paw licking, 360° circling, writhing, apneas, cardiovascular (pressor and tachycardia) responses, hypothermia, and body weight loss. NLX elicited substantially reduced withdrawal syndrome in rats that received an infusion of L-CYSee (20.8 μmol/kg/h, IV) for 36 h. NLX precipitated a marked withdrawal syndrome in rats that had received subcutaneous depots of morphine (150 mg/kg) for 48 h) and a co-infusion of vehicle. However, the NLX-precipitated withdrawal signs were markedly reduced in morphine (150 mg/kg for 48 h)-treated rats that began receiving an infusion of L-CYSee (20.8 μmol/kg/h, IV) at 36 h. In similar studies to those described previously, neither L-cysteine nor L-serine ethyl ester (both at 20.8 μmol/kg/h, IV) mimicked the effects of L-CYSee. This study demonstrates that 1) L-CYSee attenuates the development of physical dependence on morphine in male rats and 2) prior administration of L-CYSee reverses morphine dependence, most likely by intracellular actions within the brain. The lack of the effect of L-serine ethyl ester (oxygen atom instead of sulfur atom) strongly implicates thiol biochemistry in the efficacy of L-CYSee. Accordingly, L-CYSee and analogs may be a novel class of therapeutics that ameliorate the development of physical dependence on opioids in humans.

Short-chain-fatty acid valerate reduces voluntary alcohol intake in male mice

Background

Despite serious health and social consequences, effective intervention strategies for habitual alcohol binge drinking are lacking. Development of novel therapeutic and preventative approaches is highly desirable. Accumulating evidence in the past several years has established associations between the gut microbiome and microbial metabolites with drinking behavior, but druggable targets and their underlying mechanism of action are understudied.

Results

Here, using a drink-in-the-dark mouse model, we identified a microbiome metabolite-based novel treatment (sodium valerate) that can reduce excessive alcohol drinking. Sodium valerate is a sodium salt of valeric acidshort-chain-fatty-acid with similar structure as γ-aminobutyric acid (GABA). Ten days of oral sodium valerate supplementation attenuates excessive alcohol drinking by 40%, reduces blood ethanol concentration by 53%, and improves anxiety-like or approach-avoidance behavior in male mice, without affecting overall food and water intake. Mechanistically, sodium valerate supplementation increases GABA levels across stool, blood, and amygdala. It also significantly increases H4 acetylation in the amygdala of mice. Transcriptomics analysis of the amygdala revealed that sodium valerate supplementation led to changes in gene expression associated with functional pathways including potassium voltage-gated channels, inflammation, glutamate degradation, L-DOPA degradation, and psychological behaviors. 16S microbiome profiling showed that sodium valerate supplementation shifts the gut microbiome composition and decreases microbiome-derived neuroactive compounds through GABA degradation in the gut microbiome.

Conclusion

Our findings suggest that the sodium valerate holds promise as an innovative therapeutic avenue for the reduction of habitual binge drinking, potentially through multifaceted mechanisms.

Genetics and Genomics of Addiction

Substance use disorders (SUD), like many neuropsychiatric conditions, are a heterogeneous group of disorders with similar symptomatology but often different pathoetiology [...].

An allelic series of spontaneous Rorb mutant mice exhibit a gait phenotype, changes in retina morphology and behavior, and gene expression signatures associated with the unfolded protein response

The Retinoid-related orphan receptor beta (RORβ) gene encodes a developmental transcription factor and has 2 predominant isoforms created through alternative first exon usage; one specific to the retina and another present more broadly in the central nervous system, particularly regions involved in sensory processing. RORβ belongs to the nuclear receptor family and plays important roles in cell fate specification in the retina and cortical layer formation. In mice, loss of RORβ causes disorganized retina layers, postnatal degeneration, and production of immature cone photoreceptors. Hyperflexion or "high-stepping" of rear limbs caused by reduced presynaptic inhibition by Rorb-expressing inhibitory interneurons of the spinal cord is evident in RORβ-deficient mice. RORβ variants in patients are associated with susceptibility to various neurodevelopmental conditions, primarily generalized epilepsies, but including intellectual disability, bipolar, and autism spectrum disorders. The mechanisms by which RORβ variants confer susceptibility to these neurodevelopmental disorders are unknown but may involve aberrant neural circuit formation and hyperexcitability during development. Here we report an allelic series in 5 strains of spontaneous Rorb mutant mice with a high-stepping gait phenotype. We show retinal abnormalities in a subset of these mutants and demonstrate significant differences in various behavioral phenotypes related to cognition. Gene expression analyses in all 5 mutants reveal a shared over-representation of the unfolded protein response and pathways related to endoplasmic reticulum stress, suggesting a possible mechanism of susceptibility relevant to patients.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

DISCOVERY AND VALIDATION OF GENES DRIVING DRUG-INTAKE AND RELATED BEHAVIORAL TRAITS IN MICE

Substance use disorders (SUDs) are heritable disorders characterized by compulsive drug use, but the biological mechanisms driving addiction remain largely unknown. Genetic correlations reveal that predisposing drug-naïve phenotypes, including anxiety, depression, novelty preference, and sensation seeking, are predictive of drug-use phenotypes, implicating shared genetic mechanisms. Because of this relationship, high-throughput behavioral screening of predictive phenotypes in knockout (KO) mice allows efficient discovery of genes likely to be involved in drug use. We used this strategy in two rounds of screening in which we identified 33 drug-use candidate genes and ultimately validated the perturbation of 22 of these genes as causal drivers of substance intake. In our initial round of screening, we employed the two-bottle-choice paradigms to assess alcohol, methamphetamine, and nicotine intake. We identified 19 KO strains that were extreme responders on at least one predictive phenotype. Thirteen of the 19 gene deletions (68%) significantly affected alcohol use three methamphetamine use, and two both. In the second round of screening, we employed a multivariate approach to identify outliers and performed validation using methamphetamine two-bottle choice and ethanol drinking-in-the-dark protocols. We identified 15 KO strains that were extreme responders across the predisposing drug-naïve phenotypes. Eight of the 15 gene deletions (53%) significantly affected intake or preference for three alcohol, eight methamphetamine or three both (3). We observed multiple relations between predisposing behaviors and drug intake, revealing many distinct biobehavioral processes underlying these relationships. The set of mouse models identified in this study can be used to characterize these addiction-related processes further.

The microbial community dynamics of cocaine sensitization in two behaviorally divergent strains of collaborative cross mice

The gut-brain axis is increasingly recognized as an important pathway involved in cocaine use disorder. Microbial products of the murine gut have been shown to affect striatal gene expression, and depletion of the microbiome by antibiotic treatment alters cocaine-induced behavioral sensitization in C57BL/6J male mice. Some reports suggest that cocaine-induced behavioral sensitization is correlated with drug self-administration behavior in mice. Here, we profile the composition of the naïve microbiome and its response to cocaine sensitization in two collaborative cross (CC) strains. These strains display extremely divergent behavioral responses to cocaine sensitization. A high-responding strain, CC004/TauUncJ (CC04), has a gut microbiome that contains a greater amount of Lactobacillus than the cocaine-nonresponsive strain CC041/TauUncJ (CC41). The gut microbiome of CC41 is characterized by an abundance of Eisenbergella, Robinsonella and Ruminococcus. In response to cocaine, CC04 has an increased Barnsiella population, while the gut microbiome of CC41 displays no significant changes. PICRUSt functional analysis of the functional potential of the gut microbiome in CC04 shows a significant number of potential gut-brain modules altered after exposure to cocaine, specifically those encoding for tryptophan synthesis, glutamine metabolism, and menaquinone synthesis (vitamin K2). Depletion of the microbiome by antibiotic treatment revealed an altered cocaine-sensitization response following antibiotics in female CC04 mice. Depleting the microbiome by antibiotic treatment in males revealed increased infusions for CC04 during a cocaine intravenous self-administration dose-response curve. Together these data suggest that genetic differences in cocaine-related behaviors may involve the microbiome.

Gene expression genetics of the striatum of Diversity Outbred mice

Brain transcriptional variation is a heritable trait that mediates complex behaviors, including addiction. Expression quantitative trait locus (eQTL) mapping reveals genomic regions harboring genetic variants that influence transcript abundance. In this study, we profiled transcript abundance in the striatum of 386 Diversity Outbred (J:DO) mice of both sexes using RNA-Seq. All mice were characterized using a behavioral battery of widely-used exploratory and risk-taking assays prior to transcriptional profiling. We performed eQTL mapping, incorporated the results into a browser-based eQTL viewer, and deposited co-expression network members in GeneWeaver. The eQTL viewer allows researchers to query specific genes to obtain allelic effect plots, analyze SNP associations, assess gene expression correlations, and apply mediation analysis to evaluate whether the regulatory variant is acting through the expression of another gene. GeneWeaver allows multi-species comparison of gene sets using statistical and combinatorial tools. This data resource allows users to find genetic variants that regulate differentially expressed transcripts and place them in the context of other studies of striatal gene expression and function in addiction-related behavior.

Microbial glutamate metabolism predicts intravenous cocaine self-administration in diversity outbred mice

The gut microbiome is thought to play a critical role in the onset and development of psychiatric disorders, including depression and substance use disorder (SUD). To test the hypothesis that the microbiome affects addiction predisposing behaviors and cocaine intravenous self-administration (IVSA) and to identify specific microbes involved in the relationship, we performed 16S rRNA gene sequencing on feces from 228 diversity outbred mice. Twelve open field measures, two light-dark assay measures, one hole board and novelty place preference measure significantly differed between mice that acquired cocaine IVSA (ACQ) and those that failed to acquire IVSA (FACQ). We found that ACQ mice are more active and exploratory and display decreased fear than FACQ mice. The microbial abundances that differentiated ACQ from FACQ mice were an increased abundance of Barnesiella, Ruminococcus, and Robinsoniella and decreased Clostridium IV in ACQ mice. There was a sex-specific correlation between ACQ and microbial abundance, a reduced Lactobacillus abundance in ACQ male mice, and a decreased Blautia abundance in female ACQ mice. The abundance of Robinsoniella was correlated, and Clostridium IV inversely correlated with the number of doses of cocaine self-administered during acquisition. Functional analysis of the microbiome composition of a subset of mice suggested that gut-brain modules encoding glutamate metabolism genes are associated with the propensity to self-administer cocaine. These findings establish associations between the microbiome composition and glutamate metabolic potential and the ability to acquire cocaine IVSA thus indicating the potential translational impact of targeting the gut microbiome or microbial metabolites for treatment of SUD. This article is part of the Special Issue on "Microbiome & the Brain: Mechanisms & Maladies".

Identifying genetic determinants of inflammatory pain in mice using a large-scale gene-targeted screen

ABSTRACT

Identifying the genetic determinants of pain is a scientific imperative given the magnitude of the global health burden that pain causes. Here, we report a genetic screen for nociception, performed under the auspices of the International Mouse Phenotyping Consortium. A biased set of 110 single-gene knockout mouse strains was screened for 1 or more nociception and hypersensitivity assays, including chemical nociception (formalin) and mechanical and thermal nociception (von Frey filaments and Hargreaves tests, respectively), with or without an inflammatory agent (complete Freund's adjuvant). We identified 13 single-gene knockout strains with altered nocifensive behavior in 1 or more assays. All these novel mouse models are openly available to the scientific community to study gene function. Two of the 13 genes (Gria1 and Htr3a) have been previously reported with nociception-related phenotypes in genetically engineered mouse strains and represent useful benchmarking standards. One of the 13 genes (Cnrip1) is known from human studies to play a role in pain modulation and the knockout mouse reported herein can be used to explore this function further. The remaining 10 genes (Abhd13, Alg6, BC048562, Cgnl1, Cp, Mmp16, Oxa1l, Tecpr2, Trim14, and Trim2) reveal novel pathways involved in nociception and may provide new knowledge to better understand genetic mechanisms of inflammatory pain and to serve as models for therapeutic target validation and drug development.

Genes identified in rodent studies of alcohol intake are enriched for heritability of human substance use

BACKGROUND

Rodent paradigms and human genome-wide association studies (GWAS) on drug use have the potential to provide biological insight into the pathophysiology of addiction.

METHODS

Using GeneWeaver, we created rodent alcohol and nicotine gene-sets derived from 19 gene expression studies on alcohol and nicotine outcomes. We partitioned the SNP-heritability of these gene-sets using four large human GWAS: 1) alcoholic drinks per week, 2) problematic alcohol use, 3) cigarettes per day and 4) smoking cessation. We benchmarked our findings with curated human alcoholism and nicotine addiction gene-sets and performed specificity analyses using other rodent gene-sets (e.g., locomotor behavior) and other human GWAS (e.g., height).

RESULTS

The rodent alcohol gene-set was enriched for heritability of drinks per week, cigarettes per day, and smoking cessation, but not problematic alcohol use. However, the rodent nicotine gene-set was not significantly associated with any of these traits. Both rodent gene-sets showed enrichment for several non-substance use GWAS, and the extent of this relationship tended to increase as a function of trait heritability. In general, larger gene-sets demonstrated more significant enrichment. Finally, when evaluating human traits with similar heritabilities, both rodent gene-sets showed greater enrichment for substance use traits.

CONCLUSION

Our results suggest that rodent gene expression studies can help to identify genes that contribute to heritability of some substance use traits in humans, yet there was less specificity than expected. We outline various limitations, interpretations and considerations for future research.

The Gut Microbiome and Substance Use Disorder

Substance use disorders (SUDs) remain a significant public health challenge, affecting tens of millions of individuals worldwide each year. Often comorbid with other psychiatric disorders, SUD can be poly-drug and involve several different substances including cocaine, opiates, nicotine, and alcohol. SUD has a strong genetic component. Much of SUD research has focused on the neurologic and genetic facets of consumption behavior. There is now interest in the role of the gut microbiome in the pathogenesis of SUD. In this review, we summarize current animal and clinical evidence that the gut microbiome is involved in SUD, then address the underlying mechanisms by which the gut microbiome interacts with SUD through metabolomic, immune, neurological, and epigenetic mechanisms. Lastly, we discuss methods using various inbred and outbred mice models to gain an integrative understanding of the microbiome and host genetic controls in SUD.

Host genetic control of gut microbiome composition

The gut microbiome plays a significant role in health and disease, and there is mounting evidence indicating that the microbial composition is regulated in part by host genetics. Heritability estimates for microbial abundance in mice and humans range from (0.05-0.45), indicating that 5-45% of inter-individual variation can be explained by genetics. Through twin studies, genetic association studies, systems genetics, and genome-wide association studies (GWAS), hundreds of specific host genetic loci have been shown to associate with the abundance of discrete gut microbes. Using genetically engineered knock-out mice, at least 30 specific genes have now been validated as having specific effects on the microbiome. The relationships among of host genetics, microbiome composition, and abundance, and disease is now beginning to be unraveled through experiments designed to test causality. The genetic control of disease and its relationship to the microbiome can manifest in multiple ways. First, a genetic variant may directly cause the disease phenotype, resulting in an altered microbiome as a consequence of the disease phenotype. Second, a genetic variant may alter gene expression in the host, which in turn alters the microbiome, producing the disease phenotype. Finally, the genetic variant may alter the microbiome directly, which can result in the disease phenotype. In order to understand the processes that underlie the onset and progression of certain diseases, future research must take into account the relationship among host genetics, microbiome, and disease phenotype, and the resources needed to study these relationships.

Multi-omic and multi-species meta-analyses of nicotine consumption

Cross-species translational approaches to human genomic analyses are lacking. The present study uses an integrative framework to investigate how genes associated with nicotine use in model organisms contribute to the genetic architecture of human tobacco consumption. First, we created a model organism geneset by collecting results from five animal models of nicotine exposure (RNA expression changes in brain) and then tested the relevance of these genes and flanking genetic variation using genetic data from human cigarettes per day (UK BioBank N = 123,844; all European Ancestry). We tested three hypotheses: (1) DNA variation in, or around, the 'model organism geneset' will contribute to the heritability to human tobacco consumption, (2) that the model organism genes will be enriched for genes associated with human tobacco consumption, and (3) that a polygenic score based off our model organism geneset will predict tobacco consumption in the AddHealth sample (N = 1667; all European Ancestry). Our results suggested that: (1) model organism genes accounted for ~5-36% of the observed SNP-heritability in human tobacco consumption (enrichment: 1.60-31.45), (2) model organism genes, but not negative control genes, were enriched for the gene-based associations (MAGMA, H-MAGMA, SMultiXcan) for human cigarettes per day, and (3) polygenic scores based on our model organism geneset predicted cigarettes per day in an independent sample. Altogether, these findings highlight the advantages of using multiple species evidence to isolate genetic factors to better understand the etiological complexity of tobacco and other nicotine consumption.

Interpretation of psychiatric genome-wide association studies with multispecies heterogeneous functional genomic data integration

Genome-wide association studies and other discovery genetics methods provide a means to identify previously unknown biological mechanisms underlying behavioral disorders that may point to new therapeutic avenues, augment diagnostic tools, and yield a deeper understanding of the biology of psychiatric conditions. Recent advances in psychiatric genetics have been made possible through large-scale collaborative efforts. These studies have begun to unearth many novel genetic variants associated with psychiatric disorders and behavioral traits in human populations. Significant challenges remain in characterizing the resulting disease-associated genetic variants and prioritizing functional follow-up to make them useful for mechanistic understanding and development of therapeutics. Model organism research has generated extensive genomic data that can provide insight into the neurobiological mechanisms of variant action, but a cohesive effort must be made to establish which aspects of the biological modulation of behavioral traits are evolutionarily conserved across species. Scalable computing, new data integration strategies, and advanced analysis methods outlined in this review provide a framework to efficiently harness model organism data in support of clinically relevant psychiatric phenotypes.

Genetic variation regulates opioid-induced respiratory depression in mice

In the U.S., opioid prescription for treatment of pain nearly quadrupled from 1999 to 2014. The diversion and misuse of prescription opioids along with increased use of drugs like heroin and fentanyl, has led to an epidemic in addiction and overdose deaths. The most common cause of opioid overdose and death is opioid-induced respiratory depression (OIRD), a life-threatening depression in respiratory rate thought to be caused by stimulation of opioid receptors in the inspiratory-generating regions of the brain. Studies in mice have revealed that variation in opiate lethality is associated with strain differences, suggesting that sensitivity to OIRD is genetically determined. We first tested the hypothesis that genetic variation in inbred strains of mice influences the innate variability in opioid-induced responses in respiratory depression, recovery time and survival time. Using the founders of the advanced, high-diversity mouse population, the Diversity Outbred (DO), we found substantial sex and genetic effects on respiratory sensitivity and opiate lethality. We used DO mice treated with morphine to map quantitative trait loci for respiratory depression, recovery time and survival time. Trait mapping and integrative functional genomic analysis in GeneWeaver has allowed us to implicate Galnt11, an N-acetylgalactosaminyltransferase, as a gene that regulates OIRD.

Cross-Species Analyses Identify Dlgap2 as a Regulator of Age-Related Cognitive Decline and Alzheimer's Dementia

Genetic mechanisms underlying age-related cognitive decline and dementia remain poorly understood. Here, we take advantage of the Diversity Outbred mouse population to utilize quantitative trait loci mapping and identify Dlgap2 as a positional candidate responsible for modifying working memory decline. To evaluate the translational relevance of this finding, we utilize longitudinal cognitive measures from human patients, RNA expression from post-mortem brain tissue, data from a genome-wide association study (GWAS) of Alzheimer's dementia (AD), and GWAS results in African Americans. We find an association between Dlgap2 and AD phenotypes at the variant, gene and protein expression, and methylation levels. Lower cortical DLGAP2 expression is observed in AD and is associated with more plaques and tangles at autopsy and faster cognitive decline. Results will inform future studies aimed at investigating the cross-species role of Dlgap2 in regulating cognitive decline and highlight the benefit of using genetically diverse mice to prioritize novel candidates.

Do gene expression findings from mouse models of cocaine use recapitulate human cocaine use disorder in reward circuitry?

Animal models of drug use have investigated possible mechanisms governing human substance use traits for over 100 years. Most cross-species research on drug use/addiction examines behavioral overlap, but studies assessing neuromolecular (e.g. RNA) correspondence are lacking. Our study utilized transcriptome-wide data from the hippocampus and ventral tegmental area (VTA)/midbrain from a total of 35 human males with cocaine use disorder/controls and 49 male C57BL/6J cocaine/saline administering/exposed mice. We hypothesized differential expressed genes and systems of co-expressed genes (gene networks) would show appreciable overlap across mouse cocaine self-administration and human cocaine use disorder. We found modest, but significant relationships between differentially expressed genes associated with cocaine self-administration (short access) and cocaine use disorder within reward circuitry. Differentially expressed genes underlying models of acute cocaine exposure (cocaine), context re-exposure and cocaine + context re-exposure were not consistently associated with human CUD across brain regions. Investigating systems of co-expressed genes, we found several validated gene networks with weak to moderate conservation between cocaine/saline self-administering mice and disordered cocaine users/controls. The most conserved hippocampal and VTA gene networks demonstrated substantial overlap (2029 common genes) and included both novel and previously implicated targets for cocaine use/addiction. Lastly, we conducted (expression-based) phenome-wide association studies of the nine common hub genes across conserved gene networks. Common hub genes were associated with dopamine/serotonin function, cocaine self-administration and other relevant mouse traits. Overall, our study pinpointed and characterized conserved brain-related RNA patterns across mouse cocaine self-administration and human cocaine use disorder. We offer recommendations for future research and add to the dialogue surrounding pre-clinical animal research for human disease.

Discovery of a Role for Rab3b in Habituation and Cocaine Induced Locomotor Activation in Mice Using Heterogeneous Functional Genomic Analysis

Substance use disorders are prevalent and present a tremendous societal cost but the mechanisms underlying addiction behavior are poorly understood and few biological treatments exist. One strategy to identify novel molecular mechanisms of addiction is through functional genomic experimentation. However, results from individual experiments are often noisy. To address this problem, the convergent analysis of multiple genomic experiments can discern signal from these studies. In the present study, we examine genetic loci that modulate the locomotor response to cocaine identified in the recombinant inbred (BXD RI) genetic reference population. We then applied the GeneWeaver software system for heterogeneous functional genomic analysis to integrate and aggregate multiple studies of addiction genomics, resulting in the identification of Rab3b as a functional correlate of the locomotor response to cocaine in rodents. This gene encodes a member of the RAB family of Ras-like GTPases known to be involved in trafficking of secretory and endocytic vesicles in eukaryotic cells. The convergent evidence for a role of Rab3b includes co-occurrence in previously published genetic mapping studies of cocaine related behaviors; methamphetamine response and cocaine- and amphetamine-regulated transcript prepropeptide (Cartpt) transcript abundance; evidence related to other addictive substances; density of polymorphisms; and its expression pattern in reward pathways. To evaluate this finding, we examined the effect of RAB3 complex perturbation in cocaine response. B6;129-Rab3btm1Sud Rab3ctm1sud Rab3dtm1sud triple null mice (Rab3bcd -/-) exhibited significant deficits in habituation, and increased acute and repeated cocaine responses. This previously unidentified mechanism of the behavioral predisposition and response to cocaine is an example of many that can be identified and validated using aggregate genomic studies.

A Microbe Associated with Sleep Revealed by a Novel Systems Genetic Analysis of the Microbiome in Collaborative Cross Mice

The microbiome influences health and disease through complex networks of host genetics, genomics, microbes, and environment. Identifying the mechanisms of these interactions has remained challenging. Systems genetics in laboratory mice (Mus musculus) enables data-driven discovery of biological network components and mechanisms of host-microbial interactions underlying disease phenotypes. To examine the interplay among the whole host genome, transcriptome, and microbiome, we mapped QTL and correlated the abundance of cecal messenger RNA, luminal microflora, physiology, and behavior in a highly diverse Collaborative Cross breeding population. One such relationship, regulated by a variant on chromosome 7, was the association of Odoribacter (Bacteroidales) abundance and sleep phenotypes. In a test of this association in the BKS.Cg-Dock7m +/+ Leprdb/J mouse model of obesity and diabetes, known to have abnormal sleep and colonization by Odoribacter, treatment with antibiotics altered sleep in a genotype-dependent fashion. The many other relationships extracted from this study can be used to interrogate other diseases, microbes, and mechanisms.

Prospects for finding the mechanisms of sex differences in addiction with human and model organism genetic analysis

Despite substantial evidence for sex differences in addiction epidemiology, addiction-relevant behaviors and associated neurobiological phenomena, the mechanisms and implications of these differences remain unknown. Genetic analysis in model organism is a potentially powerful and effective means of discovering the mechanisms that underlie sex differences in addiction. Human genetic studies are beginning to show precise risk variants that influence the mechanisms of addiction but typically lack sufficient power or neurobiological mechanistic access, particularly for the discovery of the mechanisms that underlie sex differences. Our thesis in this review is that genetic variation in model organisms are a promising approach that can complement these investigations to show the biological mechanisms that underlie sex differences in addiction.

Genetic mapping in Diversity Outbred mice identifies a Trpa1 variant influencing late-phase formalin response

Identification of genetic variants that influence susceptibility to pain is key to identifying molecular mechanisms and targets for effective and safe therapeutic alternatives to opioids. To identify genes and variants associated with persistent pain, we measured late-phase response to formalin injection in 275 male and female Diversity Outbred mice genotyped for over 70,000 single nucleotide polymorphisms. One quantitative trait locus reached genome-wide significance on chromosome 1 with a support interval of 3.1 Mb. This locus, Nociq4 (nociceptive sensitivity quantitative trait locus 4; MGI: 5661503), harbors the well-known pain gene Trpa1 (transient receptor potential cation channel, subfamily A, member 1). Trpa1 is a cation channel known to play an important role in acute and chronic pain in both humans and mice. Analysis of Diversity Outbred founder strain allele effects revealed a significant effect of the CAST/EiJ allele at Trpa1, with CAST/EiJ carrier mice showing an early, but not late, response to formalin relative to carriers of the 7 other inbred founder alleles (A/J, C57BL/6J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, PWK/PhJ, and WSB/EiJ). We characterized possible functional consequences of sequence variants in Trpa1 by assessing channel conductance, TRPA1-TRPV1 interactions, and isoform expression. The phenotypic differences observed in CAST/EiJ relative to C57BL/6J carriers were best explained by Trpa1 isoform expression differences, implicating a splice junction variant as the causal functional variant. This study demonstrates the utility of advanced, high-precision genetic mapping populations in resolving specific molecular mechanisms of variation in pain sensitivity.

Integration of heterogeneous functional genomics data in gerontology research to find genes and pathway underlying aging across species

Understanding the biological mechanisms behind aging, lifespan and healthspan is becoming increasingly important as the proportion of the world's population over the age of 65 grows, along with the cost and complexity of their care. BigData oriented approaches and analysis methods enable current and future bio-gerontologists to synthesize, distill and interpret vast, heterogeneous data from functional genomics studies of aging. GeneWeaver is an analysis system for integration of data that allows investigators to store, search, and analyze immense amounts of data including user-submitted experimental data, data from primary publications, and data in other databases. Aging related genome-wide gene sets from primary publications were curated into this system in concert with data from other model-organism and aging-specific databases, and applied to several questions in genrontology using. For example, we identified Cd63 as a frequently represented gene among aging-related genome-wide results. To evaluate the role of Cd63 in aging, we performed RNAi knockdown of the C. elegans ortholog, tsp-7, demonstrating that this manipulation is capable of extending lifespan. The tools in GeneWeaver enable aging researchers to make new discoveries into the associations between the genes, normal biological processes, and diseases that affect aging, healthspan, and lifespan.

On Finding and Enumerating Maximal and Maximum k-Partite Cliques in k-Partite Graphs

Let k denote an integer greater than 2, let G denote a k-partite graph, and let S denote the set of all maximal k-partite cliques in G. Several open questions concerning the computation of S are resolved. A straightforward and highly-scalable modification to the classic recursive backtracking approach of Bron and Kerbosch is first described and shown to run in O(3 ^n/3) time. A series of novel graph constructions is then used to prove that this bound is best possible in the sense that it matches an asymptotically tight upper limit on |S|. The task of identifying a vertex-maximum element of S is also considered and, in contrast with the k = 2 case, shown to be NP-hard for every k ≥ 3. A special class of k-partite graphs that arises in the context of functional genomics and other problem domains is studied as well and shown to be more readily solvable via a polynomial-time transformation to bipartite graphs. Applications, limitations, potentials for faster methods, heuristic approaches, and alternate formulations are also addressed.

QTL and systems genetics analysis of mouse grooming and behavioral responses to novelty in an open field

The open field is a classic test used to assess exploratory behavior, anxiety and locomotor activity in rodents. Here, we mapped quantitative trait loci (QTLs) underlying behaviors displayed in an open field, using a panel of 53 BXD recombinant inbred mouse strains with deep replication (10 per strain and sex). The use of these strains permits the integration and comparison of data obtained in different laboratories, and also offers the possibility to study trait covariance by exploiting powerful bioinformatics tools and resources. We quantified behavioral traits during 20-min test sessions including (1) percent time spent and distance traveled near the wall (thigmotaxis), (2) leaning against the wall, (3) rearing, (4) jumping, (5) grooming duration, (6) grooming frequency, (7) locomotion and (8) defecation. All traits exhibit moderate heritability making them amenable to genetic analysis. We identified a significant QTL on chromosome M.m. 4 at approximately 104 Mb that modulates grooming duration in both males and females (likelihood ratio statistic values of approximately 18, explaining 25% and 14% of the variance, respectively) and a suggestive QTL modulating locomotion that maps to the same locus. Bioinformatic analysis indicates Disabled 1 (Dab1, a key protein in the reelin signaling pathway) as a particularly strong candidate gene modulating these behaviors. We also found 2 highly suggestive QTLs for a sex by strain interaction for grooming duration on chromosomes 13 and 17. In addition, we identified a pairwise epistatic interaction between loci on chromosomes 12 at 36-37 Mb and 14 at 34-36 Mb that influences rearing frequency in males.

Cross-Species Integrative Functional Genomics in GeneWeaver Reveals a Role for Pafah1b1 in Altered Response to Alcohol

Identifying the biological substrates of complex neurobehavioral traits such as alcohol dependency pose a tremendous challenge given the diverse model systems and phenotypic assessments used. To address this problem we have developed a platform for integrated analysis of high-throughput or genome-wide functional genomics studies. A wealth of such data exists, but it is often found in disparate, non-computable forms. Our interactive web-based software system, Gene Weaver (http://www.geneweaver.org), couples curated results from genomic studies to graph-theoretical tools for combinatorial analysis. Using this system we identified a gene underlying multiple alcohol-related phenotypes in four species. A search of over 60,000 gene sets in GeneWeaver's database revealed alcohol-related experimental results including genes identified in mouse genetic mapping studies, alcohol selected Drosophila lines, Rattus differential expression, and human alcoholic brains. We identified highly connected genes and compared these to genes currently annotated to alcohol-related behaviors and processes. The most highly connected gene not annotated to alcohol was Pafah1b1. Experimental validation using a Pafah1b1 conditional knock-out mouse confirmed that this gene is associated with an increased preference for alcohol and an altered thermoregulatory response to alcohol. Although this gene has not been previously implicated in alcohol-related behaviors, its function in various neural mechanisms makes a role in alcohol-related phenomena plausible. By making diverse cross-species functional genomics data readily computable, we were able to identify and confirm a novel alcohol-related gene that may have implications for alcohol use disorders and other effects of alcohol.

GeneWeaver: data driven alignment of cross-species genomics in biology and disease

The GeneWeaver data and analytics website (www.geneweaver.org) is a publically available resource for storing, curating and analyzing sets of genes from heterogeneous data sources. The system enables discovery of relationships among genes, variants, traits, drugs, environments, anatomical structures and diseases implicitly found through gene set intersections. Since the previous review in the 2012 Nucleic Acids Research Database issue, GeneWeaver's underlying analytics platform has been enhanced, its number and variety of publically available gene set data sources has been increased, and its advanced search mechanisms have been expanded. In addition, its interface has been redesigned to take advantage of flexible web services, programmatic data access, and a refined data model for handling gene network data in addition to its original emphasis on gene set data. By enumerating the common and distinct biological molecules associated with all subsets of curated or user submitted groups of gene sets and gene networks, GeneWeaver empowers users with the ability to construct data driven descriptions of shared and unique biological processes, diseases and traits within and across species.

GeneWeaver: finding consilience in heterogeneous cross-species functional genomics data

A persistent challenge lies in the interpretation of consensus and discord from functional genomics experimentation. Harmonizing and analyzing this data will enable investigators to discover relations of many genes to many diseases, and from many phenotypes and experimental paradigms to many diseases through their genomic substrates. The GeneWeaver.org system provides a platform for cross-species integration and interrogation of heterogeneous curated and experimentally derived functional genomics data. GeneWeaver enables researchers to store, share, analyze, and compare results of their own genome-wide functional genomics experiments in an environment containing rich companion data obtained from major curated repositories, including the Mouse Genome Database and other model organism databases, along with derived data from highly specialized resources, publications, and user submissions. The data, largely consisting of gene sets and putative biological networks, are mapped onto one another through gene identifiers and homology across species. A versatile suite of interactive tools enables investigators to perform a variety of set analysis operations to find consilience among these often noisy experimental results. Fast algorithms enable real-time analysis of large queries. Specific applications include prioritizing candidate genes for quantitative trait loci, identifying biologically valid mouse models and phenotypic assays for human disease, finding the common biological substrates of related diseases, classifying experiments and the biological concepts they represent from empirical data, and applying patterns of genomic evidence to implicate novel genes in disease. These results illustrate an alternative to strict emphasis on replicability, whereby researchers classify experimental results to identify the conditions that lead to their similarity.

Molecular identification of collagen 17a1 as a major genetic modifier of laminin gamma 2 mutation-induced junctional epidermolysis bullosa in mice

Epidermolysis Bullosa (EB) encompasses a spectrum of mechanobullous disorders caused by rare mutations that result in structural weakening of the skin and mucous membranes. While gene mutated and types of mutations present are broadly predictive of the range of disease to be expected, a remarkable amount of phenotypic variability remains unaccounted for in all but the most deleterious cases. This unexplained variance raises the possibility of genetic modifier effects. We tested this hypothesis using a mouse model that recapitulates a non-Herlitz form of junctional EB (JEB) owing to the hypomorphic jeb allele of laminin gamma 2 (Lamc2). By varying normally asymptomatic background genetics, we document the potent impact of genetic modifiers on the strength of dermal-epidermal adhesion and on the clinical severity of JEB in the context of the Lamc2(jeb) mutation. Through an unbiased genetic approach involving a combination of QTL mapping and positional cloning, we demonstrate that Col17a1 is a strong genetic modifier of the non-Herlitz JEB that develops in Lamc2(jeb) mice. This modifier is defined by variations in 1-3 neighboring amino acids in the non-collagenous 4 domain of the collagen XVII protein. These allelic variants alter the strength of dermal-epidermal adhesion in the context of the Lamc2(jeb) mutation and, consequentially, broadly impact the clinical severity of JEB. Overall the results provide an explanation for how normally innocuous allelic variants can act epistatically with a disease causing mutation to impact the severity of a rare, heritable mechanobullous disorder.

IL-21 is a double-edged sword in the systemic lupus erythematosus-like disease of BXSB.Yaa mice

The pleiotropic cytokine IL-21 is implicated in the pathogenesis of human systemic lupus erythematosus by polymorphisms in the molecule and its receptor (IL-21R). The systemic lupus erythematosus-like autoimmune disease of BXSB.Yaa mice is critically dependent on IL-21 signaling, providing a model for understanding IL-21/IL-21R signaling in lupus pathogenesis. In this study, we generated BXSB.Yaa mice selectively deficient in IL-21R on B cells, on all T cells, or on CD8(+) T cells alone and examined the effects on disease. We found that IL-21 signaling to B cells is essential for the development of all classical disease manifestations, but that IL-21 signaling also supports the expansion of central memory, CD8(+) suppressor cells and broadly represses the cytokine activity of CD4(+) T cells. These results indicate that IL-21 has both disease-promoting and disease-suppressive effects in the autoimmune disease of BXSB.Yaa mice.

High-precision genetic mapping of behavioral traits in the diversity outbred mouse population

Historically our ability to identify genetic variants underlying complex behavioral traits in mice has been limited by low mapping resolution of conventional mouse crosses. The newly developed Diversity Outbred (DO) population promises to deliver improved resolution that will circumvent costly fine-mapping studies. The DO is derived from the same founder strains as the Collaborative Cross (CC), including three wild-derived strains. Thus the DO provides more allelic diversity and greater potential for discovery compared to crosses involving standard mouse strains. We have characterized 283 male and female DO mice using open-field, light–dark box, tail-suspension and visual-cliff avoidance tests to generate 38 behavioral measures. We identified several quantitative trait loci (QTL) for these traits with support intervals ranging from 1 to 3 Mb in size. These intervals contain relatively few genes (ranging from 5 to 96). For a majority of QTL, using the founder allelic effects together with whole genome sequence data, we could further narrow the positional candidates. Several QTL replicate previously published loci. Novel loci were also identified for anxiety- and activity-related traits. Half of the QTLs are associated with wild-derived alleles, confirming the value to behavioral genetics of added genetic diversity in the DO. In the presence of wild-alleles we sometimes observe behaviors that are qualitatively different from the expected response. Our results demonstrate that high-precision mapping of behavioral traits can be achieved with moderate numbers of DO animals, representing a significant advance in our ability to leverage the mouse as a tool for behavioral genetics

Accelerating discovery for complex neurological and behavioral disorders through systems genetics and integrative genomics in the laboratory mouse

Recent advances in systems genetics and integrative functional genomics have greatly improved the study of complex neurological and behavioral traits. The methods developed for the integrated characterization of new, high-resolution mouse genetic reference populations and systems genetics enable behavioral geneticists an unprecedented opportunity to address questions of the molecular basis of neurological and psychiatric disorders and their comorbidities. Integrative genomics augment these strategies by enabling rapid informatics-assisted candidate gene prioritization, cross-species translation, and mechanistic comparison across related disorders from a wealth of existing data in mouse and other model organisms. Ultimately, through these complementary approaches, finding the mechanisms and sources of genetic variation underlying complex neurobehavioral disease related traits is becoming tractable. Furthermore, these methods enable categorization of neurobehavioral disorders through their underlying biological basis. Together, these model organism-based approaches can lead to a refinement of diagnostic categories and targeted treatment of neurological and psychiatric disease.

Development, structure, and keratin expression in C57BL/6J mouse eccrine glands

Eccrine sweat glands in the mouse are found only on the footpads and, when mature, resemble human eccrine glands. Eccrine gland anlagen were first apparent at 16.5 days postconception (DPC) in mouse embryos as small accumulations of cells in the mesenchymal tissue beneath the developing epidermis resembling hair follicle placodes. These cells extended into the dermis where significant cell organization, duct development, and evidence of the acrosyringium were observed in 6- to 7-postpartum day (PPD) mice. Mouse-specific keratin 1 (K1) and 10 (K10) expression was confined to the strata spinosum and granulosum. In 16.5 and 18.5 DPC embryos, K14 and K17 were both expressed in the stratum basale and diffusely in the gland anlagen. K5 expression closely mimicked K17 throughout gland development. K6 expression was not observed in the developing glands of the embryo but was apparent in the luminal cell layer of the duct by 6 to 7 PPD. By 21 PPD, the gland apertures appeared as depressions in the surface surrounded by cornified squames, and the footpad surface lacked the organized ridge and crease system seen in human fingers. These data serve as a valuable reference for investigators who use genetically engineered mice for skin research.

GeneWeaver: a web-based system for integrative functional genomics

High-throughput genome technologies have produced a wealth of data on the association of genes and gene products to biological functions. Investigators have discovered value in combining their experimental results with published genome-wide association studies, quantitative trait locus, microarray, RNA-sequencing and mutant phenotyping studies to identify gene-function associations across diverse experiments, species, conditions, behaviors or biological processes. These experimental results are typically derived from disparate data repositories, publication supplements or reconstructions from primary data stores. This leaves bench biologists with the complex and unscalable task of integrating data by identifying and gathering relevant studies, reanalyzing primary data, unifying gene identifiers and applying ad hoc computational analysis to the integrated set. The freely available GeneWeaver (http://www.GeneWeaver.org) powered by the Ontological Discovery Environment is a curated repository of genomic experimental results with an accompanying tool set for dynamic integration of these data sets, enabling users to interactively address questions about sets of biological functions and their relations to sets of genes. Thus, large numbers of independently published genomic results can be organized into new conceptual frameworks driven by the underlying, inferred biological relationships rather than a pre-existing semantic framework. An empirical ‘ontology’ is discovered from the aggregate of experimental knowledge around user-defined areas of biological inquiry.

MHC class I family proteins retard systemic lupus erythematosus autoimmunity and B cell lymphomagenesis

Dysregulation of the T cell-dependent Ab response can lead to numerous immunological disorders, ranging from systemic lupus erythematosus to B cell lymphomas. Cellular processes governed by MHC class II proteins play a major role in this response and its dysregulation. The extent to which processes controlled by the diverse family of MHC class I proteins impact such autoimmune and neoplastic disorders, however, is less clear. In this study, we genetically dissect the contributions of individual MHC class I family members and the pathological processes under their control in the systemic lupus erythematosus-like disease of BXSB.Yaa mice and B cell lymphomagenesis of SJL mice. This study reveals a powerful repressive regulatory axis comprised of MHC class I-dependent CD8(+) T cells and NK cells. These results indicate that the predominant role of the MHC class I protein family in such immunological disorders is to protect from more aggressive diseases.

Molecular mechanisms of induction of antigen-specific allograft tolerance by intranasal peptide administration

We have previously shown that intranasal (i.n.) administration of a single MHC class II-restricted HY peptide to female mice induces tolerance to up to five additional epitopes expressed on test male grafts, a phenomenon known as linked suppression. In this study, we investigated the molecular mechanisms involved both in the induction phase following peptide administration and during linked suppression after grafting. We report that following initial i.n. administration, peptide is widely disseminated and is presented by functionally immature dendritic cells. These fail to cause optimal stimulation of the responding HY-specific CD4(+) T cells that express genes characteristic of regulatory T cells. Following i.n. peptide plus LPS administration, causing immunization, HY-specific CD4(+) T cells express genes characteristic of activated T cells. We further find that following male skin grafting, HY-specific CD8(+) T cells from peptide-treated tolerant mice display both quantitative and qualitative differences compared with similar cells from untreated mice that reject their grafts. In tolerant mice there are fewer HY-specific CD8(+) cells and they express several genes characteristic of exhausted T cells. Furthermore, associated with specific chemokine receptor and integrin expression, HY-specific CD8(+) T cells show more limited migration from the graft draining lymph node into other tissues.

IL-21 Receptor signaling by B-cells is essential for BXSB-Yaa pathogenesis (99.39)

IL-21 is a member of the IL-2 cytokine family that exerts potent biological effects on multiple hemopoietic cell types. It is produced by activated CD4 T cells, while its receptor, IL-21R, is found on cells including NK cells, dendritic, T cells and B cells. To investigate the role of IL-21 signaling in autoimmune disease, we addressed whether this cytokine contributes to the severe systemic lupus erythematosus (SLE)-like syndrome developed by BXSB-Yaa mice. We show that BXSB-Yaa CD4 T cells produce IL-21 but lack discriminative T follicular helper cell markers. We also show that BXSB-Yaa mice completely lacking the IL-21R or Igh6 are remarkably resistant to multiple aspects of the prototypic BXSB-Yaa autoimmune disease. Finally, we show that mixed bone marrow chimeric BXSB-Yaa mice lacking IL-21R specifically in the B-cell compartment are similarly resistant. These results suggest that IL-21 produced by a novel population of BXSB-Yaa CD4 T-cells acts through B-cells to evoke multiple features characteristic of severe SLE. The results indicate that IL-21 signaling through B-cells is an essential step in the pathogenesis of this SLE-like autoimmune disease and potentially others. Funded by an Arthritis Foundation grant to JB.

Cataracts in transgenic mice caused by a human papillomavirus type 18 E7 oncogene driven by KRT1-14

Human papillomavirus type 18 (HPV18) is a common cause of cervical cancer. To create a mouse model for this common neoplastic disease, we used a human keratin 14 promoter to drive the HPV18 E7 oncogene to create transgenic mice. No mice up to a year of age developed cervical cancer. However, all transgenic mice and none of the controls developed progressive bilateral cortical cataracts. By 6 months of age, the cortex liquefied leaving the lens nucleus. Proliferation of lens epithelium formed multifocal nodules and free floating lens epithelial cells within the liquefied cortex. These cells were hyperplastic not neoplastic. Other HPV transgenic stocks develop cataracts suggesting this virus may have a broad cellular tropism.

Treatment of BXSB-Yaa Mice with IL-21R-Fc Fusion Protein Minimally Attenuates Systemic Lupus Erythematosus

Interleukin-21 (IL-21) is a pleiotropic cytokine whose function is only now being unraveled. Abundant evidence indicates that activated CD4 T cells are the primary, if not the only, source of IL-21. While it is clear that IL-21 is actively transcribed by naïve activated T cells, recent studies have shown that IL-21 potentially promotes a developmental shift of naïve T cells toward the Th2 phenotype. BXSB-Yaa mice develop an autoimmune syndrome similar to systemic lupus erythematosus (SLE), affecting males earlier than females on account of the presence of the Yaa (Y-linked autoimmune acceleration) locus. Previous results indicate the elevation of IL-21 expression by BXSB-Yaa mice at an age when the early characteristics of autoimmune processes first become evident. We set out to determine whether IL-21 was necessary for disease progression in BXSB-Yaa mice. Mice were treated for 24 weeks with soluble IL-21R-Fc in order to therapeutically neutralize the IL-21 present. The results overall suggest a biphasic effect of IL-21, negatively influencing survival early on and positively influencing survival at later stages. We propose that IL-21 exerts a pleiotropic effect in which it promotes the protective effects of CD8+ suppressor cells in the early disease phase and then promotes the humoral components of SLE in the later disease stages. This experiment provides preliminary evidence for a role of IL-21 in modulating the severity of SLE in BXSB-Yaa mice.

Aged B6.SB-Yaa/J Mice develop an SLE-like Disease (130.11)

Mouse models of SLE are considered pivotal in understanding the genetics of SLE along with identifying targets for therapy. Male mice of the BXSB/MpJ recombinant inbred strain, first described by Murphy and Roths in 1978, develop a spontaneous and severe SLE-like syndrome resulting in death averaging ~6 months of age, while females show only a very chronic syndrome. A mutant locus found on the Y chromosome is, therefore, appropriately named the Y-linked autoimmune accelerator (Yaa). Previous analysis have suggested that Yaa fails to induce an SLE-like disease in B6 mice, the caveat to this being that the SLE criteria analyzed were limited to serological analysis, histopathology, and mortality on mice up to 8–12 months of age. We report here that when mice of a significantly older age, 15–24 months, are compared to age-matched B6 controls, the B6.SB-Yaa/J mice clearly developed an SLE-like disease. This is characterized by significant splenomegaly, monocytosis, immune cells with an activation phenotype, the presence of homogeneous anti-nuclear antibody staining, glomerulonephritis, decreased weight and survival as well as a gene expression signature similar to BXSB/MpJ mice. These results suggest that Yaa alone is sufficient to induce an SLE-like syndrome in aged B6 mice.

Funded by NIH R21DK074463 and DK56597 to DCR, Arthritis Foundation to JAB.

The Yaa locus Confers Apoptosis, NF-κB and Rapamycin Signaling Resistance Through Multiple Receptors in C57BL/6J Mice (129.12)

BXSB/MpJ mice develop a spontaneous SLE-like disease that is remarkably accelerated in males, due to the presence of the Y-linked-autoimmune accelerator locus (Yaa). Recent insights into the molecular nature of the Yaa mutation were provided when an additional copy of 17 distally located X chromosome genes were discovered on the Y-chromosome of Yaa mice. Included within this region is the gene TLR7, a toll receptor which binds ssRNA and is functionally expressed on B-cells and monocyte derived lineages. Previous reports have suggested that the SLE-like disease phenotype conveyed by Yaa is due to the duplication of TLR7, which causes a shift toward anti-nucleolar antigens in B6.FcγRIIB(−/−). We show here that B6, BXSB and MRL mice that carry Yaa rarely display the anti-nucleolar antibody pattern. We found that B cells of Yaa mice display novel phenotypes consistent with increased rapamycin and NF-κB signaling through receptors other than TLR7, including TLR4, TLR9, and the BCR. In addition, Yaa B cells are resistant to anti-IgM induced apoptosis. These results demonstrate

that the copy number of a single receptor can influence the strength of signaling through other receptors and/orthat additional genes within the duplicated region contribute toward the Yaa phenotype.

Funded by NIH R21DK074463 and DK56597 to DCR, Arthritis Foundation to JAB.