Forward genetic screen of homeostatic antibody levels in the Collaborative Cross identifies MBD1 as a novel regulator of B cell homeostasis.

Variation in immune homeostasis, the state in which the immune system is maintained in the absence of stimulation, is highly variable across populations. This variation is attributed to both genetic and environmental factors. However, the identity and function of specific regulators have been difficult to identify in humans. We evaluated homeostatic antibody levels in the serum of the Collaborative Cross (CC) mouse genetic reference population. We found heritable variation in all antibody isotypes and subtypes measured. We identified 4 quantitative trait loci (QTL) associated with 3 IgG subtypes: IgG1, IgG2b, and IgG2c. While 3 of these QTL map to genome regions of known immunological significance (major histocompatibility and immunoglobulin heavy chain locus), Qih1 (associated with variation in IgG1) mapped to a novel locus on Chromosome 18. We further associated this locus with B cell proportions in the spleen and identify Methyl-CpG binding domain protein 1 under this locus as a novel regulator of homeostatic IgG1 levels in the serum and marginal zone B cells (MZB) in the spleen, consistent with a role in MZB differentiation to antibody secreting cells.

An interaction of inorganic arsenic exposure with body weight and composition on type 2 diabetes indicators in Diversity Outbred mice.

Type 2 diabetes (T2D) is a complex metabolic disorder with no cure and high morbidity. Exposure to inorganic arsenic (iAs), a ubiquitous environmental contaminant, is associated with increased T2D risk. Despite growing evidence linking iAs exposure to T2D, the factors underlying inter-individual differences in susceptibility remain unclear. This study examined the interaction between chronic iAs exposure and body composition in a cohort of 75 Diversity Outbred mice. The study design mimics that of an exposed human population where the genetic diversity of the mice provides the variation in response, in contrast to a design that includes untreated mice. Male mice were exposed to iAs in drinking water (100 ppb) for 26 weeks. Metabolic indicators used as diabetes surrogates included fasting blood glucose and plasma insulin (FBG, FPI), blood glucose and plasma insulin 15 min after glucose challenge (BG15, PI15), homeostatic model assessment for [Formula: see text]-cell function and insulin resistance (HOMA-B, HOMA-IR), and insulinogenic index. Body composition was determined using magnetic resonance imaging, and the concentrations of iAs and its methylated metabolites were measured in liver and urine. Associations between cumulative iAs consumption and FPI, PI15, HOMA-B, and HOMA-IR manifested as significant interactions between iAs and body weight/composition. Arsenic speciation analyses in liver and urine suggest little variation in the mice's ability to metabolize iAs. The observed interactions accord with current research aiming to disentangle the effects of multiple complex factors on T2D risk, highlighting the need for further research on iAs metabolism and its consequences in genetically diverse mouse strains.

Collaborative Cross mice reveal extreme epilepsy phenotypes and genetic loci for seizure susceptibility.

OBJECTIVE: Animal studies remain essential for understanding mechanisms of epilepsy and identifying new therapeutic targets. However, existing animal models of epilepsy do not reflect the high level of genetic diversity found in the human population. The Collaborative Cross (CC) population is a genetically diverse recombinant inbred panel of mice. The CC offers large genotypic and phenotypic diversity, inbred strains with stable genomes that allow for repeated phenotypic measurements, and genomic tools including whole genome sequence to identify candidate genes and candidate variants. METHODS: We evaluated multiple complex epileptic traits in a sampling of 35 CC inbred strains using the flurothyl-induced seizure and kindling paradigm. We created an F2 population of 297 mice with extreme seizure susceptibility and performed quantitative trait loci (QTL) mapping to identify genomic regions associated with seizure sensitivity. We used quantitative RNA sequencing from CC hippocampal tissue to identify candidate genes and whole genome sequence to identify genetic variants likely affecting gene expression. RESULTS: We identified new mouse models with extreme seizure susceptibility, seizure propagation, epileptogenesis, and SUDEP (sudden unexpected death in epilepsy). We performed QTL mapping and identified one known and seven novel loci associated with seizure sensitivity. We combined whole genome sequencing and hippocampal gene expression to pinpoint biologically plausible candidate genes (eg, Gabra2) and variants associated with seizure sensitivity. SIGNIFICANCE: New mouse models of epilepsy are needed to better understand the complex genetic architecture of seizures and to identify therapeutics. We performed a phenotypic screen utilizing a novel genetic reference population of CC mice. The data we provide enable the identification of protective/risk genes and novel molecular mechanisms linked to complex seizure traits that are currently challenging to study and treat.

