B cells in perivascular and peribronchiolar granuloma-associated lymphoid tissue and B-cell signatures identify asymptomatic Mycobacterium tuberculosis lung infection in Diversity Outbred mice.

Because most humans resist Mycobacterium tuberculosis infection, there is a paucity of lung samples to study. To address this gap, we infected Diversity Outbred mice with M. tuberculosis and studied the lungs of mice in different disease states. After a low-dose aerosol infection, progressors succumbed to acute, inflammatory lung disease within 60 days, while controllers maintained asymptomatic infection for at least 60 days, and then developed chronic pulmonary tuberculosis (TB) lasting months to more than 1 year. Here, we identified features of asymptomatic M. tuberculosis infection by applying computational and statistical approaches to multimodal data sets. Cytokines and anti-M. tuberculosis cell wall antibodies discriminated progressors vs controllers with chronic pulmonary TB but could not classify mice with asymptomatic infection. However, a novel deep-learning neural network trained on lung granuloma images was able to accurately classify asymptomatically infected lungs vs acute pulmonary TB in progressors vs chronic pulmonary TB in controllers, and discrimination was based on perivascular and peribronchiolar lymphocytes. Because the discriminatory lesion was rich in lymphocytes and CD4 T cell-mediated immunity is required for resistance, we expected CD4 T-cell genes would be elevated in asymptomatic infection. However, the significantly different, highly expressed genes were from B-cell pathways (e.g., Bank1, Cd19, Cd79, Fcmr, Ms4a1, Pax5, and H2-Ob), and CD20+ B cells were enriched in the perivascular and peribronchiolar regions of mice with asymptomatic M. tuberculosis infection. Together, these results indicate that genetically controlled B-cell responses are important for establishing asymptomatic M. tuberculosis lung infection.

Mouse Genetic Reference Populations: Cellular Platforms for Integrative Systems Genetics.

Interrogation of disease-relevant cellular and molecular traits exhibited by genetically diverse cell populations enables in vitro systems genetics approaches for uncovering the basic properties of cellular function and identity. Primary cells, stem cells, and organoids derived from genetically diverse mouse strains, such as Collaborative Cross and Diversity Outbred populations, offer the opportunity for parallel in vitro/in vivo screening. These panels provide genetic resolution for variant discovery and functional characterization, as well as disease modeling and in vivo validation capabilities. Here we review mouse cellular systems genetics approaches for characterizing the influence of genetic variation on signaling networks and phenotypic diversity, and we discuss approaches for data integration and cross-species validation.

Intravenous BCG Vaccination of Diversity Outbred Mice Results in Moderately Enhanced Protection against Challenge with Mycobacterium tuberculosis Compared to Intradermal Vaccination.

Tuberculosis is still the leading cause of death globally from any infectious disease, despite the widespread use of the live attenuated vaccine Bacille Calmette Guerin (BCG). While BCG has some efficacy against disseminated TB disease in children, protection wanes into adulthood resulting in over 1.8 million TB deaths per year. This has led to efforts to develop novel vaccine candidates that either replace or boost BCG, as well as to test novel delivery mechanisms to enhance BCG's efficacy. Traditional BCG vaccination is performed as an intradermal (ID) injection but delivering BCG by an alternate route may enhance the depth and breadth of protection. Previously, we demonstrated that phenotypically and genotypically disparate Diversity Outbred (DO) mice have heterogenous responses to M. tuberculosis challenge following intradermal BCG vaccination. Here, we utilize DO mice to examine BCG-induced protection when BCG is delivered systemically via intravenous (IV) administration. We find that DO mice vaccinated with IV BCG had a greater distribution of BCG throughout their organs compared to ID-vaccinated animals. However, compared to ID-vaccinated mice, M. tuberculosis burdens in lungs and spleens were not significantly reduced in animals vaccinated with BCG IV, nor was lung inflammation significantly altered. Nonetheless, DO mice that received BCG IV had increased survival over those vaccinated by the traditional ID route. Thus, our results suggest that delivering BCG by the alternate IV route enhances protection as detected in this diverse small animal model.

