R/qtl2: Software for Mapping Quantitative Trait Loci with High-Dimensional Data and Multiparent Populations

Genetic and transcriptomic analysis of lentil seed imbibition and dormancy in relation to its domestication

Seed dormancy is an adaptation that delays germination to prevent the start of this process during unsuitable conditions. It is crucial in wild species but its loss was selected during crop domestication to ensure a fast and uniform germination. Water uptake, or imbibition, is the first step of germination. In the Fabaceae family, seeds have physical dormancy, in which seed coats are impermeable to water. We used an interspecific cross between an elite lentil line (Lens culinaris) and a wild lentil (L. orientalis) to investigate the genetic basis of imbibition capacity through quantitative trait locus (QTL) mapping and by using RNA from embryos and seed coats at different development stages, and phenotypic data of seed coat thickness (SCT) and proportion of imbibed seeds (PIS). Both characteristics were consistent throughout different years and locations, suggesting a hereditary component. QTL results suggest that they are each controlled by relatively few loci. Differentially expressed genes (DEGs) within the QTL were considered candidate genes. Two glycosyl-hydrolase genes (a β-glucosidase and a β-galactosidase), which degrade complex polysaccharides in the cell wall, were found among the candidate genes, and one of them had a positive correlation (β-glucosidase) between gene expression and imbibition capacity, and the other gene (β-galactosidase) presented a negative correlation between gene expression and SCT.

A deficient CP24 allele defines variation for dynamic nonphotochemical quenching and photosystem II efficiency in maize

Maize (Zea mays L.) is a global crop species in which CO2 assimilation occurs via the C4 pathway. C4 photosynthesis is typically more efficient than C3 photosynthesis under warm and dry conditions; however, despite this inherent advantage, considerable variation remains in photosynthetic efficiency for C4 species that could be leveraged to benefit crop performance. Here, we investigate the genetic architecture of nonphotochemical quenching (NPQ) and photosystem II (PSII) efficiency using a combination of high-throughput phenotyping and quantitative trait loci (QTL) mapping in a field-grown Multi-parent Advanced Generation Inter-Cross (MAGIC) mapping population. QTL mapping was followed by the identification of putative candidate genes using a combination of genomics, transcriptomics, protein biochemistry, and targeted physiological phenotyping. We identified four genes with a putative causal role in the observed QTL effects. The highest confidence causal gene was found for a large effect QTL for photosynthetic efficiency on chromosome 10, which was underpinned by allelic variation in the expression of the minor PSII antenna protein light harvesting complex photosystem II subunit (LHCB6 or CP24), mainly driven by poor expression associated with the haplotype of the F7 founder line. The historical role of this line in breeding for early flowering time may suggest that the presence of this deficient allele could be enriched in temperate maize germplasm. These findings advance our understanding of the genetic basis of NPQ and PSII efficiency in C4 plants and highlight the potential for breeding strategies aimed at optimizing photosynthetic efficiency in maize.

Gut metagenomes reveal interactions between dietary restriction, ageing and the microbiome in genetically diverse mice

The gut microbiome changes with age and has been proposed to mediate the benefit of lifespan-extending interventions such as dietary restriction. However, the causes and consequences of microbiome ageing and the potential of such interventions remain unclear. Here we analysed 2,997 metagenomes collected longitudinally from 913 deeply phenotyped, genetically diverse mice to investigate interactions between the microbiome, ageing, dietary restriction (caloric restriction and fasting), host genetics and a range of health parameters. Among the numerous age-associated microbiome changes that we find in this cohort, increased microbiome uniqueness is the most consistent parameter across a second longitudinal mouse experiment that we performed on inbred mice and a compendium of 4,101 human metagenomes. Furthermore, cohousing experiments show that age-associated microbiome changes may be caused by an accumulation of stochastic environmental exposures (neutral theory) rather than by the influence of an ageing host (selection theory). Unexpectedly, the majority of taxonomic and functional microbiome features show small but significant heritability, and the amount of variation explained by host genetics is similar to ageing and dietary restriction. We also find that more intense dietary interventions lead to larger microbiome changes and that dietary restriction does not rejuvenate the microbiome. Lastly, we find that the microbiome is associated with multiple health parameters, including body composition, immune components and frailty, but not lifespan. Overall, this study sheds light on the factors influencing microbiome ageing and aspects of host physiology modulated by the microbiome.

Genetic architecture of the red blood cell proteome in genetically diverse mice reveals central role of hemoglobin beta cysteine redox status in maintaining circulating glutathione pools

Red blood cells (RBCs) transport oxygen but accumulate oxidative damage over time, reducing function in vivo and during storage-critical for transfusions. To explore genetic influences on RBC resilience, we profiled proteins, metabolites, and lipids from fresh and stored RBCs obtained from 350 genetically diverse mice. Our analysis identified over 6,000 quantitative trait loci (QTL). Compared to other tissues, prevalence of trans genetic effects over cis reflects the absence of de novo protein synthesis in anucleated RBCs. QTL hotspots at Hbb, Hba, Mon1a, and storage-specific Steap3 linked ferroptosis to hemolysis. Proteasome components clustered at multiple loci, underscoring the importance of degrading oxidized proteins. Post-translational modifications (PTMs) mapped predominantly to hemoglobins, particularly cysteine residues. Loss of reactive C93 in humanized mice (HBB C93A) disrupted redox balance, affecting glutathione pools, protein glutathionylation, and redox PTMs. These findings highlight genetic regulation of RBC oxidation, with implications for transfusion biology and oxidative stress-dependent hemolytic disorders.

Uncovering the genetic basis of competitiveness and the potential for cooperation in plant groups

Crop productivity was transformed by incorporating dwarfing genes that made plants smaller and less competitive (more cooperative). Beyond such major shifts in plant size, however, it is not clear how much variation in competitiveness remains and how to find its genetic basis. We performed plant density experiments, using 484 lines of the Arabidopsis thaliana multi-parent advanced generation inter-cross population, to compare methods for mapping the genetic basis of plant competitiveness. We first found that a major dwarfing gene, the erecta allele, caused reduced competitiveness and higher group productivity. Then, measuring competitiveness more generally, we found: (i) extensive variation in generic measures of competitiveness that extended beyond the effects of the erecta allele; (ii) a novel genomic region underlying variation in competitiveness; and (iii) that some measures of competitiveness were more useful than others. Our results show how modern genomic resources, including multi-parent populations, could uncover hidden genes for more cooperative crop plants.

Ferroptosis regulates hemolysis in stored murine and human red blood cells

ABSTRACT

