The Collaborative Cross, a community resource for the genetic analysis of complex traits

Preclinical murine models for the testing of antimicrobials against Mycobacterium abscessus pulmonary infections: Current practices and recommendations

Mycobacterium abscessus, a rapidly growing nontuberculous mycobacterium, is increasingly recognized as an important pathogen of the human lung, disproportionally affecting people with cystic fibrosis (CF) and other susceptible individuals with non-CF bronchiectasis and compromised immune functions. M. abscessus infections are extremely difficult to treat due to intrinsic resistance to many antibiotics, including most anti-tuberculous drugs. Current standard-of-care chemotherapy is long, includes multiple oral and parenteral repurposed drugs, and is associated with significant toxicity. The development of more effective oral antibiotics to treat M. abscessus infections has thus emerged as a high priority. While murine models have proven instrumental in predicting the efficacy of therapeutic treatments for M. tuberculosis infections, the preclinical evaluation of drugs against M. abscessus infections has proven more challenging due to the difficulty of establishing a progressive, sustained, pulmonary infection with this pathogen in mice. To address this issue, a series of three workshops were hosted in 2023 by the Cystic Fibrosis Foundation (CFF) and the National Institute of Allergy and Infectious Diseases (NIAID) to review the current murine models of M. abscessus infections, discuss current challenges and identify priorities toward establishing validated and globally harmonized preclinical models. This paper summarizes the key points from these workshops. The hope is that the recommendations that emerged from this exercise will facilitate the implementation of informative murine models of therapeutic efficacy testing across laboratories, improve reproducibility from lab-to-lab and accelerate preclinical-to-clinical translation.

Quantitative trait loci mapping provides insights into the genetic regulation of dendritic cell numbers in mouse tissues

To explore the influence of genetics on homeostatic regulation of dendritic cell (DC) numbers, we present a screen of DCs and their progenitors in lymphoid and non-lymphoid tissues in Collaborative Cross (CC) and Diversity Outbred (DO) mice. We report 30 and 71 loci with logarithm of the odds (LOD) scores >8.18 and ranging from 6.67 to 8.19, respectively. The analysis reveals the highly polygenic and pleiotropic architecture of this complex trait, including many of the previously identified genetic regulators of DC development and maturation. Two SNPs in genes potentially underlying variation in DC homeostasis, a splice variant in Gramd4 (rs235532740) and a missense variant in Orai3 (rs216659754), are confirmed by gene editing using CRISPR-Cas9. Gramd4 is a central regulator of DC homeostasis that impacts the entire DC lineage, and Orai3 regulates cDC2 numbers in tissues. Overall, the data reveal a large number of candidate genes regulating DC homeostasis in vivo.

Studying the Effect of the Host Genetic Background of Juvenile Polyposis Development Using Collaborative Cross and Smad4 Knock-Out Mouse Models

Juvenile polyposis syndrome (JPS) is a rare autosomal dominant disorder characterized by multiple juvenile polyps in the gastrointestinal tract, often associated with mutations in genes such as Smad4 and BMPR1A. This study explores the impact of Smad4 knock-out on the development of intestinal polyps using collaborative cross (CC) mice, a genetically diverse model. Our results reveal a significant increase in intestinal polyps in Smad4 knock-out mice across the entire population, emphasizing the broad influence of Smad4 on polyposis. Sex-specific analyses demonstrate higher polyp counts in knock-out males and females compared to their WT counterparts, with distinct correlation patterns. Line-specific effects highlight the nuanced response to Smad4 knock-out, underscoring the importance of genetic variability. Multimorbidity heat maps offer insights into complex relationships between polyp counts, locations, and sizes. Heritability analysis reveals a significant genetic basis for polyp counts and sizes, while machine learning models, including k-nearest neighbors and linear regression, identify key predictors, enhancing our understanding of juvenile polyposis genetics. Overall, this study provides new information on understanding the intricate genetic interplay in the context of Smad4 knock-out, offering valuable insights that could inform the identification of potential therapeutic targets for juvenile polyposis and related diseases.

Identification of Arginine-Vasopressin Receptor 1a (Avpr1a/Avpr1a) as a Novel Candidate Gene for Chronic Visceral Pain Sheds Light on the Potential Role of Enteric Neurons in the Development of Visceral Hypersensitivity

Chronic abdominal pain in the absence of ongoing disease is the hallmark of disorders of gut-brain interaction (DGBIs), including irritable bowel syndrome (IBS). While the etiology of DGBIs remains poorly understood, there is evidence that both genetic and environmental factors play a role. In this study, we report the identification and validation of arginine-vasopressin receptor 1A (Avpr1a) as a novel candidate gene for visceral hypersensitivity (VH), a primary peripheral mechanism underlying abdominal pain in DGBI/IBS. Comparing 2 C57BL/6 (BL/6) substrains (C57BL/6NTac and C57BL/6J) revealed differential susceptibility to the development of chronic VH following intrarectal zymosan instillation, a validated preclinical model for postinflammatory IBS. Using whole-genome sequencing, we identified a single-nucleotide polymorphism differentiating the 2 strains in the 5' intergenic region upstream of Avpr1a, encoding the protein Avpr1a. We used behavioral, histological, and molecular approaches to identify distal colon-specific gene expression and neuronal hyperresponsiveness covarying with Avpr1a genotype and VH susceptibility. While the 2 BL/6 substrains did not differ across other gastrointestinal phenotypes (eg, fecal water retention), VH-susceptible BL/6NTac mice had higher colonic Avpr1a mRNA and protein expression. These results parallel findings that patients' colonic Avpr1a mRNA expression corresponded to higher pain ratings. Moreover, neurons of the enteric nervous system were hyperresponsive to the Avpr1a agonist arginine-vasopressin, suggesting a role for enteric neurons in the pathology underlying VH. Taken together, these findings implicate differential regulation of Avpr1a as a novel mechanism of VH susceptibility as well as a potential therapeutic target specific to VH. PERSPECTIVE: This article presents evidence of Avpr1a as a novel candidate gene for VH in a mouse model of IBS. Avpr1a genotype and/or tissue-specific expression represents a potential biomarker for chronic abdominal pain susceptibility.

Collaborative Cross mice have diverse phenotypic responses to infection with Methicillin-resistant Staphylococcus aureus USA300

