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.

Developmental coordination disorder: What can we learn from RI mice using motor learning tasks and QTL analysis

Developmental Coordination Disorder (DCD) is a neurodevelopmental disorder of unknown etiology that affects one in 20 children. There is an indication that DCD has an underlying genetic component due to its high heritability. Therefore, we explored the use of a recombinant inbred family of mice known as the BXD panel to understand the genetic basis of complex traits (i.e., motor learning) through identification of quantitative trait loci (QTLs). The overall aim of this study was to utilize the QTL approach to evaluate the genome-to-phenome correlation in BXD strains of mice in order to better understand the human presentation of DCD. Results of this current study confirm differences in motor learning in selected BXD strains and strains with altered cerebellar volume. Five strains - BXD15, BXD27, BXD28, BXD75, and BXD86 - exhibited the most DCD-like phenotype when compared with other BXD strains of interest. Results indicate that BXD15 and BXD75 struggled primarily with gross motor skills, BXD28 primarily had difficulties with fine motor skills, and BXD27 and BXD86 strains struggled with both fine and gross motor skills. The functional roles of genes within significant QTLs were assessed in relation to DCD-like behavior. Only Rab3a (Ras-related protein Rab-3A) emerged as a high likelihood candidate gene for the horizontal ladder rung task. This gene is associated with brain and skeletal muscle development, but lacked nonsynonymous polymorphisms. This study along with Gill et al. (same issue) is the first studies to specifically examine the genetic linkage of DCD using BXD strains of mice.

Increased development of T-bet+CD11c+ B cells predisposes to lupus in females: Analysis in BXD2 mouse and genetic crosses

Cardinal features of lupus include elevated B cell activation and autoantibody production with a female sex preponderance. We quantified interactions of sex and genetic variation on the development of autoimmune B-cell phenotypes and autoantibodies in the BXD2 murine model of lupus using a cohort of backcrossed progeny (BXD2 x C57BL/6J) x BXD2. Sex was the key factor leading to increased total IgG, IgG2b, and autoantibodies. The percentage of T-bet+CD11c+ IgD+ activated naive B cells (aNAV) was higher in females and was associated with increased T-bet+CD11c+ IgD- age-related B cells, Fas+GL7+ germinal center B cells, Cxcr5-Icos+ peripheral T-helper cells, and Cxcr5+Icos+ follicular T-helper cells. IFN-β was elevated in females. Variation in aNAV cells was mapped to Chr 7 in a locus that shows significant interactions between the female sex and heterozygous B/D variant. Our results suggest that activation of naive B cells forms the basis for the female-predominant development of autoantibodies in lupus-susceptible BXD2 mice.

Polyglucosan body density in the aged mouse hippocampus is controlled by a novel modifier locus on chromosome 1

Aging can be associated with the accumulation of hypobranched glycogen molecules (polyglucosan bodies, PGBs), particularly in astrocytes of the hippocampus. While PGBs have a detrimental effect on cognition in diseases such as adult polyglucosan body disease and Lafora disease, the underlying mechanism and clinical relevance of age-related PGB accumulation remains unknown. Here, we have investigated the genetic basis and functional impact of age-related PGB accumulation in 32 fully sequenced BXD-type strains of mice which exhibit a 400-fold variation in PGB burden in 16-18 month old females. We mapped a major locus controlling PGB density in the hippocampus to chromosome 1 at 72-75 Mb (linkage of 4.9 -logP), which we defined as the Pgb1 locus. To identify potentially causal gene variants within Pgb1, we generated extensive hippocampal transcriptome datasets and identified two strong candidate genes for which mRNA correlates with PGB density-Smarcal1 and Usp37. In addition, both Smarcal1 and Usp37 contain non-synonymous allele variations likely to impact protein function. A phenome-wide association analysis highlighted a trans-regulatory effect of the Pgb1 locus on expression of Hp1bp3, a gene known to play a role in age-related changes in learning and memory. To investigate the potential impact of PGBs on cognition, we performed conditioned fear memory testing on strains displaying varying degrees of PGB burden, and a phenome-wide association scan of ~12,000 traits. Importantly, we did not find any evidence suggesting a negative impact of PGB burden on cognitive capacity. Taken together, we have identified a major modifier locus controlling PGB burden in the hippocampus and shed light on the genetic architecture and clinical relevance of this strikingly heterogeneous hippocampal phenotype.

Genome-wide screen identifies host loci that modulate Mycobacterium tuberculosis fitness in immunodivergent mice

Genetic differences among mammalian hosts and among strains of Mycobacterium tuberculosis (Mtb) are well-established determinants of tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host-pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the highly diverse BXD family of strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb-resistant C57BL/6J (B6 or B) and Mtb-susceptible DBA/2J (D2 or D) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters of "endophenotypes," each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted quantitative trait loci (QTL) mapping of these bacterial fitness endophenotypes and identified 140 host-pathogen QTL (hpQTL). We located a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes: Rv0127 (mak), Rv0359 (rip2), Rv0955 (perM), and Rv3849 (espR). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Expression Levels of the Tnni3k Gene in the Heart Are Highly Associated with Cardiac and Glucose Metabolism-Related Phenotypes and Functional Pathways

BACKGROUND

Troponin-I interacting kinase encoded by the TNNI3K gene is expressed in nuclei and Z-discs of cardiomyocytes. Mutations in TNNI3K were identified in patients with cardiac conduction diseases, arrhythmias, and cardiomyopathy.