Red blood cell (RBC) metabolism regulates hemolysis during aging in vivo and in the blood bank. However, the genetic underpinnings of RBC metabolic heterogeneity and extravascular hemolysis at population scale are incompletely understood. On the basis of the breeding of 8 founder strains with extreme genetic diversity, the Jackson Laboratory diversity outbred population can capture the impact of genetic heterogeneity in like manner to population-based studies. RBCs from 350 outbred mice, either fresh or stored for 7 days, were tested for posttransfusion recovery, as well as metabolomics and lipidomics analyses. Metabolite and lipid quantitative trait loci (QTL) mapped >400 gene-metabolite associations, which we collated into an online interactive portal. Relevant to RBC storage, we identified a QTL hotspot on chromosome 1, mapping on the region coding for the ferrireductase 6-transmembrane epithelial antigen of the prostate 3 (Steap3), a transcriptional target to p53. Steap3 regulated posttransfusion recovery, contributing to a ferroptosis-like process of lipid peroxidation, as validated via genetic manipulation in mice. Translational validation of murine findings in humans, STEAP3 polymorphisms were associated with RBC iron content, lipid peroxidation, and in vitro hemolysis in 13 091 blood donors from the Recipient Epidemiology and Donor Evaluation Study. QTL analyses in humans identified a network of gene products (fatty acid desaturases 1 and 2, epoxide hydrolase 2, lysophosphatidylcholine acetyl-transferase 3, solute carrier family 22 member 16, glucose 6-phosphate dehydrogenase, very long chain fatty acid elongase, and phospholipase A2 group VI) associated with altered levels of oxylipins. These polymorphisms were prevalent in donors of African descent and were linked to allele frequency of hemolysis-linked polymorphisms for Steap3 or p53. These genetic variants were also associated with lower hemoglobin increments in thousands of single-unit transfusion recipients from the vein-to-vein database.

Arginine metabolism is a biomarker of red blood cell and human aging

Increasing global life expectancy motivates investigations of molecular mechanisms of aging and age-related diseases. This study examines age-associated changes in red blood cells (RBCs), the most numerous host cell in humans. Four cohorts, including healthy individuals and patients with sickle cell disease, were analyzed to define age-dependent changes in RBC metabolism. Over 15,700 specimens from 13,757 humans were examined, a major expansion over previous studies of RBCs in aging. Multi-omics approaches identified chronological age-related alterations in the arginine pathway with increased arginine utilization in RBCs from older individuals. These changes were consistent across healthy and sickle cell disease cohorts and were influenced by genetic variation, sex, and body mass index. Integrating multi-omics data and metabolite quantitative trait loci (mQTL) in humans and 525 diversity outbred mice functionally linked metabolism of arginine during RBC storage to increased vesiculation-a hallmark of RBC aging-and lower post-transfusion hemoglobin increments. Thus, arginine metabolism is a biomarker of RBC and organismal aging, suggesting potential new targets for addressing sequelae of aging.

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1 and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis and other bioinformatics inquiries. The combined behavioural and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6. Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.

Cross cell-type systems genetics reveals the influence of eQTL at multiple points in the developmental trajectory of mouse neural progenitor cells

Genetic variation leads to phenotypic variability in pluripotent stem cells that presents challenges for regenerative medicine. Although recent studies have investigated the impact of genetic variation on pluripotency maintenance and differentiation capacity, less is known about how genetic variants affecting the pluripotent state influence gene regulation in later stages of development. Here, we characterized expression of more than 12,000 genes in 127 donor-matched Diversity Outbred (DO) mouse embryonic stem cell (mESC) and neural progenitor cell (mNPC) lines. Quantitative trait locus (QTL) mapping identified 2,947 expression QTL (eQTL) unique to DO mNPCs and 1,113 eQTL observed in both mNPCs and mESCs with highly concordant allele effects. We mapped three eQTL hotspots on Chromosomes (Chrs) 1, 10, and 11 that were unique to mNPCs. Target genes of the Chr 1 hotspot were overrepresented for those involved in mRNA processing, DNA repair, chromatin organization, protein degradation, and cell cycle. Mediation analysis of the Chr 1 hotspot identified Rnf152 as the best candidate mediator expressed in mNPCs, while cross-cell type mediation using mESC gene expression along with partial correlation analysis strongly implicated genetic variant(s) affecting Pign expression in the mESC state as regulating the mNPC Chr 1 eQTL hotspot. Together these findings highlight that many local eQTL confer similar effects on gene expression in multiple cell states; distant eQTL in DO mNPCs are numerous and largely unique to that cell state, with many co-localizing to mNPC-specific hotspots; and mediation analysis across cell types suggests that expression of Pign early in development (mESCs) shapes the transcriptome of the more specialized mNPC state.

Diurnal Regulation of SOS Pathway and Sodium Excretion Underlying Salinity Tolerance of Vigna marina

Vigna marina (Barm.) Merr. is adapted to tropical marine beaches and has an outstanding tolerance to salt stress. Given there are growing demands for cultivating crops in saline soil or with saline water, it is important to understand how halophytic species are adapted to the saline environments. Here we revealed by positron-emitting tracer imaging system (PETIS) that V. marina actively excretes sodium from the root during the light period but not the dark period. The following whole genome sequencing accompanied with forward genetic study identified a QTL region harbouring SOS1, encoding plasma membrane Na+/H+ antiporter, which was associated with not only salt tolerance but also the ability of sodium excretion. We also found the QTL region contained a large structural rearrangement that suppressed recombination across ~14 Mbp, fixing multiple gene loci potentially involved in salt tolerance. RNA-seq and promoter analyses revealed SOS1 in V. marina was highly expressed even without salt stress and its promoter shared common cis-regulatory motifs with those exhibiting similar expression profiles. Interestingly, the cis-regulatory motifs seemed installed by a transposable element (TE) insertion. Though not identified by genetic analysis, the transcriptome data also revealed SOS2 transcription was under diurnal regulation, explaining the pattern of sodium excretion together with upregulated expression of SOS1. Altogether, the study elucidated one aspect of the strategy adopted by V. marina to adapt to marine beach, which is highly saline and transpiring.

Unlocking the power of AI for phenotyping fruit morphology in Arabidopsis

Deep learning can revolutionise high-throughput image-based phenotyping by automating the measurement of complex traits, a task that is often labour-intensive, time-consuming, and prone to human error. However, its precision and adaptability in accurately phenotyping organ-level traits, such as fruit morphology, remain to be fully evaluated. Establishing the links between phenotypic and genotypic variation is essential for uncovering the genetic basis of traits and can also provide an orthologous test of pipeline effectiveness. In this study, we assess the efficacy of deep learning for measuring variation in fruit morphology in Arabidopsis using images from a multiparent advanced generation intercross (MAGIC) mapping family. We trained an instance segmentation model and developed a pipeline to phenotype Arabidopsis fruit morphology, based on the model outputs. Our model achieved strong performance with an average precision of 88.0% for detection and 55.9% for segmentation. Quantitative trait locus analysis of the derived phenotypic metrics of the MAGIC population identified significant loci associated with fruit morphology. This analysis, based on automated phenotyping of 332,194 individual fruits, underscores the capability of deep learning as a robust tool for phenotyping large populations. Our pipeline for quantifying pod morphological traits is scalable and provides high-quality phenotype data, facilitating genetic analysis and gene discovery, as well as advancing crop breeding research.

