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

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.

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.

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.

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.

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.

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.

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.

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.

QTL mapping and candidate gene analysis reveal two major loci regulating green leaf color in non-heading Chinese cabbage

KEY MESSAGE

Two major-effect QTL GlcA07.1 and GlcA09.1 for green leaf color were fine mapped into 170.25 kb and 191.41 kb intervals on chromosomes A07 and A09, respectively, and were validated by transcriptome analysis. Non-heading Chinese cabbage (NHCC) is a leafy vegetable with a wide range of green colors. Understanding the genetic mechanism behind broad spectrum of green may facilitate the breeding of high-quality NHCC. Here, we used F2 and F7:8 recombination inbred line (RIL) population from a cross between Wutacai (dark-green) and Erqing (lime-green) to undertake the genetic analysis and quantitative trait locus (QTL) mapping in NHCC. The genetic investigation of the F2 population revealed that the variation of green leaf color was controlled by two recessive genes. Six pigments associated with green leaf color, including total chlorophyll, chlorophyll a, chlorophyll b, total carotenoids, lutein, and carotene were quantified and applied for QTL mapping in the RIL population. A total of 7 QTL were detected across the whole genome. Among them, two major-effect QTL were mapped on chromosomes A07 (GlcA07.1) and A09 (GlcA09.1) corresponding to two QTL identified in the F2 population. The QTL GlcA07.1 and GlcA09.1 were further fine mapped into 170.25 kb and 191.41 kb genomic regions, respectively. By comparing gene expression level and gene annotation, BraC07g023810 and BraC07g023970 were proposed as the best candidates for GlcA07.1, while BraC09g052220 and BraC09g052270 were suggested for GlcA09.1. Two InDel molecular markers (GlcA07.1-BcGUN4 and GlcA09.1-BcSG1) associated with BraC07gA023810 and BraC09g052220 were developed and could effectively identify leaf color in natural NHCC accessions, suggesting their potential for marker-assisted leaf color selection in NHCC breeding.

Gene expression variation underlying tissue-specific responses to copper stress in Drosophila melanogaster

Copper is one of a handful of biologically necessary heavy metals that is also a common environmental pollutant. Under normal conditions, copper ions are required for many key physiological processes. However, in excess, copper results in cell and tissue damage ranging in severity from temporary injury to permanent neurological damage. Because of its biological relevance, and because many conserved copper-responsive genes respond to nonessential heavy metal pollutants, copper resistance in Drosophila melanogaster is a useful model system with which to investigate the genetic control of the heavy metal stress response. Because heavy metal toxicity has the potential to differently impact specific tissues, we genetically characterized the control of the gene expression response to copper stress in a tissue-specific manner in this study. We assessed the copper stress response in head and gut tissue of 96 inbred strains from the Drosophila Synthetic Population Resource using a combination of differential expression analysis and expression quantitative trait locus mapping. Differential expression analysis revealed clear patterns of tissue-specific expression. Tissue and treatment specific responses to copper stress were also detected using expression quantitative trait locus mapping. Expression quantitative trait locus associated with MtnA, Mdr49, Mdr50, and Sod3 exhibited both genotype-by-tissue and genotype-by-treatment effects on gene expression under copper stress, illuminating tissue- and treatment-specific patterns of gene expression control. Together, our data build a nuanced description of the roles and interactions between allelic and expression variation in copper-responsive genes, provide valuable insight into the genomic architecture of susceptibility to metal toxicity, and highlight candidate genes for future functional characterization.

The regulation of methylation on the Z chromosome and the identification of multiple novel Male Hyper-Methylated regions in the chicken

DNA methylation is a key regulator of eukaryote genomes, and is of particular relevance in the regulation of gene expression on the sex chromosomes, with a key role in dosage compensation in mammalian XY systems. In the case of birds, dosage compensation is largely absent, with it being restricted to two small Male Hyper-Methylated (MHM) regions on the Z chromosome. To investigate how variation in DNA methylation is regulated on the Z chromosome we utilised a wild x domestic advanced intercross in the chicken, with both hypothalamic methylomes and transcriptomes assayed in 124 individuals. The relatively large numbers of individuals allowed us to identify additional genomic MHM regions on the Z chromosome that were significantly differentially methylated between the sexes. These regions appear to down-regulate local gene expression in males, but not remove it entirely (unlike the lncRNAs identified in the initial MHM regions). These MHM regions were further tested and the most balanced genes appear to show decreased expression in males, whilst methylation appeared to be far more correlated with gene expression in the less balanced, as compared to the most balanced genes. In addition, quantitative trait loci (QTL) that regulate variation in methylation on the Z chromosome, and those loci that regulate methylation on the autosomes that derive from the Z chromosome were mapped. Trans-effect hotspots were also identified that were based on the autosomes but affected the Z, and also one that was based on the Z chromosome but that affected both autosomal and sex chromosome DNA methylation regulation. We show that both cis and trans loci that originate from the Z chromosome never exhibit an interaction with sex, whereas trans loci originating from the autosomes but affecting the Z chromosome always display such an interaction. Our results highlight how additional MHM regions are actually present on the Z chromosome, and they appear to have smaller-scale effects on gene expression in males. Quantitative variation in methylation is also regulated both from the autosomes to the Z chromosome, and from the Z chromosome to the autosomes.