METHODS

We performed cardiac gene expression, whole genome sequencing (WGS), and cardiac function analysis in 40 strains of BXD recombinant inbred mice derived from C57BL/6J (B6) and DBA/2J (D2) strains. Expression quantitative trait loci (eQTLs) mapping and gene enrichment analysis was performed, followed by validation of candidate Tnni3k-regulatory genes.

RESULTS

WGS identified compound splicing and missense T659I Tnni3k variants in the D2 parent and some BXD strains (D allele) and these strains had significantly lower Tnni3k expression than those carrying wild-type Tnni3k (B allele). Expression levels of Tnni3k significantly correlated with multiple cardiac (heart rate, wall thickness, PR duration, and T amplitude) and metabolic (glucose levels and insulin resistance) phenotypes in BXDs. A significant cis-eQTL on chromosome 3 was identified for the regulation of Tnni3k expression. Furthermore, Tnni3k-correlated genes were primarily involved in cardiac and glucose metabolism-related functions and pathways. Genes Nodal, Gnas, Nfkb1, Bmpr2, Bmp7, Smad7, Acvr1b, Acvr2b, Chrd, Tgfb3, Irs1, and Ppp1cb were differentially expressed between the B and D alleles.

CONCLUSIONS

Compound splicing and T659I Tnni3k variants reduce cardiac Tnni3k expression and Tnni3k levels are associated with cardiac and glucose metabolism-related phenotypes.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize 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 the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. 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.

Epigenetic analysis in a murine genetic model of Gulf War illness

Of the nearly 1 million military personnel who participated in the 1990-1991 Gulf War, between 25% and 35% became ill with what now is referred to as Gulf War Illness (GWI) by the Department of Defense. Symptoms varied from gastrointestinal distress to lethargy, memory loss, inability to concentrate, depression, respiratory, and reproductive problems. The symptoms have persisted for 30 years in those afflicted but the basis of the illness remains largely unknown. Nerve agents and other chemical exposures in the war zone have been implicated but the long-term effects of these acute exposures have left few if any identifiable signatures. The major aim of this study is to elucidate the possible genomic basis for the persistence of symptoms, especially of the neurological and behavioral effects. To address this, we performed a whole genome epigenetic analysis of the proposed cause of GWI, viz., exposure to organophosphate neurotoxicants combined with high circulating glucocorticoids in two inbred mouse strains, C57BL/6J and DBA/2J. The animals received corticosterone in their drinking water for 7 days followed by injection of diisopropylfluorophosphate, a nerve agent surrogate. Six weeks after DFP injection, the animals were euthanized and medial prefrontal cortex harvested for genome-wide DNA methylation analysis using high-throughput sequencing. We observed 67 differentially methylated genes, notably among them, Ttll7, Akr1c14, Slc44a4, and Rusc2, all related to different symptoms of GWI. Our results support proof of principle of genetic differences in the chronic effects of GWI-related exposures and may reveal why the disease has persisted in many of the now aging Gulf War veterans.

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Merged magnetic resonance and light sheet microscopy of the whole mouse brain

We have developed workflows to align 3D magnetic resonance histology (MRH) of the mouse brain with light sheet microscopy (LSM) and 3D delineations of the same specimen. We start with MRH of the brain in the skull with gradient echo and diffusion tensor imaging (DTI) at 15 μm isotropic resolution which is ~ 1,000 times higher than that of most preclinical MRI. Connectomes are generated with superresolution tract density images of ~5 μm. Brains are cleared, stained for selected proteins, and imaged by LSM at 1.8 μm/pixel. LSM data are registered into the reference MRH space with labels derived from the ABA common coordinate framework. The result is a high-dimensional integrated volume with registration (HiDiver) with alignment precision better than 50 µm. Throughput is sufficiently high that HiDiver is being used in quantitative studies of the impact of gene variants and aging on mouse brain cytoarchitecture and connectomics.

Building pangenome graphs

Pangenome graphs can represent all variation between multiple genomes, but existing methods for constructing them are biased due to reference-guided approaches. In response, we have developed PanGenome Graph Builder (PGGB), a reference-free pipeline for constructing unbi-ased pangenome graphs. PGGB uses all-to-all whole-genome alignments and learned graph embeddings to build and iteratively refine a model in which we can identify variation, measure conservation, detect recombination events, and infer phylogenetic relationships.

Genome-wide screen identifies host loci that modulate M. tuberculosis fitness in immunodivergent mice

Genetic differences among mammalian hosts and Mycobacterium tuberculosis ( Mtb ) strains determine diverse tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host- pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the BXD family of mouse strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb -resistant C57BL/6J (B6 or B ) and Mtb -susceptible DBA/2J (D2 or D ) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters for "endophenotypes", each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted QTL mapping of these bacterial fitness endophenotypes and identified 140 h ost- p athogen quantitative trait loci ( hp QTL). We identified a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes; Rv0127 ( mak ), Rv0359 ( rip2 ), Rv0955 ( perM ), and Rv3849 ( espR ). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Weighted gene co-expression network analysis identifies key hub genes and pathways in acute myeloid leukemia