PITX2 expression and Neanderthal introgression in HS3ST3A1 contribute to variation in tooth dimensions in modern humans

Dental morphology varies greatly throughout evolution, including in the human lineage, but little is known about the biology of this variation. Here, we use multiomics analyses to examine the genetics of variation in tooth crown dimensions. In a human cohort with mixed continental ancestry, we detected genome-wide significant associations at 18 genome regions. One region includes EDAR, a gene known to impact dental features in East Asians. Furthermore, we find that EDAR variants increase the mesiodistal diameter of all teeth, following an anterior-posterior gradient of decreasing strength. Among the 17 novel-associated regions, we replicate 7/13 in an independent human cohort and find that 4/12 orthologous regions affect molar size in mice. Two association signals point to compelling candidate genes. One is ∼61 kb from PITX2, a major determinant of tooth development. Another overlaps HS3ST3A1, a paralogous neighbor of HS3ST3B1, a tooth enamel knot factor. We document the expression of Pitx2 and Hs3st3a1 in enamel knot and dental epithelial cells of developing mouse incisors. Furthermore, associated SNPs in PITX2 and HS3ST3A1 overlap enhancers active in these cells, suggesting a role for these SNPs in gene regulation during dental development. In addition, we document that Pitx2 and Hs3st3a1/Hs3st3b1 knockout mice show alterations in dental morphology. Finally, we find that associated SNPs in HS3ST3A1 are in a DNA tract introgressed from Neanderthals, consistent with an involvement of HS3ST3A1 in tooth size variation during human evolution.

Addressing Sample Mix-Ups: Tools and Approaches for Large-Scale Multi-Omics Studies

Advances in high-throughput omics technologies have enabled system-wide characterization of biological samples across multiple molecular levels, such as the genome, transcriptome, and proteome. However, as sample sizes rapidly increase in large-scale multi-omics studies, sample mix-ups have become a prevalent issue, compromising data integrity and leading to erroneous conclusions. The interconnected nature of multi-omics data presents an opportunity to identify and correct these errors. This review examines the potential sources of sample mix-ups and evaluates the methodologies and tools developed for detecting and correcting these errors, with an emphasis on approaches applicable to proteomics data. We categorize existing tools into three main groups: expression/protein quantitative trait loci-based, genotype concordance-based, and gene/protein expression correlation-based approaches. Notably, only a handful of tools currently utilize the proteogenomics approach for correcting sample mix-ups at the proteomics level. Integrating the strengths of current tools across diverse data types could enable the development of more versatile and comprehensive solutions. In conclusion, verifying sample identity is a critical first step to reduce bias and increase precision in subsequent analyses for large-scale multi-omics studies. By leveraging these tools for identifying and correcting sample mix-ups, researchers can significantly improve the reliability and reproducibility of biomedical research.

Chromosome-scale Reference Genome and RAD-based Genetic Map of Yellow Starthistle (Centaurea solstitialis) Reveal Putative Structural Variation and QTL Associated With Invader Traits

Invasive species offer outstanding opportunities to identify the genomic sources of variation that contribute to rapid adaptation, as well as the genetic mechanisms facilitating invasions. The Eurasian plant yellow starthistle (Centaurea solstitialis) is highly invasive in North and South American grasslands and known to have evolved increased growth and reproduction during invasion. Here, we develop new genomic resources for C. solstitialis and map the genetic basis of invasiveness traits. We present a chromosome-scale (1N = 8) reference genome using PacBio CLR and Dovetail Omni-C technologies, and functional gene annotation using RNAseq. We find repeat structure typical of the family Asteraceae, with over 25% of gene content derived from ancestral whole-genome duplications (paleologs). Using an F2 mapping population derived from a cross between native and invading parents, with a restriction site-associated DNA (RAD)-based genetic map, we validate the assembly and identify 13 quantitative trait loci underpinning size traits that have evolved during invasion. We find evidence that large effects of quantitative trait loci may be associated with structural variants between native and invading genotypes, including a variant with an overdominant and pleiotropic effect on key invader traits. We also find evidence of significant paleolog enrichment under two quantitative trait loci. Our results add to growing evidence of the importance of structural variants in evolution, and to understanding of the rapid evolution of invaders.

Analysis of lifespan across Diversity Outbred mouse studies identifies multiple longevity-associated loci

Lifespan is an integrative phenotype whose genetic architecture is likely to highlight multiple processes with high impact on health and aging. Here, we conduct a genetic meta-analysis of longevity in Diversity Outbred (DO) mice that includes 2,444 animals from three independently conducted lifespan studies. We identify six loci that contribute significantly to lifespan independently of diet and drug treatment, one of which also influences lifespan in a sex-dependent manner, as well as an additional locus with a diet-specific effect on lifespan. Collectively, these loci explain over half of the estimated heritable variation in lifespan across these studies and provide insight into the genetic architecture of lifespan in DO mice.

Analysis of CNS autoimmunity in genetically diverse mice reveals unique phenotypes and mechanisms

Multiple sclerosis (MS) is a complex disease with significant heterogeneity in disease course and progression. Genetic studies have identified numerous loci associated with MS risk, but the genetic basis of disease progression remains elusive. To address this, we leveraged the Collaborative Cross (CC), a genetically diverse mouse strain panel, and experimental autoimmune encephalomyelitis (EAE). The 32 CC strains studied captured a wide spectrum of EAE severity, trajectory, and presentation, including severe-progressive, monophasic, relapsing remitting, and axial rotary-EAE (AR-EAE), accompanied by distinct immunopathology. Sex differences in EAE severity were observed in 6 strains. Quantitative trait locus analysis revealed distinct genetic linkage patterns for different EAE phenotypes, including EAE severity and incidence of AR-EAE. Machine learning-based approaches prioritized candidate genes for loci underlying EAE severity (Abcc4 and Gpc6) and AR-EAE (Yap1 and Dync2h1). This work expands the EAE phenotypic repertoire and identifies potentially novel loci controlling unique EAE phenotypes, supporting the hypothesis that heterogeneity in MS disease course is driven by genetic variation.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Identification of novel genetic loci and candidate genes for progressive ethanol consumption in diversity outbred mice