Effects of Preconception and in Utero Inorganic Arsenic Exposure on the Metabolic Phenotype of Genetically Diverse Collaborative Cross Mice.

In humans and mice, in utero exposure to inorganic arsenic (iAs) is associated with adverse health outcomes later in life. The contribution of preconception exposure to the adverse outcomes in offspring has never been studied. Here combined in utero and postnatal exposures produce insulin resistance in two collaborative cross strains. Furthermore, combined preconception and in utero exposure resulted in increased birth weight and developed insulin resistance in one strain. Thus, preconception exposure to arsenic may contribute to the metabolic disorders later in life, but the susceptibility to the effects of this exposure is determined, at least in part, by genetics.

Genome Wide Identification of SARS-CoV Susceptibility Loci Using the Collaborative Cross.

New systems genetics approaches are needed to rapidly identify host genes and genetic networks that regulate complex disease outcomes. Using genetically diverse animals from incipient lines of the Collaborative Cross mouse panel, we demonstrate a greatly expanded range of phenotypes relative to classical mouse models of SARS-CoV infection including lung pathology, weight loss and viral titer. Genetic mapping revealed several loci contributing to differential disease responses, including an 8.5Mb locus associated with vascular cuffing on chromosome 3 that contained 23 genes and 13 noncoding RNAs. Integrating phenotypic and genetic data narrowed this region to a single gene, Trim55, an E3 ubiquitin ligase with a role in muscle fiber maintenance. Lung pathology and transcriptomic data from mice genetically deficient in Trim55 were used to validate its role in SARS-CoV-induced vascular cuffing and inflammation. These data establish the Collaborative Cross platform as a powerful genetic resource for uncovering genetic contributions of complex traits in microbial disease severity, inflammation and virus replication in models of outbred populations.

A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2.

Significant departures from expected Mendelian inheritance ratios (transmission ratio distortion, TRD) are frequently observed in both experimental crosses and natural populations. TRD on mouse Chromosome (Chr) 2 has been reported in multiple experimental crosses, including the Collaborative Cross (CC). Among the eight CC founder inbred strains, we found that Chr 2 TRD was exclusive to females that were heterozygous for the WSB/EiJ allele within a 9.3 Mb region (Chr 2 76.9 - 86.2 Mb). A copy number gain of a 127 kb-long DNA segment (designated as responder to drive, R2d) emerged as the strongest candidate for the causative allele. We mapped R2d sequences to two loci within the candidate interval. R2d1 is located near the proximal boundary, and contains a single copy of R2d in all strains tested. R2d2 maps to a 900 kb interval, and the number of R2d copies varies from zero in classical strains (including the mouse reference genome) to more than 30 in wild-derived strains. Using real-time PCR assays for the copy number, we identified a mutation (R2d2WSBdel1) that eliminates the majority of the R2d2WSB copies without apparent alterations of the surrounding WSB/EiJ haplotype. In a three-generation pedigree segregating for R2d2WSBdel1, the mutation is transmitted to the progeny and Mendelian segregation is restored in females heterozygous for R2d2WSBdel1, thus providing direct evidence that the copy number gain is causal for maternal TRD. We found that transmission ratios in R2d2WSB heterozygous females vary between Mendelian segregation and complete distortion depending on the genetic background, and that TRD is under genetic control of unlinked distorter loci. Although the R2d2WSB transmission ratio was inversely correlated with average litter size, several independent lines of evidence support the contention that female meiotic drive is the cause of the distortion. We discuss the implications and potential applications of this novel meiotic drive system.

R2d2 Drives Selfish Sweeps in the House Mouse.

A selective sweep is the result of strong positive selection driving newly occurring or standing genetic variants to fixation, and can dramatically alter the pattern and distribution of allelic diversity in a population. Population-level sequencing data have enabled discoveries of selective sweeps associated with genes involved in recent adaptations in many species. In contrast, much debate but little evidence addresses whether "selfish" genes are capable of fixation-thereby leaving signatures identical to classical selective sweeps-despite being neutral or deleterious to organismal fitness. We previously described R2d2, a large copy-number variant that causes nonrandom segregation of mouse Chromosome 2 in females due to meiotic drive. Here we show population-genetic data consistent with a selfish sweep driven by alleles of R2d2 with high copy number (R2d2(HC)) in natural populations. We replicate this finding in multiple closed breeding populations from six outbred backgrounds segregating for R2d2 alleles. We find that R2d2(HC) rapidly increases in frequency, and in most cases becomes fixed in significantly fewer generations than can be explained by genetic drift. R2d2(HC) is also associated with significantly reduced litter sizes in heterozygous mothers, making it a true selfish allele. Our data provide direct evidence of populations actively undergoing selfish sweeps, and demonstrate that meiotic drive can rapidly alter the genomic landscape in favor of mutations with neutral or even negative effects on overall Darwinian fitness. Further study will reveal the incidence of selfish sweeps, and will elucidate the relative contributions of selfish genes, adaptation and genetic drift to evolution.

