Into the Wild: A novel wild-derived inbred strain resource expands the genomic and phenotypic diversity of laboratory mouse models.

The laboratory mouse has served as the premier animal model system for both basic and preclinical investigations for over a century. However, laboratory mice capture only a subset of the genetic variation found in wild mouse populations, ultimately limiting the potential of classical inbred strains to uncover phenotype-associated variants and pathways. Wild mouse populations are reservoirs of genetic diversity that could facilitate the discovery of new functional and disease-associated alleles, but the scarcity of commercially available, well-characterized wild mouse strains limits their broader adoption in biomedical research. To overcome this barrier, we have recently developed, sequenced, and phenotyped a set of 11 inbred strains derived from wild-caught Mus musculus domesticus. Each of these "Nachman strains" immortalizes a unique wild haplotype sampled from one of five environmentally distinct locations across North and South America. Whole genome sequence analysis reveals that each strain carries between 4.73-6.54 million single nucleotide differences relative to the GRCm39 mouse reference, with 42.5% of variants in the Nachman strain genomes absent from current classical inbred mouse strain panels. We phenotyped the Nachman strains on a customized pipeline to assess the scope of disease-relevant neurobehavioral, biochemical, physiological, metabolic, and morphological trait variation. The Nachman strains exhibit significant inter-strain variation in >90% of 1119 surveyed traits and expand the range of phenotypic diversity captured in classical inbred strain panels. These novel wild-derived inbred mouse strain resources are set to empower new discoveries in both basic and preclinical research.

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.

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.

B cells in perivascular and peribronchiolar granuloma-associated lymphoid tissue and B-cell signatures identify asymptomatic Mycobacterium tuberculosis lung infection in Diversity Outbred mice.

Because most humans resist Mycobacterium tuberculosis infection, there is a paucity of lung samples to study. To address this gap, we infected Diversity Outbred mice with M. tuberculosis and studied the lungs of mice in different disease states. After a low-dose aerosol infection, progressors succumbed to acute, inflammatory lung disease within 60 days, while controllers maintained asymptomatic infection for at least 60 days, and then developed chronic pulmonary tuberculosis (TB) lasting months to more than 1 year. Here, we identified features of asymptomatic M. tuberculosis infection by applying computational and statistical approaches to multimodal data sets. Cytokines and anti-M. tuberculosis cell wall antibodies discriminated progressors vs controllers with chronic pulmonary TB but could not classify mice with asymptomatic infection. However, a novel deep-learning neural network trained on lung granuloma images was able to accurately classify asymptomatically infected lungs vs acute pulmonary TB in progressors vs chronic pulmonary TB in controllers, and discrimination was based on perivascular and peribronchiolar lymphocytes. Because the discriminatory lesion was rich in lymphocytes and CD4 T cell-mediated immunity is required for resistance, we expected CD4 T-cell genes would be elevated in asymptomatic infection. However, the significantly different, highly expressed genes were from B-cell pathways (e.g., Bank1, Cd19, Cd79, Fcmr, Ms4a1, Pax5, and H2-Ob), and CD20+ B cells were enriched in the perivascular and peribronchiolar regions of mice with asymptomatic M. tuberculosis infection. Together, these results indicate that genetically controlled B-cell responses are important for establishing asymptomatic M. tuberculosis lung infection.

CXCL1: A new diagnostic biomarker for human tuberculosis discovered using Diversity Outbred mice.

More humans have died of tuberculosis (TB) than any other infectious disease and millions still die each year. Experts advocate for blood-based, serum protein biomarkers to help diagnose TB, which afflicts millions of people in high-burden countries. However, the protein biomarker pipeline is small. Here, we used the Diversity Outbred (DO) mouse population to address this gap, identifying five protein biomarker candidates. One protein biomarker, serum CXCL1, met the World Health Organization's Targeted Product Profile for a triage test to diagnose active TB from latent M.tb infection (LTBI), non-TB lung disease, and normal sera in HIV-negative, adults from South Africa and Vietnam. To find the biomarker candidates, we quantified seven immune cytokines and four inflammatory proteins corresponding to highly expressed genes unique to progressor DO mice. Next, we applied statistical and machine learning methods to the data, i.e., 11 proteins in lungs from 453 infected and 29 non-infected mice. After searching all combinations of five algorithms and 239 protein subsets, validating, and testing the findings on independent data, two combinations accurately diagnosed progressor DO mice: Logistic Regression using MMP8; and Gradient Tree Boosting using a panel of 4: CXCL1, CXCL2, TNF, IL-10. Of those five protein biomarker candidates, two (MMP8 and CXCL1) were crucial for classifying DO mice; were above the limit of detection in most human serum samples; and had not been widely assessed for diagnostic performance in humans before. In patient sera, CXCL1 exceeded the triage diagnostic test criteria (>90% sensitivity; >70% specificity), while MMP8 did not. Using Area Under the Curve analyses, CXCL1 averaged 94.5% sensitivity and 88.8% specificity for active pulmonary TB (ATB) vs LTBI; 90.9% sensitivity and 71.4% specificity for ATB vs non-TB; and 100.0% sensitivity and 98.4% specificity for ATB vs normal sera. Our findings overall show that the DO mouse population can discover diagnostic-quality, serum protein biomarkers of human TB.

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.