Mouse behavioral genetic mapping studies can identify genomic intervals modulating complex traits under well-controlled environmental conditions and have been used to study ethanol behaviors to aid in understanding genetic risk and the neurobiology of alcohol use disorder (AUD). However, historically such studies have produced large confidence intervals, thus complicating identification of potential causal candidate genes. Diversity Outbred (DO) mice offer the ability to perform high-resolution quantitative trait loci (QTL) mapping on a very genetically diverse background, thus facilitating identification of candidate genes. Here, we studied a population of 636 male DO mice with four weeks of intermittent ethanol access via a three-bottle choice procedure, producing a progressive ethanol consumption phenotype. QTL analysis identified 3 significant (Chrs 3, 4, and 12) and 13 suggestive loci for ethanol-drinking behaviors with narrow confidence intervals (1-4 Mbp for significant QTLs). Results suggested that genetic influences on initial versus progressive ethanol consumption were localized to different genomic intervals. A defined set of positional candidate genes were prioritized using haplotype analysis, identified coding polymorphisms, prefrontal cortex transcriptomics data, human GWAS data and prior rodent gene set data for ethanol or other misused substances. These candidates included Car8, the lone gene with a significant cis-eQTL within a Chr 4 QTL for week four ethanol consumption. These results represent the highest-resolution genetic mapping of ethanol consumption behaviors in mice to date, providing identification of novel loci and candidate genes for study in relation to the neurobiology of AUD.

High-content phenotypic analysis of a C. elegans recombinant inbred population identifies genetic and molecular regulators of lifespan

Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 C. elegans recombinant inbred advanced intercross lines (RIAILs). We assessed molecular profiles-transcriptome, proteome, and lipidome-and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which correlated positively with developmental time. We validated three longevity modulators, including rict-1, gfm-1, and mltn-1, among the top candidates obtained from multiomics data integration and quantitative trait locus (QTL) mapping. We translated their relevance to humans using UK Biobank data and showed that variants in GFM1 are associated with an elevated risk of age-related heart failure. We organized our dataset as a resource that allows interactive explorations for new longevity targets.

Dietary restriction impacts health and lifespan of genetically diverse mice

Caloric restriction extends healthy lifespan in multiple species1. Intermittent fasting, an alternative form of dietary restriction, is potentially more sustainable in humans, but its effectiveness remains largely unexplored2-8. Identifying the most efficacious forms of dietary restriction is key for developing interventions to improve human health and longevity9. Here we performed an extensive assessment of graded levels of caloric restriction (20% and 40%) and intermittent fasting (1 and 2 days fasting per week) on the health and survival of 960 genetically diverse female mice. We show that caloric restriction and intermittent fasting both resulted in lifespan extension in proportion to the degree of restriction. Lifespan was heritable and genetics had a larger influence on lifespan than dietary restriction. The strongest trait associations with lifespan included retention of body weight through periods of handling-an indicator of stress resilience, high lymphocyte proportion, low red blood cell distribution width and high adiposity in late life. Health effects differed between interventions and exhibited inconsistent relationships with lifespan extension. 40% caloric restriction had the strongest lifespan extension effect but led to a loss of lean mass and changes in the immune repertoire that could confer susceptibility to infections. Intermittent fasting did not extend the lifespan of mice with high pre-intervention body weight, and two-day intermittent fasting was associated with disruption of erythroid cell populations. Metabolic responses to dietary restriction, including reduced adiposity and lower fasting glucose, were not associated with increased lifespan, suggesting that dietary restriction does more than just counteract the negative effects of obesity. Our findings indicate that improving health and extending lifespan are not synonymous and raise questions about which end points are the most relevant for evaluating aging interventions in preclinical models and clinical trials.

Transcripts with high distal heritability mediate genetic effects on complex metabolic traits

Although many genes are subject to local regulation, recent evidence suggests that complex distal regulation may be more important in mediating phenotypic variability. To assess the role of distal gene regulation in complex traits, we combined multi-tissue transcriptomes with physiological outcomes to model diet-induced obesity and metabolic disease in a population of Diversity Outbred mice. Using a novel high-dimensional mediation analysis, we identified a composite transcriptome signature that summarized genetic effects on gene expression and explained 30% of the variation across all metabolic traits. The signature was heritable, interpretable in biological terms, and predicted obesity status from gene expression in an independently derived mouse cohort and multiple human studies. Transcripts contributing most strongly to this composite mediator frequently had complex, distal regulation distributed throughout the genome. These results suggest that trait-relevant variation in transcription is largely distally regulated, but is nonetheless identifiable, interpretable, and translatable across species.

Integration of QTL and transcriptome approaches for the identification of genes involved in tomato response to nitrogen deficiency

Optimizing plant nitrogen (N) usage and inhibiting N leaching loss in the soil-crop system is crucial to maintaining crop yield and reducing environmental pollution. This study aimed at identifying quantitative trait loci (QTLs) and differentially expressed genes (DEGs) between two N treatments in order to list candidate genes related to nitrogen-related contrasting traits in tomato varieties. We characterized a genetic diversity core-collection (CC) and a multi-parental advanced generation intercross (MAGIC) tomato population grown in a greenhouse under two nitrogen levels and assessed several N-related traits and mapped QTLs. Transcriptome response under the two N conditions was also investigated through RNA sequencing of fruit and leaves in four parents of the MAGIC population. Significant differences in response to N input reduction were observed at the phenotypic level for biomass and N-related traits. Twenty-seven QTLs were detected for three target traits (leaf N content, leaf nitrogen balance index, and petiole NO3- content), 10 and six in the low and high N condition, respectively, while 19 QTLs were identified for plasticity traits. At the transcriptome level, 4752 and 2405 DEGs were detected between the two N conditions in leaves and fruits, respectively, among which 3628 (50.6%) in leaves and 1717 (71.4%) in fruit were genotype specific. When considering all the genotypes, 1677 DEGs were shared between organs or tissues. Finally, we integrated DEG and QTL analyses to identify the most promising candidate genes. The results highlighted a complex genetic architecture of N homeostasis in tomato and novel putative genes useful for breeding tomato varieties requiring less N input.

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.

Enhlink infers distal and context-specific enhancer-promoter linkages

Enhlink is a computational tool for scATAC-seq data analysis, facilitating precise interrogation of enhancer function at the single-cell level. It employs an ensemble approach incorporating technical and biological covariates to infer condition-specific regulatory DNA linkages. Enhlink can integrate multi-omic data for enhanced specificity, when available. Evaluation with simulated and real data, including multi-omic datasets from the mouse striatum and novel promoter capture Hi-C data, demonstrate that Enhlink outperfoms alternative methods. Coupled with eQTL analysis, it identified a putative super-enhancer in striatal neurons. Overall, Enhlink offers accuracy, power, and potential for revealing novel biological insights in gene regulation.

