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

A novel quantitative trait locus implicates Msh3 in the propensity for genome-wide short tandem repeat expansions in mice

Short tandem repeats (STRs) are a class of rapidly mutating genetic elements typically characterized by repeated units of 1-6 bp. We leveraged whole-genome sequencing data for 152 recombinant inbred (RI) strains from the BXD family of mice to map loci that modulate genome-wide patterns of new mutations arising during parent-to-offspring transmission at STRs. We defined quantitative phenotypes describing the numbers and types of germline STR mutations in each strain and performed quantitative trait locus (QTL) analyses for each of these phenotypes. We identified a locus on Chromosome 13 at which strains inheriting the C57BL/6J (B) haplotype have a higher rate of STR expansions than those inheriting the DBA/2J (D) haplotype. The strongest candidate gene in this locus is Msh3, a known modifier of STR stability in cancer and at pathogenic repeat expansions in mice and humans, as well as a current drug target against Huntington's disease. The D haplotype at this locus harbors a cluster of variants near the 5' end of Msh3, including multiple missense variants near the DNA mismatch recognition domain. In contrast, the B haplotype contains a unique retrotransposon insertion. The rate of expansion covaries positively with Msh3 expression-with higher expression from the B haplotype. Finally, detailed analysis of mutation patterns showed that strains carrying the B allele have higher expansion rates, but slightly lower overall total mutation rates, compared with those with the D allele, particularly at tetranucleotide repeats. Our results suggest an important role for inherited variants in Msh3 in modulating genome-wide patterns of germline mutations at STRs.

Identification of actionable targets for breast cancer intervention using a diversity outbred mouse model

HER2-targeted therapy has improved breast cancer survival, but treatment resistance and disease prevention remain major challenges. Genes that enable HER2/Neu oncogenesis are the next intervention targets. A bioinformatics discovery platform of HER2/Neu-expressing Diversity Outbred (DO) F1 Mice was established to identify cancer-enabling genes. Quantitative Trait Loci (QTL) associated with onset ages and growth rates of spontaneous mammary tumors were sought. Twenty-six genes in 3 QTL contain sequence variations unique to the genetic backgrounds that are linked to aggressive tumors and 21 genes are associated with human breast cancer survival. Concurrent identification of TSC22D3, a transcription factor, and its target gene LILRB4, a myeloid cell checkpoint receptor, suggests an immune axis for regulation, or intervention, of disease. We also investigated TIEG1 gene that impedes tumor immunity but suppresses tumor growth. Although not an actionable target, TIEG1 study revealed genetic regulation of tumor progression, forming the basis of the genetics-based discovery platform.

Genetic dissection of the pluripotent proteome through multi-omics data integration

Genetic background drives phenotypic variability in pluripotent stem cells (PSCs). Most studies to date have used transcript abundance as the primary molecular readout of cell state in PSCs. We performed a comprehensive proteogenomics analysis of 190 genetically diverse mouse embryonic stem cell (mESC) lines. The quantitative proteome is highly variable across lines, and we identified pluripotency-associated pathways that were differentially activated in the proteomics data that were not evident in transcriptome data from the same lines. Integration of protein abundance to transcript levels and chromatin accessibility revealed broad co-variation across molecular layers as well as shared and unique drivers of quantitative variation in pluripotency-associated pathways. Quantitative trait locus (QTL) mapping localized the drivers of these multi-omic signatures to genomic hotspots. This study reveals post-transcriptional mechanisms and genetic interactions that underlie quantitative variability in the pluripotent proteome and provides a regulatory map for mESCs that can provide a basis for future mechanistic studies.

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

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

Functional genomic effects of indels using Bayesian genome-phenome wide association studies in sorghum

High-throughput genomic and phenomic data have enhanced the ability to detect genotype-to-phenotype associations that can resolve broad pleiotropic effects of mutations on plant phenotypes. As the scale of genotyping and phenotyping has advanced, rigorous methodologies have been developed to accommodate larger datasets and maintain statistical precision. However, determining the functional effects of associated genes/loci is expensive and limited due to the complexity associated with cloning and subsequent characterization. Here, we utilized phenomic imputation of a multi-year, multi-environment dataset using PHENIX which imputes missing data using kinship and correlated traits, and we screened insertions and deletions (InDels) from the recently whole-genome sequenced Sorghum Association Panel for putative loss-of-function effects. Candidate loci from genome-wide association results were screened for potential loss of function using a Bayesian Genome-Phenome Wide Association Study (BGPWAS) model across both functionally characterized and uncharacterized loci. Our approach is designed to facilitate in silico validation of associations beyond traditional candidate gene and literature-search approaches and to facilitate the identification of putative variants for functional analysis and reduce the incidence of false-positive candidates in current functional validation methods. Using this Bayesian GPWAS model, we identified associations for previously characterized genes with known loss-of-function alleles, specific genes falling within known quantitative trait loci, and genes without any previous genome-wide associations while additionally detecting putative pleiotropic effects. In particular, we were able to identify the major tannin haplotypes at the Tan1 locus and effects of InDels on the protein folding. Depending on the haplotype present, heterodimer formation with Tan2 was significantly affected. We also identified major effect InDels in Dw2 and Ma1, where proteins were truncated due to frameshift mutations that resulted in early stop codons. These truncated proteins also lost most of their functional domains, suggesting that these indels likely result in loss of function. Here, we show that the Bayesian GPWAS model is able to identify loss-of-function alleles that can have significant effects upon protein structure and folding as well as multimer formation. Our approach to characterize loss-of-function mutations and their functional repercussions will facilitate precision genomics and breeding by identifying key targets for gene editing and trait integration.

Lipidomic QTL in Diversity Outbred mice identifies a novel function for α/β hydrolase domain 2 ( Abhd2 ) as an enzyme that metabolizes phosphatidylcholine and cardiolipin

We and others have previously shown that genetic association can be used to make causal connections between gene loci and small molecules measured by mass spectrometry in the bloodstream and in tissues. We identified a locus on mouse chromosome 7 where several phospholipids in liver showed strong genetic association to distinct gene loci. In this study, we integrated gene expression data with genetic association data to identify a single gene at the chromosome 7 locus as the driver of the phospholipid phenotypes. The gene encodes α/β-hydrolase domain 2 ( Abhd2 ), one of 23 members of the ABHD gene family. We validated this observation by measuring lipids in a mouse with a whole-body deletion of Abhd2 . The Abhd2 KO mice had a significant increase in liver levels of phosphatidylcholine and phosphatidylethanolamine. Unexpectedly, we also found a decrease in two key mitochondrial lipids, cardiolipin and phosphatidylglycerol, in male Abhd2 KO mice. These data suggest that Abhd2 plays a role in the synthesis, turnover, or remodeling of liver phospholipids.

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.

Genome-wide screen identifies host loci that modulate M. tuberculosis fitness in immunodivergent mice