Genetic dissection of domestication traits in interspecific chickpea populations

Chickpea (Cicer arietinum) is a pulse crop that provides an integral source of nutrition for human consumption. The close wild relatives Cicer reticulatum and Cicer echinospermum harbor untapped genetic diversity that can be exploited by chickpea breeders to improve domestic varieties. Knowledge of genomic loci that control important chickpea domestication traits will expedite the development of improved chickpea varieties derived from interspecific crosses. Therefore, we set out to identify genomic loci underlying key chickpea domestication traits by both association and quantitative trait locus (QTL) mapping using interspecific F2 populations. Diverse phenotypes were recorded for various agronomic traits. A total of 11 high-confidence markers were detected on chromosomes 1, 3, and 7 by both association and QTL mapping; these were associated with growth habit, flowering time, and seed traits. Furthermore, we identified candidate genes linked to these markers, which advanced our understanding of the genetic basis of domestication traits and validated known genes such as the FLOWERING LOCUS gene cluster that regulates flowering time. Collectively, this study has elucidated the genetic basis of chickpea domestication traits, which can facilitate the development of superior chickpea varieties.

Multiple genetic loci influence vaccine-induced protection against Mycobacterium tuberculosis in genetically diverse mice

Mycobacterium tuberculosis (M.tb.) infection leads to over 1.5 million deaths annually, despite widespread vaccination with BCG at birth. Causes for the ongoing tuberculosis endemic are complex and include the failure of BCG to protect many against progressive pulmonary disease. Host genetics is one of the known factors implicated in susceptibility to primary tuberculosis, but less is known about the role that host genetics plays in controlling host responses to vaccination against M.tb. Here, we addressed this gap by utilizing Diversity Outbred (DO) mice as a small animal model to query genetic drivers of vaccine-induced protection against M.tb. DO mice are a highly genetically and phenotypically diverse outbred population that is well suited for fine genetic mapping. Similar to outcomes in people, our previous studies demonstrated that DO mice have a wide range of disease outcomes following BCG vaccination and M.tb. challenge. In the current study, we used a large population of BCG-vaccinated/M.tb.-challenged mice to perform quantitative trait loci mapping of complex infection traits; these included lung and spleen M.tb. burdens, as well as lung cytokines measured at necropsy. We found sixteen chromosomal loci associated with complex infection traits and cytokine production. QTL associated with bacterial burdens included a region encoding major histocompatibility antigens that are known to affect susceptibility to tuberculosis, supporting validity of the approach. Most of the other QTL represent novel associations with immune responses to M.tb. and novel pathways of cytokine regulation. Most importantly, we discovered that protection induced by BCG is a multigenic trait, in which genetic loci harboring functionally-distinct candidate genes influence different aspects of immune responses that are crucial collectively for successful protection. These data provide exciting new avenues to explore and exploit in developing new vaccines against M.tb.

Genetic mapping identifies Homer1 as a developmental modifier of attention

Attention is required for most higher-order cognitive functions. Prior studies have revealed functional roles for the prefrontal cortex and its extended circuits to enabling attention, but the underlying molecular processes and their impacts on cellular and circuit function remain poorly understood. To develop insights, we here took an unbiased forward genetics approach to identify single genes of large effect on attention. We studied 200 genetically diverse mice on measures of pre-attentive processing and through genetic mapping identified a small locus on chromosome 13 (95%CI: 92.22-94.09 Mb) driving substantial variation (19%) in this trait. Further characterization of the locus revealed a causative gene, Homer1, encoding a synaptic protein, where down-regulation of its short isoforms in prefrontal cortex (PFC) during early postnatal development led to improvements in multiple measures of attention in the adult. Subsequent mechanistic studies revealed that prefrontal Homer1 down-regulation is associated with GABAergic receptor up-regulation in those same cells. This enhanced inhibitory influence, together with dynamic neuromodulatory coupling, led to strikingly low PFC activity at baseline periods of the task but targeted elevations at cue onset, predicting short-latency correct choices. Notably high-Homer1, low-attentional performers, exhibited uniformly elevated PFC activity throughout the task. We thus identify a single gene of large effect on attention - Homer1 - and find that it improves prefrontal inhibitory tone and signal-to-noise (SNR) to enhance attentional performance. A therapeutic strategy focused on reducing prefrontal activity and increasing SNR, rather than uniformly elevating PFC activity, may complement the use of stimulants to improve attention.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