Modeling human heterogeneity of obesity with diversity outbred mice reveals a fat mass-dependent therapeutic window for resolvin E1.

Resolvin E1 (RvE1), a specialized pro-resolving mediator (SPM), improves glucose homeostasis in inbred mouse models of obesity. However, an impediment toward translation is that obesity is a highly heterogenous disease in which individuals will respond very differently to interventions such as RvE1. Thus, there is a need to study SPMs in the context of modeling the heterogeneity of obesity that is observed in humans. We investigated how RvE1 controls the concentration of key circulating metabolic biomarkers using diversity outbred (DO) mice, which mimic human heterogeneity. We first demonstrate that weights of DO mice can be classified into distinct distributions of fat mass (i.e., modeling differing classes of obesity) in response to a high-fat diet and in the human population when examining body composition. Next, we show RvE1 administration based on body weight for four consecutive days after giving mice a high-fat diet led to approximately half of the mice responding positively for serum total gastric inhibitory polypeptide (GIP), glucagon, insulin, glucose, leptin, and resistin. Interestingly, RvE1 improved hyperleptinemia most effectively in the lowest class of fat mass despite adjusting the dose of RvE1 with increasing adiposity. Furthermore, leptin levels after RvE1 treatment were the lowest in those mice that were also RvE1 positive responders for insulin and resistin. Collectively, these results suggest a therapeutic fat mass-dependent window for RvE1, which should be considered in future clinical trials. Moreover, the data underscore the importance of studying SPMs with heterogenous mice as a step toward precision SPM administration in humans.

The genetic basis of neurocranial size and shape across varied lab mouse populations.

Brain and skull tissues interact through molecular signalling and mechanical forces during head development, leading to a strong correlation between the neurocranium and the external brain surface. Therefore, when brain tissue is unavailable, neurocranial endocasts are often used to approximate brain size and shape. Evolutionary changes in brain morphology may have resulted in secondary changes to neurocranial morphology, but the developmental and genetic processes underlying this relationship are not well understood. Using automated phenotyping methods, we quantified the genetic basis of endocast variation across large genetically varied populations of laboratory mice in two ways: (1) to determine the contributions of various genetic factors to neurocranial form and (2) to help clarify whether a neurocranial variation is based on genetic variation that primarily impacts bone development or on genetic variation that primarily impacts brain development, leading to secondary changes in bone morphology. Our results indicate that endocast size is highly heritable and is primarily determined by additive genetic factors. In addition, a non-additive inbreeding effect led to founder strains with lower neurocranial size, but relatively large brains compared to skull size; suggesting stronger canalization of brain size and/or a general allometric effect. Within an outbred sample of mice, we identified a locus on mouse chromosome 1 that is significantly associated with variation in several positively correlated endocast size measures. Because the protein-coding genes at this locus have been previously associated with brain development and not with bone development, we propose that genetic variation at this locus leads primarily to variation in brain volume that secondarily leads to changes in neurocranial globularity. We identify a strain-specific missense mutation within Akt3 that is a strong causal candidate for this genetic effect. Whilst it is not appropriate to generalize our hypothesis for this single locus to all other loci that also contribute to the complex trait of neurocranial skull morphology, our results further reveal the genetic basis of neurocranial variation and highlight the importance of the mechanical influence of brain growth in determining skull morphology.

Genome-wide association mapping of ethanol sensitivity in the Diversity Outbred mouse population.