Introduction: Acute myeloid leukemia (AML) is the most common type of leukemia in adults. However, there is a gap in understanding the molecular basis of the disease, partly because key genes associated with AML have not been extensively explored. In the current study, we aimed to identify genes that have strong association with AML based on a cross-species integrative approach. Methods: We used Weighted Gene Co-Expression Network Analysis (WGCNA) to identify co-expressed gene modules significantly correlated with human AML, and further selected the genes exhibiting a significant difference in expression between AML and healthy mouse. Protein-protein interactions, transcription factors, gene function, genetic regulation, and coding sequence variants were integrated to identify key hub genes in AML. Results: The cross-species approach identified a total of 412 genes associated with both human and mouse AML. Enrichment analysis confirmed an association of these genes with hematopoietic and immune-related functions, phenotypes, processes, and pathways. Further, the integrated analysis approach identified a set of important module genes including Nfe2, Trim27, Mef2c, Ets1, Tal1, Foxo1, and Gata1 in AML. Six of these genes (except ETS1) showed significant differential expression between human AML and healthy samples in an independent microarray dataset. All of these genes are known to be involved in immune/hematopoietic functions, and in transcriptional regulation. In addition, Nfe2, Trim27, Mef2c, and Ets1 harbor coding sequence variants, whereas Nfe2 and Trim27 are cis-regulated, making them attractive candidates for validation. Furthermore, subtype-specific analysis of the hub genes in human AML indicated high expression of NFE2 across all the subtypes (M0 through M7) and enriched expression of ETS1, LEF1, GATA1, and TAL1 in M6 and M7 subtypes. A significant correlation between methylation status and expression level was observed for most of these genes in AML patients. Conclusion: Findings from the current study highlight the importance of our cross-species approach in the identification of multiple key candidate genes in AML, which can be further studied to explore their detailed role in leukemia/AML.

Genomic Basis for Individual Differences in Susceptibility to the Neurotoxic Effects of Diesel Exhaust

Air pollution is a known environmental health hazard. A major source of air pollution includes diesel exhaust (DE). Initially, research on DE focused on respiratory morbidities; however, more recently, exposures to DE have been associated with neurological developmental disorders and neurodegeneration. In this study, we investigated the effects of sub-chronic inhalation exposure to DE on neuroinflammatory markers in two inbred mouse strains and both sexes, including whole transcriptome examination of the medial prefrontal cortex. We exposed aged male and female C57BL/6J (B6) and DBA/2J (D2) mice to DE, which was cooled and diluted with HEPA-filtered compressed air for 2 h per day, 5 days a week, for 4 weeks. Control animals were exposed to HEPA-filtered air on the same schedule as DE-exposed animals. The prefrontal cortex was harvested and analyzed for proinflammatory cytokine gene expression (Il1β, Il6, Tnfα) and transcriptome-wide response by RNA-seq. We observed differential cytokine gene expression between strains and sexes in the DE-exposed vs. control-exposed groups for Il1β, Tnfα, and Il6. For RNA-seq, we identified 150 differentially expressed genes between air and DE treatment related to natural killer cell-mediated cytotoxicity per Kyoto Encyclopedia of Genes and Genomes pathways. Overall, our data show differential strain-related effects of DE on neuroinflammation and neurotoxicity and demonstrate that B6 are more susceptible than D2 to gene expression changes due to DE exposures than D2. These results are important because B6 mice are often used as the default mouse model for DE studies and strain-related effects of DE neurotoxicity warrant expanded studies.

Identification of cyclin D1 as a major modulator of 3-nitropropionic acid-induced striatal neurodegeneration

Mitochondria dysfunction occurs in the aging brain as well as in several neurodegenerative disorders and predisposes neuronal cells to enhanced sensitivity to neurotoxins. 3-nitropropionic acid (3-NP) is a naturally occurring plant and fungal neurotoxin that causes neurodegeneration predominantly in the striatum by irreversibly inhibiting the tricarboxylic acid respiratory chain enzyme, succinate dehydrogenase (SDH), the main constituent of the mitochondria respiratory chain complex II. Significantly, although 3-NP-induced inhibition of SDH occurs in all brain regions, neurodegeneration occurs primarily and almost exclusively in the striatum for reasons still not understood. In rodents, 3-NP-induced striatal neurodegeneration depends on the strain background suggesting that genetic differences among genotypes modulate toxicant variability and mechanisms that underlie 3-NP-induced neuronal cell death. Using the large BXD family of recombinant inbred (RI) strains we demonstrate that variants in Ccnd1 - the gene encoding cyclin D1 - of the DBA/2 J parent underlie the resistance to 3-NP-induced striatal neurodegeneration. In contrast, the Ccnd1 variant inherited from the widely used C57BL/6 J parental strain confers sensitivity. Given that cellular stress triggers induction of cyclin D1 expression followed by cell-cycle re-entry and consequent neuronal cell death, we sought to determine if the C57BL/6 J and DBA/2 J Ccnd1 variants are differentially modulated in response to 3-NP. We confirm that 3-NP induces cyclin D1 expression in striatal neuronal cells of C57BL/6 J, but this response is blunted in the DBA/2 J. We further show that striatal-specific alternative processing of a highly conserved 3'UTR negative regulatory region of Ccnd1 co-segregates with the C57BL/6 J parental Ccnd1 allele in BXD strains and that its differential processing accounts for sensitivity or resistance to 3-NP. Our results indicate that naturally occurring Ccnd1 variants may play a role in the variability observed in neurodegenerative disorders involving mitochondria complex II dysfunction and point to cyclin D1 as a possible therapeutic target.