The type 1 diabetes susceptibility locus Idd5 favours robust neonatal development of highly autoreactive regulatory T cells in the NOD mouse

Regulatory T lymphocytes expressing the transcription factor Foxp3 (Tregs) play an important role in the prevention of autoimmune diseases and other immunopathologies. Aberrations in Treg-mediated immunosuppression are therefore thought to be involved in the development of autoimmune pathologies, but few have been documented. Recent reports indicated a central role for Tregs developing during the neonatal period in the prevention of autoimmune pathology. We therefore investigated the development of Tregs in neonatal NOD mice, an important animal model for autoimmune type 1 diabetes. Surprisingly, we found that, as compared with seven other commonly studied inbred mouse strains, in neonatal NOD mice, exceptionally large proportions of developing Tregs express high levels of GITR and PD-1. The latter phenotype was previously associated with high Treg autoreactivity in C57BL/6 mice, which we here confirm for NOD animals. The proportions of newly developing GITRhighPD-1+ Tregs rapidly drop during the first week of age. A genome-wide genetic screen indicated the involvement of several diabetes susceptibility loci in this trait. Analysis of a congenic mouse strain confirmed that Idd5 contributes to the genetic control of GITRhighPD-1+ Treg development in neonates. Our data thus demonstrate an intriguing and paradoxical correlation between an idiosyncrasy in Treg development in NOD mice and their susceptibility to type 1 diabetes.

Dissection of the E8 locus in two early maturing Canadian soybean populations

Soybean [Glycine max (L.) Merr.] is a short-day crop for which breeders want to expand the cultivation range to more northern agro-environments by introgressing alleles involved in early reproductive traits. To do so, we investigated quantitative trait loci (QTL) and expression quantitative trait loci (eQTL) regions comprised within the E8 locus, a large undeciphered region (~7.0 Mbp to 44.5 Mbp) associated with early maturity located on chromosome GM04. We used a combination of two mapping algorithms, (i) inclusive composite interval mapping (ICIM) and (ii) genome-wide composite interval mapping (GCIM), to identify major and minor regions in two soybean populations (QS15524F2:F3 and QS15544RIL) having fixed E1, E2, E3, and E4 alleles. Using this approach, we identified three main QTL regions with high logarithm of the odds (LODs), phenotypic variation explained (PVE), and additive effects for maturity and pod-filling within the E8 region: GM04:16,974,874-17,152,230 (E8-r1); GM04:35,168,111-37,664,017 (E8-r2); and GM04:41,808,599-42,376,237 (E8-r3). Using a five-step variant analysis pipeline, we identified Protein far-red elongated hypocotyl 3 (Glyma.04G124300; E8-r1), E1-like-a (Glyma.04G156400; E8-r2), Light-harvesting chlorophyll-protein complex I subunit A4 (Glyma.04G167900; E8-r3), and Cycling dof factor 3 (Glyma.04G168300; E8-r3) as the most promising candidate genes for these regions. A combinatorial eQTL mapping approach identified significant regulatory interactions for 13 expression traits (e-traits), including Glyma.04G050200 (Early flowering 3/E6 locus), with the E8-r3 region. Four other important QTL regions close to or encompassing major flowering genes were also detected on chromosomes GM07, GM08, and GM16. In GM07:5,256,305-5,404,971, a missense polymorphism was detected in the candidate gene Glyma.07G058200 (Protein suppressor of PHYA-105). These findings demonstrate that the locus known as E8 is regulated by at least three distinct genomic regions, all of which comprise major flowering genes.

The chromosome-scale genome and the genetic resistance machinery against insect herbivores of the Mexican toloache, Datura stramonium

Plant resistance refers to the heritable ability of plants to reduce damage caused by natural enemies, such as herbivores and pathogens, either through constitutive or induced traits like chemical compounds or trichomes. However, the genetic architecture-the number and genome location of genes that affect plant defense and the magnitude of their effects-of plant resistance to arthropod herbivores in natural populations remains poorly understood. In this study, we aimed to unveil the genetic architecture of plant resistance to insect herbivores in the annual herb Datura stramonium (Solanaceae) through quantitative trait loci mapping. We achieved this by assembling the species' genome and constructing a linkage map using an F2 progeny transplanted into natural habitats. Furthermore, we conducted differential gene expression analysis between undamaged and damaged plants caused by the primary folivore, Lema daturaphila larvae. Our genome assembly resulted in 6,109 scaffolds distributed across 12 haploid chromosomes. A single quantitative trait loci region on chromosome 3 was associated with plant resistance, spanning 0 to 5.17 cM. The explained variance by the quantitative trait loci was 8.44%. Our findings imply that the resistance mechanisms of D. stramonium are shaped by the complex interplay of multiple genes with minor effects. Protein-protein interaction networks involving genes within the quantitative trait loci region and overexpressed genes uncovered the key role of receptor-like cytoplasmic kinases in signaling and regulating tropane alkaloids and terpenoids, which serve as powerful chemical defenses against D. stramonium herbivores. The data generated in our study constitute important resources for delving into the evolution and ecology of secondary compounds mediating plant-insect interactions.

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Methods to study xenografted human cancer in genetically diverse mice