BACKGROUND: A strong predictor for the development of alcohol use disorder (AUD) is altered sensitivity to the intoxicating effects of alcohol. Individual differences in the initial sensitivity to alcohol are controlled in part by genetic factors. Mice offer a powerful tool to elucidate the genetic basis of behavioral and physiological traits relevant to AUD, but conventional experimental crosses have only been able to identify large chromosomal regions rather than specific genes. Genetically diverse, highly recombinant mouse populations make it possible to observe a wider range of phenotypic variation, offer greater mapping precision, and thus increase the potential for efficient gene identification. METHODS: We have taken advantage of the Diversity Outbred (DO) mouse population to identify and precisely map quantitative trait loci (QTL) associated with ethanol sensitivity. We phenotyped 798 male J:DO mice for three measures of ethanol sensitivity: ataxia, hypothermia, and loss of the righting response. We used high-density MegaMUGA and GigaMUGA to obtain genotypes ranging from 77,808 to 143,259 SNPs. We also performed RNA sequencing in striatum to map expression QTLs and identify gene expression-trait correlations. We then applied a systems genetic strategy to identify narrow QTLs and construct the network of correlations that exists between DNA sequence, gene expression values, and ethanol-related phenotypes to prioritize our list of positional candidate genes. RESULTS: We observed large amounts of phenotypic variation with the DO population and identified suggestive and significant QTLs associated with ethanol sensitivity on chromosomes 1, 2, and 16. The implicated regions were narrow (4.5-6.9 Mb in size) and each QTL explained ~4-5% of the variance. CONCLUSIONS: Our results can be used to identify alleles that contribute to AUD in humans, elucidate causative biological mechanisms, or assist in the development of novel therapeutic interventions.

Genes, other than Muc5b, play a role in bleomycin-induced lung fibrosis.

Idiopathic pulmonary fibrosis (IPF) is an incurable genetic disease that affects 5 million people worldwide. The gain-of-function MUC5B promoter variant rs35705950 is the dominant genetic risk factor for IPF, yet has a low penetrance. This raises the possibility that other genes and transcripts affect the penetrance of MUC5B. Previously, we have shown that the concentration of Muc5b in bronchoalveolar epithelia is directly associated with the extent and persistence of bleomycin-induced lung fibrosis in mice. In this study, we investigated whether bleomycin-induced lung injury is Muc5b dependent in genetically divergent strains of mice. Specifically, mice from the eight Diversity Outbred (DO) founders were phenotyped for Muc5b expression and lung fibrosis 3 wk after intratracheal bleomycin administration. Although we identified strains with low Muc5b expression and minimal lung fibrosis (CAST/EiJ and PWK/PhJ) and strains with high Muc5b expression and extensive lung fibrosis (NZO/H1LtJ and WSB/EiJ), there also were strains that did not demonstrate a clear relationship between Muc5b expression and lung fibrosis (129S1/SvlmJ, NOD/ShiLtJ, and C57BL/6J, A/J). Hierarchical clustering suggests that other factors may work in concert with or potentially independent of Muc5b to promote bleomycin-induced lung injury and fibrosis. This study suggests that these strains and their recombinant inbred crosses may prove helpful in identifying the genes and transcripts that interact with Muc5b and cause lung fibrosis.

Resolvin E1 derived from eicosapentaenoic acid prevents hyperinsulinemia and hyperglycemia in a host genetic manner.

Eicosapentaenoic acid (EPA) has garnered attention after the success of the REDUCE-IT trial, which contradicted previous conclusions on EPA for cardiovascular disease risk. Here we first investigated EPA's preventative role on hyperglycemia and hyperinsulinemia. EPA ethyl esters prevented obesity-induced glucose intolerance, hyperinsulinemia, and hyperglycemia in C57BL/6J mice. Supporting NHANES analyses showed that fasting glucose levels of obese adults were inversely related to EPA intake. We next investigated how EPA improved murine hyperinsulinemia and hyperglycemia. EPA overturned the obesity-driven decrement in the concentration of 18-hydroxyeicosapentaenoic acid (18-HEPE) in white adipose tissue and liver. Treatment of obese inbred mice with RvE1, the downstream immunoresolvant metabolite of 18-HEPE, but not 18-HEPE itself, reversed hyperinsulinemia and hyperglycemia through the G-protein coupled receptor ERV1/ChemR23. To translate the findings, we determined if the effects of RvE1 were dependent on host genetics. RvE1's effects on hyperinsulinemia and hyperglycemia were divergent in diversity outbred mice that model human genetic variation. Secondary SNP analyses further confirmed extensive genetic variation in human RvE1/EPA-metabolizing genes. Collectively, the data suggest EPA prevents hyperinsulinemia and hyperglycemia, in part, through RvE1's activation of ERV1/ChemR23 in a host genetic manner. The studies underscore the need for personalized administration of RvE1 based on genetic/metabolic enzyme profiles.