Gene-by-environment modulation of lifespan and weight gain in the murine BXD family

How lifespan and body weight vary as a function of diet and genetic differences is not well understood. Here we quantify the impact of differences in diet on lifespan in a genetically diverse family of female mice, split into matched isogenic cohorts fed a low-fat chow diet (CD, n = 663) or a high-fat diet (HFD, n = 685). We further generate key metabolic data in a parallel cohort euthanized at four time points. HFD feeding shortens lifespan by 12%: equivalent to a decade in humans. Initial body weight and early weight gains account for longevity differences of roughly 4-6 days per gram. At 500 days, animals on a HFD typically gain four times as much weight as control, but variation in weight gain does not correlate with lifespan. Classic serum metabolites, often regarded as health biomarkers, are not necessarily strong predictors of longevity. Our data indicate that responses to a HFD are substantially modulated by gene-by-environment interactions, highlighting the importance of genetic variation in making accurate individualized dietary recommendations.

Abstract 2919: Novel pre-clinical model to identify genetic modifiers of triple negative breast cancer

Triple negative breast cancer (TNBC) is an aggressive breast cancer subtype with poor outcomes. This is a grave clinical challenge for the ~30,000 patients diagnosed with this disease every year. Discovering genetic modifiers of differential TNBC vulnerability and disease progression is critical to improving predictive and personalized treatments. We hypothesized that using a well-established recombinant inbred strain, novel genetic modifiers of TNBC risk and aggression will be identified. The C3(1)-T antigen (C3Tag) genetically engineered mouse model (GEMM) recapitulates many facets of human basal-like TNBC to demonstrate promoting effects of exposures on tumor phenotypes. However, GEMM are highly constrained by their inbred genotype and do not allow a robust interrogation of the manner in which individual genetic variation might impact tumor initiation, progression, and response to therapy. Therefore, we developed a novel murine model of TNBC in the background of the largest and best characterized genetic reference population. Systems genetics is used to identify gene candidates. Cross-species comparison of our findings with publicly available human GWAS and genomic databases is an effective approach to validate conserved biologically relevant and targetable pathways. To our knowledge, this is the first study to explore modifier genes for TNBC phenotypes using a systems genetics approach in a GEMM for TNBC. Our results will contribute to significant advances in understanding risk and improving outcomes for breast cancer.

Citation Format: Laura M. Sipe, Emily B. Korba, Lu Lu, Robert W. Williams, David G. Ashbrook, Liza Makowski. Novel pre-clinical model to identify genetic modifiers of triple negative breast cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2919.

The role of interindividual licking received and dopamine genotype on later-life licking provisioning in female rat offspring

INTRODUCTION

Rat mothers exhibit natural variations in care that propagate between generations of female offspring. However, there is limited information on genetic variation that could influence this propagation.

METHODS

We assessed early-life maternal care received by individual female rat offspring, later-life maternal care provisioning, and dopaminergic activity in the maternal brain in relation to naturally occurring genetic polymorphisms linked to the dopaminergic system. We also conducted a systematic analysis of other genetic variants potentially related to maternal behavior in our Long-Evans rat population.

RESULTS

While we did not find a direct relationship between early-life licking received and later-life licking provisioning, this relationship was indirectly affected by dopamine levels in the nucleus accumbens and dependent on variation in the dopamine receptor 2 gene (rs107017253). More specifically, female rat offspring with the A/G genotype showed a positive relationship between average licking received and dopamine levels in the nucleus accumbens of the maternal brain; there was no relationship with female rat offspring with the A/A genotype. The higher dopamine levels in the nucleus accumbens corresponded with higher maternal licking provisioning from postnatal days 2-9. We also discovered and validated several new variants that were predicted by our systematic analysis.

CONCLUSION

Our findings suggest that genetic variation influences the relationship between early-life maternal care received and the dopaminergic system of the maternal brain, which can indirectly influence later-life maternal care provisioning.

A platform for experimental precision medicine: The extended BXD mouse family

The challenge of precision medicine is to model complex interactions among DNA variants, phenotypes, development, environments, and treatments. We address this challenge by expanding the BXD family of mice to 140 fully isogenic strains, creating a uniquely powerful model for precision medicine. This family segregates for 6 million common DNA variants-a level that exceeds many human populations. Because each member can be replicated, heritable traits can be mapped with high power and precision. Current BXD phenomes are unsurpassed in coverage and include much omics data and thousands of quantitative traits. BXDs can be extended by a single-generation cross to as many as 19,460 isogenic F1 progeny, and this extended BXD family is an effective platform for testing causal modeling and for predictive validation. BXDs are a unique core resource for the field of experimental precision medicine.

Genetic Dissection of the Regulatory Mechanisms of Ace2 in the Infected Mouse Lung