On the subspecific origin of the laboratory mouse.

The genome of the laboratory mouse is thought to be a mosaic of regions with distinct subspecific origins. We have developed a high-resolution map of the origin of the laboratory mouse by generating 25,400 phylogenetic trees at 100-kb intervals spanning the genome. On average, 92% of the genome is of Mus musculus domesticus origin, and the distribution of diversity is markedly nonrandom among the chromosomes. There are large regions of extremely low diversity, which represent blind spots for studies of natural variation and complex traits, and hot spots of diversity. In contrast with the mosaic model, we found that most of the genome has intermediate levels of variation of intrasubspecific origin. Finally, mouse strains derived from the wild that are supposed to represent different mouse subspecies show substantial intersubspecific introgression, which has strong implications for evolutionary studies that assume these are pure representatives of a given subspecies.

Modeling host genetic regulation of influenza pathogenesis in the collaborative cross.

Genetic variation contributes to host responses and outcomes following infection by influenza A virus or other viral infections. Yet narrow windows of disease symptoms and confounding environmental factors have made it difficult to identify polymorphic genes that contribute to differential disease outcomes in human populations. Therefore, to control for these confounding environmental variables in a system that models the levels of genetic diversity found in outbred populations such as humans, we used incipient lines of the highly genetically diverse Collaborative Cross (CC) recombinant inbred (RI) panel (the pre-CC population) to study how genetic variation impacts influenza associated disease across a genetically diverse population. A wide range of variation in influenza disease related phenotypes including virus replication, virus-induced inflammation, and weight loss was observed. Many of the disease associated phenotypes were correlated, with viral replication and virus-induced inflammation being predictors of virus-induced weight loss. Despite these correlations, pre-CC mice with unique and novel disease phenotype combinations were observed. We also identified sets of transcripts (modules) that were correlated with aspects of disease. In order to identify how host genetic polymorphisms contribute to the observed variation in disease, we conducted quantitative trait loci (QTL) mapping. We identified several QTL contributing to specific aspects of the host response including virus-induced weight loss, titer, pulmonary edema, neutrophil recruitment to the airways, and transcriptional expression. Existing whole-genome sequence data was applied to identify high priority candidate genes within QTL regions. A key host response QTL was located at the site of the known anti-influenza Mx1 gene. We sequenced the coding regions of Mx1 in the eight CC founder strains, and identified a novel Mx1 allele that showed reduced ability to inhibit viral replication, while maintaining protection from weight loss.

Genetic analysis of complex traits in the emerging Collaborative Cross.

The Collaborative Cross (CC) is a mouse recombinant inbred strain panel that is being developed as a resource for mammalian systems genetics. Here we describe an experiment that uses partially inbred CC lines to evaluate the genetic properties and utility of this emerging resource. Genome-wide analysis of the incipient strains reveals high genetic diversity, balanced allele frequencies, and dense, evenly distributed recombination sites-all ideal qualities for a systems genetics resource. We map discrete, complex, and biomolecular traits and contrast two quantitative trait locus (QTL) mapping approaches. Analysis based on inferred haplotypes improves power, reduces false discovery, and provides information to identify and prioritize candidate genes that is unique to multifounder crosses like the CC. The number of expression QTLs discovered here exceeds all previous efforts at eQTL mapping in mice, and we map local eQTL at 1-Mb resolution. We demonstrate that the genetic diversity of the CC, which derives from random mixing of eight founder strains, results in high phenotypic diversity and enhances our ability to map causative loci underlying complex disease-related traits.

The Collaborative Cross as a resource for modeling human disease: CC011/Unc, a new mouse model for spontaneous colitis.