Staphylococcus aureus (S. aureus) is an opportunistic pathogen causing diseases ranging from mild skin infections to life threatening conditions, including endocarditis, pneumonia, and sepsis. To identify host genes modulating this host-pathogen interaction, we infected 25 Collaborative Cross (CC) mouse strains with methicillin-resistant S. aureus (MRSA) and monitored disease progression for seven days using a surgically implanted telemetry system. CC strains varied widely in their response to intravenous MRSA infection. We identified eight 'susceptible' CC strains with high bacterial load, tissue damage, and reduced survival. Among the surviving strains, six with minimal colonization were classified as 'resistant', while the remaining six tolerated higher organ colonization ('tolerant'). The kidney was the most heavily colonized organ, but liver, spleen and lung colonization were better correlated with reduced survival. Resistant strains had higher pre-infection circulating neutrophils and lower post-infection tissue damage compared to susceptible and tolerant strains. We identified four CC strains with sexual dimorphism: all females survived the study period while all males met our euthanasia criteria earlier. In these CC strains, males had more baseline circulating monocytes and red blood cells. We identified several CC strains that may be useful as new models for endocarditis, myocarditis, pneumonia, and resistance to MRSA infection. Quantitative Trait Locus (QTL) analysis identified two significant loci, on Chromosomes 18 and 3, involved in early susceptibility and late survival after infection. We prioritized Npc1 and Ifi44l genes as the strongest candidates influencing survival using variant analysis and mRNA expression data from kidneys within these intervals.

A mutation in Themis contributes to anaphylaxis severity following oral peanut challenge in CC027 mice

BACKGROUND

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized Collaborative Cross strain CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, in contrast to C3H/HeJ (C3H) mice.

OBJECTIVE

This study aimed to determine the genetic basis of orally induced anaphylaxis to peanut in CC027 mice.

METHODS

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 mice and 5 additional Collaborative Cross strains.

RESULTS

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis and 4% having severe anaphylaxis. There were 8 genetic loci associated with variation in response to peanut challenge-6 associated with anaphylaxis (temperature decrease) and 2 associated with peanut-specific IgE levels. There were 2 major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis gene. Consistent with described functions of Themis, we found that CC027 mice have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

CONCLUSIONS

Our results demonstrate a key role for Themis in the orally reactive CC027 mouse model of peanut allergy.

Genetic context modulates aging and degeneration in the murine retina

Background

Age is the principal risk factor for neurodegeneration in both the retina and brain. The retina and brain share many biological properties; thus, insights into retinal aging and degeneration may shed light onto similar processes in the brain. Genetic makeup strongly influences susceptibility to age-related retinal disease. However, studies investigating retinal aging have not sufficiently accounted for genetic diversity. Therefore, examining molecular aging in the retina across different genetic backgrounds will enhance our understanding of human-relevant aging and degeneration in both the retina and brain-potentially improving therapeutic approaches to these debilitating conditions.

Methods

Transcriptomics and proteomics were employed to elucidate retinal aging signatures in nine genetically diverse mouse strains (C57BL/6J, 129S1/SvlmJ, NZO/HlLtJ, WSB/EiJ, CAST/EiJ, PWK/PhK, NOD/ShiLtJ, A/J, and BALB/cJ) across lifespan. These data predicted human disease-relevant changes in WSB and NZO strains. Accordingly, B6, WSB and NZO mice were subjected to human-relevant in vivo examinations at 4, 8, 12, and/or 18M, including: slit lamp, fundus imaging, optical coherence tomography, fluorescein angiography, and pattern/full-field electroretinography. Retinal morphology, vascular structure, and cell counts were assessed ex vivo.

Results

We identified common molecular aging signatures across the nine mouse strains, which included genes associated with photoreceptor function and immune activation. Genetic background strongly modulated these aging signatures. Analysis of cell type-specific marker genes predicted age-related loss of photoreceptors and retinal ganglion cells (RGCs) in WSB and NZO, respectively. Fundus exams revealed retinitis pigmentosa-relevant pigmentary abnormalities in WSB retinas and diabetic retinopathy (DR)-relevant cotton wool spots and exudates in NZO retinas. Profound photoreceptor dysfunction and loss were confirmed in WSB. Molecular analyses indicated changes in photoreceptor-specific proteins prior to loss, suggesting photoreceptor-intrinsic dysfunction in WSB. In addition, age-associated RGC dysfunction, loss, and concomitant microvascular dysfunction was observed in NZO mice. Proteomic analyses revealed an early reduction in protective antioxidant processes, which may underlie increased susceptibility to DR-relevant pathology in NZO.

Conclusions

Genetic context is a strong determinant of retinal aging, and our multi-omics resource can aid in understanding age-related diseases of the eye and brain. Our investigations identified and validated WSB and NZO mice as improved preclinical models relevant to common retinal neurodegenerative diseases.

The MexTAg collaborative cross: host genetics affects asbestos related disease latency, but has little influence once tumours develop

Objectives: This study combines two innovative mouse models in a major gene discovery project to assess the influence of host genetics on asbestos related disease (ARD). Conventional genetics studies provided evidence that some susceptibility to mesothelioma is genetic. However, the identification of host modifier genes, the roles they may play, and whether they contribute to disease susceptibility remain unknown. Here we report a study designed to rapidly identify genes associated with mesothelioma susceptibility by combining the Collaborative Cross (CC) resource with the well-characterised MexTAg mesothelioma mouse model. Methods: The CC is a powerful mouse resource that harnesses over 90% of common genetic variation in the mouse species, allowing rapid identification of genes mediating complex traits. MexTAg mice rapidly, uniformly, and predictably develop mesothelioma, but only after asbestos exposure. To assess the influence of host genetics on ARD, we crossed 72 genetically distinct CC mouse strains with MexTAg mice and exposed the resulting CC-MexTAg (CCMT) progeny to asbestos and monitored them for traits including overall survival, the time to ARD onset (latency), the time between ARD onset and euthanasia (disease progression) and ascites volume. We identified phenotype-specific modifier genes associated with these traits and we validated the role of human orthologues in asbestos-induced carcinogenesis using human mesothelioma datasets. Results: We generated 72 genetically distinct CCMT strains and exposed their progeny (2,562 in total) to asbestos. Reflecting the genetic diversity of the CC, there was considerable variation in overall survival and disease latency. Surprisingly, however, there was no variation in disease progression, demonstrating that host genetic factors do have a significant influence during disease latency but have a limited role once disease is established. Quantitative trait loci (QTL) affecting ARD survival/latency were identified on chromosomes 6, 12 and X. Of the 97-protein coding candidate modifier genes that spanned these QTL, eight genes (CPED1, ORS1, NDUFA1, HS1BP3, IL13RA1, LSM8, TES and TSPAN12) were found to significantly affect outcome in both CCMT and human mesothelioma datasets. Conclusion: Host genetic factors affect susceptibility to development of asbestos associated disease. However, following mesothelioma establishment, genetic variation in molecular or immunological mechanisms did not affect disease progression. Identification of multiple candidate modifier genes and their human homologues with known associations in other advanced stage or metastatic cancers highlights the complexity of ARD and may provide a pathway to identify novel therapeutic targets.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

Embracing cancer complexity: Hallmarks of systemic disease