Defining the consequences of genetic variation on a proteome-wide scale.

Genetic variation modulates protein expression through both transcriptional and post-transcriptional mechanisms. To characterize the consequences of natural genetic diversity on the proteome, here we combine a multiplexed, mass spectrometry-based method for protein quantification with an emerging outbred mouse model containing extensive genetic variation from eight inbred founder strains. By measuring genome-wide transcript and protein expression in livers from 192 Diversity outbred mice, we identify 2,866 protein quantitative trait loci (pQTL) with twice as many local as distant genetic variants. These data support distinct transcriptional and post-transcriptional models underlying the observed pQTL effects. Using a sensitive approach to mediation analysis, we often identified a second protein or transcript as the causal mediator of distant pQTL. Our analysis reveals an extensive network of direct protein-protein interactions. Finally, we show that local genotype can provide accurate predictions of protein abundance in an independent cohort of collaborative cross mice.

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome 3@52.5 Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass.

Genetic determinants of gut microbiota composition and bile acid profiles in mice.

The microbial communities that inhabit the distal gut of humans and other mammals exhibit large inter-individual variation. While host genetics is a known factor that influences gut microbiota composition, the mechanisms underlying this variation remain largely unknown. Bile acids (BAs) are hormones that are produced by the host and chemically modified by gut bacteria. BAs serve as environmental cues and nutrients to microbes, but they can also have antibacterial effects. We hypothesized that host genetic variation in BA metabolism and homeostasis influence gut microbiota composition. To address this, we used the Diversity Outbred (DO) stock, a population of genetically distinct mice derived from eight founder strains. We characterized the fecal microbiota composition and plasma and cecal BA profiles from 400 DO mice maintained on a high-fat high-sucrose diet for ~22 weeks. Using quantitative trait locus (QTL) analysis, we identified several genomic regions associated with variations in both bacterial and BA profiles. Notably, we found overlapping QTL for Turicibacter sp. and plasma cholic acid, which mapped to a locus containing the gene for the ileal bile acid transporter, Slc10a2. Mediation analysis and subsequent follow-up validation experiments suggest that differences in Slc10a2 gene expression associated with the different strains influences levels of both traits and revealed novel interactions between Turicibacter and BAs. This work illustrates how systems genetics can be utilized to generate testable hypotheses and provide insight into host-microbe interactions.

Uncovering Modifier Genes of X-Linked Alport Syndrome Using a Novel Multiparent Mouse Model.

BACKGROUND: Mutations in COL4A5 are responsible for 80% of cases of X-linked Alport Syndrome (XLAS). Although genes that cause AS are well characterized, people with AS who have similar genetic mutations present with a wide variation in the extent of kidney impairment and age of onset, suggesting the activities of modifier genes. METHODS: We created a cohort of genetically diverse XLAS male and female mice using the Diversity Outbred mouse resource and measured albuminuria, GFR, and gene expression. Using a quantitative trait locus approach, we mapped modifier genes that can best explain the underlying phenotypic variation measured in our diverse population. RESULTS: Genetic analysis identified several loci associated with the variation in albuminuria and GFR, including a locus on the X chromosome associated with X inactivation and a locus on chromosome 2 containing Fmn1. Subsequent analysis of genetically reduced Fmn1 expression in Col4a5 knockout mice showed a decrease in albuminuria, podocyte effacement, and podocyte protrusions in the glomerular basement membrane, which support the candidacy of Fmn1 as a modifier gene for AS. CONCLUSION: With this novel approach, we emulated the variability in the severity of kidney phenotypes found in human patients with Alport Syndrome through albuminuria and GFR measurements. This approach can identify modifier genes in kidney disease that can be used as novel therapeutic targets.