Genetic differences among mammalian hosts and Mycobacterium tuberculosis ( Mtb ) strains determine diverse tuberculosis (TB) patient outcomes. The advent of recombinant inbred mouse panels and next-generation transposon mutagenesis and sequencing approaches has enabled dissection of complex host- pathogen interactions. To identify host and pathogen genetic determinants of Mtb pathogenesis, we infected members of the BXD family of mouse strains with a comprehensive library of Mtb transposon mutants (TnSeq). Members of the BXD family segregate for Mtb -resistant C57BL/6J (B6 or B ) and Mtb -susceptible DBA/2J (D2 or D ) haplotypes. The survival of each bacterial mutant was quantified within each BXD host, and we identified those bacterial genes that were differentially required for Mtb fitness across BXD genotypes. Mutants that varied in survival among the host family of strains were leveraged as reporters for "endophenotypes", each bacterial fitness profile directly probing specific components of the infection microenvironment. We conducted QTL mapping of these bacterial fitness endophenotypes and identified 140 h ost- p athogen quantitative trait loci ( hp QTL). We identified a QTL hotspot on chromosome 6 (75.97-88.58 Mb) associated with the genetic requirement of multiple Mtb genes; Rv0127 ( mak ), Rv0359 ( rip2 ), Rv0955 ( perM ), and Rv3849 ( espR ). Together, this screen reinforces the utility of bacterial mutant libraries as precise reporters of the host immunological microenvironment during infection and highlights specific host-pathogen genetic interactions for further investigation. To enable downstream follow-up for both bacterial and mammalian genetic research communities, all bacterial fitness profiles have been deposited into GeneNetwork.org and added into the comprehensive collection of TnSeq libraries in MtbTnDB.

Genetic mapping of microbial and host traits reveals production of immunomodulatory lipids by Akkermansia muciniphila in the murine gut

The molecular bases of how host genetic variation impacts the gut microbiome remain largely unknown. Here we used a genetically diverse mouse population and applied systems genetics strategies to identify interactions between host and microbe phenotypes including microbial functions, using faecal metagenomics, small intestinal transcripts and caecal lipids that influence microbe-host dynamics. Quantitative trait locus (QTL) mapping identified murine genomic regions associated with variations in bacterial taxa; bacterial functions including motility, sporulation and lipopolysaccharide production and levels of bacterial- and host-derived lipids. We found overlapping QTL for the abundance of Akkermansia muciniphila and caecal levels of ornithine lipids. Follow-up in vitro and in vivo studies revealed that A. muciniphila is a major source of these lipids in the gut, provided evidence that ornithine lipids have immunomodulatory effects and identified intestinal transcripts co-regulated with these traits including Atf3, which encodes for a transcription factor that plays vital roles in modulating metabolism and immunity. Collectively, these results suggest that ornithine lipids are potentially important for A. muciniphila-host interactions and support the role of host genetics as a determinant of responses to gut microbes.

Rate of tau propagation is a heritable disease trait in genetically diverse mouse strains

The speed and scope of cognitive deterioration in Alzheimer's disease is highly associated with the advancement of tau neurofibrillary lesions across brain networks. We tested whether the rate of tau propagation is a heritable disease trait in a large, well-characterized cohort of genetically divergent mouse strains. Using an AAV-based model system, P301L-mutant human tau (hTau) was introduced into the entorhinal cortex of mice derived from 18 distinct lines. The extent of tau propagation was measured by distinguishing hTau-producing cells from neurons that were recipients of tau transfer. Heritability calculation revealed that 43% of the variability in tau spread was due to genetic variants segregating across background strains. Strain differences in glial markers were also observed, but did not correlate with tau propagation. Identifying unique genetic variants that influence the progression of pathological tau may uncover novel molecular targets to prevent or slow the pace of tau spread and cognitive decline.

Sample multiplexing-based targeted pathway proteomics with real-time analytics reveals the impact of genetic variation on protein expression

Targeted proteomics enables hypothesis-driven research by measuring the cellular expression of protein cohorts related by function, disease, or class after perturbation. Here, we present a pathway-centric approach and an assay builder resource for targeting entire pathways of up to 200 proteins selected from >10,000 expressed proteins to directly measure their abundances, exploiting sample multiplexing to increase throughput by 16-fold. The strategy, termed GoDig, requires only a single-shot LC-MS analysis, ~1 µg combined peptide material, a list of up to 200 proteins, and real-time analytics to trigger simultaneous quantification of up to 16 samples for hundreds of analytes. We apply GoDig to quantify the impact of genetic variation on protein expression in mice fed a high-fat diet. We create several GoDig assays to quantify the expression of multiple protein families (kinases, lipid metabolism- and lipid droplet-associated proteins) across 480 fully-genotyped Diversity Outbred mice, revealing protein quantitative trait loci and establishing potential linkages between specific proteins and lipid homeostasis.

Mouse genomic associations with in vitro sensitivity to simulated space radiation

Exposure to ionizing radiation is considered by NASA to be a major health hazard for deep space exploration missions. Ionizing radiation sensitivity is modulated by both genomic and environmental factors. Understanding their contributions is crucial for designing experiments in model organisms, evaluating the risk of deep space (i.e. high-linear energy transfer, or LET, particle) radiation exposure in astronauts, and also selecting therapeutic irradiation regimes for cancer patients. We identified single nucleotide polymorphisms in 15 strains of mice, including 10 collaborative cross model strains and 5 founder strains, associated with spontaneous and ionizing radiation-induced in vitro DNA damage quantified based on immunofluorescent tumor protein p53 binding protein (53BP1) positive nuclear foci. Statistical analysis suggested an association with pathways primarily related to cellular signaling, metabolism, tumorigenesis and nervous system damage. We observed different genomic associations in early (4 and 8 h) responses to different LET radiation, while later (24 hour) DNA damage responses showed a stronger overlap across all LETs. Furthermore, a subset of pathways was associated with spontaneous DNA damage, suggesting 53BP1 positive foci as a potential biomarker for DNA integrity in mouse models. Our results suggest several mouse strains as new models to further study the impact of ionizing radiation and validate the identified genetic loci. We also highlight the importance of future human in vitro studies to refine the association of genes and pathways with the DNA damage response to ionizing radiation and identify targets for space travel countermeasures.

Inferring causal structures of gut microbiota diversity and feed efficiency traits in poultry using Bayesian learning and genomic structural equation models

Feed and phosphorus (P) efficiency are of increasing importance in poultry breeding. It has been shown recently that these efficiency traits are influenced by the gut microbiota composition of the birds. The efficiency traits and the gut microbiota composition are partly under control of the host genome. Thus, the gut microbiota composition can be seen as a mediator trait between the host genome and the efficiency traits. The present study used data from 749 individuals of a Japanese quail F2 cross. The birds were genotyped for 4k single-nucleotide polymorphism (SNP) and trait recorded for P utilization (PU) and P retention (PR), body weight gain (BWG), and feed per gain ratio (F:G). The gut microbiota composition was characterized by targeted amplicon sequencing. The alpha diversity was calculated as the Pielou's evenness index (J'). A stable Bayesian network was established using a Hill-Climbing learning algorithm. Pielou's evenness index was placed as the most upstream trait and BWG as the most downstream trait, with direct and indirect links via PR, PU, and F:G. The direct and indirect effects between J', PU, and PR were quantified with structural equation models (SEM), which revealed a causal link from J' to PU and from PU to PR. Quantitative trait loci (QTL) linkage mapping revealed three genome-wide significant QTL regions for these traits with in total 49 trait-associated SNP within the QTL regions. SEM association mapping separated the total SNP effect for a trait into a direct effect and indirect effects mediated by upstream traits. Although the indirect effects were in general small, they contributed to the total SNP effect in some cases. This enabled us to detect some shared genetic effects. The method applied allows for the detection of shared genetic architecture of quantitative traits and microbiota compositions.

Genetic Mapping of Multiple Traits Identifies Novel Genes for Adiposity, Lipids, and Insulin Secretory Capacity in Outbred Rats