Acute lung injury (ALI) is an important cause of morbidity and mortality after viral infections, including influenza A virus H1N1, SARS-CoV, MERS-CoV, and SARS-CoV-2. The angiotensin I converting enzyme 2 (ACE2) is a key host membrane-bound protein that modulates ALI induced by viral infection, pulmonary acid aspiration, and sepsis. However, the contributions of ACE2 sequence variants to individual differences in disease risk and severity after viral infection are not understood. In this study, we quantified H1N1 influenza-infected lung transcriptomes across a family of 41 BXD recombinant inbred strains of mice and both parents—C57BL/6J and DBA/2J. In response to infection Ace2 mRNA levels decreased significantly for both parental strains and the expression levels was associated with disease severity (body weight loss) and viral load (expression levels of viral NA segment) across the BXD family members. Pulmonary RNA-seq for 43 lines was analyzed using weighted gene co-expression network analysis (WGCNA) and Bayesian network approaches. Ace2 not only participated in virus-induced ALI by interacting with TNF, MAPK, and NOTCH signaling pathways, but was also linked with high confidence to gene products that have important functions in the pulmonary epithelium, including Rnf128, Muc5b, and Tmprss2. Comparable sets of transcripts were also highlighted in parallel studies of human SARS-CoV-infected primary human airway epithelial cells. Using conventional mapping methods, we determined that weight loss at two and three days after viral infection maps to chromosome X—the location of Ace2. This finding motivated the hierarchical Bayesian network analysis, which defined molecular endophenotypes of lung infection linked to Ace2 expression and to a key disease outcome. Core members of this Bayesian network include Ace2, Atf4, Csf2, Cxcl2, Lif, Maml3, Muc5b, Reg3g, Ripk3, and Traf3. Collectively, these findings define a causally-rooted Ace2 modulatory network relevant to host response to viral infection and identify potential therapeutic targets for virus-induced respiratory diseases, including those caused by influenza and coronaviruses.

Variability and heritability of mouse brain structure: Microscopic MRI atlases and connectomes for diverse strains

Genome-wide association studies have demonstrated significant links between human brain structure and common DNA variants. Similar studies with rodents have been challenging because of smaller brain volumes. Using high field MRI (9.4 T) and compressed sensing, we have achieved microscopic resolution and sufficiently high throughput for rodent population studies. We generated whole brain structural MRI and diffusion connectomes for four diverse isogenic lines of mice (C57BL/6J, DBA/2J, CAST/EiJ, and BTBR) at spatial resolution 20,000 times higher than human connectomes. We measured narrow sense heritability (h2) I.e. the fraction of variance explained by strains in a simple ANOVA model for volumes and scalar diffusion metrics, and estimates of residual technical error for 166 regions in each hemisphere and connectivity between the regions. Volumes of discrete brain regions had the highest mean heritability (0.71 ± 0.23 SD, n = 332), followed by fractional anisotropy (0.54 ± 0.26), radial diffusivity (0.34 ± 0.022), and axial diffusivity (0.28 ± 0.19). Connection profiles were statistically different in 280 of 322 nodes across all four strains. Nearly 150 of the connection profiles were statistically different between the C57BL/6J, DBA/2J, and CAST/EiJ lines. Microscopic whole brain MRI/DTI has allowed us to identify significant heritable phenotypes in brain volume, scalar DTI metrics, and quantitative connectomes.

Genome-wide transcriptome architecture in a mouse model of Gulf War Illness

Gulf War Illness (GWI) is thought to be a chronic neuroimmune disorder caused by in-theater exposure during the 1990-1991 Gulf War. There is a consensus that the illness is caused by exposure to insecticides and nerve agent toxicants. However, the heterogeneity in both development of disease and clinical outcomes strongly suggests a genetic contribution. Here, we modeled GWI in 30 BXD recombinant inbred mouse strains with a combined treatment of corticosterone (CORT) and diisopropyl fluorophosphate (DFP). We quantified transcriptomes from 409 prefrontal cortex samples. Compared to the untreated and DFP treated controls, the combined treatment significantly activated pathways such as cytokine-cytokine receptor interaction and TNF signaling pathway. Protein-protein interaction analysis defined 6 subnetworks for CORT + DFP, with the key regulators being Cxcl1, Il6, Ccnb1, Tnf, Agt, and Itgam. We also identified 21 differentially expressed genes having significant QTLs related to CORT + DFP, but without evidence for untreated and DFP treated controls, suggesting regions of the genome specifically involved in the response to CORT + DFP. We identified Adamts9 as a potential contributor to response to CORT + DFP and found links to symptoms of GWI. Furthermore, we observed a significant effect of CORT + DFP treatment on the relative proportion of myelinating oligodendrocytes, with a QTL on Chromosome 5. We highlight three candidates, Magi2, Sema3c, and Gnai1, based on their high expression in the brain and oligodendrocyte. In summary, our results show significant genetic effects of the CORT + DFP treatment, which mirrors gene and protein expression changes seen in GWI sufferers, providing insight into the disease and a testbed for future interventions.

Body weight and high-fat diet are associated with epigenetic aging in female members of the BXD murine family

DNA methylation (DNAm) is shaped by genetic and environmental factors and modulated by aging. Here, we examine interrelations between epigenetic aging, body weight (BW), and life span in 12 isogenic strains from the BXD family of mice that exhibit over twofold variation in longevity. Genome-wide DNAm was assayed in 70 liver specimens from predominantly female cases, 6-25 months old, that were maintained on normal chow or high-fat diet (HFD). We defined subsets of CpG regions associated with age, BW at young adulthood, and strain-by-diet-dependent life span. These age-associated differentially methylated CpG regions (age-DMRs) featured distinct genomic characteristics, with DNAm gains over time occurring in sites such as promoters and exons that have high CpG density and low average methylation. CpG regions associated with BW were enriched in introns, tended to have lower methylation in mice with higher BW, and were inversely correlated with gene expression (i.e., higher mRNA levels in mice with higher BW). CpG regions associated with life span were linked to genes involved in life span modulation, including the telomerase reverse transcriptase gene, Tert, which had both lower methylation and higher expression in long-lived strains. An epigenetic clock defined from age-DMRs revealed accelerated aging in mice belonging to strains with shorter life spans. Both higher BW and the HFD were associated with accelerated epigenetic aging. Our results highlight the age-accelerating effect of heavier BW. Furthermore, we demonstrate that the measure of epigenetic aging derived from age-DMRs can predict genotype and diet-induced differences in life span among female BXD members.