microRNA-146a-5p association with the cardiometabolic disease risk factor TMAO.

Trimethylamine-N-oxide (TMAO), a microbial choline metabolism byproduct that is processed in the liver and excreted into circulation, is associated with increased atherosclerotic lesion formation and cardiovascular disease risk. Genetic regulators of TMAO levels are largely unknown. In the present study, we used 288 mice from a genetically heterogeneous mouse population [Diversity Outbred (DO)] to determine hepatic microRNA associations with TMAO in the context of an atherogenic diet. We also validated findings in two additional animal models of atherosclerosis: liver-specific insulin receptor knockout mice fed a chow diet (LIRKO) and African green monkeys fed high-fat/high-cholesterol diet. Small RNA-sequencing analysis in DO mice, LIRKO mice, and African green monkeys identified only one hepatic microRNA (miR-146a-5p) that is aberrantly expressed across all three models. Moreover, miR-146a-5p levels are associated with circulating TMAO after atherogenic diet in each of these models. We also performed high-resolution genetic mapping and identified a novel quantitative trait locus on Chromosome 12 for TMAO levels. This interval includes two genes, Numb and Dlst, which are inversely correlated with both miR-146a and TMAO and are predicted targets of miR-146a. Both of these genes have been validated as direct targets of miR-146a, though in other cellular contexts. This is the first report to our knowledge of a link between miR-146 and TMAO. Our findings suggest that miR-146-5p, as well as one or more genes at the Chromosome 12 QTL (possibly Numb or Dlst), is strongly linked to TMAO levels and likely involved in the control of atherosclerosis.

Systems genetics identifies a co-regulated module of liver microRNAs associated with plasma LDL cholesterol in murine diet-induced dyslipidemia.

Chronically altered levels of circulating lipids, termed dyslipidemia, is a significant risk factor for a number of metabolic and cardiovascular morbidities. MicroRNAs (miRNAs) have emerged as important regulators of lipid balance, have been implicated in dyslipidemia, and have been proposed as candidate therapeutic targets in lipid-related disorders including atherosclerosis. A major limitation of most murine studies of miRNAs in lipid metabolic disorders is that they have been performed in just one (or very few) inbred strains, such as C57BL/6. Moreover, although individual miRNAs have been associated with lipid phenotypes, it is well understood that miRNAs likely work together in functional modules. To address these limitations, we implemented a systems genetics strategy using the Diversity Outbred (DO) mouse population. Specifically, we performed gene and miRNA expression profiling in the livers from ~300 genetically distinct DO mice after 18 wk on either a high-fat/high-cholesterol diet or a high-protein diet. Large-scale correlative analysis of these data with a wide range of cardio-metabolic end points revealed a co-regulated module of miRNAs significantly associated with circulating low-density lipoprotein cholesterol (LDL-C) levels. The hubs of this module were identified as miR-199a, miR-181b, miR-27a, miR-21_-_1, and miR-24. In sum, we demonstrate that a high-fat/high-cholesterol diet robustly rewires the miRNA regulatory network, and we identify a small group of co-regulated miRNAs that may exert coordinated effects to control circulating LDL-C.

QTL mapping in outbred populations: successes and challenges.

Quantitative trait locus (QTL) mapping in animal populations has been a successful strategy for identifying genomic regions that play a role in complex diseases and traits. When conducted in an F2 intercross or backcross population, the resulting QTL is frequently large, often encompassing 30 Mb or more and containing hundreds of genes. To narrow the locus and identify candidate genes, additional strategies are needed. Congenic strains have proven useful but work less well when there are multiple tightly linked loci, frequently resulting in loss of phenotype. As an alternative, we discuss the use of highly recombinant outbred models for directly fine-mapping QTL to only a few megabases. We discuss the use of several currently available models such as the advanced intercross (AI), heterogeneous stocks (HS), the diversity outbred (DO), and commercially available outbred stocks (CO). Once a QTL has been fine-mapped, founder sequence and expression QTL mapping can be used to identify candidate genes. In this regard, the large number of alleles found in outbred stocks can be leveraged to identify causative genes and variants. We end this review by discussing some important statistical considerations when analyzing outbred populations. Fine-resolution mapping in outbred models, coupled with full genome sequence, has already led to the identification of several underlying causative genes for many complex traits and diseases. These resources will likely lead to additional successes in the coming years.