The last 50 years have witnessed extraordinary developments in understanding mechanisms of carcinogenesis, synthesized as the hallmarks of cancer. Despite this logical framework, our understanding of the molecular basis of systemic manifestations and the underlying causes of cancer-related death remains incomplete. Looking forward, elucidating how tumors interact with distant organs and how multifaceted environmental and physiological parameters impinge on tumors and their hosts will be crucial for advances in preventing and more effectively treating human cancers. In this perspective, we discuss complexities of cancer as a systemic disease, including tumor initiation and promotion, tumor micro- and immune macro-environments, aging, metabolism and obesity, cancer cachexia, circadian rhythms, nervous system interactions, tumor-related thrombosis, and the microbiome. Model systems incorporating human genetic variation will be essential to decipher the mechanistic basis of these phenomena and unravel gene-environment interactions, providing a modern synthesis of molecular oncology that is primed to prevent cancers and improve patient quality of life and cancer outcomes.

Unique immune profiles in collaborative cross mice linked to survival and viral clearance upon infection

The response to infection is generally heterogeneous and diverse, with some individuals remaining asymptomatic while others present with severe disease or a diverse range of symptoms. Here, we address the role of host genetics on immune phenotypes and clinical outcomes following viral infection by studying genetically diverse mice from the Collaborative Cross (CC), allowing for use of a small animal model with controlled genetic diversity while maintaining genetic replicates. We demonstrate variation by deeply profiling a broad range of innate and adaptive immune cell phenotypes at steady-state in 63 genetically distinct CC mouse strains and link baseline immune signatures with virologic and clinical disease outcomes following infection of mice with herpes simplex virus 2 (HSV-2) or severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This work serves as a resource for CC strain selection based on steady-state immune phenotypes or disease presentation upon viral infection, and further, points to possible pre-infection immune correlates of survival and early viral clearance upon infection.

Antigen identification strategies and preclinical evaluation models for advancing tuberculosis vaccine development

In its myriad devastating forms, Tuberculosis (TB) has existed for centuries, and humanity is still affected by it. Mycobacterium tuberculosis (M. tuberculosis), the causative agent of TB, was the foremost killer among infectious agents until the COVID-19 pandemic. One of the key healthcare strategies available to reduce the risk of TB is immunization with bacilli Calmette-Guerin (BCG). Although BCG has been widely used to protect against TB, reports show that BCG confers highly variable efficacy (0-80%) against adult pulmonary TB. Unwavering efforts have been made over the past 20 years to develop and evaluate new TB vaccine candidates. The failure of conventional preclinical animal models to fully recapitulate human response to TB, as also seen for the failure of MVA85A in clinical trials, signifies the need to develop better preclinical models for TB vaccine evaluation. In the present review article, we outline various approaches used to identify protective mycobacterial antigens and recent advancements in preclinical models for assessing the efficacy of candidate TB vaccines.

Cross-sectional association between blood cholesterol and calcium levels in genetically diverse strains of mice

Genetically diverse outbred mice allow for the study of genetic variation in the context of high dietary and environmental control. Using a machine learning approach, we investigated clinical and morphometric factors that associate with serum cholesterol levels in 840 genetically unique Diversity Outbred mice of both sexes (n = 417 male and 423 female), and on both a control chow (% kcals in diet: protein 22%, carbohydrate 62%, fat 16%, no cholesterol) and high fat high sucrose (% kcals in diet: protein 15%, carbohydrate 41%, fat 45%, 0.05% cholesterol). We find expected elevations of cholesterol in male mice, as well as in mice with elevated serum triglycerides and/or fed a high fat high sucrose diet. The third strongest predictor was serum calcium which correlated with serum cholesterol across both diets and sexes (r = 0.39-0.48) in both Diversity Outbred (P = 3.0 × 10-43 ) and BXD (P = 0.005) mice. This is in-line with several human cohort studies which show associations between calcium and cholesterol, and calcium as an independent predictor of cardiovascular events.

Genetic variation influences the skeletal response to hindlimb unloading in the eight founder strains of the diversity outbred mouse population

During disuse, mechanical unloading causes extensive bone loss, decreasing bone volume and strength. Variations in bone mass and risk of osteoporosis are influenced by genetics; however, it remains unclear how genetic variation affects the skeletal response to unloading. We previously found that genetic variation affects the musculoskeletal response to 3 weeks of immobilization in the 8 Jackson Laboratory J:DO founder strains: C57Bl/6J, A/J, 129S1/SvImJ, NOD/ShiLtJ, NZO/HlLtJ, CAST/EiJ, PWK/PhJ, and WSB/EiJ. Hindlimb unloading (HLU) is the best model for simulating local and systemic contributors of disuse and therefore may have a greater impact on bones than immobilization. We hypothesized that genetic variation would affect the response to HLU across the eight founder strains. Mice of each founder strain were placed in HLU for 3 weeks, and the femurs and tibias were analyzed. There were significant HLU and mouse strain interactions on body weight, femur trabecular BV/TV, and femur ultimate force. This indicates that unloading only caused significant catabolic effects in some mouse strains. C57BL/6 J mice were most affected by unloading while other strains were more protected. There were significant HLU and mouse strain interactions on gene expression of genes encoding bone metabolism genes in the tibia. This indicates that unloading only caused significant effects on bone metabolism genes in some mouse strains. Different mouse strains respond to HLU differently, and this can be explained by genetic differences. These results suggest the outbred J:DO mice will be a powerful model for examining the effects of genetics on the skeletal response to HLU.

Genetic and microbial determinants of azoxymethane-induced colorectal tumor susceptibility in Collaborative Cross mice and their implication in human cancer

The insights into interactions between host genetics and gut microbiome (GM) in colorectal tumor susceptibility (CTS) remains lacking. We used Collaborative Cross mouse population model to identify genetic and microbial determinants of Azoxymethane-induced CTS. We identified 4417 CTS-associated single nucleotide polymorphisms (SNPs) containing 334 genes that were transcriptionally altered in human colorectal cancers (CRCs) and consistently clustered independent human CRC cohorts into two subgroups with different prognosis. We discovered a set of genera in early-life associated with CTS and defined a 16-genus signature that accurately predicted CTS, the majority of which were correlated with human CRCs. We identified 547 SNPs associated with abundances of these genera. Mediation analysis revealed GM as mediators partially exerting the effect of SNP UNC3869242 within Duox2 on CTS. Intestine cell-specific depletion of Duox2 altered GM composition and contribution of Duox2 depletion to CTS was significantly influenced by GM. Our findings provide potential novel targets for personalized CRC prevention and treatment.

Identification of arginine-vasopressin receptor 1a (Avpr1a/AVPR1A) as a novel candidate gene for chronic visceral pain