Modeling the Genetic Basis of Individual Differences in Susceptibility to Gulf War Illness

Between 25% and 30% of the nearly one million military personnel who participated in the 1991 Persian Gulf War became ill with chronic symptoms ranging from gastrointestinal to nervous system dysfunction. This disorder is now referred to as Gulf War Illness (GWI) and the underlying pathophysiology has been linked to exposure-based neuroinflammation caused by organophosphorous (OP) compounds coupled with high circulating glucocorticoids. In a mouse model of GWI we developed, corticosterone was shown to act synergistically with an OP (diisopropylflurophosphate) to dramatically increase proinflammatory cytokine gene expression in the brain. Because not all Gulf War participants became sick, the question arises as to whether differential genetic constitution might underlie individual differences in susceptibility. To address this question of genetic liability, we tested the impact of OP and glucocorticoid exposure in a genetic reference population of 30 inbred mouse strains. We also studied both sexes. The results showed wide differences among strains and overall that females were less sensitive to the combined treatment than males. Furthermore, we identified one OP-glucocorticoid locus and nominated a candidate gene-Spon1-that may underlie the marked differences in response.

A Cross-Species Systems Genetics Analysis Links APBB1IP as a Candidate for Schizophrenia and Prepulse Inhibition

Background: Prepulse inhibition (PPI) of the startle response is a highly conserved form of sensorimotor gating, disruption of which is found in schizophrenia patients and their unaffected first-degree relatives. PPI can be measured in many species, and shows considerable phenotypic variation between and within rodent models. This makes PPI a useful endophenotype. Genome-wide association studies (GWAS) have been carried out to identify genetic variants underlying schizophrenia, and these suggest that schizophrenia is highly polygenic. GWAS have been unable to account for the high heritability of schizophrenia seen in family studies, partly because of the low power of GWAS due to multiple comparisons. By contrast, complementary mouse model linkage studies often have high statistical power to detect variants for behavioral traits but lower resolution, producing loci that include tens or hundreds of genes. To capitalize on the advantages of both GWAS and genetic mouse models, our study uses a cross-species approach to identify novel genes associated with PPI regulation, which thus may contribute to the PPI deficits seen in schizophrenia.

Results: Using experimental data from the recombinant inbred (RI) mouse panel BXD, we identified two significant loci affecting PPI. These genomic regions contain genetic variants which influence PPI in mice and are therefore candidates that may be influencing aspects of schizophrenia in humans. We next investigated these regions in whole-genome data from the Psychiatric Genomics Consortium (PGC) schizophrenia GWAS and identify one novel candidate gene (ABPP1IP) that was significantly associated with PPI in mice and risk of schizophrenia in humans. A systems genetics approach demonstrates that APBB1IP coexpresses with several other genes related to schizophrenia in several brain regions. Gene coexpression and enrichment analysis shows clear links between APBB1IP and the immune system.

Conclusion: The combination of human GWAS and mouse quantitative trait loci (QTL) from some of the largest study systems available has enabled us to identify a novel gene, APBB1IP, which influences schizophrenia in humans and PPI in mice.

Offspring genetic effects on maternal care

Parental care is found widely across animal taxa and is manifest in a range of behaviours from basic provisioning in cockroaches to highly complex behaviours seen in mammals. The evolution of parental care is viewed as the outcome of an evolutionary cost/benefit trade-off between investing in current and future offspring, leading to the selection of traits in offspring that influence parental behaviour. Thus, level and quality of parental care are affected by both parental and offspring genetic differences that directly and indirectly influence parental care behaviour. While significant research effort has gone into understanding how parental genomes affect parental, and mostly maternal, behaviour, few studies have investigated how offspring genomes affect parental care. In this review, we bring together recent findings across different fields focussing on the mechanism and genetics of offspring effects on maternal care in mammals.

Epigenetic impacts of stress priming of the neuroinflammatory response to sarin surrogate in mice: a model of Gulf War illness

Background

Gulf War illness (GWI) is an archetypal, medically unexplained, chronic condition characterised by persistent sickness behaviour and neuroimmune and neuroinflammatory components. An estimated 25–32% of the over 900,000 veterans of the 1991 Gulf War fulfil the requirements of a GWI diagnosis. It has been hypothesised that the high physical and psychological stress of combat may have increased vulnerability to irreversible acetylcholinesterase (AChE) inhibitors leading to a priming of the neuroimmune system. A number of studies have linked high levels of psychophysiological stress and toxicant exposures to epigenetic modifications that regulate gene expression. Recent research in a mouse model of GWI has shown that pre-exposure with the stress hormone corticosterone (CORT) causes an increase in expression of specific chemokines and cytokines in response to diisopropyl fluorophosphate (DFP), a sarin surrogate and irreversible AChE inhibitor.