Despite the successes of human genome-wide association studies, the causal genes underlying most metabolic traits remain unclear. We used outbred heterogeneous stock (HS) rats, coupled with expression data and mediation analysis, to identify quantitative trait loci (QTLs) and candidate gene mediators for adiposity, glucose tolerance, serum lipids, and other metabolic traits. Physiological traits were measured in 1,519 male HS rats, with liver and adipose transcriptomes measured in >410 rats. Genotypes were imputed from low-coverage whole-genome sequencing. Linear mixed models were used to detect physiological and expression QTLs (pQTLs and eQTLs, respectively), using both single nucleotide polymorphism (SNP)- and haplotype-based models for pQTL mapping. Genes with cis-eQTLs that overlapped pQTLs were assessed as causal candidates through mediation analysis. We identified 14 SNP-based pQTLs and 19 haplotype-based pQTLs, of which 10 were in common. Using mediation, we identified the following genes as candidate mediators of pQTLs: Grk5 for fat pad weight and serum triglyceride pQTLs on Chr1, Krtcap3 for fat pad weight and serum triglyceride pQTLs on Chr6, Ilrun for a fat pad weight pQTL on Chr20, and Rfx6 for a whole pancreatic insulin content pQTL on Chr20. Furthermore, we verified Grk5 and Ktrcap3 using gene knockdown/out models, thereby shedding light on novel regulators of obesity.

Covariation of brain and skull shapes as a model to understand the role of crosstalk in development and evolution

Covariation among discrete phenotypes can arise due to selection for shared functions, and/or shared genetic and developmental underpinnings. The consequences of such phenotypic integration are far-reaching and can act to either facilitate or limit morphological variation. The vertebrate brain is known to act as an "organizer" of craniofacial development, secreting morphogens that can affect the shape of the growing neurocranium, consistent with roles for pleiotropy in brain-neurocranium covariation. Here, we test this hypothesis in cichlid fishes by first examining the degree of shape integration between the brain and the neurocranium using three-dimensional geometric morphometrics in an F5 hybrid population, and then genetically mapping trait covariation using quantitative trait loci (QTL) analysis. We observe shape associations between the brain and the neurocranium, a pattern that holds even when we assess associations between the brain and constituent parts of the neurocranium: the rostrum and braincase. We also recover robust genetic signals for both hard- and soft-tissue traits and identify a genomic region where QTL for the brain and braincase overlap, implicating a role for pleiotropy in patterning trait covariation. Fine mapping of the overlapping genomic region identifies a candidate gene, notch1a, which is known to be involved in patterning skeletal and neural tissues during development. Taken together, these data offer a genetic hypothesis for brain-neurocranium covariation, as well as a potential mechanism by which behavioral shifts may simultaneously drive rapid change in neuroanatomy and craniofacial morphology.

Genetic control of the pluripotency epigenome determines differentiation bias in mouse embryonic stem cells

Genetically diverse pluripotent stem cells display varied, heritable responses to differentiation cues. Here, we harnessed these disparities through derivation of mouse embryonic stem cells from the BXD genetic reference panel, along with C57BL/6J (B6) and DBA/2J (D2) parental strains, to identify loci regulating cell state transitions. Upon transition to formative pluripotency, B6 stem cells quickly dissolved naïve networks adopting gene expression modules indicative of neuroectoderm lineages, whereas D2 retained aspects of naïve pluripotency. Spontaneous formation of embryoid bodies identified divergent differentiation where B6 showed a propensity toward neuroectoderm and D2 toward definitive endoderm. Genetic mapping identified major trans-acting loci co-regulating chromatin accessibility and gene expression in both naïve and formative pluripotency. These loci distally modulated occupancy of pluripotency factors at hundreds of regulatory elements. One trans-acting locus on Chr 12 primarily impacted chromatin accessibility in embryonic stem cells, while in epiblast-like cells, the same locus subsequently influenced expression of genes enriched for neurogenesis, suggesting early chromatin priming. These results demonstrate genetically determined biases in lineage commitment and identify major regulators of the pluripotency epigenome.

Behavioral phenotypes revealed during reversal learning are linked with novel genetic loci in diversity outbred mice

Impulsive behavior and impulsivity are heritable phenotypes that are strongly associated with risk for substance use disorders. Identifying the neurogenetic mechanisms that influence impulsivity may also reveal novel biological insights into addiction vulnerability. Our past studies using the BXD and Collaborative Cross (CC) recombinant inbred mouse panels have revealed that behavioral indicators of impulsivity measured in a reversal-learning task are heritable and are genetically correlated with aspects of intravenous cocaine self-administration. Genome-wide linkage studies in the BXD panel revealed a quantitative trait locus (QTL) on chromosome 10, but we expect to identify additional QTL by testing in a population with more genetic diversity. To this end, we turned to Diversity Outbred (DO) mice; 392 DO mice (156 males, 236 females) were phenotyped using the same reversal learning test utilized previously. Our primary indicator of impulsive responding, a measure that isolates the relative difficulty mice have with reaching performance criteria under reversal conditions, revealed a genome-wide significant QTL on chromosome 7 (max LOD score = 8.73, genome-wide corrected p<0.05). A measure of premature responding akin to that implemented in the 5-choice serial reaction time task yielded a suggestive QTL on chromosome 17 (max LOD score = 9.14, genome-wide corrected <0.1). Candidate genes were prioritized (2900076A07Rik, Wdr73 and Zscan2) based upon expression QTL data we collected in DO and CC mice and analyses using publicly available gene expression and phenotype databases. These findings may advance understanding of the genetics that drive impulsive behavior and enhance risk for substance use disorders.

Architecture of the outbred brown fat proteome defines regulators of metabolic physiology

Brown adipose tissue (BAT) regulates metabolic physiology. However, nearly all mechanistic studies of BAT protein function occur in a single inbred mouse strain, which has limited the understanding of generalizable mechanisms of BAT regulation over physiology. Here, we perform deep quantitative proteomics of BAT across a cohort of 163 genetically defined diversity outbred mice, a model that parallels the genetic and phenotypic variation found in humans. We leverage this diversity to define the functional architecture of the outbred BAT proteome, comprising 10,479 proteins. We assign co-operative functions to 2,578 proteins, enabling systematic discovery of regulators of BAT. We also identify 638 proteins that correlate with protection from, or sensitivity to, at least one parameter of metabolic disease. We use these findings to uncover SFXN5, LETMD1, and ATP1A2 as modulators of BAT thermogenesis or adiposity, and provide OPABAT as a resource for understanding the conserved mechanisms of BAT regulation over metabolic physiology.

Intramembranous bone regeneration in diversity outbred mice is heritable

There are over one million cases of failed bone repair in the U.S. annually, resulting in substantial patient morbidity and societal costs. Multiple candidate genes affecting bone traits such as bone mineral density have been identified in human subjects and animal models using genome-wide association studies (GWAS). This approach for understanding the genetic factors affecting bone repair is impractical in human subjects but could be performed in a model organism if there is sufficient variability and heritability in the bone regeneration response. Diversity Outbred (DO) mice, which have significant genetic diversity and have been used to examine multiple intact bone traits, would be an excellent possibility. Thus, we sought to evaluate the phenotypic distribution of bone regeneration, sex effects and heritability of intramembranous bone regeneration on day 7 following femoral marrow ablation in 47 12-week old DO mice (23 males, 24 females). Compared to a previous study using 4 inbred mouse strains, we found similar levels of variability in the amount of regenerated bone (coefficient of variation of 86 % v. 88 %) with approximately the same degree of heritability (0.42 v. 0.49). There was a trend toward more bone regeneration in males than females. The amount of regenerated bone was either weakly or not correlated with bone mass at intact sites, suggesting that the genetic factors responsible for bone regeneration and intact bone phenotypes are at least partially independent. In conclusion, we demonstrate that DO mice exhibit variation and heritability of intramembranous bone regeneration that will be suitable for future GWAS.