Chronic abdominal pain in the absence of ongoing disease is the hallmark of disorders of gut-brain interaction (DGBIs), including irritable bowel syndrome (IBS). While the etiology of DGBIs remains poorly understood, there is evidence that both genetic and environmental factors play a role. In this study, we report the identification and validation of Avpr1a as a novel candidate gene for visceral hypersensitivity (VH), a primary peripheral mechanism underlying abdominal pain in DGBI/IBS. Comparing two C57BL/6 (BL/6) substrains (C57BL/6NTac and C57BL/6J) revealed differential susceptibility to the development of chronic VH following intrarectal zymosan (ZYM) instillation, a validated preclinical model for post-inflammatory IBS. Using whole genome sequencing, we identified a SNP differentiating the two strains in the 5' intergenic region upstream of Avpr1a, encoding the protein arginine-vasopressin receptor 1A (AVPR1A). We used behavioral, histological, and molecular approaches to identify distal colon-specific gene expression differences and neuronal hyperresponsiveness covarying with Avpr1a genotype and VH susceptibility. While the two BL/6 substrains did not differ across other gastrointestinal (GI) phenotypes (e.g., GI motility), VH-susceptible BL/6NTac mice had higher colonic Avpr1a mRNA and protein expression. Moreover, neurons of the enteric nervous system were hyperresponsive to the AVPR1A agonist AVP, suggesting a role for enteric neurons in the pathology underlying VH. These results parallel our findings that patients' colonic Avpr1a mRNA expression was higher in patients with higher pain ratings. Taken together, these findings implicate differential regulation of Avpr1a as a novel mechanism of VH-susceptibility as well as a potential therapeutic target specific to VH.

Towards system genetics analysis of head and neck squamous cell carcinoma using the mouse model, cellular platform, and clinical human data

Head and neck squamous cell cancer (HNSCC) is a leading global malignancy. Every year, More than 830 000 people are diagnosed with HNSCC globally, with more than 430 000 fatalities. HNSCC is a deadly diverse malignancy with many tumor locations and biological characteristics. It originates from the squamous epithelium of the oral cavity, oropharynx, nasopharynx, larynx, and hypopharynx. The most frequently impacted regions are the tongue and larynx. Previous investigations have demonstrated the critical role of host genetic susceptibility in the progression of HNSCC. Despite the advances in our knowledge, the improved survival rate of HNSCC patients over the last 40 years has been limited. Failure to identify the molecular origins of development of HNSCC and the genetic basis of the disease and its biological heterogeneity impedes the development of new therapeutic methods. These results indicate a need to identify more genetic factors underlying this complex disease, which can be better used in early detection and prevention strategies. The lack of reliable animal models to investigate the underlying molecular processes is one of the most significant barriers to understanding HNSCC tumors. In this report, we explore and discuss potential research prospects utilizing the Collaborative Cross mouse model and crossing it to mice carrying single or double knockout genes (e.g. Smad4 and P53 genes) to identify genetic factors affecting the development of this complex disease using genome-wide association studies, epigenetics, microRNA, long noncoding RNA, lncRNA, histone modifications, methylation, phosphorylation, and proteomics.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes

OBJECTIVES

Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome.

METHODS

We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci.

RESULTS

We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation.

CONCLUSIONS

The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.

Intrinsic factors and CD1d1 but not CD1d2 expression levels control invariant natural killer T cell subset differentiation

Invariant natural killer T (NKT) cell subsets are defined based on their cytokine-production profiles and transcription factors. Their distribution is different in C57BL/6 (B6) and BALB/c mice, with a bias for NKT1 and NKT2/NKT17 subsets, respectively. Here, we show that the non-classical class I-like major histocompatibility complex CD1 molecules CD1d2, expressed in BALB/c and not in B6 mice, could not account for this difference. We find however that NKT cell subset distribution is intrinsic to bone marrow derived NKT cells, regardless of syngeneic CD1d-ligand recognition, and that multiple intrinsic factors are likely involved. Finally, we find that CD1d expression levels in combination with T cell antigen receptor signal strength could also influence NKT cell distribution and function. Overall, this study indicates that CD1d-mediated TCR signals and other intrinsic signals integrate to influence strain-specific NKT cell differentiation programs and subset distributions.

Modeling Brain Gene Expression in Alcohol Use Disorder with Genetic Animal Models

Animal genetic models have and will continue to provide important new information about the behavioral and physiological adaptations associated with alcohol use disorder (AUD). This chapter focuses on two models, ethanol preference and drinking in the dark (DID), their usefulness in interrogating brain gene expression data and the relevance of the data obtained to interpret AUD-related GWAS and TWAS studies. Both the animal and human data point to the importance for AUD of changes in synaptic transmission (particularly glutamate and GABA transmission), of changes in the extracellular matrix (specifically including collagens, cadherins and protocadherins) and of changes in neuroimmune processes. The implementation of new technologies (e.g., cell type-specific gene expression) is expected to further enhance the value of genetic animal models in understanding AUD.

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

Host Genetic Background Effect on Body Weight Changes Influenced by Heterozygous Smad4 Knockout Using Collaborative Cross Mouse Population

Obesity and its attendant conditions have become major health problems worldwide, and obesity is currently ranked as the fifth most common cause of death globally. Complex environmental and genetic factors are causes of the current obesity epidemic. Diet, lifestyle, chemical exposure, and other confounding factors are difficult to manage in humans. The mice model is helpful in researching genetic BW gain because genetic and environmental risk factors can be controlled in mice. Studies in mouse strains with various genetic backgrounds and established genetic structures provide unparalleled opportunities to find and analyze trait-related genomic loci. In this study, we used the Collaborative Cross (CC), a large panel of recombinant inbred mouse strains, to present a predictive study using heterozygous Smad4 knockout profiles of CC mice to understand and effectively identify predispositions to body weight gain. Male C57Bl/6J Smad4+/- mice were mated with female mice from 10 different CC lines to create F1 mice (Smad4+/-x CC). Body weight (BW) was measured weekly until week 16 and then monthly until the end of the study (week 48). The heritability (H2) of the assessed traits was estimated and presented. Comparative analysis of various machine learning algorithms for predicting the BW changes and genotype of mice was conducted. Our data showed that the body weight records of F1 mice with different CC lines differed between wild-type and mutant Smad4 mice during the experiment. Genetic background affects weight gain and some lines gained more weight in the presence of heterozygous Smad4 knockout, while others gained less, but, in general, the mutation caused overweight mice, except for a few lines. In both control and mutant groups, female %BW had a higher heritability (H2) value than males. Additionally, both sexes with wild-type genotypes showed higher heritability values than the mutant group. Logistic regression provides the most accurate mouse genotype predictions using machine learning. We plan to validate the proposed method on more CC lines and mice per line to expand the literature on machine learning for BW prediction.

Genetic architecture of the acute and persistent immune cell response after radiation exposure

Hematologic toxicity is a common side effect of multimodal cancer therapy. Nearly all animal studies investigating the causes of radiotherapy-induced hematologic toxicity use inbred strains with limited genetic diversity and do not reflect the diverse responses observed in humans. We used the population-based Collaborative Cross (CC) mouse resource to investigate the genetic architecture of the acute and persistent immune response after radiation exposure by measuring 22 immune parameters in 1,720 CC mice representing 35 strains. We determined relative acute and persistent radiation resistance scores at the individual strain level considering contributions from all immune parameters. Genome-wide association analysis identified quantitative trait loci associated with baseline and radiation responses. A cross-species radiation resistance score predicted recurrence-free survival in medulloblastoma patients. We present a community resource of immune parameters and genome-wide association analyses before and after radiation exposure for future investigations of the contributions of host genetics on radiosensitivity.

