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

Mapping the genetic landscape establishing a tumor immune microenvironment favorable for anti-PD-1 response in mice and humans

Identifying host genetic factors modulating immune checkpoint inhibitor (ICI) efficacy has been experimentally challenging because of variations in both host and tumor genomes, differences in the microbiome, and patient life exposures. Utilizing the Collaborative Cross (CC) multi-parent mouse genetic resource population, we developed an approach that fixes the tumor genomic configuration while varying host genetics. With this approach, we discovered that response to anti-PD-1 (aPD1) immunotherapy was significantly heritable in four distinct murine tumor models (H2 between 0.18-0.40). For the MC38 colorectal carcinoma system (H2 = 0.40), we mapped four significant ICI response quantitative trait loci (QTL) localized to mouse chromosomes (mChr) 5, 9, 15 and 17, and identified significant epistatic interactions between specific QTL pairs. Differentially expressed genes within these QTL were highly enriched for immune genes and pathways mediating allograft rejection and graft vs host disease. Using a cross species analytical approach, we found a core network of 48 genes within the four QTLs that showed significant prognostic value for overall survival in aPD1 treated human cohorts that outperformed all other existing validated immunotherapy biomarkers, especially in human tumors of the previously defined immune subtype 4. Functional blockade of two top candidate immune targets within the 48 gene network, GM-CSF and high affinity IL-2/IL-15 signaling, completely abrogated the MC38 tumor transcriptional response to aPD1 therapy in vivo. Thus, we have established a powerful cross species in vivo platform capable of uncovering host genetic factors that establish the tumor immune microenvironment configuration propitious for ICI response.

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.

Genetic regulation of carnitine metabolism controls lipid damage repair and aging RBC hemolysis in vivo and in vitro

ABSTRACT

Recent large-scale multiomics studies suggest that genetic factors influence the chemical individuality of donated blood. To examine this concept, we performed metabolomics analyses of 643 blood units from volunteers who donated units of packed red blood cells (RBCs) on 2 separate occasions. These analyses identified carnitine metabolism as the most reproducible pathway across multiple donations from the same donor. We also measured l-carnitine and acyl-carnitines in 13 091 packed RBC units from donors in the Recipient Epidemiology and Donor Evaluation study. Genome-wide association studies against 879 000 polymorphisms identified critical genetic factors contributing to interdonor heterogeneity in end-of-storage carnitine levels, including common nonsynonymous polymorphisms in genes encoding carnitine transporters (SLC22A16, SLC22A5, and SLC16A9); carnitine synthesis (FLVCR1 and MTDH) and metabolism (CPT1A, CPT2, CRAT, and ACSS2), and carnitine-dependent repair of lipids oxidized by ALOX5. Significant associations between genetic polymorphisms on SLC22 transporters and carnitine pools in stored RBCs were validated in 525 Diversity Outbred mice. Donors carrying 2 alleles of the rs12210538 SLC22A16 single-nucleotide polymorphism exhibited the lowest l-carnitine levels, significant elevations of in vitro hemolysis, and the highest degree of vesiculation, accompanied by increases in lipid peroxidation markers. Separation of RBCs by age, via in vivo biotinylation in mice, and Percoll density gradients of human RBCs, showed age-dependent depletions of l-carnitine and acyl-carnitine pools, accompanied by progressive failure of the reacylation process after chemically induced membrane lipid damage. Supplementation of stored murine RBCs with l-carnitine boosted posttransfusion recovery, suggesting this could represent a viable strategy to improve RBC storage quality.

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 the blood bank. Here we leveraged 13,029 volunteers from the Recipient Epidemiology and Donor Evaluation Study to identify an association 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, 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 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.

eTOC and Highlights

Highlights

Blood donor age and sex affect glycolysis in stored RBCs from 13,029 volunteers;Ancestry, genetic polymorphisms in PFKP, HK1, CD38/BST1 influence RBC glycolysis;Modeled PFKP effects relate to preventing loss of the total AXP pool in stored RBCs;ATP and hypoxanthine are biomarkers of hemolysis in vitro and in vivo.

Ferroptosis regulates hemolysis in stored murine and human red blood cells

Red blood cell (RBC) metabolism regulates hemolysis during aging in vivo and in the blood bank. Here, we leveraged a diversity outbred mouse population to map the genetic drivers of fresh/stored RBC metabolism and extravascular hemolysis upon storage and transfusion in 350 mice. We identify the ferrireductase Steap3 as a critical regulator of a ferroptosis-like process of lipid peroxidation. Steap3 polymorphisms were associated with RBC iron content, in vitro hemolysis, and in vivo extravascular hemolysis both in mice and 13,091 blood donors from the Recipient Epidemiology and Donor evaluation Study. Using metabolite Quantitative Trait Loci analyses, we identified a network of gene products (FADS1/2, EPHX2 and LPCAT3) - enriched in donors of African descent - associated with oxylipin metabolism in stored human RBCs and related to Steap3 or its transcriptional regulator, the tumor protein TP53. Genetic variants were associated with lower in vivo hemolysis in thousands of single-unit transfusion recipients.