Discovery of genomic loci for liver health and steatosis reveals overlap with glutathione redox genetics

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver condition in the United States, encompassing a wide spectrum of liver pathologies including steatosis, steatohepatitis, fibrosis, and cirrhosis. Despite its high prevalence, there are no medications currently approved by the Food and Drug Administration for the treatment of NAFLD. Recent work has suggested that NAFLD has a strong genetic component and identifying causative genes will improve our understanding of the molecular mechanisms contributing to NAFLD and yield targets for future therapeutic investigations. Oxidative stress is known to play an important role in NAFLD pathogenesis, yet the underlying mechanisms accounting for disturbances in redox status are not entirely understood. To better understand the relationship between the glutathione redox system and signs of NAFLD in a genetically-diverse population, we measured liver weight, serum biomarkers aspartate aminotransferase (AST) and alanine aminotransferase (ALT), and graded liver pathology in a large cohort of Diversity Outbred mice. We compared hepatic endpoints to those of the glutathione redox system previously measured in the livers and kidneys of the same mice, and we screened for statistical and genetic associations using the R/qtl2 software. We discovered several novel genetic loci associated with markers of liver health, including loci that were associated with both liver steatosis and glutathione redox status. Candidate genes within each locus point to possible new mechanisms underlying the complex relationship between NAFLD and the glutathione redox system, which could have translational implications for future studies targeting NAFLD pathology.

Genetic drivers of age-related changes in urinary magnesium excretion

Although age-dependent alterations in urinary magnesium (Mg2+) excretion have been described, the underlying mechanism remains elusive. As heritability significantly contributes to variations in urinary Mg2+ excretion, we measured urinary Mg2+ excretion at different ages in a cohort of genetically variable Diversity Outbred (DO) mice. Compared with animals aged 6 mo, an increase in Mg2+ excretion was observed at 12 and 18 mo. Quantitative trait locus (QTL) analysis revealed an association of a locus on chromosome 10 with Mg2+ excretion at 6 mo of age, with Oit3 (encoding oncoprotein-induced transcript 3; OIT3) as our primary candidate gene. To study the possible role of OIT3 in renal Mg2+ handling, we generated and characterized Oit3 knockout (Oit3-/-) mice. Although a slightly lower serum Mg2+ concentration was present in male Oit3-/- mice, this effect was not observed in female Oit3-/- mice. In addition, urinary Mg2+ excretion and the expression of renal magnesiotropic genes were unaltered in Oit3-/- mice. For animals aged 12 and 18 mo, QTL analysis revealed an association with a locus on chromosome 19, which contains the gene encoding TRPM6, a known Mg2+ channel involved in renal Mg2+ reabsorption. Comparison with RNA sequencing (RNA-Seq) data revealed that Trpm6 mRNA expression is inversely correlated with the QTL effect, implying that TRPM6 may be involved in age-dependent changes in urinary Mg2+ excretion in mice. In conclusion, we show here that variants in Oit3 and Trpm6 are associated with urinary Mg2+ excretion at distinct periods of life, although OIT3 is unlikely to affect renal Mg2+ handling.NEW & NOTEWORTHY Aging increased urinary magnesium (Mg2+) excretion in mice. We show here that variation in Oit3, a candidate gene for the locus associated with Mg2+ excretion in young mice, is unlikely to be involved as knockout of Oit3 did not affect Mg2+ excretion. Differences in the expression of the renal Mg2+ channel TRPM6 may contribute to the variation in urinary Mg2+ excretion in older mice.

Characterization of quantitative trait loci from DJ123 (aus) independently affecting panicle structure traits in indica rice cultivar IR64

The rice panicle is the principal organ to influence productivity and traits affecting panicle architecture determine sink size and yield potential. Improving panicle architecture may be effective in increasing yield under low-input conditions, but which traits are of importance under such conditions and how they are genetically controlled is not well understood. Using recombinant inbred lines (RILs) derived from a cross between a modern variety IR64 and a low fertility tolerant accession DJ123, quantitative trait locus (QTL) mapping was conducted under high soil fertility in Japan and low fertility in Madagascar. Among QTL for panicle length (PL) detected, the DJ123 allele increased rachis length at qCL1 and qPL9, while the IR64 allele increased primary branch length at qPL7. DJ123 further contributed two QTL for grain width whereas IR64 contributed two grain length QTL. Analysis of lines carrying different combinations of detected QTL indicates that rachis and primary branch lengths are independently regulated, explaining strong transgressive segregation for PL. The positive effects of PL-related QTL were further confirmed by a genome-wide analysis of allelic states in two breeding lines that had been selected repeatedly for total panicle weight per plant under low input conditions. This study provides the genetic basis for complex panicle architecture in rice and will aid in designing an ideal panicle architecture that leads to increased yield under low fertility conditions.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-024-01494-5.

Continuous variation in the shell colour of the snail Cepaea nemoralis is associated with the colour locus of the supergene

While the shell of the land snail Cepaea nemoralis is typically classed as yellow, pink, or brown, the reality is that colour variation is continuously distributed. To further understand the origin of the continuous variation, we used crosses of C. nemoralis to compare quantitative measures of the colour with the inferred genotype of the underlying supergene locus. We also used a recently developed linkage map to find quantitative trait loci that may influence colour. The results show that the colour locus of the supergene-at around 31.385 cM on linkage group 11-is involved in determining the quantitative chromatic differences that are perceptible to human vision. We also found some evidence that variation within colour classes may be due to allelic variation at or around the supergene. There are likely other unlinked loci involved in determining colour within classes, but confirmation will require greater statistical power. Although not investigated here, environmental factors, including diet, may also impact upon variation within colour types.

Heterogeneous genomic architecture of skeletal armour traits in sticklebacks

Whether populations adapt to similar selection pressures using the same underlying genetic variants depends on population history and the distribution of standing genetic variation at the metapopulation level. Studies of sticklebacks provide a case in point: when colonizing and adapting to freshwater habitats, three-spined sticklebacks (Gasterosteus aculeatus) with high gene flow tend to fix the same adaptive alleles in the same major loci, whereas nine-spined sticklebacks (Pungitius pungitius) with limited gene flow tend to utilize a more heterogeneous set of loci. In accordance with this, we report results of quantitative trait locus (QTL) analyses using a backcross design showing that lateral plate number variation in the western European nine-spined sticklebacks mapped to 3 moderate-effect QTL, contrary to the major-effect QTL in three-spined sticklebacks and different from the 4 QTL previously identified in the eastern European nine-spined sticklebacks. Furthermore, several QTL were identified associated with variation in lateral plate size, and 3 moderate-effect QTL with body size. Together, these findings indicate more heterogenous and polygenic genetic underpinnings of skeletal armour variation in nine-spined than three-spined sticklebacks, indicating limited genetic parallelism underlying armour trait evolution in the family Gasterostidae.

Genetic background affects the strength of crossover interference in house mice

Meiotic recombination is required for faithful chromosome segregation in most sexually reproducing organisms and shapes the distribution of genetic variation in populations. Both the overall rate and the spatial distribution of crossovers vary within and between species. Adjacent crossovers on the same chromosome tend to be spaced more evenly than expected at random, a phenomenon known as crossover interference. Although interference has been observed in many taxa, the factors that influence the strength of interference are not well understood. We used house mice (Mus musculus), a well-established model system for understanding recombination, to study the effects of genetics and age on recombination rate and interference in the male germline. We analyzed crossover positions in 503 progeny from reciprocal F1 hybrids between inbred strains representing the three major subspecies of house mice. Consistent with previous studies, autosomal alleles from M. m. musculus tend to increase recombination rate, while inheriting a M. m. musculus X chromosome decreases recombination rate. Old males transmit an average of 0.6 more crossovers per meiosis (5.0%) than young males, though the effect varies across genetic backgrounds. We show that the strength of crossover interference depends on genotype, providing a rare demonstration that interference evolves over short timescales. Differences between reciprocal F1s suggest that X-linked factors modulate the strength of interference. Our findings motivate additional comparisons of interference among recently diverged species and further examination of the role of paternal age in determining the number and positioning of crossovers.