The regulatory landscape of multiple brain regions in outbred heterogeneous stock rats

Heterogeneous Stock (HS) rats are a genetically diverse outbred rat population that is widely used for studying genetics of behavioral and physiological traits. Mapping Quantitative Trait Loci (QTL) associated with transcriptional changes would help to identify mechanisms underlying these traits. We generated genotype and transcriptome data for five brain regions from 88 HS rats. We identified 21 392 cis-QTLs associated with expression and splicing changes across all five brain regions and validated their effects using allele specific expression data. We identified 80 cases where eQTLs were colocalized with genome-wide association study (GWAS) results from nine physiological traits. Comparing our dataset to human data from the Genotype-Tissue Expression (GTEx) project, we found that the HS rat data yields twice as many significant eQTLs as a similarly sized human dataset. We also identified a modest but highly significant correlation between genetic regulatory variation among orthologous genes. Surprisingly, we found less genetic variation in gene regulation in HS rats relative to humans, though we still found eQTLs for the orthologs of many human genes for which eQTLs had not been found. These data are available from the RatGTEx data portal (RatGTEx.org) and will enable new discoveries of the genetic influences of complex traits.

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

Salt tolerance QTLs of an endemic rice landrace, Horkuch at seedling and reproductive stages

Salinity has a significant negative impact on production of rice. To cope with the increased soil salinity due to climate change, we need to develop salt tolerant rice varieties that can maintain their high yield. Rice landraces indigenous to coastal Bangladesh can be a great resource to study the genetic basis of salt adaptation. In this study, we implemented a QTL analysis framework with a reciprocal mapping population developed from a salt tolerant landrace Horkuch and a high yielding rice variety IR29. Our aim was to detect genetic loci that contributes to the salt adaptive responses of the two different developmental stages of rice which are very sensitive to salinity stress. We identified 14 QTLs for 9 traits and found that most are unique to specific developmental stages. In addition, we detected a significant effect of the cytoplasmic genome on the QTL model for some traits such as leaf total potassium and filled grain weight. This underscores the importance of considering cytoplasm-nuclear interaction for breeding programs. Finally, we identified QTLs co-localization for multiple traits that highlights the possible constraint of multiple QTL selection for breeding programs due to different contributions of a donor allele for different traits.

A Pleiotropic Flowering Time QTL Exhibits Gene-by-Environment Interaction for Fitness in a Perennial Grass

Appropriate flowering time is a crucial adaptation impacting fitness in natural plant populations. Although the genetic basis of flowering variation has been extensively studied, its mechanisms in nonmodel organisms and its adaptive value in the field are still poorly understood. Here, we report new insights into the genetic basis of flowering time and its effect on fitness in Panicum hallii, a native perennial grass. Genetic mapping in populations derived from inland and coastal ecotypes identified flowering time quantitative trait loci (QTL) and many exhibited extensive QTL-by-environment interactions. Patterns of segregation within recombinant hybrids provide strong support for directional selection driving ecotypic divergence in flowering time. A major QTL on chromosome 5 (q-FT5) was detected in all experiments. Fine-mapping and expression studies identified a gene with orthology to a rice FLOWERING LOCUS T-like 9 (PhFTL9) as the candidate underlying q-FT5. We used a reciprocal transplant experiment to test for local adaptation and the specific impact of q-FT5 on performance. We did not observe local adaptation in terms of fitness tradeoffs when contrasting ecotypes in home versus away habitats. However, we observed that the coastal allele of q-FT5 conferred a fitness advantage only in its local habitat but not at the inland site. Sequence analyses identified an excess of low-frequency polymorphisms at the PhFTL9 promoter in the inland lineage, suggesting a role for either selection or population expansion on promoter evolution. Together, our findings demonstrate the genetic basis of flowering variation in a perennial grass and provide evidence for conditional neutrality underlying flowering time divergence.

Integrative systems analysis identifies genetic and dietary modulators of bile acid homeostasis

Bile acids (BAs) are complex and incompletely understood enterohepatic-derived hormones that control whole-body metabolism. Here, we profiled postprandial BAs in the liver, feces, and plasma of 360 chow- or high-fat-diet-fed BXD male mice and demonstrated that both genetics and diet strongly influence BA abundance, composition, and correlation with metabolic traits. Through an integrated systems approach, we mapped hundreds of quantitative trait loci that modulate BAs and identified both known and unknown regulators of BA homeostasis. In particular, we discovered carboxylesterase 1c (Ces1c) as a genetic determinant of plasma tauroursodeoxycholic acid (TUDCA), a BA species with established disease-preventing actions. The association between Ces1c and plasma TUDCA was validated using data from independent mouse cohorts and a Ces1c knockout mouse model. Collectively, our data are a unique resource to dissect the physiological importance of BAs as determinants of metabolic traits, as underscored by the identification of CES1C as a master regulator of plasma TUDCA levels.

Life-long dietary restrictions have negligible or damaging effects on late-life cognitive performance: A key role for genetics in outcomes

Several studies report that caloric restriction (CR) or intermittent fasting (IF) can improve cognition, while others report limited or no cognitive benefits. Here, we compare the effects of 20% CR, 40% CR, 1-day IF, and 2-day IF feeding paradigms to ad libitum controls on Y-maze working memory (WM) and contextual fear memory (CFM) in a large population of Diversity Outbred mice that model the genetic diversity of humans. While CR and IF interventions improve lifespan, we observed no enhancement of working memory or CFM in mice on these feeding paradigms, and report 40% CR to be damaging to recall of CFM. Using Quantitative Trait Loci mapping, we identified the gene Slc16a7 to be associated with CFM outcomes in aged mice on lifespan promoting feeding paradigms. Limited utility of dieting and fasting on memory in mice that recapitulate genetic diversity in the human population highlights the need for anti-aging therapeutics that promote cognitive function, with the neuronal monocarboxylate transporter MCT2 encoded by Slc16a7 highlighted as novel target.

Using the Collaborative Cross and Diversity Outbred Mice in Immunology

The Collaborative Cross (CC) and the Diversity Outbred (DO) stock mouse panels are the most powerful murine genetics tools available to the genetics community. Together, they combine the strength of inbred animal models with the diversity of outbred populations. Using the 63 CC strains or a panel of DO mice, each derived from the same 8 parental mouse strains, researchers can map genetic contributions to exceptionally complex immunological and infectious disease traits that would require far greater powering if performed by genome-wide association studies (GWAS) in human populations. These tools allow genes to be studied in heterozygous and homozygous states and provide a platform to study epistasis between interacting loci. Most importantly, once a quantitative phenotype is investigated and quantitative trait loci are identified, confirmatory genetic studies can be performed, which is often problematic using the GWAS approach. In addition, novel stable mouse models for immune phenotypes are often derived from studies utilizing the DO and CC mice that can serve as stronger model systems than existing ones in the field. The CC/DO systems have contributed to the fields of cancer immunology, autoimmunity, vaccinology, infectious disease, allergy, tissue rejection, and tolerance but have thus far been greatly underutilized. In this article, we present a recent review of the field and point out key areas of immunology that are ripe for further investigation and awaiting new CC/DO research projects. We also highlight some of the strong computational tools that have been developed for analyzing CC/DO genetic and phenotypic data. Additionally, we have formed a centralized community on the CyVerse infrastructure where immunogeneticists can utilize those software tools, collaborate with groups across the world, and expand the use of the CC and DO systems for investigating immunogenetic phenomena. © 2022 Wiley Periodicals LLC.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to severe acute respiratory syndrome coronavirus (SARS-CoV) disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6, that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2, and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species. IMPORTANCE Host genetic variation is an important determinant that predicts disease outcomes following infection. In the setting of highly pathogenic coronavirus infections genetic determinants underlying host susceptibility and mortality remain unclear. To elucidate the role of host genetic variation on sarbecovirus pathogenesis and disease outcomes, we utilized the Collaborative Cross (CC) mouse genetic reference population as a model to identify susceptibility alleles to SARS-CoV and SARS-CoV-2 infections. Our findings reveal that a multitrait loci found in chromosome 9 is an important regulator of sarbecovirus pathogenesis in mice. Within this locus, we identified and validated CCR9 and CXCR6 as important regulators of host disease outcomes. Specifically, both CCR9 and CXCR6 are protective against severe SARS-CoV, SARS-CoV-2, and SARS-related HKU3 virus disease in mice. This chromosome 9 multitrait locus may be important to help identify genes that regulate coronavirus disease outcomes in humans.