Network analysis identifies strain-dependent response to tau and tau seeding-associated genes

Previous research demonstrated that genetic heterogeneity is a critical factor in modeling amyloid accumulation and other Alzheimer's disease phenotypes. However, it is unknown what mechanisms underlie these effects of genetic background on modeling tau aggregate-driven pathogenicity. In this study, we induced tau aggregation in wild-derived mice by expressing MAPT. To investigate the effect of genetic background on the action of tau aggregates, we performed RNA sequencing with brains of C57BL/6J, CAST/EiJ, PWK/PhJ, and WSB/EiJ mice (n = 64) and determined core transcriptional signature conserved in all genetic backgrounds and signature unique to wild-derived backgrounds. By measuring tau seeding activity using the cortex, we identified 19 key genes associated with tau seeding and amyloid response. Interestingly, microglial pathways were strongly associated with tau seeding activity in CAST/EiJ and PWK/PhJ backgrounds. Collectively, our study demonstrates that mouse genetic context affects tau-mediated alteration of transcriptome and tau seeding. The gene modules associated with tau seeding provide an important resource to better model tauopathy.

Key advances in vaccine development for tuberculosis-success and challenges

Breakthrough findings in the clinical and preclinical development of tuberculosis (TB) vaccines have galvanized the field and suggest, for the first time since the development of bacille Calmette-Guérin (BCG), that a novel and protective TB vaccine is on the horizon. Here we highlight the TB vaccines that are in the development pipeline and review the basis for optimism in both the clinical and preclinical space. We describe immune signatures that could act as immunological correlates of protection (CoP) to facilitate the development and comparison of vaccines. Finally, we discuss new animal models that are expected to more faithfully model the pathology and complex immune responses observed in human populations.

Narrating the Genetic Landscape of Human Class I Occlusion: A Perspective-Infused Review

This review examines a prevalent condition with multifaceted etiology encompassing genetic, environmental, and oral behavioral factors. It stands as a significant ailment impacting oral functionality, aesthetics, and quality of life. Longitudinal studies indicate that malocclusion in primary dentition may progress to permanent malocclusion. Recognizing and managing malocclusion in primary dentition is gaining prominence. The World Health Organization ranks malocclusions as the third most widespread oral health issue globally. Angle's classification system is widely used to categorize malocclusions, with Class I occlusion considered the norm. However, its prevalence varies across populations due to genetic and examination disparities. Genetic factors, including variants in genes like MSX1, PAX9, and AXIN2, have been associated with an increased risk of Class I occlusion. This review aims to provide a comprehensive overview of clinical strategies for managing Class I occlusion and consolidate genetic insights from both human and murine populations. Additionally, genomic relationships among craniofacial genes will be assessed in individuals with Class I occlusion, along with a murine model, shedding light on phenotype-genotype associations of clinical relevance. The prevalence of Class I occlusion, its impact, and treatment approaches will be discussed, emphasizing the importance of early intervention. Additionally, the role of RNA alterations in skeletal Class I occlusion will be explored, focusing on variations in expression or structure that influence craniofacial development. Mouse models will be highlighted as crucial tools for investigating mandible size and prognathism and conducting QTL analysis to gain deeper genetic insights. This review amalgamates cellular, molecular, and clinical trait data to unravel correlations between malocclusion and Class I phenotypes.

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for a century. However, laboratory mice capture a narrow subset of the genetic variation found in wild mouse populations. This consideration inherently restricts the scope of potential discovery in laboratory models and narrows the pool of potentially identified phenotype-associated variants and pathways. Wild mouse populations are reservoirs of predicted functional and disease-associated alleles, but the sparsity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently imported, sequenced, and phenotyped a set of 11 wild-derived inbred strains developed from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from five environmentally diverse locations across North and South America: Saratoga Springs, New York, USA; Gainesville, Florida, USA; Manaus, Brazil; Tucson, Arizona, USA; and Edmonton, Alberta, Canada. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the mouse reference assembly, with 42.5% of variants in the Nachman strain genomes absent from classical inbred mouse strains. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels alone. Taken together, our work introduces a novel wild-derived inbred mouse strain resource that will enable new discoveries in basic and preclinical research. These strains are currently available through The Jackson Laboratory Repository under laboratory code NachJ.

A universal method for generating knockout mice in multiple genetic backgrounds using zygote electroporation

Genetically engineered mouse models are essential tools for understanding mammalian gene functions and disease pathogenesis. Genome editing allows the generation of these models in multiple inbred strains of mice without backcrossing. Zygote electroporation dramatically removed the barrier for introducing the CRISPR-Cas9 complex in terms of cost and labour. Here, we demonstrate that the generalised zygote electroporation method is also effective for generating knockout mice in multiple inbred strains. By combining in vitro fertilisation and electroporation, we obtained founders for knockout alleles in eight common inbred strains. Long-read sequencing analysis detected not only intended mutant alleles but also differences in read frequency of intended and unintended alleles among strains. Successful germline transmission of knockout alleles demonstrated that our approach can establish mutant mice targeting the same locus in multiple inbred strains for phenotyping analysis, contributing to reverse genetics and human disease research.

A mutation in Themis contributes to peanut-induced oral anaphylaxis in CC027 mice

Background

The development of peanut allergy is due to a combination of genetic and environmental factors, although specific genes have proven difficult to identify. Previously, we reported that peanut-sensitized CC027/GeniUnc (CC027) mice develop anaphylaxis upon oral challenge to peanut, unlike C3H/HeJ (C3H) mice.

Objective

To determine the genetic basis of orally-induced anaphylaxis to peanut in CC027 mice.

Methods

A genetic mapping population between CC027 and C3H mice was designed to identify the genetic factors that drive oral anaphylaxis. A total of 356 CC027xC3H backcrossed mice were generated, sensitized to peanut, then challenged to peanut by oral gavage. Anaphylaxis and peanut-specific IgE were quantified for all mice. T-cell phenotyping was conducted on CC027 and five additional CC strains.

Results

Anaphylaxis to peanut was absent in 77% of backcrossed mice, with 19% showing moderate anaphylaxis, and 4% having severe anaphylaxis. A total of eight genetic loci were associated with variation in response to peanut challenge, six associated with anaphylaxis (temperature decrease) and two associated with peanut-specific IgE levels. There were two major loci that impacted multiple aspects of the severity of acute anaphylaxis, at which the CC027 allele was associated with worse outcome. At one of these loci, CC027 has a private genetic variant in the Themis (thymocyte-expressed molecule involved in selection) gene. Consistent with Themis' described functions, we found that CC027 have more immature T cells with fewer CD8+, CD4+, and CD4+CD25+CD127- regulatory T cells.