Identification of a genetic region linked to tolerance to MRSA infection using Collaborative Cross mice

Staphylococcus aureus (S. aureus) colonizes humans asymptomatically but can also cause opportunistic infections, ranging from mild skin infections to severe life-threatening conditions. Resistance and tolerance are two ways a host can survive an infection. Resistance is limiting the pathogen burden, while tolerance is limiting the health impact of a given pathogen burden. In previous work, we established that collaborative cross (CC) mouse line CC061 is highly susceptible to Methicillin-resistant S. aureus infection (MRSA, USA300), while CC024 is tolerant. To identify host genes involved in tolerance after S. aureus infection, we crossed CC061 mice and CC024 mice to generate F1 and F2 populations. Survival after MRSA infection in the F1 and F2 generations was 65% and 55% and followed a complex dominant inheritance pattern for the CC024 increased survival phenotype. Colonization in F2 animals was more extreme than in their parents, suggesting successful segregation of genetic factors. We identified a Quantitative Trait Locus (QTL) peak on chromosome 7 for survival and weight change after infection. In this QTL, the WSB/EiJ (WSB) allele was present in CC024 mice and contributed to their MRSA tolerant phenotype. Two genes, C5ar1 and C5ar2, have high-impact variants in this region. C5ar1 and C5ar2 are receptors for the complement factor C5a, an anaphylatoxin that can trigger a massive immune response by binding to these receptors. We hypothesize that C5a may have altered binding to variant receptors in CC024 mice, reducing damage caused by the cytokine storm and resulting in the ability to tolerate a higher pathogen burden and longer survival.

Optix and cortex/ivory/mir-193 again: the repeated use of two mimicry hotspot loci

The extent to which evolution is repeatable has been a debated topic among evolutionary biologists. Although rewinding the tape of life perhaps would not lead to the same outcome every time, repeated evolution of analogous genes for similar functions has been extensively reported. Wing phenotypes of butterflies and moths have provided a wealth of examples of gene re-use, with certain 'hotspot loci' controlling wing patterns across diverse taxa. Here, we present an example of convergent evolution in the molecular genetic basis of Batesian wing mimicry in two Hypolimnas butterfly species. We show that mimicry is controlled by variation near cortex/ivory/mir-193, a known butterfly hotspot locus. By dissecting the genetic architecture of mimicry in Hypolimnas misippus and Hypolimnas bolina, we present evidence that distinct non-coding regions control the development of white pattern elements in the forewing and hindwing of the two species, suggesting independent evolution, and that no structural variation is found at the locus. Finally, we also show that orange coloration in H. bolina is associated with optix, a well-known patterning gene. Overall, our study once again implicates variation near the hotspot loci cortex/ivory/mir-193 and optix in butterfly wing mimicry and thereby highlights the repeatability of adaptive evolution.

Genotype imputation in F2 crosses of inbred lines

Motivation

Crosses among inbred lines are a fundamental tool for the discovery of genetic loci associated with phenotypes of interest. In organisms for which large reference panels or SNP chips are not available, imputation from low-pass whole-genome sequencing is an effective method for obtaining genotype data from a large number of individuals. To date, a structured analysis of the conditions required for optimal genotype imputation has not been performed.

Results

We report a systematic exploration of the effect of several design variables on imputation performance in F2 crosses of inbred medaka lines using the imputation software STITCH. We determined that, depending on the number of samples, imputation performance reaches a plateau when increasing the per-sample sequencing coverage. We also systematically explored the trade-offs between cost, imputation accuracy, and sample numbers. We developed a computational pipeline to streamline the process, enabling other researchers to perform a similar cost-benefit analysis on their population of interest.

Availability and implementation

The source code for the pipeline is available at https://github.com/birneylab/stitchimpute. While our pipeline has been developed and tested for an F2 population, the software can also be used to analyse populations with a different structure.

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.

Quantitative trait locus analysis of gray leaf spot resistance in the maize IBM Syn10 DH population

KEY MESSAGE

The exploration and dissection of a set of QTLs and candidate genes for gray leaf spot disease resistance using two fully assembled parental genomes may help expedite maize resistance breeding. The fungal disease of maize known as gray leaf spot (GLS), caused by Cercospora zeae-maydis and Cercospora zeina, is a significant concern in China, Southern Africa, and the USA. Resistance to GLS is governed by multiple genes with an additive effect and is influenced by both genotype and environment. The most effective way to reduce the cost of production is to develop resistant hybrids. In this study, we utilized the IBM Syn 10 Doubled Haploid (IBM Syn10 DH) population to identify quantitative trait loci (QTLs) associated with resistance to gray leaf spot (GLS) in multiple locations. Analysis of seven distinct environments revealed a total of 58 QTLs, 49 of which formed 12 discrete clusters distributed across chromosomes 1, 2, 3, 4, 8 and 10. By comparing these findings with published research, we identified colocalized QTLs or GWAS loci within eleven clustering intervals. By integrating transcriptome data with genomic structural variations between parental individuals, we identified a total of 110 genes that exhibit both robust disparities in gene expression and structural alterations. Further analysis revealed 19 potential candidate genes encoding conserved resistance gene domains, including putative leucine-rich repeat receptors, NLP transcription factors, fucosyltransferases, and putative xyloglucan galactosyltransferases. Our results provide a valuable resource and linked loci for GLS marker resistance selection breeding in maize.

Parallel evolution of integrated craniofacial traits in trophic specialist pupfishes

Populations may adapt to similar environments via parallel or non-parallel genetic changes, but the frequency of these alternative mechanisms and underlying contributing factors are still poorly understood outside model systems. We used QTL mapping to investigate the genetic basis of highly divergent craniofacial traits between the scale-eater (Cyprinodon desquamator) and molluscivore (C. brontotheroides) pupfish adapting to two different hypersaline lake environments on San Salvador Island, Bahamas. We lab-reared F2 scale-eater x molluscivore intercrosses from two different lake populations, estimated linkage maps, scanned for significant QTL for 29 skeletal and craniofacial traits, female mate preference, and sex. We compared the location of QTL between lakes to quantify parallel and non-parallel genetic changes. We detected significant QTL for six craniofacial traits in at least one lake. However, nearly all shared QTL loci were associated with a different craniofacial trait within each lake. Therefore, our estimate of parallel evolution of craniofacial genetic architecture could range from one out of six identical trait QTL (low parallelism) to five out of six integrated trait QTL (high parallelism). We suggest that pleiotropy and trait integration can affect estimates of parallel evolution, particularly within rapid radiations. We also observed increased adaptive introgression in shared QTL regions, suggesting that gene flow contributed to parallel evolution. Overall, our results suggest that the same genomic regions may contribute to parallel adaptation across integrated suites of craniofacial traits, rather than specific traits, and highlight the need for a more expansive definition of parallel evolution.