Methods

C57BL/6J mice were exposed to CORT for 4 days, and exposed to DFP on day 5, before sacrifice 6 h later. The transcriptome was examined using RNA-seq, and the epigenome was examined using reduced representation bisulfite sequencing and H3K27ac ChIP-seq.

Results

We show transcriptional, histone modification (H3K27ac) and DNA methylation changes in genes related to the immune and neuronal system, potentially relevant to neuroinflammatory and cognitive symptoms of GWI. Further evidence suggests altered proportions of myelinating oligodendrocytes in the frontal cortex, perhaps connected to white matter deficits seen in GWI sufferers.

Conclusions

Our findings may reflect the early changes which occurred in GWI veterans, and we observe alterations in several pathways altered in GWI sufferers. These close links to changes seen in veterans with GWI indicates that this model reflects the environmental exposures related to GWI and may provide a model for biomarker development and testing future treatments.

Electronic supplementary material

The online version of this article (10.1186/s12974-018-1113-9) contains supplementary material, which is available to authorized users.

Post-genomic behavioral genetics: From revolution to routine

What was once expensive and revolutionary-full-genome sequence-is now affordable and routine. Costs will continue to drop, opening up new frontiers in behavioral genetics. This shift in costs from the genome to the phenome is most notable in large clinical studies of behavior and associated diseases in cohorts that exceed hundreds of thousands of subjects. Examples include the Women's Health Initiative (www.whi.org), the Million Veterans Program (www.

RESEARCH

va.gov/MVP), the 100 000 Genomes Project (genomicsengland.co.uk) and commercial efforts such as those by deCode (www.decode.com) and 23andme (www.23andme.com). The same transition is happening in experimental neuro- and behavioral genetics, and sample sizes of many hundreds of cases are becoming routine (www.genenetwork.org, www.mousephenotyping.org). There are two major consequences of this new affordability of massive omics datasets: (1) it is now far more practical to explore genetic modulation of behavioral differences and the key role of gene-by-environment interactions. Researchers are already doing the hard part-the quantitative analysis of behavior. Adding the omics component can provide powerful links to molecules, cells, circuits and even better treatment. (2) There is an acute need to highlight and train behavioral scientists in how best to exploit new omics approaches. This review addresses this second issue and highlights several new trends and opportunities that will be of interest to experts in animal and human behaviors.

Offspring genes indirectly influence sibling and maternal behavioural strategies over resource share

Family members show behavioural strategies predicted to maximize individual fitness. These behaviours depend directly on genes expressed in focal individuals but also indirectly on genes expressed in other family members. However, how sibling and parental behavioural strategies are modified by genes expressed in family members, and to what degree, remains unclear. To answer this question, we have used a split litter design in an experimental population of genetically variable mouse families, and identified loci that indirectly affected sibling and maternal behaviour simultaneously. These loci map to genomic regions that also show a direct effect on offspring behaviour. Directly and indirectly affected traits were significantly correlated at the phenotypic level, illustrating how indirect effects are caused. Genetic variants in offspring that influence solicitation also impacted their siblings' and maternal behaviour. However, in contrast to predictions from sibling competition, unrelated litter mates benefited from increased solicitation. Overall, such indirect genetic effects explained a large proportion of variation seen in behaviours, with candidate genes involved in metabolism to neuronal development. These results reveal that we need to view behavioural strategies as the result of conjoint selection on genetic variation in all interacting family members.

Social Interactions and Indirect Genetic Effects on Complex Juvenile and Adult Traits

Most animal species are social in one form or another, yet many studies in rodent model systems use either individually housed animals or ignore potential confounds caused by group housing. While such social interaction effects on developmental and behavioral traits are well established, the genetic basis of social interactions has not been researched in as much detail. Specifically, the effects of genetic variation in social partners on the phenotype of a focal individual have mostly been studied at the phenotypic level. Such indirect genetic effects (IGEs), where the genotype of one individual influences the phenotype of a second individual, can have important evolutionary and medically relevant consequences. In this chapter, we give a brief outline of social interaction effects, and how systems genetics approaches using recombinant inbred populations can be used to investigate indirect genetic effects specifically, including maternal genetic effects. We discuss experimental designs for the study of IGEs and show how indirect genetic loci can be identified that underlie social interaction effects, their mechanisms, and consequences for trait variation in focal individuals.

A cross-species genetic analysis identifies candidate genes for mouse anxiety and human bipolar disorder

Bipolar disorder (BD) is a significant neuropsychiatric disorder with a lifetime prevalence of ~1%. To identify genetic variants underlying BD genome-wide association studies (GWAS) have been carried out. While many variants of small effect associated with BD have been identified few have yet been confirmed, partly because of the low power of GWAS due to multiple comparisons being made. Complementary mapping studies using murine models have identified genetic variants for behavioral traits linked to BD, often with high power, but these identified regions often contain too many genes for clear identification of candidate genes. In the current study we have aligned human BD GWAS results and mouse linkage studies to help define and evaluate candidate genes linked to BD, seeking to use the power of the mouse mapping with the precision of GWAS. We use quantitative trait mapping for open field test and elevated zero maze data in the largest mammalian model system, the BXD recombinant inbred mouse population, to identify genomic regions associated with these BD-like phenotypes. We then investigate these regions in whole genome data from the Psychiatric Genomics Consortium's bipolar disorder GWAS to identify candidate genes associated with BD. Finally we establish the biological relevance and pathways of these genes in a comprehensive systems genetics analysis. We identify four genes associated with both mouse anxiety and human BD. While TNR is a novel candidate for BD, we can confirm previously suggested associations with CMYA5, MCTP1, and RXRG. A cross-species, systems genetics analysis shows that MCTP1, RXRG, and TNR coexpress with genes linked to psychiatric disorders and identify the striatum as a potential site of action. CMYA5, MCTP1, RXRG, and TNR are associated with mouse anxiety and human BD. We hypothesize that MCTP1, RXRG, and TNR influence intercellular signaling in the striatum.