Conclusion

Our results demonstrate a key role for Themis in the orally-reactive CC027 mouse model of peanut allergy.

Anterior chamber depth in mice is controlled by several quantitative trait loci

Anterior chamber depth (ACD) is a quantitative trait associated with primary angle closure glaucoma (PACG). Although ACD is highly heritable, known genetic variations explain a small fraction of the phenotypic variability. The purpose of this study was to identify additional ACD-influencing loci using strains of mice. Cohorts of 86 N2 and 111 F2 mice were generated from crosses between recombinant inbred BXD24/TyJ and wild-derived CAST/EiJ mice. Using anterior chamber optical coherence tomography, mice were phenotyped at 10-12 weeks of age, genotyped based on 93 genome-wide SNPs, and subjected to quantitative trait locus (QTL) analysis. In an analysis of ACD among all mice, six loci passed the significance threshold of p = 0.05 and persisted after multiple regression analysis. These were on chromosomes 6, 7, 11, 12, 15 and 17 (named Acdq6, Acdq7, Acdq11, Acdq12, Acdq15, and Acdq17, respectively). Our findings demonstrate a quantitative multi-genic pattern of ACD inheritance in mice and identify six previously unrecognized ACD-influencing loci. We have taken a unique approach to studying the anterior chamber depth phenotype by using mice as genetic tool to examine this continuously distributed trait.

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.

A study of the influence of genetic variance and sex on the density and thickness of the calvarial bone in collaborative cross mice

BACKGROUND

Bone microarchitecture is affected by multiple genes, each having a small effect on the external appearance. It is thus challenging to characterize the genes and their specific effect on bone thickness and porosity. The purpose of this study was to assess the heritability and the genetic variation effect, as well as the sex effect on the calvarial bone thickness (Ca.Th) and calvarial porosity (%PoV) using the Collaborative Cross (CC) mouse population.

METHODS

In the study we examined the parietal bones of 56 mice from 9 lines of CC mice. Morphometric parameters were evaluated using microcomputed tomography (μCT) and included Ca.Th and %PoV. We then evaluated heritability, genetic versus environmental variance and the sex effect for these parameters.

RESULTS

Our morphometric analysis showed that Ca.Th and %PoV are both significantly different among the CC lines with a broad sense heritability of 0.78 and 0.90, respectively. The sex effect within the lines was significant in line IL111 and showed higher values of Ca.Th and %PoV in females compared to males. In line IL19 there was a borderline sex effect in Ca.Th in which males showed higher values than females.

CONCLUSIONS

These results stress the complexity of sex and genotype interactions controlling Ca.Th and %PoV, as the skeletal sexual dimorphism was dependent on the genetic background. This study also shows that the CC population is a powerful tool for establishing the genetic effect on these traits.

Refining risk estimates for lead in drinking water based on the impact of genetics and diet on blood lead levels using the Collaborative Cross mouse population

Blood lead (Pb) level (BLL) is a commonly used biomarker to evaluate associations with health effects. However, interventions to reduce the adverse effects of Pb require relating BLL to external exposure. Moreover, risk mitigation actions need to ensure protection of more susceptible individuals with a greater tendency to accumulate Pb. Because little data is available to quantify inter-individual variability in biokinetics of Pb, we investigated the influence of genetics and diet on BLL in the genetically diverse Collaborative Cross (CC) mouse population. Adult female mice from 49 CC strains received either a standard mouse chow or a chow mimicking the American diet while being provided water ad libitum with 1000 ppm Pb for 4 weeks. In both arms of the study, inter-strain variability was observed; however, in American diet-fed animals, the BLL was greater and more variable. Importantly, the degree of variation in BLL among strains on the American diet was greater (2.3) than the default variability estimate (1.6) used in setting the regulatory standards. Genetic analysis identified suggestive diet-associated haplotypes that were associated with variation in BLL, largely contributed by the PWK/PhJ strain. This study quantified the variation in BLL that is due to genetic background, diet, and their interactions, and observed that it may be greater than that assumed for current regulatory standards for Pb in drinking water. Moreover, this work highlights the need of characterizing inter-individual variation in BLL to ensure adequate public health interventions aimed at reducing human health risks from Pb.

Using the collaborative cross to identify the role of host genetics in defining the murine gut microbiome

BACKGROUND

The human gut microbiota is a complex community comprised of trillions of bacteria and is critical for the digestion and absorption of nutrients. Bacterial communities of the intestinal microbiota influence the development of several conditions and diseases. We studied the effect of host genetics on gut microbial composition using Collaborative Cross (CC) mice. CC mice are a panel of mice that are genetically diverse across strains, but genetically identical within a given strain allowing repetition and deeper analysis than is possible with other collections of genetically diverse mice.

RESULTS

16S rRNA from the feces of 167 mice from 28 different CC strains was sequenced and analyzed using the Qiime2 pipeline. We observed a large variance in the bacterial composition across CC strains starting at the phylum level. Using bacterial composition data, we identified 17 significant Quantitative Trait Loci (QTL) linked to 14 genera on 9 different mouse chromosomes. Genes within these intervals were analyzed for significant association with pathways and the previously known human GWAS database using Enrichr analysis and Genecards database. Multiple host genes involved in obesity, glucose homeostasis, immunity, neurological diseases, and many other protein-coding genes located in these regions may play roles in determining the composition of the gut microbiota. A subset of these CC mice was infected with Salmonella Typhimurium. Using infection outcome data, an increase in abundance of genus Lachnospiraceae and decrease in genus Parasutterella correlated with positive health outcomes after infection. Machine learning classifiers accurately predicted the CC strain and the infection outcome using pre-infection bacterial composition data from the feces.

CONCLUSION

Our study supports the hypothesis that multiple host genes influence the gut microbiome composition and homeostasis, and that certain organisms may influence health outcomes after S. Typhimurium infection. Video Abstract.

Host genetic background is a barrier to broadly effective vaccine-mediated protection against tuberculosis

Heterogeneity in human immune responses is difficult to model in standard laboratory mice. To understand how host variation affects Bacillus Calmette Guerin-induced (BCG-induced) immunity against Mycobacterium tuberculosis, we studied 24 unique collaborative cross (CC) mouse strains, which differ primarily in the genes and alleles they inherit from founder strains. The CC strains were vaccinated with or without BCG and challenged with aerosolized M. tuberculosis. Since BCG protects only half of the CC strains tested, we concluded that host genetics has a major influence on BCG-induced immunity against M. tuberculosis infection, making it an important barrier to vaccine-mediated protection. Importantly, BCG efficacy is dissociable from inherent susceptibility to tuberculosis (TB). T cell immunity was extensively characterized to identify components associated with protection that were stimulated by BCG and recalled after M. tuberculosis infection. Although considerable diversity is observed, BCG has little impact on the composition of T cells in the lung after infection. Instead, variability is largely shaped by host genetics. BCG-elicited protection against TB correlated with changes in immune function. Thus, CC mice can be used to define correlates of protection and to identify vaccine strategies that protect a larger fraction of genetically diverse individuals instead of optimizing protection for a single genotype.