Genetic variation modulates susceptibility to aberrant DNA hypomethylation and imprint deregulation in naïve pluripotent stem cells

Naïve pluripotent stem cells (nPSC) frequently undergo pathological and not readily reversible loss of DNA methylation marks at imprinted gene loci. This abnormality poses a hurdle for using pluripotent cell lines in biomedical applications and underscores the need to identify the causes of imprint instability in these cells. We show that nPSCs from inbred mouse strains exhibit pronounced strain-specific susceptibility to locus-specific deregulation of imprinting marks during reprogramming to pluripotency and upon culture with MAP kinase inhibitors, a common approach to maintain naïve pluripotency. Analysis of genetically highly diverse nPSCs from the Diversity Outbred (DO) stock confirms that genetic variation is a major determinant of epigenome stability in pluripotent cells. We leverage the variable DNA hypomethylation in DO lines to identify several trans-acting quantitative trait loci (QTLs) that determine epigenome stability at either specific target loci or genome-wide. Candidate factors encoded by two multi-target QTLs on chromosomes 4 and 17 suggest specific transcriptional regulators that contribute to DNA methylation maintenance in nPSCs. We propose that genetic variants represent candidate biomarkers to identify pluripotent cell lines with desirable properties and might serve as entry points for the targeted engineering of nPSCs with stable epigenomes.

Highlights

Naïve pluripotent stem cells from distinct inbred mouse strains exhibit variable DNA methylation levels at imprinted gene loci.The vulnerability of pluripotent stem cells to loss of genomic imprinting caused by MAP kinase inhibition strongly differs between inbred mouse strains.Genetically diverse pluripotent stem cell lines from Diversity Outbred mouse stock allow the identification of quantitative trait loci controlling DNA methylation stability.Genetic variants may serve as biomarkers to identify naïve pluripotent stem cell lines that are epigenetically stable in specific culture conditions.

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.

Systems genetics of metabolic health in the BXD mouse genetic reference population

Susceptibility to metabolic syndrome (MetS) is dependent on genetics, environment, and gene-by-environment interactions, rendering the study of underlying mechanisms challenging. The majority of experiments in model organisms do not incorporate genetic variation and lack specific evaluation criteria for MetS. Here, we derived a continuous metric, the metabolic health score (MHS), based on standard clinical parameters and defined its molecular signatures in the liver and circulation. In human UK Biobank, the MHS associated with MetS status and was predictive of future disease incidence, even in individuals without MetS. Using quantitative trait locus analyses in mice, we found two MHS-associated genetic loci and replicated them in unrelated mouse populations. Through a prioritization scheme in mice and human genetic data, we identified TNKS and MCPH1 as candidates mediating differences in the MHS. Our findings provide insights into the molecular mechanisms sustaining metabolic health across species and uncover likely regulators. A record of this paper's transparent peer review process is included in the supplemental information.

Identification of ethanol analgesia quantitative trait loci and candidate genes in BXD recombinant inbred mouse lines

Alcohol consumption produces acute analgesic effects, and people experiencing pain conditions may drink alcohol to alleviate discomfort. However, tolerance to the analgesic properties of alcohol could prompt escalating consumption and dependence. Both nociception and alcohol-induced analgesia are under significant genetic control. Understanding the genetic architecture of these processes could inform better treatment options for people with pain conditions. This study aims to identify quantitative trait loci (QTL) driving variation in ethanol-induced analgesia across BXD recombinant inbred mouse lines. Male and female mice from 62 BXD strains received ethanol or saline oral gavage for five days and were tested for hot plate (HP) latency at baseline, Day 1, and Day 5. QTL mapping of HP phenotypes identified a significant provisional QTL on chromosome 17 for Day 1 HP latency in mice receiving ethanol. An additional highly suggestive QTL was present on chromosome 9 for the difference in pre- and post-ethanol thermal nociception. Candidate genes within QTL support intervals were provisionally identified using HP phenotypic correlations to transcriptomic database, expression QTL analysis, and other bioinformatics inquiries. The combined behavioral and bioinformatic analyses yielded strong ethanol analgesia candidate genes, specifically Myo6. Thus, the results of this genetic study of ethanol-induced analgesia in BXD mouse strains may contribute significantly to our understanding of the molecular basis for individual variation in the analgesic response to acute ethanol.

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

Fine Mapping of QTLs for Alkaline Tolerance in Crucian Carp (Carassius auratus) Using Genome-Wide SNP Markers

Crucian carp (Carassius auratus) is widely distributed in the world and has become an economically freshwater fish. The population in Lake Dali Nur can tolerate the extreme alkaline environment with alkalinity over 50 mmol/L (pH 9.6), thus providing a special model for exploring alkali-tolerant molecular markers in an extremely alkaline environment. In this study, we constructed a high-density and high-resolution linkage map with 16,224 SNP markers based on genotyping-by-sequencing (GBS) consisting of 152 progenies and conducted QTL studies for alkali-tolerant traits. The total length of the linkage map was 3918.893 cM, with an average distance of 0.241 cM. Two QTLs for the ammonia-N-tolerant trait were detected on LG27 and LG45. A QTL for the urea-N-tolerant trait was detected on LG27. Interestingly, mapping the two QTLs on LG27 revealed that the mapped genes were both located in the intron of CDC42. GO functional annotation and KEGG enrichment analysis results indicated that the biological functions might be involved in the cell cycle, cellular senescence, MAPK, and Ras signaling pathways. These findings suggest that CDC42 may play an important role in the process of dealing with extremely alkaline environments.

A PLURIPOTENT STEM CELL PLATFORM FOR IN VITRO SYSTEMS GENETICS STUDIES OF MOUSE DEVELOPMENT

The directed differentiation of pluripotent stem cells (PSCs) from panels of genetically diverse individuals is emerging as a powerful experimental system for characterizing the impact of natural genetic variation on developing cell types and tissues. Here, we establish new PSC lines and experimental approaches for modeling embryonic development in a genetically diverse, outbred mouse stock (Diversity Outbred mice). We show that a range of inbred and outbred PSC lines can be stably maintained in the primed pluripotent state (epiblast stem cells -- EpiSCs) and establish the contribution of genetic variation to phenotypic differences in gene regulation and directed differentiation. Using pooled in vitro fertilization, we generate and characterize a genetic reference panel of Diversity Outbred PSCs (n = 230). Finally, we demonstrate the feasibility of pooled culture of Diversity Outbred EpiSCs as "cell villages", which can facilitate the differentiation of large numbers of EpiSC lines for forward genetic screens. These data can complement and inform similar efforts within the stem cell biology and human genetics communities to model the impact of natural genetic variation on phenotypic variation and disease-risk.