Genetically diverse mice are novel and valuable models of age-associated susceptibility to Mycobacterium tuberculosis.

BACKGROUND: Tuberculosis, the disease due to Mycobacterium tuberculosis, is an important cause of morbidity and mortality in the elderly. Use of mouse models may accelerate insight into the disease and tests of therapies since mice age thirty times faster than humans. However, the majority of TB research relies on inbred mouse strains, and these results might not extrapolate well to the genetically diverse human population. We report here the first tests of M. tuberculosis infection in genetically heterogeneous aging mice, testing if old mice benefit from rapamycin. FINDINGS: We find that genetically diverse aging mice are much more susceptible than young mice to M. tuberculosis, as are aging human beings. We also find that rapamycin boosts immune responses during primary infection but fails to increase survival. CONCLUSIONS: Genetically diverse mouse models provide a valuable resource to study how age influences responses and susceptibility to pathogens and to test interventions. Additionally, surrogate markers such as immune measures may not predict whether interventions improve survival.

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.

Host Genetic Variation Impacts SARS-CoV-2 Vaccination Response in the Diversity Outbred Mouse Population.

The COVID-19 pandemic led to the rapid and worldwide development of highly effective vaccines against SARS-CoV-2. However, there is significant individual-to-individual variation in vaccine efficacy due to factors including viral variants, host age, immune status, environmental and host genetic factors. Understanding those determinants driving this variation may inform the development of more broadly protective vaccine strategies. While host genetic factors are known to impact vaccine efficacy for respiratory pathogens such as influenza and tuberculosis, the impact of host genetic variation on vaccine efficacy against COVID-19 is not well understood. To model the impact of host genetic variation on SARS-CoV-2 vaccine efficacy, while controlling for the impact of non-genetic factors, we used the Diversity Outbred (DO) mouse model. We found that DO mice immunized against SARS-CoV-2 exhibited high levels of variation in vaccine-induced neutralizing antibody responses. While the majority of the vaccinated mice were protected from virus-induced disease, similar to human populations, we observed vaccine breakthrough in a subset of mice. Importantly, we found that this variation in neutralizing antibody, virus-induced disease, and viral titer is heritable, indicating that the DO serves as a useful model system for studying the contribution of genetic variation of both vaccines and disease outcomes.

Mt10 Vaccine Protects Diversity Outbred Mice from CVB3 Infection by Producing Virus-Specific Neutralizing Antibodies and Diverse Antibody Isotypes.

Group B coxsackieviruses (CVBs) cause a wide range of diseases in humans, but no vaccines are currently available to prevent these infections. Previously, we had demonstrated that a live attenuated CVB3 vaccine virus, Mutant 10 (Mt10), offers protection against multiple CVB serotypes as evaluated in various inbred mouse strains; however, the applicability of these findings to the outbred human population remains uncertain. To address this issue, we used Diversity Outbred (DO) mice, whose genome is derived from eight inbred mouse strains that may capture the level of genetic diversity of the outbred human population. To determine the efficacy of the Mt10 vaccine, we established the CVB3 infection model in the DO mice. We noted that CVB3 infection resulted mainly in pancreatitis, although viral RNA was detected in both the pancreas and heart. Histologically, the pancreatic lesions comprised of necrosis, post-necrotic atrophy, and lymphocyte infiltration. In evaluating the efficacy of the Mt10 vaccine, both male and female DO mice were completely protected in challenge studies with CVB3, and viral RNA was not detected in the heart or pancreas. Likewise, vaccine recipients of both sexes showed significant levels of virus-neutralizing antibodies. Furthermore, by using the CVB3 viral protein 1, virus-reactive antibodies were found to be diverse in the order of IgG2c, followed by IgG2a, IgG2b/IgG3, and IgG1. Together, the data suggest that the Mt10 vaccine virus can offer protection against CVB infections that may have translational significance.