Inflammatory bowel disease (IBD) is an immune-mediated condition driven by improper responses to intestinal microflora in the context of environmental and genetic background. GWAS in humans have identified many loci associated with IBD, but animal models are valuable for dissecting the underlying molecular mechanisms, characterizing environmental and genetic contributions and developing treatments. Mouse models rely on interventions such as chemical treatment or introduction of an infectious agent to induce disease. Here, we describe a new model for IBD in which the disease develops spontaneously in 20-week-old mice in the absence of known murine pathogens. The model is part of the Collaborative Cross and came to our attention due to a high incidence of rectal prolapse in an incompletely inbred line. Necropsies revealed a profound proliferative colitis with variable degrees of ulceration and vasculitis, splenomegaly and enlarged mesenteric lymph nodes with no discernible anomalies of other organ systems. Phenotypic characterization of the CC011/Unc mice with homozygosity ranging from 94.1 to 99.8% suggested that the trait was fixed and acted recessively in crosses to the colitis-resistant C57BL/6J inbred strain. Using a QTL approach, we identified four loci, Ccc1, Ccc2, Ccc3 and Ccc4 on chromosomes 12, 14, 1 and 8 that collectively explain 27.7% of the phenotypic variation. Surprisingly, we also found that minute levels of residual heterozygosity in CC011/Unc have significant impact on the phenotype. This work demonstrates the utility of the CC as a source of models of human disease that arises through new combinations of alleles at susceptibility loci.

Sarbecovirus disease susceptibility is conserved across viral and host models.

Coronaviruses have caused three severe epidemics since the start of the 21st century: SARS, MERS and COVID-19. The severity of the ongoing COVID-19 pandemic and increasing likelihood of future coronavirus outbreaks motivates greater understanding of factors leading to severe coronavirus disease. We screened ten strains from the Collaborative Cross mouse genetic reference panel and identified strains CC006/TauUnc (CC006) and CC044/Unc (CC044) as coronavirus-susceptible and resistant, respectively, as indicated by variable weight loss and lung congestion scores four days post-infection. We generated a genetic mapping population of 755 CC006xCC044 F2 mice and exposed the mice to one of three genetically distinct mouse-adapted coronaviruses: clade 1a SARS-CoV MA15 (n=391), clade 1b SARS-CoV-2 MA10 (n=274), and clade 2 HKU3-CoV MA (n=90). Quantitative trait loci (QTL) mapping in SARS-CoV MA15- and SARS-CoV-2 MA10-infected F2 mice identified genetic loci associated with disease severity. Specifically, we identified seven loci associated with variation in outcome following infection with either virus, including one, HrS43, that is present in both groups. Three of these QTL, including HrS43, were also associated with HKU3-CoV MA outcome. HrS43 overlaps with a QTL previously reported by our lab that is associated with SARS-CoV MA15 outcome in CC011xCC074 F2 mice and is also syntenic with a human chromosomal region associated with severe COVID-19 outcomes in humans GWAS. The results reported here provide: (a) additional support for the involvement of this locus in SARS-CoV MA15 infection, (b) the first conclusive evidence that this locus is associated with susceptibility across the Sarbecovirus subgenus, and (c) demonstration of the relevance of mouse models in the study of coronavirus disease susceptibility in humans.

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.

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.

Genetic loci regulate Sarbecovirus pathogenesis: A comparison across mice and humans.

Coronavirus (CoV) cause considerable morbidity and mortality in humans and other mammals, as evidenced by the emergence of Severe Acute Respiratory CoV (SARS-CoV) in 2003, Middle East Respiratory CoV (MERS-CoV) in 2012, and SARS-CoV-2 in 2019. Although poorly characterized, natural genetic variation in human and other mammals modulate virus pathogenesis, as reflected by the spectrum of clinical outcomes ranging from asymptomatic infections to lethal disease. Using multiple human epidemic and zoonotic Sarbecoviruses, coupled with murine Collaborative Cross genetic reference populations, we identify several dozen quantitative trait loci that regulate SARS-like group-2B CoV pathogenesis and replication. Under a Chr4 QTL, we deleted a candidate interferon stimulated gene, Trim14 which resulted in enhanced SARS-CoV titers and clinical disease, suggesting an antiviral role during infection. Importantly, about 60 % of the murine QTL encode susceptibility genes identified as priority candidates from human genome-wide association studies (GWAS) studies after SARS-CoV-2 infection, suggesting that similar selective forces have targeted analogous genes and pathways to regulate Sarbecovirus disease across diverse mammalian hosts. These studies provide an experimental platform in rodents to investigate the molecular-genetic mechanisms by which potential cross mammalian susceptibility loci and genes regulate type-specific and cross-SARS-like group 2B CoV replication, immunity, and pathogenesis in rodent models. Our study also provides a paradigm for identifying susceptibility loci for other highly heterogeneous and virulent viruses that sporadically emerge from zoonotic reservoirs to plague human and animal populations.