Xenografting human cancer tissues into mice to test new cures against cancers is critical for understanding and treating the disease. However, only a few inbred strains of mice are used to study cancers, and derivatives of mainly one strain, mostly NOD/ShiLtJ, are used for therapy efficacy studies. As it has been demonstrated when human cancer cell lines or patient-derived tissues (PDX) are xenografted into mice, the neoplastic cells are human but the supporting cells that comprise the tumor (the stroma) are from the mouse. Therefore, results of studies of xenografted tissues are influenced by the host strain. We previously published that when the same neoplastic cells are xenografted into different mouse strains, the pattern of tumor growth, histology of the tumor, number of immune cells infiltrating the tumor, and types of circulating cytokines differ depending on the strain. Therefore, to better comprehend the behavior of cancer in vivo, one must xenograft multiple mouse strains. Here we describe and report a series of methods that we used to reveal the genes and proteins expressed when the same cancer cell line, MDA-MB-231, is xenografted in different hosts. First, using proteomic analysis, we show how to use the same cell line in vivo to reveal the protein changes in the neoplastic cell that help it adapt to its host. Then, we show how different hosts respond molecularly to the same cell line. We also find that using multiple strains can reveal a more suitable host than those traditionally used for a "difficult to xenograft" PDX. In addition, using complex trait genetics, we illustrate a feasible method for uncovering the alleles of the host that support tumor growth. Finally, we demonstrate that Diversity Outbred mice, the epitome of a model of mouse-strain genetic diversity, can be xenografted with human cell lines or PDX using 2-deoxy-D-glucose treatment.

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 genetically diverse C. elegans recombinant intercross advanced inbred 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 positively correlated with developmental time. Among the top candidates obtained from multi-omics data integration and QTL mapping, we validated known and novel longevity modulators, including rict-1, gfm-1 and mltn-1. We translated their relevance to humans using UK Biobank data and showed that variants in RICTOR and GFM1 are associated with an elevated risk of age-related heart disease, dementia, diabetes, kidney, and liver diseases. We organized our dataset as a resource (https://lisp-lms.shinyapps.io/RIAILs/) that allows interactive explorations for new longevity targets.

Hippocampus Glutathione S Reductase Potentially Confers Genetic Resilience to Cognitive Decline in the AD-BXD Mouse Population

Alzheimer's disease (AD) is a prevalent and costly age-related dementia. Heritable factors account for 58-79% of variation in late-onset AD, but substantial variation remains in age-of- onset, disease severity, and whether those with high-risk genotypes acquire AD. To emulate the diversity of human populations, we utilized the AD-BXD mouse panel. This genetically diverse resource combines AD genotypes with multiple BXD strains to discover new genetic drivers of AD resilience. Comparing AD-BXD carriers to noncarrier littermates, we computed a novel quantitative metric for resilience to cognitive decline in the AD-BXDs. Our quantitative AD resilience trait was heritable and genetic mapping identified a locus on chr8 associated with resilience to AD mutations that resulted in amyloid brain pathology. Using a hippocampus proteomics dataset, we nominated the mitochondrial glutathione S reductase protein (GR or GSHR) as a resilience factor, finding that the DBA/2J genotype was associated with substantially higher GR abundance. By mapping protein QTLs (pQTLs), we identified synaptic organization and mitochondrial proteins coregulated in trans with a cis-pQTL for GR. We found four coexpression modules correlated with the quantitative resilience score in aged 5XFAD mice using paracliques, which were related to cell structure, protein folding, and postsynaptic densities. Finally, we found significant positive associations between human GSR transcript abundance in the brain and better outcomes on AD-related cognitive and pathology traits in the Religious Orders Study/Memory and Aging project (ROSMAP). Taken together, these data support a framework for resilience in which neuronal antioxidant pathway activity provides for stability of synapses within the hippocampus.

Functional allele of a MATE gene selected during domestication modulates seed color in chickpea

Seed color is one of the key target traits of domestication and artificial selection in chickpeas due to its implications on consumer preference and market value. The complex seed color trait has been well dissected in several crop species; however, the genetic mechanism underlying seed color variation in chickpea remains poorly understood. Here, we employed an integrated genomics strategy involving QTL mapping, high-density mapping, map-based cloning, association analysis, and molecular haplotyping in an inter-specific RIL mapping population, association panel, wild accessions, and introgression lines (ILs) of Cicer gene pool. This delineated a MATE gene, CaMATE23, encoding a Transparent Testa (TT) and its natural allele (8-bp insertion) and haplotype underlying a major QTL governing seed color on chickpea chromosome 4. Signatures of selective sweep and a strong purifying selection reflected that CaMATE23, especially its 8-bp insertion natural allelic variant, underwent selection during chickpea domestication. Functional investigations revealed that the 8-bp insertion containing the third cis-regulatory RY-motif element in the CaMATE23 promoter is critical for enhanced binding of CaFUSCA3 transcription factor, a key regulator of seed development and flavonoid biosynthesis, thereby affecting CaMATE23 expression and proanthocyanidin (PA) accumulation in the seed coat to impart varied seed color in chickpea. Consequently, overexpression of CaMATE23 in Arabidopsis tt12 mutant partially restored the seed color phenotype to brown pigmentation, ascertaining its functional role in PA accumulation in the seed coat. These findings shed new light on the seed color regulation and evolutionary history, and highlight the transcriptional regulation of CaMATE23 by CaFUSCA3 in modulating seed color in chickpea. The functionally relevant InDel variation, natural allele, and haplotype from CaMATE23 are vital for translational genomic research, including marker-assisted breeding, for developing chickpea cultivars with desirable seed color that appeal to consumers and meet global market demand.

Genetic analysis of pericarp pigmentation variation in Corn Belt dent maize

The US standard for maize commercially grown for grain specifies that yellow corn can contain at maximum 5% corn of other colors. Inbred parents of commercial hybrids typically have clear pericarp, but transgressive segregants in breeding populations can display variation in pericarp pigmentation. We identified 10 doubled haploid biparental populations segregating for pigmented pericarp and evaluated qualitative genetic models using chi-square tests of observed and expected frequencies. Pigmentation ranged from light to dark brown color, and pigmentation intensity was quantitatively measured across 1,327 inbred lines using hue calculated from RGB pixel values. Genetic mapping was used to identify loci associated with pigmentation intensity. For 9 populations, pigmentation inheritance best fit a hypothesis of a 2- or 3-gene epistatic model. Significant differences in pigment intensity were observed across populations. W606S-derived inbred lines with the darkest pericarp often had clear glumes, suggesting the presence of a novel P1-rw allele, a hypothesis supported by a significant quantitative trait locus peak at P1. A separate quantitative trait locus region on chromosome 2 between 221.64 and 226.66 Mbp was identified in LH82-derived populations, and the peak near p1 was absent. A genome-wide association study using 416 inbred lines from the Wisconsin Diversity panel with full genome resequencing revealed 4 significant associations including the region near P1. This study supports that pericarp pigmentation among dent maize inbreds can arise by transgressive segregation when pigmentation in the parental generation is absent and is partially explained by functional allelic variation at the P1 locus.

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

Mycobacterium tuberculosis, the bacillus that causes tuberculosis (TB), infects 2 billion people across the globe, and results in 8-9 million new TB cases and 1-1.5 million deaths each year. Most patients have no known genetic basis that predisposes them to disease. We investigated 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 clinical indicators of disease, granuloma histopathological features, and immune response traits identified five new loci on mouse chromosomes 1, 2, 4, 16 and three previously identified loci on chromosomes 3 and 17. Quantitative trait loci (QTLs) on chromosomes 1, 16, and 17, associated with multiple correlated traits and had similar patterns of allele effects, suggesting these QTLs contain important genetic regulators of responses to M. tuberculosis. To narrow the list of candidate genes in QTLs, we used a machine learning strategy that integrated gene expression signatures from lungs of M. tuberculosis-infected Diversity Outbred mice with gene interaction networks, generating functional scores. The scores were then used to rank candidates for each mapped trait in each locus, resulting in 11 candidates: Ncf2, Fam20b, S100a8, S100a9, Itgb5, Fstl1, Zbtb20, Ddr1, Ier3, Vegfa, and Zfp318. Importantly, all 11 candidates have roles in infection, inflammation, cell migration, extracellular matrix remodeling, or intracellular signaling. Further, all candidates contain single nucleotide polymorphisms (SNPs), and some but not all SNPs were predicted to have deleterious consequences on protein functions. Multiple methods were used for validation including (i) a statistical method that showed Diversity Outbred mice carrying PWH/PhJ alleles on chromosome 17 QTL have shorter survival; (ii) quantification of S100A8 protein levels, confirming predicted allele effects; and (iii) infection of C57BL/6 mice deficient for the S100a8 gene. Overall, this work demonstrates that systems genetics using Diversity Outbred mice can identify new (and known) QTLs and new functionally relevant gene candidates that may be major regulators of granuloma necrosis and acute inflammation in pulmonary TB.

BulkLMM: Real-time genome scans for multiple quantitative traits using linear mixed models

Genetic studies often collect data using high-throughput phenotyping. That has led to the need for fast genomewide scans for large number of traits using linear mixed models (LMMs). Computing the scans one by one on each trait is time consuming. We have developed new algorithms for performing genome scans on a large number of quantitative traits using LMMs, BulkLMM, that speeds up the computation by orders of magnitude compared to one trait at a time scans. On a mouse BXD Liver Proteome data with more than 35,000 traits and 7,000 markers, BulkLMM completed in a few seconds. We use vectorized, multi-threaded operations and regularization to improve optimization, and numerical approximations to speed up the computations. Our software implementation in the Julia programming language also provides permutation testing for LMMs and is available at https://github.com/senresearch/BulkLMM.jl.

Fine mapping of the Cepaea nemoralis shell colour and mid-banded loci using a high-density linkage map

Molluscs are a highly speciose phylum that exhibits an astonishing array of colours and patterns, yet relatively little progress has been made in identifying the underlying genes that determine phenotypic variation. One prominent example is the land snail Cepaea nemoralis for which classical genetic studies have shown that around nine loci, several physically linked and inherited together as a 'supergene', control the shell colour and banding polymorphism. As a first step towards identifying the genes involved, we used whole-genome resequencing of individuals from a laboratory cross to construct a high-density linkage map, and then trait mapping to identify 95% confidence intervals for the chromosomal region that contains the supergene, specifically the colour locus (C), and the unlinked mid-banded locus (U). The linkage map is made up of 215,593 markers, ordered into 22 linkage groups, with one large group making up ~27% of the genome. The C locus was mapped to a ~1.3 cM region on linkage group 11, and the U locus was mapped to a ~0.7 cM region on linkage group 15. The linkage map will serve as an important resource for further evolutionary and population genomic studies of C. nemoralis and related species, as well as the identification of candidate genes within the supergene and for the mid-banding phenotype.

Increased development of T-bet+CD11c+ B cells predisposes to lupus in females: Analysis in BXD2 mouse and genetic crosses

Cardinal features of lupus include elevated B cell activation and autoantibody production with a female sex preponderance. We quantified interactions of sex and genetic variation on the development of autoimmune B-cell phenotypes and autoantibodies in the BXD2 murine model of lupus using a cohort of backcrossed progeny (BXD2 x C57BL/6J) x BXD2. Sex was the key factor leading to increased total IgG, IgG2b, and autoantibodies. The percentage of T-bet+CD11c+ IgD+ activated naive B cells (aNAV) was higher in females and was associated with increased T-bet+CD11c+ IgD- age-related B cells, Fas+GL7+ germinal center B cells, Cxcr5-Icos+ peripheral T-helper cells, and Cxcr5+Icos+ follicular T-helper cells. IFN-β was elevated in females. Variation in aNAV cells was mapped to Chr 7 in a locus that shows significant interactions between the female sex and heterozygous B/D variant. Our results suggest that activation of naive B cells forms the basis for the female-predominant development of autoantibodies in lupus-susceptible BXD2 mice.

Genetic control of grain amino acid composition in a UK soft wheat mapping population

Wheat (Triticum aestivum L.) is a major source of nutrients for populations across the globe, but the amino acid composition of wheat grain does not provide optimal nutrition. The nutritional value of wheat grain is limited by low concentrations of lysine (the most limiting essential amino acid) and high concentrations of free asparagine (precursor to the processing contaminant acrylamide). There are currently few available solutions for asparagine reduction and lysine biofortification through breeding. In this study, we investigated the genetic architecture controlling grain free amino acid composition and its relationship to other traits in a Robigus × Claire doubled haploid population. Multivariate analysis of amino acids and other traits showed that the two groups are largely independent of one another, with the largest effect on amino acids being from the environment. Linkage analysis of the population allowed identification of quantitative trait loci (QTL) controlling free amino acids and other traits, and this was compared against genomic prediction methods. Following identification of a QTL controlling free lysine content, wheat pangenome resources facilitated analysis of candidate genes in this region of the genome. These findings can be used to select appropriate strategies for lysine biofortification and free asparagine reduction in wheat breeding programs.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

Epistasis between mutator alleles contributes to germline mutation spectra variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair [1], mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs [2,3]. In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh [4]. Its effect depended on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci had greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Mapping the genomic landscape of multidrug resistance in Plasmodium falciparum and its impact on parasite fitness

Drug-resistant Plasmodium falciparum parasites have swept across Southeast Asia and now threaten Africa. By implementing a P. falciparum genetic cross using humanized mice, we report the identification of key determinants of resistance to artemisinin (ART) and piperaquine (PPQ) in the dominant Asian KEL1/PLA1 lineage. We mapped k13 as the central mediator of ART resistance in vitro and identified secondary markers. Applying bulk segregant analysis, quantitative trait loci mapping using 34 recombinant haplotypes, and gene editing, our data reveal an epistatic interaction between mutant PfCRT and multicopy plasmepsins 2/3 in mediating high-grade PPQ resistance. Susceptibility and parasite fitness assays implicate PPQ as a driver of selection for KEL1/PLA1 parasites. Mutant PfCRT enhanced susceptibility to lumefantrine, the first-line partner drug in Africa, highlighting a potential benefit of opposing selective pressures with this drug and PPQ. We also identified that the ABCI3 transporter can operate in concert with PfCRT and plasmepsins 2/3 in mediating multigenic resistance to antimalarial agents.

Development of genomic markers associated to growth-related traits and sex determination in lumpfish (Cyclopterus lumpus)

Cleaner fish species have gained great importance in the control of sea lice, among them, lumpfish (Cyclopterus lumpus) has become one of the most popular. Lumpfish life cycle has been closed, and hatchery reproduction is now possible, however, current production is reliant on wild caught broodstock to meet the increasing demand. Selective breeding practices are called to play an important role in the successful breeding of most aquaculture species, including lumpfish. In this study we analysed a lumpfish population for the identification of genomic markers linked to production traits. Sequencing of RAD libraries allowed us to identify, 7193 informative markers within the sampled individuals. Genome wide association analysis for sex, weight, condition factor and standard length was performed. One single major QTL region was identified for sex, while nine QTL regions were detected for weight, and three QTL regions for standard length. A total of 177 SNP markers of interest (from QTL regions) and 399 high Fst SNP markers were combined in a low-density panel, useful to obtain relevant genetic information from lumpfish populations. Moreover, a robust combined subset of 29 SNP markers (10 associated to sex, 14 to weight and 18 to standard length) provided over 90% accuracy in predicting the animal's phenotype by machine learning. Overall, our findings provide significant insights into the genetic control of important traits in lumpfish and deliver important genomic resources that will facilitate the establishment of selective breeding programmes in lumpfish.

A polygenic explanation for Haldane's rule in butterflies

Two robust rules have been discovered about animal hybrids: Heterogametic hybrids are more unfit (Haldane's rule), and sex chromosomes are disproportionately involved in hybrid incompatibility (the large-X/Z effect). The exact mechanisms causing these rules in female heterogametic taxa such as butterflies are unknown but are suggested by theory to involve dominance on the sex chromosome. We investigate hybrid incompatibilities adhering to both rules in Papilio and Heliconius butterflies and show that dominance theory cannot explain our data. Instead, many defects coincide with unbalanced multilocus introgression between the Z chromosome and all autosomes. Our polygenic explanation predicts both rules because the imbalance is likely greater in heterogametic females, and the proportion of introgressed ancestry is more variable on the Z chromosome. We also show that mapping traits polygenic on a single chromosome in backcrosses can generate spurious large-effect QTLs. This mirage is caused by statistical linkage among polygenes that inflates estimated effect sizes. By controlling for statistical linkage, most incompatibility QTLs in our hybrid crosses are consistent with a polygenic basis. Since the two genera are very distantly related, polygenic hybrid incompatibilities are likely common in butterflies.

Exploring horticultural traits and disease resistance in Capsicum baccatum through segmental introgression lines

KEY MESSAGE

Segmental introgression and advanced backcross lines were developed and validated as important tools for improving agronomically important traits in pepper, offering improved sensitivity in detecting quantitative trait loci for breeding. Segmental introgression lines (SILs) and advanced backcross lines (ABs) can accelerate genetics and genomics research and breeding in crop plants. This study presents the development of a complete collection of SILs and ABs in pepper using Capsicum annuum cv. 'CM334' as the recipient parent and Capsicum baccatum 'PBC81', which displays various agronomically important traits including powdery mildew and anthracnose resistance, as donor parent. Using embryo rescue to overcome abortion in interspecific crosses, and marker-assisted selection with genotyping-in-thousands by sequencing (GT-seq) to develop SILs and ABs containing different segments of the C. baccatum genome, we obtained 63 SILs and 44 ABs, covering 94.8% of the C. baccatum genome. We characterized them for traits including powdery mildew resistance, anthracnose resistance, anthocyanin accumulation, trichome density, plant architecture, and fruit morphology. We validated previously known loci for these traits and discovered new sources of variation and quantitative trait loci (QTLs). A total of 15 QTLs were identified, including four for anthracnose resistance with three novel loci, seven for plant architecture, and four for fruit morphology. This is the first complete collection of pepper SILs and ABs validated for agronomic traits and will enhance QTL detection and serve as valuable breeding resources. Further, these SILs and ABs will be useful for comparative genomics and to better understand the genetic mechanisms underlying important agronomic traits in pepper, ultimately leading to improved crop productivity and sustainability.

Genetic Variability in Seed Longevity and Germination Traits in a Tomato MAGIC Population in Contrasting Environments

The stable production of high vigorous seeds is pivotal to crop yield. Also, a high longevity is essential to avoid progressive loss of seed vigour during storage. Both seed traits are strongly influenced by the environment during seed development. Here, we investigated the impact of heat stress (HS) during fruit ripening on tomato seed lifespan during storage at moderate relative humidity, speed (t50) and homogeneity of germination, using a MAGIC population that was produced under optimal and HS conditions. A plasticity index was used to assess the extent of the impact of HS for each trait. HS reduced the average longevity and germination homogeneity by 50% within the parents and MAGIC population. However, there was a high genetic variability in the seed response to heat stress. A total of 39 QTLs were identified, including six longevity QTLs for seeds from control (3) and HS (3) conditions, and six plasticity QTLs for longevity, with only one overlapping with a longevity QTL under HS. Four out of the six longevity QTL co-located with t50 QTL, revealing hotspots for seed quality traits. Twenty-one QTLs with intervals below 3 cM were analyzed using previous transcriptome and gene network data to propose candidate genes for seed vigour and longevity traits.

Genetic and dietary modulators of the inflammatory response in the gastrointestinal tract of the BXD mouse genetic reference population

Inflammatory gut disorders, including inflammatory bowel disease (IBD), can be impacted by dietary, environmental, and genetic factors. While the incidence of IBD is increasing worldwide, we still lack a complete understanding of the gene-by-environment interactions underlying inflammation and IBD. Here, we profiled the colon transcriptome of 52 BXD mouse strains fed with a chow or high-fat diet (HFD) and identified a subset of BXD strains that exhibit an IBD-like transcriptome signature on HFD, indicating that an interplay of genetics and diet can significantly affect intestinal inflammation. Using gene co-expression analyses, we identified modules that are enriched for IBD-dysregulated genes and found that these IBD-related modules share cis-regulatory elements that are responsive to the STAT2, SMAD3, and REL transcription factors. We used module quantitative trait locus analyses to identify genetic loci associated with the expression of these modules. Through a prioritization scheme involving systems genetics in the mouse and integration with external human datasets, we identified Muc4 and Epha6 as the top candidates mediating differences in HFD-driven intestinal inflammation. This work provides insights into the contribution of genetics and diet to IBD risk and identifies two candidate genes, MUC4 and EPHA6, that may mediate IBD susceptibility in humans.

Identifying QTLs involved in hybrid performance and heterotic group complementarity: new GWAS models applied to factorial and admixed diallel maize hybrid panels

KEY MESSAGE

An original GWAS model integrating the ancestry of alleles was proposed and allowed the detection of background specific additive and dominance QTLs involved in heterotic group complementarity and hybrid performance. Maize genetic diversity is structured into genetic groups selected and improved relative to each other. This process increases group complementarity and differentiation over time and ensures that the hybrids produced from inter-group crosses exhibit high performances and heterosis. To identify loci involved in hybrid performance and heterotic group complementarity, we introduced an original association study model that disentangles allelic effects from the heterotic group origin of the alleles and compared it with a conventional additive/dominance model. This new model was applied on a factorial between Dent and Flint lines and a diallel between Dent-Flint admixed lines with two different layers of analysis: within each environment and in a multiple-environment context. We identified several strong additive QTLs for all traits, including some well-known additive QTLs for flowering time (in the region of Vgt1/2 on chromosome 8). Yield trait displayed significant non-additive effects in the diallel panel. Most of the detected Yield QTLs exhibited overdominance or, more likely, pseudo-overdominance effects. Apparent overdominance at these QTLs contributed to a part of the genetic group complementarity. The comparison between environments revealed a higher stability of additive QTL effects than non-additive ones. Several QTLs showed variations of effects according to the local heterotic group origin. We also revealed large chromosomic regions that display genetic group origin effects. Altogether, our results illustrate how admixed panels combined with dedicated GWAS modeling allow the identification of new QTLs that could not be revealed by a classical hybrid panel analyzed with traditional modeling.