Immune resilience despite inflammatory stress promotes longevity and favorable health outcomes including resistance to infection

Some people remain healthier throughout life than others but the underlying reasons are poorly understood. Here we hypothesize this advantage is attributable in part to optimal immune resilience (IR), defined as the capacity to preserve and/or rapidly restore immune functions that promote disease resistance (immunocompetence) and control inflammation in infectious diseases as well as other causes of inflammatory stress. We gauge IR levels with two distinct peripheral blood metrics that quantify the balance between (i) CD8+ and CD4+ T-cell levels and (ii) gene expression signatures tracking longevity-associated immunocompetence and mortality-associated inflammation. Profiles of IR metrics in ~48,500 individuals collectively indicate that some persons resist degradation of IR both during aging and when challenged with varied inflammatory stressors. With this resistance, preservation of optimal IR tracked (i) a lower risk of HIV acquisition, AIDS development, symptomatic influenza infection, and recurrent skin cancer; (ii) survival during COVID-19 and sepsis; and (iii) longevity. IR degradation is potentially reversible by decreasing inflammatory stress. Overall, we show that optimal IR is a trait observed across the age spectrum, more common in females, and aligned with a specific immunocompetence-inflammation balance linked to favorable immunity-dependent health outcomes. IR metrics and mechanisms have utility both as biomarkers for measuring immune health and for improving health outcomes.

Lowering fasting blood glucose with non-dialyzable material of cranberry extract is dependent on host genetic background, sex and diet

BACKGROUND

Type 2 diabetes (T2D) is a polygenic metabolic disease, characterized by high fasting blood glucose (FBG). The ability of cranberry (CRN) fruit to regulate glycemia in T2D patients is well known. Here, a cohort of 13 lines of the genetically diverse Collaborative Cross (CC) mouse model was assessed for the effect of non-dialyzable material (NDM) of cranberry extract in lowering fasting blood glucose.

METHODS

Eight-week-old mice were maintained on either a standard chow diet (control group) or a high-fat diet (HFD) for 12 weeks, followed by injections of intraperitoneal (IP) NDM (50 mg/kg) per mouse, three times a week for the next 6 weeks. Absolute FBG (mg/dl) was measured bi-weekly and percentage changes in FBG (%FBG) between weeks 0 and 12 were calculated.

RESULTS

Statistical analysis showed a significant decrease in FBG between weeks 0 and 12 in male and female mice maintained on CHD. However, a non-significant increase in FBG values was observed in male and female mice maintained on HFD during the same period. Following administration of NDM during the following 6 weeks, the results show a variation in significant levels of FBG lowering between lines, male and female mice and under the different diets.

CONCLUSION

The results suggest that the efficacy of NDM treatment in lowering FGB depends on host genetic background (pharmacogenetics), sex of the mouse (pharmacosex), and diet (pharmacodiet). All these results support the need for follow-up research to better understand and implement a personalized medicine approach/utilization of NDM for reducing FBG.

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.

Resolution of structural variation in diverse mouse genomes reveals chromatin remodeling due to transposable elements

Diverse inbred mouse strains are important biomedical research models, yet genome characterization of many strains is fundamentally lacking in comparison with humans. In particular, catalogs of structural variants (SVs) (variants ≥ 50 bp) are incomplete, limiting the discovery of causative alleles for phenotypic variation. Here, we resolve genome-wide SVs in 20 genetically distinct inbred mice with long-read sequencing. We report 413,758 site-specific SVs affecting 13% (356 Mbp) of the mouse reference assembly, including 510 previously unannotated coding variants. We substantially improve the Mus musculus transposable element (TE) callset, and we find that TEs comprise 39% of SVs and account for 75% of altered bases. We further utilize this callset to investigate how TE heterogeneity affects mouse embryonic stem cells and find multiple TE classes that influence chromatin accessibility. Our work provides a comprehensive analysis of SVs found in diverse mouse genomes and illustrates the role of TEs in epigenetic differences.

Emerging role of microRNAs and long non-coding RNAs in COVID-19 with implications to therapeutics

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection which is commonly known as COVID-19 (COronaVIrus Disease 2019) has creeped into the human population taking tolls of life and causing tremendous economic crisis. It is indeed crucial to gain knowledge about their characteristics and interactions with human host cells. It has been shown that the majority of our genome consists of non-coding RNAs. Non-coding RNAs including micro RNAs (miRNAs) and long non-coding RNAs (lncRNAs) display significant roles in regulating gene expression in almost all cancers and viral diseases. It is intriguing that miRNAs and lncRNAs remarkably regulate the function and expression of major immune components of SARS-CoV-2. MiRNAs act via RNA interference mechanism in which they bind to the complementary sequences of the viral RNA strand, inducing the formation of silencing complex that eventually degrades or inhibits the viral RNA and viral protein expression. LncRNAs have been extensively shown to regulate gene expression in cytokine storm and thus emerges as a critical target for COVID-19 treatment. These lncRNAs also act as competing endogenous RNAs (ceRNAs) by sponging miRNAs and thus affecting the expression of downstream targets during SARS-CoV-2 infection. In this review, we extensively discuss the role of miRNAs and lncRNAs, describe their mechanism of action and their different interacting human targets cells during SARS-CoV-2 infection. Finally, we discuss possible ways how an interference with their molecular function could be exploited for new therapies against SARS-CoV-2.

Which mouse multiparental population is right for your study? The Collaborative Cross inbred strains, their F1 hybrids, or the Diversity Outbred population

Multiparental populations (MPPs) encompass greater genetic diversity than traditional experimental crosses of two inbred strains, enabling broader surveys of genetic variation underlying complex traits. Two such mouse MPPs are the Collaborative Cross (CC) inbred panel and the Diversity Outbred (DO) population, which are descended from the same eight inbred strains. Additionally, the F1 intercrosses of CC strains (CC-RIX) have been used and enable study designs with replicate outbred mice. Genetic analyses commonly used by researchers to investigate complex traits in these populations include characterizing how heritable a trait is, i.e. its heritability, and mapping its underlying genetic loci, i.e. its quantitative trait loci (QTLs). Here we evaluate the relative merits of these populations for these tasks through simulation, as well as provide recommendations for performing the quantitative genetic analyses. We find that sample populations that include replicate animals, as possible with the CC and CC-RIX, provide more efficient and precise estimates of heritability. We report QTL mapping power curves for the CC, CC-RIX, and DO across a range of QTL effect sizes and polygenic backgrounds for samples of 174 and 500 mice. The utility of replicate animals in the CC and CC-RIX for mapping QTLs rapidly decreased as traits became more polygenic. Only large sample populations of 500 DO mice were well-powered to detect smaller effect loci (7.5-10%) for highly complex traits (80% polygenic background). All results were generated with our R package musppr, which we developed to simulate data from these MPPs and evaluate genetic analyses from user-provided genotypes.