Highlights

Steap3 regulates lipid peroxidation and extravascular hemolysis in 350 diversity outbred miceSteap3 SNPs are linked to RBC iron, hemolysis, vesiculation in 13,091 blood donorsmQTL analyses of oxylipins identified ferroptosis-related gene products FADS1/2, EPHX2, LPCAT3Ferroptosis markers are linked to hemoglobin increments in transfusion recipients.

Graphical abstract

The role of gonadal hormones and sex chromosomes in sex-dependent effects of early nutrition on metabolic health

Early-life conditions such as prenatal nutrition can have long-term effects on metabolic health, and these effects may differ between males and females. Understanding the biological mechanisms underlying sex differences in the response to early-life environment will improve interventions, but few such mechanisms have been identified, and there is no overall framework for understanding sex differences. Biological sex differences may be due to chromosomal sex, gonadal sex, or interactions between the two. This review describes approaches to distinguish between the roles of chromosomal and gonadal sex, and summarizes findings regarding sex differences in metabolism. The Four Core Genotypes (FCG) mouse model allows dissociation of the sex chromosome genotype from gonadal type, whereas the XY* mouse model can be used to distinguish effects of X chromosome dosage vs the presence of the Y chromosome. Gonadectomy can be used to distinguish between organizational (permanent) and activational (reversible) effects of sex hormones. Baseline sex differences in a variety of metabolic traits are influenced by both activational and organizational effects of gonadal hormones, as well as sex chromosome complement. Thus far, these approaches have not been widely applied to examine sex-dependent effects of prenatal conditions, although a number of studies have found activational effects of estradiol to be protective against the development of hypertension following early-life adversity. Genes that escape X chromosome inactivation (XCI), such as Kdm5c, contribute to baseline sex-differences in metabolism, while Ogt, another XCI escapee, leads to sex-dependent responses to prenatal maternal stress. Genome-wide approaches to the study of sex differences include mapping genetic loci influencing metabolic traits in a sex-dependent manner. Seeking enrichment for binding sites of hormone receptors among genes showing sexually-dimorphic expression can elucidate the relative roles of hormones. Using the approaches described herein to identify mechanisms underlying sex-dependent effects of early nutrition on metabolic health may enable the identification of fundamental mechanisms and potential interventions.

Cell morphology QTL reveal gene by environment interactions in a genetically diverse cell population

Genetically heterogenous cell lines from laboratory mice are promising tools for population-based screening as they offer power for genetic mapping, and potentially, predictive value for in vivo experimentation in genetically matched individuals. To explore this further, we derived a panel of fibroblast lines from a genetic reference population of laboratory mice (the Diversity Outbred, DO). We then used high-content imaging to capture hundreds of cell morphology traits in cells exposed to the oxidative stress-inducing arsenic metabolite monomethylarsonous acid (MMAIII). We employed dose-response modeling to capture latent parameters of response and we then used these parameters to identify several hundred cell morphology quantitative trait loci (cmQTL). Response cmQTL encompass genes with established associations with cellular responses to arsenic exposure, including Abcc4 and Txnrd1, as well as novel gene candidates like Xrcc2. Moreover, baseline trait cmQTL highlight the influence of natural variation on fundamental aspects of nuclear morphology. We show that the natural variants influencing response include both coding and non-coding variation, and that cmQTL haplotypes can be used to predict response in orthogonal cell lines. Our study sheds light on the major molecular initiating events of oxidative stress that are under genetic regulation, including the NRF2-mediated antioxidant response, cellular detoxification pathways, DNA damage repair response, and cell death trajectories.

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.

Heterogeneity in the response to n-3 polyunsaturated fatty acids

PURPOSE OF REVIEW

A central goal in the study of long chain n-3 polyunsaturated fatty acids (PUFA) is to translate findings from the basic sciences to the population level to improve human health and prevent chronic diseases. A tenet of this vision is to think in terms of precision medicine and nutrition, that is, stratification of individuals into differing groups that will have different needs across the lifespan for n-3 PUFAs. Therefore, there is a critical need to identify the sources of heterogeneity in the human population in the dietary response to n-3 PUFA intervention.

RECENT FINDINGS

We briefly review key sources of heterogeneity in the response to intake of long chain n-3 PUFAs. These include background diet, host genome, composition of the gut microbiome, and sex. We also discuss the need to integrate data from newer rodent models (e.g. population-based approaches), multi -omics, and analyses of big data using machine learning and data-driven cluster analyses.

SUMMARY

Accounting for vast heterogeneity in the human population, particularly with the use of big data integrated with preclinical evidence, will drive the next generation of precision nutrition studies and randomized clinical trials with long-chain n-3 PUFAs.

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.