stuart: an R package for the curation of SNP genotypes from experimental crosses

Genetic mapping in 2-generation crosses requires genotyping, usually performed with single nucleotide polymorphism markers arrays which provide high-density genetic information. However, genetic analysis on raw genotypes can lead to spurious or unreliable results due to defective single nucleotide polymorphism assays or wrong genotype interpretation. Here, we introduce stuart, an open-source R package, which analyzes raw genotyping data to filter single nucleotide polymorphism markers based on informativeness, Mendelian inheritance pattern, and consistency with parental genotypes. The functions of this package provide a curation pipeline and formatting adequate for genetic analysis with the R/qtl package. stuart is available with detailed documentation from https://gitlab.pasteur.fr/mouselab/stuart/.

Genetic and transcriptome analyses reveal the candidate genes and pathways involved in the inactive shade-avoidance response enabling high-density planting of soybean

High-density planting is a major way to improve crop yields. However, shade-avoidance syndrome (SAS) is a major factor limiting increased planting density. First Green Revolution addressed grass lodging problem by using dwarf/semi-dwarf genes. However, it is not suitable for soybean, which bear seeds on stalk and whose seed yield depends on plant height. Hence, mining shade-tolerant germplasms and elucidating the underlying mechanism could provide meaningful resources and information for high-yield breeding. Here, we report a high-plant density-tolerant soybean cultivar, JiDou 17, which exhibited an inactive SAS (iSAS) phenotype under high-plant density or low-light conditions at the seedling stage. A quantitative trait locus (QTL) mapping analysis using a recombinant inbred line (RIL) population showed that this iSAS phenotype is related to a major QTL, named shade-avoidance response 1 (qSAR1), which was detected. The mapping region was narrowed by a haplotype analysis into a 554 kb interval harboring 44 genes, including 4 known to be key regulators of the SAS network and 4 with a variance response to low-light conditions between near isogenic line (NIL) stems. Via RNA-seq, we identified iSAS-specific genes based on one pair of near isogenic lines (NILs) and their parents. The iSAS-specific genes expressed in the stems were significantly enriched in the "proteasomal protein catabolic" process and the proteasome pathway, which were recently suggested to promote the shade-avoidance response by enhancing PIF7 stability. Most iSAS-specific proteasome-related genes were downregulated under low-light conditions. The expression of genes related to ABA, CK, and GA significantly varied between the low- and normal-light conditions. This finding is meaningful for the cloning of genes that harbor beneficial variation(s) conferring the iSAS phenotype fixed in domestication and breeding practice.

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

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

QTLViewer: an interactive webtool for genetic analysis in the Collaborative Cross and Diversity Outbred mouse populations

The Collaborative Cross and the Diversity Outbred mouse populations are related multiparental populations, derived from the same 8 isogenic founder strains. They carry >50 M known genetic variants, which makes them ideal tools for mapping genetic loci that regulate phenotypes, including physiological and molecular traits. Mapping quantitative trait loci requires statistical and computational training, which can present a barrier to access for some researchers. The QTLViewer software allows users to graphically explore Collaborative Cross and Diversity Outbred quantitative trait locus mapping and related analyses performed through the R/qtl2 package. Additionally, the QTLViewer website serves as a repository for published Collaborative Cross and Diversity Outbred studies, increasing the accessibility of these genetic resources to the broader scientific community.

Age and diet shape the genetic architecture of body weight in diversity outbred mice

Understanding how genetic variation shapes a complex trait relies on accurately quantifying both the additive genetic and genotype–environment interaction effects in an age-dependent manner. We used a linear mixed model to quantify diet-dependent genetic contributions to body weight measured through adulthood in diversity outbred female mice under five diets. We observed that heritability of body weight declined with age under all diets, except the 40% calorie restriction diet. We identified 14 loci with age-dependent associations and 19 loci with age- and diet-dependent associations, with many diet-dependent loci previously linked to neurological function and behavior in mice or humans. We found their allelic effects to be dynamic with respect to genomic background, age, and diet, identifying several loci where distinct alleles affect body weight at different ages. These results enable us to more fully understand and predict the effectiveness of dietary intervention on overall health throughout age in distinct genetic backgrounds.

Body weight is one trait influenced by genes, age and environmental factors. Both internal and external environmental pressures are known to affect genetic variation over time. However, it is largely unknown how all factors – including age – interact to shape metabolism and bodyweight.

Wright et al. set out to quantify the interactions between genes and diet in ageing mice and found that the effect of genetics on mouse body weight changes with age. In the experiments, Wright et al. weighed 960 female mice with diverse genetic backgrounds, starting at two months of age into adulthood. The animals were randomized to different diets at six months of age. Some mice had unlimited food access, others received 20% or 40% less calories than a typical mouse diet, and some fasted one or two days per week.

Variations in their genetic background explained about 80% of differences in mice’s weight, but the influence of genetics relative to non-genetic factors decreased as they aged. Mice on the 40% calorie restriction diet were an exception to this rule and genetics accounted for 80% of their weight throughout adulthood, likely due to reduced influence from diet and reduced interactions between diet and genes. Several genes involved in metabolism, neurological function, or behavior, were associated with mouse weight.

The experiments highlight the importance of considering interactions between genetics, environment, and age in determining complex traits like body weight. The results and the approaches used by Wright et al. may help other scientists learn more about how the genetic predisposition to disease changes with environmental stimuli and age.

Accounting for diet and age

The diet and age of mice can modulate how different genetic variants impact body weight, demonstrating the need to take context into account when performing genetic studies.

SoyMAGIC: An Unprecedented Platform for Genetic Studies and Breeding Activities in Soybean