Mapping of susceptibility loci for Ebola virus pathogenesis in mice.

Ebola virus (EBOV), a major global health concern, causes severe, often fatal EBOV disease (EVD) in humans. Host genetic variation plays a critical role, yet the identity of host susceptibility loci in mammals remains unknown. Using genetic reference populations, we generate an F2 mapping cohort to identify host susceptibility loci that regulate EVD. While disease-resistant mice display minimal pathogenesis, susceptible mice display severe liver pathology consistent with EVD-like disease and transcriptional signatures associated with inflammatory and liver metabolic processes. A significant quantitative trait locus (QTL) for virus RNA load in blood is identified in chromosome (chr)8, and a severe clinical disease and mortality QTL is mapped to chr7, which includes the Trim5 locus. Using knockout mice, we validate the Trim5 locus as one potential driver of liver failure and mortality after infection. The identification of susceptibility loci provides insight into molecular genetic mechanisms regulating EVD progression and severity, potentially informing therapeutics and vaccination strategies.

The Updated Mouse Universal Genotyping Array Bioinformatic Pipeline Improves Genetic QC in Laboratory Mice.

The MiniMUGA genotyping array is a popular tool for genetic QC of laboratory mice and genotyping of samples from most types of experimental crosses involving laboratory strains, particularly for reduced complexity crosses. The content of the production version of the MiniMUGA array is fixed; however, there is the opportunity to improve array's performance and the associated report's usefulness by leveraging thousands of samples genotyped since the initial description of MiniMUGA in 2020. Here we report our efforts to update and improve marker annotation, increase the number and the reliability of the consensus genotypes for inbred strains and increase the number of constructs that can reliably be detected with MiniMUGA. In addition, we have implemented key changes in the informatics pipeline to identify and quantify the contribution of specific genetic backgrounds to the makeup of a given sample, remove arbitrary thresholds, include the Y Chromosome and mitochondrial genome in the ideogram, and improve robust detection of the presence of commercially available substrains based on diagnostic alleles. Finally, we have made changes to the layout of the report, to simplify the interpretation and completeness of the analysis and added a table summarizing the ideogram. We believe that these changes will be of general interest to the mouse research community and will be instrumental in our goal of improving the rigor and reproducibility of mouse-based biomedical research.

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants.

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

Large differences in collateral blood vessel abundance among individuals arise from multiple genetic variants.

Collateral blood flow varies greatly among humans for reasons that remain unclear, resulting in significant differences in ischemic tissue damage. A similarly large variation has also been found in mice that is caused by genetic background-dependent differences in the extent of collateral formation, termed collaterogenesis-a unique angiogenic process that occurs during development and determines collateral number and diameter in the adult. Previous studies have identified several quantitative trait loci (QTL) linked to this variation. However, understanding has been hampered by the use of closely related inbred strains that do not model the wide genetic variation present in the "outbred" human population. The Collaborative Cross (CC) multiparent mouse genetic reference panel was developed to address this limitation. Herein we measured the number and average diameter of cerebral collaterals in 60 CC strains, their 8 founder strains, 8 F1 crosses of CC strains selected for abundant versus sparse collaterals, and 2 intercross populations created from the latter. Collateral number evidenced 47-fold variation among the 60 CC strains, with 14% having poor, 25% poor-to-intermediate, 47% intermediate-to-good, and 13% good collateral abundance, that was associated with large differences in post-stroke infarct volume. Collateral number in skeletal muscle and intestine of selected high- and low-collateral strains evidenced the same relative abundance as in brain. Genome-wide mapping demonstrated that collateral abundance is a highly polymorphic trait. Subsequent analysis identified: 6 novel QTL circumscribing 28 high-priority candidate genes harboring putative loss-of-function polymorphisms (SNPs) associated with low collateral number; 335 predicted-deleterious SNPs present in their human orthologs; and 32 genes associated with vascular development but lacking protein coding variants. Six additional suggestive QTL (LOD > 4.5) were also identified in CC-wide QTL mapping. This study provides a comprehensive set of candidate genes for future investigations aimed at identifying signaling proteins within the collaterogenesis pathway whose variants potentially underlie genetic-dependent collateral insufficiency in brain and other tissues.

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.