Nine quick tips for efficient bioinformatics curriculum development and training.

Biomedical research is becoming increasingly data driven. New technologies that generate large-scale, complex data are continually emerging and evolving. As a result, there is a concurrent need for training researchers to use and understand new computational tools. Here we describe an efficient and effective approach to developing curriculum materials that can be deployed in a research environment to meet this need.

Correction to: Quantitative trait mapping in Diversity Outbred mice identifies two genomic regions associated with heart size.

The original article has been published with an incorrect text in Materials and Methods section. The corrected text should read as.

Identification of a novel microRNA profile in pediatric patients with cancer treated with anthracycline chemotherapy.

Anthracycline chemotherapy (AC) is associated with decline in left ventricular ejection fraction (LVEF), yet the mechanisms remain unclear. Although changes in microRNAs (miRs) have been identified in adult cardiovascular disease, miR profiles in pediatric patients with AC have not been well studied. The goal of this study was to examine miR profiles (unbiased array) in pediatric patients with AC compared with age-matched referent normal patients. We hypothesize that pediatric patients with AC will express a unique miR profile at the initiation and completion of therapy and will be related to LVEF. Serum was collected in pediatric patients (10-22 yr, n = 12) with newly diagnosed malignancy requiring AC within 24-48 h after the initiation of therapy (30-60 mg/m2) and ~1 yr after completing therapy. A custom microarray of 84 miRs associated with cardiovascular disease was used (quantitative RT-PCR) and indexed to referent normal profiles (13-17 yr, n = 17). LVEF was computed by cardiac MRI. LVEF fell from AC initiation at ~1 yr after AC completion (64.28 ± 1.78% vs. 57.53 ± 0.95%, respectively, P = 0.004). Of the 84 miRs profiled, significant shifts in 17 miRs occurred relative to referent normal ( P ≤ 0.05). Moreover, the functional domain of miRs associated with myocardial differentiation and development fell over threefold at the completion of AC ( P ≤ 0.05). Moreover, eight miRs were significantly downregulated after AC completion in those patients with the greatest decline in LVEF (≥10%, P < 0.05). This study demonstrates, for the first time, that changes in miR expression occur in pediatric patients with AC. These findings suggest that miRs are a potential strategy for the early identification of patients with AC susceptible to left ventricular dysfunction. NEW & NOTEWORTHY Although anthracycline chemotherapy (AC) is effective for a number of pediatric cancers, an all too often consequence of AC is the development of left ventricular failure. The present study identified that specific shifts in the pattern of microRNAs, which regulate myocardial growth, function, and viability, occurred during and after AC in pediatric patients, whereby the magnitude of this shift was associated with the degree of left ventricular failure.

Quantitative trait mapping in Diversity Outbred mice identifies two genomic regions associated with heart size.

Heart size is an important factor in cardiac health and disease. In particular, increased heart weight is predictive of adverse cardiovascular outcomes in multiple large community-based studies. We use two cohorts of Diversity Outbred (DO) mice to investigate the role of genetics, sex, age, and diet on heart size. DO mice (n = 289) of both sexes from generation 10 were fed a standard chow diet, and analyzed at 12-15 weeks of age. Another cohort of female DO mice (n = 258) from generation 11 were fed either a high-fat, cholesterol-containing (HFC) diet or a low-fat, high-protein diet, and analyzed at 24-25 weeks. We did not observe an effect of diet on body or heart weight in generation 11 mice, although we previously reported an effect on other cardiovascular risk factors, including cholesterol, triglycerides, and insulin. We do observe a significant genetic effect on heart weight in this population. We identified two quantitative trait loci for heart weight, one (Hwtf1) at a genome-wide significance level of p ≤ 0.05 on MMU15 and one (Hwtf2) at a genome-wide suggestive level of p ≤ 0.1 on MMU10, that together explain 13.3% of the phenotypic variance. Hwtf1 contained collagen type XXII alpha 1 chain (Col22a1), and the NZO/HlLtJ and WSB/EiJ haplotypes were associated with larger hearts. This is consistent with heart tissue Col22a1 expression in DO founders and SNP patterns within Hwtf1 for Col22a1. Col22a1 has been previously associated with cardiac fibrosis in mice, suggesting that Col22a1 may be involved in pathological cardiac hypertrophy.