Multi-Parent Advanced Generation Inter-Cross (MAGIC) populations are emerging genetic platforms for high-resolution and fine mapping of quantitative traits, such as agronomic and seed composition traits in soybean (Glycine max L.). We have established an eight-parent MAGIC population, comprising 721 recombinant inbred lines (RILs), through conical inter-mating of eight soybean lines. The parental lines were genetically diverse elite cultivars carrying different agronomic and seed composition characteristics, including amino acids and fatty acids, as well as oil and protein concentrations. This study aimed to introduce soybean MAGIC (SoyMAGIC) population as an unprecedented platform for genotypic and phenotypic investigation of agronomic and seed quality traits in soybean. The RILs were evaluated for important seed composition traits using replicated field trials during 2020 and 2021. To measure the seed composition traits, near-infrared reflectance (NIR) was employed. The RILs were genotyped using genotyping-by-sequencing (GBS) method to decipher the genome and discover single-nucleotide polymorphic (SNP) markers among the RILs. A high-density linkage map was constructed through inclusive composite interval mapping (ICIM). The linkage map was 3,770.75 cM in length and contained 12,007 SNP markers. Chromosomes 11 and 18 were recorded as the shortest and longest linkage groups with 71.01 and 341.15 cM in length, respectively. Observed transgressive segregation of the selected traits and higher recombination frequency across the genome confirmed the capability of MAGIC population in reshuffling the diversity in the soybean genome among the RILs. The assessment of haplotype blocks indicated an uneven distribution of the parents' genomes in RILs, suggesting cryptic influence against or in favor of certain parental genomes. The SoyMAGIC population is a recombined genetic material that will accelerate further genomic studies and the development of soybean cultivars with improved seed quality traits through the development and implementation of reliable molecular-based toolkits.

Disentangling the genetic basis of rhizosphere microbiome assembly in tomato

Microbiomes play a pivotal role in plant growth and health, but the genetic factors involved in microbiome assembly remain largely elusive. Here, we map the molecular features of the rhizosphere microbiome as quantitative traits of a diverse hybrid population of wild and domesticated tomato. Gene content analysis of prioritized tomato quantitative trait loci suggests a genetic basis for differential recruitment of various rhizobacterial lineages, including a Streptomyces-associated 6.31 Mbp region harboring tomato domestication sweeps and encoding, among others, the iron regulator FIT and the water channel aquaporin SlTIP2.3. Within metagenome-assembled genomes of root-associated Streptomyces and Cellvibrio, we identify bacterial genes involved in metabolism of plant polysaccharides, iron, sulfur, trehalose, and vitamins, whose genetic variation associates with specific tomato QTLs. By integrating 'microbiomics' and quantitative plant genetics, we pinpoint putative plant and reciprocal rhizobacterial traits underlying microbiome assembly, thereby providing a first step towards plant-microbiome breeding programs.

A Multitrait Locus Regulates Sarbecovirus Pathogenesis

Infectious diseases have shaped the human population genetic structure, and genetic variation influences the susceptibility to many viral diseases. However, a variety of challenges have made the implementation of traditional human Genome-wide Association Studies (GWAS) approaches to study these infectious outcomes challenging. In contrast, mouse models of infectious diseases provide an experimental control and precision, which facilitates analyses and mechanistic studies of the role of genetic variation on infection. Here we use a genetic mapping cross between two distinct Collaborative Cross mouse strains with respect to SARS-CoV disease outcomes. We find several loci control differential disease outcome for a variety of traits in the context of SARS-CoV infection. Importantly, we identify a locus on mouse Chromosome 9 that shows conserved synteny with a human GWAS locus for SARS-CoV-2 severe disease. We follow-up and confirm a role for this locus, and identify two candidate genes, CCR9 and CXCR6 that both play a key role in regulating the severity of SARS-CoV, SARS-CoV-2 and a distantly related bat sarbecovirus disease outcomes. As such we provide a template for using experimental mouse crosses to identify and characterize multitrait loci that regulate pathogenic infectious outcomes across species.

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

BACKGROUND

A strong predictor for the development of alcohol use disorder (AUD) is altered sensitivity to the intoxicating effects of alcohol. Individual differences in the initial sensitivity to alcohol are controlled in part by genetic factors. Mice offer a powerful tool to elucidate the genetic basis of behavioral and physiological traits relevant to AUD, but conventional experimental crosses have only been able to identify large chromosomal regions rather than specific genes. Genetically diverse, highly recombinant mouse populations make it possible to observe a wider range of phenotypic variation, offer greater mapping precision, and thus increase the potential for efficient gene identification.

METHODS

We have taken advantage of the Diversity Outbred (DO) mouse population to identify and precisely map quantitative trait loci (QTL) associated with ethanol sensitivity. We phenotyped 798 male J:DO mice for three measures of ethanol sensitivity: ataxia, hypothermia, and loss of the righting response. We used high-density MegaMUGA and GigaMUGA to obtain genotypes ranging from 77,808 to 143,259 SNPs. We also performed RNA sequencing in striatum to map expression QTLs and identify gene expression-trait correlations. We then applied a systems genetic strategy to identify narrow QTLs and construct the network of correlations that exists between DNA sequence, gene expression values, and ethanol-related phenotypes to prioritize our list of positional candidate genes.

RESULTS

We observed large amounts of phenotypic variation with the DO population and identified suggestive and significant QTLs associated with ethanol sensitivity on chromosomes 1, 2, and 16. The implicated regions were narrow (4.5-6.9 Mb in size) and each QTL explained ~4-5% of the variance.

CONCLUSIONS

Our results can be used to identify alleles that contribute to AUD in humans, elucidate causative biological mechanisms, or assist in the development of novel therapeutic interventions.

The genetic basis of structural colour variation in mimetic Heliconius butterflies

Structural colours, produced by the reflection of light from ultrastructures, have evolved multiple times in butterflies. Unlike pigmentary colours and patterns, little is known about the genetic basis of these colours. Reflective structures on wing-scale ridges are responsible for iridescent structural colour in many butterflies, including the Müllerian mimics Heliconius erato and Heliconius melpomene. Here, we quantify aspects of scale ultrastructure variation and colour in crosses between iridescent and non-iridescent subspecies of both of these species and perform quantitative trait locus (QTL) mapping. We show that iridescent structural colour has a complex genetic basis in both species, with offspring from crosses having a wide variation in blue colour (both hue and brightness) and scale structure measurements. We detect two different genomic regions in each species that explain modest amounts of this variation, with a sex-linked QTL in H. erato but not H. melpomene. We also find differences between species in the relationships between structure and colour, overall suggesting that these species have followed different evolutionary trajectories in their evolution of structural colour. We then identify genes within the QTL intervals that are differentially expressed between subspecies and/or wing regions, revealing likely candidates for genes controlling structural colour formation.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

Using breeding and quantitative genetics to understand the C4 pathway

Reducing photorespiration in C3 crops could significantly increase rates of photosynthesis and yield. One method to achieve this would be to integrate C4 photosynthesis into C3 species. This objective is challenging as it involves engineering incompletely understood traits into C3 leaves, including complex changes to their biochemistry, cell biology, and anatomy. Quantitative genetics and selective breeding offer underexplored routes to identify regulators of these processes. We first review examples of natural intraspecific variation in C4 photosynthesis as well as the potential for hybridization between C3 and C4 species. We then discuss how quantitative genetic approaches including artificial selection and genome-wide association could be used to better understand the C4 syndrome and in so doing guide the engineering of the C4 pathway into C3 crops.

Dissecting the Genetic Architecture of Carbon Partitioning in Sorghum Using Multiscale Phenotypes