Common genetic variants influence human subcortical brain structures

The highly complex structure of the human brain is strongly shaped by genetic influences. Subcortical brain regions form circuits with cortical areas to coordinate movement, learning, memory and motivation, and altered circuits can lead to abnormal behaviour and disease. To investigate how common genetic variants affect the structure of these brain regions, here we conduct genome-wide association studies of the volumes of seven subcortical regions and the intracranial volume derived from magnetic resonance images of 30,717 individuals from 50 cohorts. We identify five novel genetic variants influencing the volumes of the putamen and caudate nucleus. We also find stronger evidence for three loci with previously established influences on hippocampal volume and intracranial volume. These variants show specific volumetric effects on brain structures rather than global effects across structures. The strongest effects were found for the putamen, where a novel intergenic locus with replicable influence on volume (rs945270; P = 1.08 × 10(-33); 0.52% variance explained) showed evidence of altering the expression of the KTN1 gene in both brain and blood tissue. Variants influencing putamen volume clustered near developmental genes that regulate apoptosis, axon guidance and vesicle transport. Identification of these genetic variants provides insight into the causes of variability in human brain development, and may help to determine mechanisms of neuropsychiatric dysfunction.

Transcript co-variance with Nestin in two mouse genetic reference populations identifies Lef1 as a novel candidate regulator of neural precursor cell proliferation in the adult hippocampus

Adult neurogenesis, the lifelong production of new neurons in the adult brain, is under complex genetic control but many of the genes involved remain to be identified. In this study, we have integrated publicly available gene expression data from the BXD and CXB recombinant inbred mouse lines to discover genes co-expressed in the adult hippocampus with Nestin, a common marker of the neural precursor cell population. In addition, we incorporated spatial expression information to restrict candidates to genes with high differential gene expression in the hippocampal dentate gyrus. Incorporating data from curated protein-protein interaction databases revealed interactions between our candidate genes and those already known to be involved in adult neurogenesis. Enrichment analysis suggested a link to the Wnt/β-catenin pathway, known to be involved in adult neurogenesis. In particular, our candidates were enriched in targets of Lef1, a modulator of the Wnt pathway. In conclusion, our combination of bioinformatics approaches identified six novel candidate genes involved in adult neurogenesis; Amer3, Eya3, Mtdh, Nr4a3, Polr2a, and Tbkbp1. Further, we propose a role for Lef1 transcriptional control in the regulation of adult hippocampal precursor cell proliferation.

Joint genetic analysis of hippocampal size in mouse and human identifies a novel gene linked to neurodegenerative disease

BACKGROUND

Variation in hippocampal volume has been linked to significant differences in memory, behavior, and cognition among individuals. To identify genetic variants underlying such differences and associated disease phenotypes, multinational consortia such as ENIGMA have used large magnetic resonance imaging (MRI) data sets in human GWAS studies. In addition, mapping studies in mouse model systems have identified genetic variants for brain structure variation with great power. A key challenge is to understand how genetically based differences in brain structure lead to the propensity to develop specific neurological disorders.

RESULTS

We combine the largest human GWAS of brain structure with the largest mammalian model system, the BXD recombinant inbred mouse population, to identify novel genetic targets influencing brain structure variation that are linked to increased risk for neurological disorders. We first use a novel cross-species, comparative analysis using mouse and human genetic data to identify a candidate gene, MGST3, associated with adult hippocampus size in both systems. We then establish the coregulation and function of this gene in a comprehensive systems-analysis.

CONCLUSIONS

We find that MGST3 is associated with hippocampus size and is linked to a group of neurodegenerative disorders, such as Alzheimer's.

Empirical testing of hypotheses about the evolution of genomic imprinting in mammals

The close interaction between mother and offspring in mammals is thought to contribute to the evolution of genomic imprinting or parent-of-origin dependent gene expression. Empirical tests of theories about the evolution of imprinting have been scant for several reasons. Models make different assumptions about the traits affected by imprinted genes and the scenarios in which imprinting is predicted to have been selected for. Thus, competing hypotheses cannot readily be tested against each other. Further, it is far from clear how predictions about expression patterns of genes with specific phenotypic effects can be tested given current methodology of assaying gene expression levels, be it in the brain or in other tissues. We first set out a scenario for testing competing hypotheses and delineate the different assumptions and predictions of models. We then outline how predictions may be tested using mouse models such as intercrosses or recombinant inbred (RI) systems that can be phenotyped for traits relevant to imprinting theories. Further, we briefly discuss different molecular approaches that may be used in conjunction with experiments to ascertain expression patterns of imprinted genes and thus the testing of predictions.