RSPO4 is a potential risk gene of stages III-IV, grade C periodontitis through effects on innate immune response and oral barrier integrity

AIM

R-spondin 4 (RSPO4) is a suggestive risk gene of stage III-IV, grade C periodontitis and upregulated in gingiva of mice resistant to bacteria-induced alveolar bone loss. We aimed to replicate the association, identify and characterize the putative causal variant(s) and molecular effects, and understand the downstream effects of RSPO4 upregulation.

MATERIALS AND METHODS

We performed a two-step association study for RSPO4 with imputed genotypes of a German-Dutch (896 stage III-IV, grade C periodontitis cases, 7104 controls) and Spanish sample (441 cases and 1141 controls). We analysed the allelic effects on transcription factor binding sites with reporter gene and antibody electrophoretic mobility shift assays. We used CRISPR/dCas9 activation and RNA sequencing to pinpoint RSPO4 as the target gene and to analyse downstream effects.

RESULTS

RSPO4 was associated with periodontitis (rs6056178, p_meta  = 4.6 × 10^-5 ). rs6056178 contains a GATA-binding motif. The rs6056178 T-allele abolished reporter activity (p = .004) and reduced GATA binding (-14.5%). CRISPRa of the associated region increased RSPO4 expression (25.8 ± 6.5-fold, p = .003). RSPO4 activation showed strongest induction of Gliomedin (439-fold) and Mucin 21 (178-fold) and of the gene set "response to interferon-alpha" (area under the curve [AUC] = 0.8, p < 5 × 10^-6 ). The most repressed gene set was "extracellular matrix interactions" (AUC = 0.8, p_adj  = .00016).

CONCLUSION

RSPO4 is a potential periodontitis risk gene and modifies host defence and barrier integrity.

Thirdhand tobacco smoke exposure increases the genetic background-dependent risk of pan-tumor development in Collaborative Cross mice

Increasing evidence has shown that thirdhand smoke (THS) exposure is likely to induce adverse health effects. An important knowledge gap remains in our understanding of THS exposure related to cancer risk in the human population. Population-based animal models are useful and powerful in investigating the interplay between host genetics and THS exposure on cancer risk. Here, we used the Collaborative Cross (CC) mouse population-based model system, which recapitulates the genetic and phenotypic diversity observed in the human population, to assess cancer risk after a short period of exposure, between 4 and 9 weeks of age. Eight CC strains (CC001, CC019, CC026, CC036, CC037, CC041, CC042 and CC051) were included in our study. We quantified pan-tumor incidence, tumor burden per mouse, organ tumor spectrum and tumor-free survival until 18 months of age. At the population level, we observed a significantly increased pan-tumor incidence and tumor burden per mouse in THS-treated mice as compared to the control (p = 3.04E-06). Lung and liver tissues exhibited the largest risk of undergoing tumorigenesis after THS exposure. Tumor-free survival was significantly reduced in THS-treated mice compared to control (p = 0.044). At the individual strain level, we observed a large variation in tumor incidence across the 8 CC strains. CC036 and CC041 exhibited a significant increase in pan-tumor incidence (p = 0.0084 and p = 0.000066, respectively) after THS exposure compared to control. We conclude that early-life THS exposure increases tumor development in CC mice and that host genetic background plays an important role in individual susceptibility to THS-induced tumorigenesis. Genetic background is an important factor that should be taken into account when determining human cancer risk of THS exposure.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Multi-omics analysis identifies drivers of protein phosphorylation

BACKGROUND

Phosphorylation of proteins is a key step in the regulation of many cellular processes including activation of enzymes and signaling cascades. The abundance of a phosphorylated peptide (phosphopeptide) is determined by the abundance of its parent protein and the proportion of target sites that are phosphorylated.

RESULTS

We quantified phosphopeptides, proteins, and transcripts in heart, liver, and kidney tissue samples of mice from 58 strains of the Collaborative Cross strain panel. We mapped ~700 phosphorylation quantitative trait loci (phQTL) across the three tissues and applied genetic mediation analysis to identify causal drivers of phosphorylation. We identified kinases, phosphatases, cytokines, and other factors, including both known and potentially novel interactions between target proteins and genes that regulate site-specific phosphorylation. Our analysis highlights multiple targets of pyruvate dehydrogenase kinase 1 (PDK1), a regulator of mitochondrial function that shows reduced activity in the NZO/HILtJ mouse, a polygenic model of obesity and type 2 diabetes.

CONCLUSIONS

Together, this integrative multi-omics analysis in genetically diverse CC strains provides a powerful tool to identify regulators of protein phosphorylation. The data generated in this study provides a resource for further exploration.

DNA methylation entropy is associated with DNA sequence features and developmental epigenetic divergence

Epigenetic information defines tissue identity and is largely inherited in development through DNA methylation. While studied mostly for mean differences, methylation also encodes stochastic change, defined as entropy in information theory. Analyzing allele-specific methylation in 49 human tissue sample datasets, we find that methylation entropy is associated with specific DNA binding motifs, regulatory DNA, and CpG density. Then applying information theory to 42 mouse embryo methylation datasets, we find that the contribution of methylation entropy to time- and tissue-specific patterns of development is comparable to the contribution of methylation mean, and methylation entropy is associated with sequence and chromatin features conserved with human. Moreover, methylation entropy is directly related to gene expression variability in development, suggesting a role for epigenetic entropy in developmental plasticity.

The delayed effects of acute radiation exposure (DEARE): characteristics, mechanisms, animal models, and promising medical countermeasures

PURPOSE

Terrorist use of nuclear weapons and radiation accidents put the human population at risk for exposure to life-threatening levels of radiation. Victims of lethal radiation exposure face potentially lethal acute injury, while survivors of the acute phase are plagued with chronic debilitating multi-organ injuries for years after exposure. Developing effective medical countermeasures (MCM) for the treatment of radiation exposure is an urgent need that relies heavily on studies conducted in reliable and well-characterized animal models according to the FDA Animal Rule. Although relevant animal models have been developed in several species and four MCM for treatment of the acute radiation syndrome are now FDA-approved, animal models for the delayed effects of acute radiation exposure (DEARE) have only recently been developed, and there are no licensed MCM for DEARE. Herein, we provide a review of the DEARE including key characteristics of the DEARE gleaned from human data as well as animal, mechanisms common to multi-organ DEARE, small and large animal models used to study the DEARE, and promising new or repurposed MCM under development for alleviation of the DEARE.

CONCLUSIONS

Intensification of research efforts and support focused on better understanding of mechanisms and natural history of DEARE are urgently needed. Such knowledge provides the necessary first steps toward the design and development of MCM that effectively alleviate the life-debilitating consequences of the DEARE for the benefit of humankind worldwide.

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".