Genome-wide association for testis weight in the diversity outbred mouse population.

Testis weight is a genetically mediated trait associated with reproductive efficiency across numerous species. We sought to evaluate the genetically diverse, highly recombinant Diversity Outbred (DO) mouse population as a tool to identify and map quantitative trait loci (QTLs) associated with testis weight. Testis weights were recorded for 502 male DO mice and the mice were genotyped on the GIGAMuga array at ~ 143,000 SNPs. We performed a genome-wide association analysis and identified one significant and two suggestive QTLs associated with testis weight. Using bioinformatic approaches, we developed a list of candidate genes and identified those with known roles in testicular size and development. Candidates of particular interest include the RNA demethylase gene Alkbh5, the cyclin-dependent kinase inhibitor gene Cdkn2c, the dynein axonemal heavy chain gene Dnah11, the phospholipase D gene Pld6, the trans-acting transcription factor gene Sp4, and the spermatogenesis-associated gene Spata6, each of which has a human ortholog. Our results demonstrate the utility of DO mice in high-resolution genetic mapping of complex traits, enabling us to identify developmentally important genes in adult mice. Understanding how genetic variation in these genes influence testis weight could aid in the understanding of mechanisms of mammalian reproductive function.

CloudNeo: a cloud pipeline for identifying patient-specific tumor neoantigens.

SUMMARY: We present CloudNeo, a cloud-based computational workflow for identifying patient-specific tumor neoantigens from next generation sequencing data. Tumor-specific mutant peptides can be detected by the immune system through their interactions with the human leukocyte antigen complex, and neoantigen presence has recently been shown to correlate with anti T-cell immunity and efficacy of checkpoint inhibitor therapy. However computing capabilities to identify neoantigens from genomic sequencing data are a limiting factor for understanding their role. This challenge has grown as cancer datasets become increasingly abundant, making them cumbersome to store and analyze on local servers. Our cloud-based pipeline provides scalable computation capabilities for neoantigen identification while eliminating the need to invest in local infrastructure for data transfer, storage or compute. The pipeline is a Common Workflow Language (CWL) implementation of human leukocyte antigen (HLA) typing using Polysolver or HLAminer combined with custom scripts for mutant peptide identification and NetMHCpan for neoantigen prediction. We have demonstrated the efficacy of these pipelines on Amazon cloud instances through the Seven Bridges Genomics implementation of the NCI Cancer Genomics Cloud, which provides graphical interfaces for running and editing, infrastructure for workflow sharing and version tracking, and access to TCGA data. AVAILABILITY AND IMPLEMENTATION: The CWL implementation is at: https://github.com/TheJacksonLaboratory/CloudNeo. For users who have obtained licenses for all internal software, integrated versions in CWL and on the Seven Bridges Cancer Genomics Cloud platform (https://cgc.sbgenomics.com/, recommended version) can be obtained by contacting the authors. CONTACT: jeff.chuang@jax.org. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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

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

R2d2 Drives Selfish Sweeps in the House Mouse.

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

Block network mapping approach to quantitative trait locus analysis.

BACKGROUND: Advances in experimental biology have enabled the collection of enormous troves of data on genomic variation in living organisms. The interpretation of this data to extract actionable information is one of the keys to developing novel therapeutic strategies to treat complex diseases. Network organization of biological data overcomes measurement noise in several biological contexts. Does a network approach, combining information about the linear organization of genomic markers with correlative information on these markers in a Bayesian formulation, lead to an analytic method with higher power for detecting quantitative trait loci? RESULTS: Block Network Mapping, combining Similarity Network Fusion (Wang et al., NM 11:333-337, 2014) with a Bayesian locus likelihood evaluation, leads to large improvements in area under the receiver operating characteristic and power over interval mapping with expectation maximization. The method has a monotonically decreasing false discovery rate as a function of effect size, unlike interval mapping. CONCLUSIONS: Block Network Mapping provides an alternative data-driven approach to mapping quantitative trait loci that leverages correlations in the sampled genotypes. The evaluation methodology can be combined with existing approaches such as Interval Mapping. Python scripts are available at http://lbm.niddk.nih.gov/vipulp/ . Genotype data is available at http://churchill-lab.jax.org/website/GattiDOQTL .