Carbon partitioning in plants may be viewed as a dynamic process composed of the many interactions between sources and sinks. The accumulation and distribution of fixed carbon is not dictated simply by the sink strength and number but is dependent upon the source, pathways, and interactions of the system. As such, the study of carbon partitioning through perturbations to the system or through focus on individual traits may fail to produce actionable developments or a comprehensive understanding of the mechanisms underlying this complex process. Using the recently published sorghum carbon-partitioning panel, we collected both macroscale phenotypic characteristics such as plant height, above-ground biomass, and dry weight along with microscale compositional traits to deconvolute the carbon-partitioning pathways in this multipurpose crop. Multivariate analyses of traits resulted in the identification of numerous loci associated with several distinct carbon-partitioning traits, which putatively regulate sugar content, manganese homeostasis, and nitrate transportation. Using a multivariate adaptive shrinkage approach, we identified several loci associated with multiple traits suggesting that pleiotropic and/or interactive effects may positively influence multiple carbon-partitioning traits, or these overlaps may represent molecular switches mediating basal carbon allocating or partitioning networks. Conversely, we also identify a carbon tradeoff where reduced lignin content is associated with increased sugar content. The results presented here support previous studies demonstrating the convoluted nature of carbon partitioning in sorghum and emphasize the importance of taking a holistic approach to the study of carbon partitioning by utilizing multiscale phenotypes.

Further evidence from common garden rearing experiments of heritable traits separating lean and siscowet lake charr (Salvelinus namaycush) ecotypes

Genetic evidence of selection for complex and polygenically regulated phenotypes can easily become masked by neutral population genetic structure and phenotypic plasticity. Without direct evidence of genotype‐phenotype associations it can be difficult to conclude to what degree a phenotype is heritable or a product of environment. Common garden laboratory studies control for environmental stochasticity and help to determine the mechanism that regulate traits. Here we assess lipid content, growth, weight, and length variation in full and hybrid F₁ crosses of deep and shallow water sympatric lake charr ecotypes reared for nine years in a common garden experiment. Redundancy analysis (RDA) and quantitative‐trait‐loci (QTL) genomic scans are used to identify associations between genotypes at 19,714 single nucleotide polymorphisms (SNPs) aligned to the lake charr genome and individual phenotypes to determine the role that genetic inheritance plays in ecotype phenotypic diversity. Lipid content, growth, length, and weight differed significantly among lake charr crosses throughout the experiment suggesting that pedigree plays a large roll in lake charr development. Polygenic scores of 15 SNPs putatively associated with lipid content and/or condition factor indicated that ecotype distinguishing traits are polygenically regulated and additive. A QTL identified on chromosome 38 contained >200 genes, some of which were associated with lipid metabolism and growth, demonstrating the complex nature of ecotype diversity. The results of our common garden study further indicate that lake charr ecotypes observed in nature are predetermined at birth and that ecotypes differ fundamentally in lipid metabolism and growth.

Genome-wide transcript and protein analysis highlights the role of protein homeostasis in the aging mouse heart

Investigation of the molecular mechanisms of aging in the human heart is challenging because of confounding factors, such as diet and medications, as well as limited access to tissues from healthy aging individuals. The laboratory mouse provides an ideal model to study aging in healthy individuals in a controlled environment. However, previous mouse studies have examined only a narrow range of the genetic variation that shapes individual differences during aging. Here, we analyze transcriptome and proteome data from 185 genetically diverse male and female mice at ages 6, 12, and 18 mo to characterize molecular changes that occur in the aging heart. Transcripts and proteins reveal activation of pathways related to exocytosis and cellular transport with age, whereas processes involved in protein folding decrease with age. Additional changes are apparent only in the protein data including reduced fatty acid oxidation and increased autophagy. For proteins that form complexes, we see a decline in correlation between their component subunits with age, suggesting age-related loss of stoichiometry. The most affected complexes are themselves involved in protein homeostasis, which potentially contributes to a cycle of progressive breakdown in protein quality control with age. Our findings highlight the important role of post-transcriptional regulation in aging. In addition, we identify genetic loci that modulate age-related changes in protein homeostasis, suggesting that genetic variation can alter the molecular aging process.

A Bayesian model selection approach to mediation analysis

Genetic studies often seek to establish a causal chain of events originating from genetic variation through to molecular and clinical phenotypes. When multiple phenotypes share a common genetic association, one phenotype may act as an intermediate for the genetic effects on the other. Alternatively, the phenotypes may be causally unrelated but share genetic loci. Mediation analysis represents a class of causal inference approaches used to determine which of these scenarios is most plausible. We have developed a general approach to mediation analysis based on Bayesian model selection and have implemented it in an R package, bmediatR. Bayesian model selection provides a flexible framework that can be tailored to different analyses. Our approach can incorporate prior information about the likelihood of models and the strength of causal effects. It can also accommodate multiple genetic variants or multi-state haplotypes. Our approach reports posterior probabilities that can be useful in interpreting uncertainty among competing models. We compared bmediatR with other popular methods, including the Sobel test, Mendelian randomization, and Bayesian network analysis using simulated data. We found that bmediatR performed as well or better than these alternatives in most scenarios. We applied bmediatR to proteome data from Diversity Outbred (DO) mice, a multi-parent population, and demonstrate the power of mediation with multi-state haplotypes. We also applied bmediatR to data from human cell lines to identify transcripts that are mediated through or are expressed independently from local chromatin accessibility. We demonstrate that Bayesian model selection provides a powerful and versatile approach to identify causal relationships in genetic studies using model organism or human data.

A natural mutator allele shapes mutation spectrum variation in mice

Although germline mutation rates and spectra can vary within and between species, common genetic modifiers of the mutation rate have not been identified in nature1. Here we searched for loci that influence germline mutagenesis using a uniquely powerful resource: a panel of recombinant inbred mouse lines known as the BXD, descended from the laboratory strains C57BL/6J (B haplotype) and DBA/2J (D haplotype). Each BXD lineage has been maintained by brother-sister mating in the near absence of natural selection, accumulating de novo mutations for up to 50 years on a known genetic background that is a unique linear mosaic of B and D haplotypes2. We show that mice inheriting D haplotypes at a quantitative trait locus on chromosome 4 accumulate C>A germline mutations at a 50% higher rate than those inheriting B haplotypes, primarily owing to the activity of a C>A-dominated mutational signature known as SBS18. The B and D quantitative trait locus haplotypes encode different alleles of Mutyh, a DNA repair gene that underlies the heritable cancer predisposition syndrome that causes colorectal tumors with a high SBS18 mutation load3,4. Both B and D Mutyh alleles are present in wild populations of Mus musculus domesticus, providing evidence that common genetic variation modulates germline mutagenesis in a model mammalian species.

A diversity outbred F1 mouse model identifies host-intrinsic genetic regulators of response to immune checkpoint inhibitors

Immune checkpoint inhibitors (ICI) have improved outcomes for a variety of malignancies; however, many patients fail to benefit. While tumor-intrinsic mechanisms are likely involved in therapy resistance, it is unclear to what extent host genetic background influences response. To investigate this, we utilized the Diversity Outbred (DO) and Collaborative Cross (CC) mouse models. DO mice are an outbred stock generated by crossbreeding eight inbred founder strains, and CC mice are recombinant inbred mice generated from the same eight founders. We generated 207 DOB6F1 mice representing 48 DO dams and demonstrated that these mice reliably accept the C57BL/6-syngeneic B16F0 tumor and that host genetic background influences response to ICI. Genetic linkage analysis from 142 mice identified multiple regions including one within chromosome 13 that associated with therapeutic response. We utilized 6 CC strains bearing the positive (NZO) or negative (C57BL/6) driver genotype in this locus. We found that 2/3 of predicted responder CCB6F1 crosses show reproducible ICI response. The chromosome 13 locus contains the murine prolactin family, which is a known immunomodulating cytokine associated with various autoimmune disorders. To directly test whether prolactin influences ICI response rates, we implanted inbred C57BL/6 mice with subcutaneous slow-release prolactin pellets to induce mild hyperprolactinemia. Prolactin augmented ICI response against B16F0, with increased CD8 infiltration and 5/8 mice exhibiting slowed tumor growth relative to controls. This study highlights the role of host genetics in ICI response and supports the use of F1 crosses in the DO and CC mouse populations as powerful cancer immunotherapy models.