Host Genetic Background Influences BCG-Induced Antibodies Cross-Reactive to SARS-CoV-2 Spike Protein.

Mycobacterium bovis Bacillus Calmette-Guérin (BCG) protects against childhood tuberculosis; and unlike most vaccines, BCG broadly impacts immunity to other pathogens and even some cancers. Early in the COVID-19 pandemic, epidemiological studies identified a protective association between BCG vaccination and outcomes of SARS-CoV-2, but the associations in later studies were inconsistent. We sought possible reasons and noticed the study populations often lived in the same country. Since individuals from the same regions can share common ancestors, we hypothesized that genetic background could influence associations between BCG and SARS-CoV-2. To explore this hypothesis in a controlled environment, we performed a pilot study using Diversity Outbred mice. First, we identified amino acid sequences shared by BCG and SARS-CoV-2 spike protein. Next, we tested for IgG reactive to spike protein from BCG-vaccinated mice. Sera from some, but not all, BCG-vaccinated Diversity Outbred mice contained higher levels of IgG cross-reactive to SARS-CoV-2 spike protein than sera from BCG-vaccinated C57BL/6J inbred mice and unvaccinated mice. Although larger experimental studies are needed to obtain mechanistic insight, these findings suggest that genetic background may be an important variable contributing to different associations observed in human randomized clinical trials evaluating BCG vaccination on SARS-CoV-2 and COVID-19.

Biological and genetic determinants of glycolysis: Phosphofructokinase isoforms boost energy status of stored red blood cells and transfusion outcomes.

Mature red blood cells (RBCs) lack mitochondria and thus exclusively rely on glycolysis to generate adenosine triphosphate (ATP) during aging in vivo or storage in blood banks. Here, we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify associations between end-of-storage levels of glycolytic metabolites and donor age, sex, and ancestry-specific genetic polymorphisms in regions encoding phosphofructokinase 1, platelet (detected in mature RBCs); hexokinase 1 (HK1); and ADP-ribosyl cyclase 1 and 2 (CD38/BST1). Gene-metabolite associations were validated in fresh and stored RBCs from 525 Diversity Outbred mice and via multi-omics characterization of 1,929 samples from 643 human RBC units during storage. ATP and hypoxanthine (HYPX) levels-and the genetic traits linked to them-were associated with hemolysis in vitro and in vivo, both in healthy autologous transfusion recipients and in 5,816 critically ill patients receiving heterologous transfusions, suggesting their potential as markers to improve transfusion outcomes.

Voluntary nicotine consumption and reward in a subset of diversity outbred founder strains.

BACKGROUND: Nicotine is largely responsible for the initiation and maintenance of tobacco dependence and contributes to a global health problem. AIMS: This study characterizes nicotine oral consumption and preference in male and female mice of several Diversity Outbred (DO) founder strains: C57BL/6J, A/J, 129S1/SvImJ, PWK/PhJ, NOD/ShiLtJ, and CAST/EiJ. It assesses the impact of nicotine concentration on intake and preference, the potential interaction of strain with sex, and estimates the degree of heritable variation in nicotine consumption. METHODS: Two-bottle choice oral self-administration paradigm was used to assess nicotine intake, nicotine preference, and total fluid intake in male and female mice of each strain in a concentration-response manner. A conditioned place preference (CPP) test was performed to evaluate the rewarding and aversive effects of nicotine in certain strains after systemic administration of the drug. RESULTS: The highest nicotine-consuming strain was found to be 129S1/SvlmJ, and the lowest nicotine-consuming strain was A/J. Strain differences in nicotine intake were not due to differences in bitter and sweet tastes as shown in the saccharine and quinine two-bottle choice tests. A/J strain showed no significant CPP for nicotine while the 129S1/SvImJ strain showed a significant CPP for nicotine and a higher preference when compared to the C57BL/6J strain. Heritability estimates of nicotine intake were sex dependent and concentration dependent. CONCLUSIONS: Data support that nicotine consumption patterns are heritable with an influence of genotype in a voluntary oral self-administration paradigm. Results pave the way for future studies with the highly recombinant DO mice that might lead to the identification of novel genetic loci and genes influencing nicotine consumption.

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.