Probing the basis of disease heterogeneity in multiple sclerosis using genetically diverse mice

Multiple sclerosis (MS) is a complex disease with significant heterogeneity in disease course and progression. Genetic studies have identified numerous loci associated with MS risk, but the genetic basis of disease progression remains elusive. To address this, we leveraged the Collaborative Cross (CC), a genetically diverse mouse strain panel, and experimental autoimmune encephalomyelitis (EAE). The thirty-two CC strains studied captured a wide spectrum of EAE severity, trajectory, and presentation, including severe-progressive, monophasic, relapsing remitting, and axial rotary (AR)-EAE, accompanied by distinct immunopathology. Sex differences in EAE severity were observed in six strains. Quantitative trait locus analysis revealed distinct genetic linkage patterns for different EAE phenotypes, including EAE severity and incidence of AR-EAE. Machine learning-based approaches prioritized candidate genes for loci underlying EAE severity ( Abcc4 and Gpc6 ) and AR-EAE ( Yap1 and Dync2h1 ). This work expands the EAE phenotypic repertoire and identifies novel loci controlling unique EAE phenotypes, supporting the hypothesis that heterogeneity in MS disease course is driven by genetic variation.

Summary

The genetic basis of disease heterogeneity in multiple sclerosis (MS) remains elusive. We leveraged the Collaborative Cross to expand the phenotypic repertoire of the experimental autoimmune encephalomyelitis (EAE) model of MS and identify loci controlling EAE severity, trajectory, and presentation.

Systems genetics uncover new loci containing functional gene candidates in Mycobacterium tuberculosis-infected Diversity Outbred mice

Mycobacterium tuberculosis infects two billion people across the globe, and results in 8-9 million new tuberculosis (TB) cases and 1-1.5 million deaths each year. Most patients have no known genetic basis that predisposes them to disease. Here, we investigate the complex genetic basis of pulmonary TB by modelling human genetic diversity with the Diversity Outbred mouse population. When infected with M. tuberculosis, one-third develop early onset, rapidly progressive, necrotizing granulomas and succumb within 60 days. The remaining develop non-necrotizing granulomas and survive longer than 60 days. Genetic mapping using immune and inflammatory mediators; and clinical, microbiological, and granuloma correlates of disease identified five new loci on mouse chromosomes 1, 2, 4, 16; and three known loci on chromosomes 3 and 17. Further, multiple positively correlated traits shared loci on chromosomes 1, 16, and 17 and had similar patterns of allele effects, suggesting these loci contain critical genetic regulators of inflammatory responses to M. tuberculosis. To narrow the list of candidate genes, we used a machine learning strategy that integrated gene expression signatures from lungs of M. tuberculosis-infected Diversity Outbred mice with gene interaction networks to generate scores representing functional relationships. The scores were used to rank candidates for each mapped trait, resulting in 11 candidate genes: Ncf2, Fam20b, S100a8, S100a9, Itgb5, Fstl1, Zbtb20, Ddr1, Ier3, Vegfa, and Zfp318. Although all candidates have roles in infection, inflammation, cell migration, extracellular matrix remodeling, or intracellular signaling, and all contain single nucleotide polymorphisms (SNPs), SNPs in only four genes (S100a8, Itgb5, Fstl1, Zfp318) are predicted to have deleterious effects on protein functions. We performed methodological and candidate validations to (i) assess biological relevance of predicted allele effects by showing that Diversity Outbred mice carrying PWK/PhJ alleles at the H-2 locus on chromosome 17 QTL have shorter survival; (ii) confirm accuracy of predicted allele effects by quantifying S100A8 protein in inbred founder strains; and (iii) infection of C57BL/6 mice deficient for the S100a8 gene. Overall, this body of work demonstrates that systems genetics using Diversity Outbred mice can identify new (and known) QTLs and functionally relevant gene candidates that may be major regulators of complex host-pathogens interactions contributing to granuloma necrosis and acute inflammation in pulmonary TB.

Cell type-specific network analysis in Diversity Outbred mice identifies genes potentially responsible for human bone mineral density GWAS associations

Genome-wide association studies (GWASs) have identified many sources of genetic variation associated with bone mineral density (BMD), a clinical predictor of fracture risk and osteoporosis. Aside from the identification of causal genes, other difficult challenges to informing GWAS include characterizing the roles of predicted causal genes in disease and providing additional functional context, such as the cell type predictions or biological pathways in which causal genes operate. Leveraging single-cell transcriptomics (scRNA-seq) can assist in informing BMD GWAS by linking disease-associated variants to genes and providing a cell type context for which these causal genes drive disease. Here, we use large-scale scRNA-seq data from bone marrow-derived stromal cells cultured under osteogenic conditions (BMSC-OBs) from Diversity Outbred (DO) mice to generate cell type-specific networks and contextualize BMD GWAS-implicated genes. Using trajectories inferred from the scRNA-seq data, we identify networks enriched with genes that exhibit the most dynamic changes in expression across trajectories. We discover 21 network driver genes, which are likely to be causal for human BMD GWAS associations that colocalize with expression/splicing quantitative trait loci (eQTL/sQTL). These driver genes, including Fgfrl1 and Tpx2, along with their associated networks, are predicted to be novel regulators of BMD via their roles in the differentiation of mesenchymal lineage cells. In this work, we showcase the use of single-cell transcriptomics from mouse bone-relevant cells to inform human BMD GWAS and prioritize genetic targets with potential causal roles in the development of osteoporosis.

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.

QTL for the Kinematic Traits That Define the Arabidopsis Root Elongation Zone and Their Relationship to Gravitropism

Cell expansion in a discrete region called the elongation zone drives root elongation. Analyzing time lapse images can quantify the expansion in kinematic terms as if it were fluid flow. We used horizontal microscopes to collect images from which custom software extracted the length of the elongation zone, the peak relative elemental growth rate (REGR) within it, the axial position of the REGR peak, and the root elongation rate. Automation enabled these kinematic traits to be measured in 1575 Arabidopsis seedlings representing 162 recombinant inbred lines (RILs) derived from a cross of Cvi and Ler ecotypes. We mapped ten quantitative trait loci (QTL), affecting the four kinematic traits. Three QTL affected two or more traits in these vertically oriented seedlings. We compared this genetic architecture with that previously determined for gravitropism using the same RIL population. The major QTL peaks for the kinematic traits did not overlap with the gravitropism QTL. Furthermore, no single kinematic trait correlated with quantitative descriptors of the gravitropism response curve across this population. In addition to mapping QTL for growth zone traits, this study showed that the size and shape of the elongation zone may vary widely without affecting the differential growth induced by gravity.