Genetic mapping and prediction of flowering time and plant height in a maize Stiff Stalk MAGIC population

The Stiff Stalk heterotic pool is a foundation of US maize seed parent germplasm and has been heavily utilized by both public and private maize breeders since its inception in the 1930's. Flowering time and plant height are critical characteristics for both inbred parents and their test crossed hybrid progeny. To study these traits, a six parent multiparent advanced generation intercross (MAGIC) population was developed including maize inbred lines B73, B84, PHB47 (B37 type), LH145 (B14 type), PHJ40 (novel early Stiff Stalk), and NKH8431 (B73/B14 type). A set of 779 doubled haploid lines were evaluated for flowering time and plant height in two field replicates in 2016 and 2017, and a subset of 689 and 561 doubled haploid lines were crossed to two testers, respectively, and evaluated as hybrids in two locations in 2018 and 2019 using an incomplete block design. Markers were derived from a Practical Haplotype Graph built from the founder whole genome assemblies and genotype-by-sequencing and exome capture-based sequencing of the population. Genetic mapping utilizing an update to R/qtl2 revealed differing profiles of significant loci for both traits between 635 of the DH lines and two sets of 570 and 471 derived hybrids. Genomic prediction was used to test the feasibility of predicting hybrid phenotypes based on the per se data. Predictive abilities were highest on direct models trained using the data they would predict (0.55 to 0.63), and indirect models trained using per se data to predict hybrid traits had slightly lower predictive abilities (0.49 to 0.55). Overall, this finding is consistent with the overlapping and non-overlapping significant QTL found within the per se and hybrid populations and suggests that selections for phenology traits can be made effectively on doubled haploid lines before hybrid data is available.

Automated, high-dimensional evaluation of physiological aging and resilience in outbred mice

Behavior and physiology are essential readouts in many studies but have not benefited from the high-dimensional data revolution that has transformed molecular and cellular phenotyping. To address this, we developed an approach that combines commercially available automated phenotyping hardware with a systems biology analysis pipeline to generate a high-dimensional readout of mouse behavior/physiology, as well as intuitive and health-relevant summary statistics (resilience and biological age). We used this platform to longitudinally evaluate aging in hundreds of outbred mice across an age range from 3 months to 3.4 years. In contrast to the assumption that aging can only be measured at the limits of animal ability via challenge-based tasks, we observed widespread physiological and behavioral aging starting in early life. Using network connectivity analysis, we found that organism-level resilience exhibited an accelerating decline with age that was distinct from the trajectory of individual phenotypes. We developed a method, Combined Aging and Survival Prediction of Aging Rate (CASPAR), for jointly predicting chronological age and survival time and showed that the resulting model is able to predict both variables simultaneously, a behavior that is not captured by separate age and mortality prediction models. This study provides a uniquely high-resolution view of physiological aging in mice and demonstrates that systems-level analysis of physiology provides insights not captured by individual phenotypes. The approach described here allows aging, and other processes that affect behavior and physiology, to be studied with improved throughput, resolution, and phenotypic scope.

Prediction performance of linear models and gradient boosting machine on complex phenotypes in outbred mice

We compared the performance of linear (GBLUP, BayesB, and elastic net) methods to a nonparametric tree-based ensemble (gradient boosting machine) method for genomic prediction of complex traits in mice. The dataset used contained genotypes for 50,112 SNP markers and phenotypes for 835 animals from 6 generations. Traits analyzed were bone mineral density, body weight at 10, 15, and 20 weeks, fat percentage, circulating cholesterol, glucose, insulin, triglycerides, and urine creatinine. The youngest generation was used as a validation subset, and predictions were based on all older generations. Model performance was evaluated by comparing predictions for animals in the validation subset against their adjusted phenotypes. Linear models outperformed gradient boosting machine for 7 out of 10 traits. For bone mineral density, cholesterol, and glucose, the gradient boosting machine model showed better prediction accuracy and lower relative root mean squared error than the linear models. Interestingly, for these 3 traits, there is evidence of a relevant portion of phenotypic variance being explained by epistatic effects. Using a subset of top markers selected from a gradient boosting machine model helped for some of the traits to improve the accuracy of prediction when these were fitted into linear and gradient boosting machine models. Our results indicate that gradient boosting machine is more strongly affected by data size and decreased connectedness between reference and validation sets than the linear models. Although the linear models outperformed gradient boosting machine for the polygenic traits, our results suggest that gradient boosting machine is a competitive method to predict complex traits with assumed epistatic effects.

Linkage mapping, comparative genome analysis, and QTL detection for growth in a non-model teleost, the meagre Argyrosomus regius, using ddRAD sequencing

Meagre (Argyrosomus regius), is a benthopelagic species rapidly emerging in aquaculture, due to its low food to biomass conversion rate, good fillet yield and ease of production. Tracing a species genomic background along with describing the genetic basis of important traits can greatly influence both conservation strategies and production perspectives. In this study, we employed ddRAD sequencing of 266 fish from six F1 meagre families, to construct a high-density genetic map comprising 4529 polymorphic SNP markers. The QTL mapping analysis provided a genomic appreciation for the weight trait identifying a statistically significant QTL on linkage group 15 (LG15). The comparative genomics analysis with six teleost species revealed an evolutionarily conserved karyotype structure. The synteny observed, verified the already well-known fusion events of the three-spine stickleback genome, reinforced the evidence of reduced evolutionary distance of Sciaenids with the Sparidae family, reflected the evolutionary proximity with Dicentrarchus labrax, traced several putative chromosomal rearrangements and a prominent putative fusion event in meagre’s LG17. This study presents novel elements concerning the genome evolutionary history of a non-model teleost species recently adopted in aquaculture, starts to unravel the genetic basis of the species growth-related traits, and provides a high-density genetic map as a tool that can help to further establish meagre as a valuable resource for research and production.

New Insights on Gene by Environmental Effects of Drugs of Abuse in Animal Models Using GeneNetwork

Gene-by-environment interactions are important for all facets of biology, especially behaviour. Families of isogenic strains of mice, such as the BXD strains, are excellently placed to study these interactions, as the same genome can be tested in multiple environments. BXD strains are recombinant inbred mouse strains derived from crossing two inbred strains—C57BL/6J and DBA/2J mice. Many reproducible genometypes can be leveraged, and old data can be reanalysed with new tools to produce novel insights. We obtained drug and behavioural phenotypes from Philip et al. Genes, Brain and Behaviour 2010, and reanalysed their data with new genotypes from sequencing, as well as new models (Genome-wide Efficient Mixed Model Association (GEMMA) and R/qtl2). We discovered QTLs on chromosomes 3, 5, 9, 11, and 14, not found in the original study. We reduced the candidate genes based on their ability to alter gene expression or protein function. Candidate genes included Slitrk6 and Cdk14. Slitrk6, in a Chromosome14 QTL for locomotion, was found to be part of a co-expression network involved in voluntary movement and associated with neuropsychiatric phenotypes. Cdk14, one of only three genes in a Chromosome5 QTL, is associated with handling induced convulsions after ethanol treatment, that is regulated by the anticonvulsant drug valproic acid. By using families of isogenic strains, we can reanalyse data to discover novel candidate genes involved in response to drugs of abuse.