Advertisement

  • Loading metrics

Open Access

Review

Mouse protein coding diversity: What’s left to discover?

Abstract

For over a century, mice have been used to model human disease, leading to many fundamental discoveries about mammalian biology and the development of new therapies. Mouse genetics research has been further catalysed by a plethora of genomic resources developed in the last 20 years, including the genome sequence of C57BL/6J and more recently the first draft reference genomes for 16 additional laboratory strains. Collectively, the comparison of these genomes highlights the extreme diversity that exists at loci associated with the immune system, pathogen response, and key sensory functions, which form the foundation for dissecting phenotypic traits in vivo. We review the current status of the mouse genome across the diversity of the mouse lineage and discuss the value of mice to understanding human disease.

Author summary

For decades, the laboratory mouse has been widely used to make fundamental discoveries about human biology, model human disease, and develop new treatments. The mouse reference genome is based on the C57BL/6J; however, researchers use a variety of strains to model human disease. Recent genome analysis has identified that the most highly variable regions of the mouse genome are enriched with genes relevant to disease and infection response. In this review, we discuss what is currently known about these regions, why they are important for human disease modelling, and what is known about their ancestral origins.

Citation: Lilue J, Shivalikanjli A, Adams DJ, Keane TM (2019) Mouse protein coding diversity: What’s left to discover? PLoS Genet 15(11): e1008446. https://doi.org/10.1371/journal.pgen.1008446

Editor: Elizabeth M. C. Fisher, University College London, UNITED KINGDOM

Published: November 14, 2019

Copyright: © 2019 Lilue et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Funding: This work was supported by the Medical Research Council (MR/R017565/1). D.J.A. is supported by Cancer Research-UK and the Wellcome Trust (100956/Z/13/Z). The funders had no role in the preparation of the article.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Although mice and humans have coexisted for many millennia, modern mouse genetics was initiated in the early 20th century [1]. The first genetically homozygous mouse strain was DBA [2], developed to study coat color inheritance and cancer susceptibility. Subsequently, hundreds of genetically defined strains for modelling human diseases and biological processes (e.g., behaviour, carcinogenesis, and immune response against pathogens) were developed. As one of the most important model organisms in biomedical research, the mouse was the second mammalian genome to be sequenced after the human genome [3]. The C57BL/6J reference genome has enabled the creation of detailed molecular maps of mouse diversity [4,5], generation of null alleles, and phenotyping across thousands of genes; and enabled genetic screens at an unprecedented rate [6].

Modern-day laboratory mouse strains are comprised of classical and wild-derived strains. The classical inbred strains have ‘fancy mice’ as their founding ancestors and are largely Mus musculus domesticus derived. Other subspecies—M. musculus musculus and M. musculus castaneus—contribute approximately 4%–14% to the classical strains [7]. A genome-wide haplotype map of 100 classical mouse strains showed that over 97% of their genome is a mosaic of less than 10 haplotypes [8]. Nevertheless, there are many loci in classical mouse strains, like the major histocompatibility complex (MHC), which have extensive haplotypic diversity [9,10]. Wild-derived inbred mouse strains are recent progenies of wild-caught individuals of M. musculus musculus, M. musculus castaneus, and M. spretus origin and therefore contain many divergent haplotypes not shared with classical inbred strains. Wild-derived strains are increasingly employed as mouse models to study phenotypes—such as resistance against Orthomyxovirus [11], virulent Toxoplasma gondii strains (CIM, CAST/EiJ, and PWK/PhJ) [12,13], resistance to anticoagulant rodenticides (SPRET/EiJ) [14], and resistance to cerebral Plasmodium berghei (WLA/Pas) [15].

Limitations of a single mouse reference genome

Wild-derived mouse strains have hundreds of thousands of structural differences and novel haplotypes compared to C57BL/6J [4,5]. Most SNP discoveries, if not all, are based on high-density genotyping or short-read sequencing. Paired-end reads are aligned to the C57BL/6J reference genome to identify SNPs, indels, and structure variations (SVs) [16,17]. This means that using the C57BL/6J reference genome to study these strains is blind to many nonreference loci [16]. In these strain-specific diverse regions (SSDRs), next-generation sequencing (NGS) reads are forced to map incorrectly to other paralogous loci in the reference and are often represented as dense regions of heterozygous SNPs (hSNPs) that disrupt the collinearity between the genome of a mouse strain and the reference [13,16]. SSDRs are enriched for genes associated with immunity, sensory, sexual reproduction, and behaviour [16]. In this review, we will introduce the SSDRs among 16 mouse strains and their potential importance in human biomedical research.

Individual SSDRs associated with phenotypes in mouse inbred strains have long been studied (Fig 1). In 2004, a high resolution whole genome Bacterial Artificial Chromosome (BAC) array analysis reported ‘segmental polymorphisms’ between mouse strains C57BL/6J and 129/Sv [18]. Subsequent work found similar patterns by comparative genomic hybridization analysis and reported 2,094 ‘copy number variations’ (CNVs) in 41 inbred strains [19]. In 2015, an analysis of 351 high-density microarray data for mouse tail samples highlighted 9,634 putative autosomal CNVs affecting 6.87% of the mouse genome [20]. In 2016, Morgan and colleagues performed a genome-wide subsequence diversity test in seven mouse strains and two wild mice samples and reported at least 0.8% of the mouse genome is a ‘genomic revolving door’ with high mutation and recombination rates [21]. In 2018, the first draft de novo assemblies of 16 mouse strains successfully assembled some of these regions, reporting a total of 2,567 SSDRs that encompass 0.5%–2.8% of the mouse genome (Fig 2A and S1 Data), encoding 1,828 coding genes. These genes can be classified into 468 gene families (S1 Table), and 318 (67.7%) have previously been studied in detail. Only 3.1% of gene families have complete sequences (introns and intergenic regions) for multiple mouse strains, 9.8% have coding regions from multiple mouse strains, and most (87.1%) studies draw scientific conclusions based on a single laboratory mouse strain, typically C57BL/6J or a 129 substrain (Fig 2B). SSDRs are enriched for recently transposed long interspersed nuclear elements (LINEs) and long-terminal repeat (LTR) elements, posing a challenge for genome assembly [16] and consequently are often incomplete in the current mouse reference genome.

thumbnail
Download:
Fig 1. Discovery of SSDRs in the mouse genome.

Loci for chromosome 1 and 6 are shown. Li and colleagues [18] defined segmental polymorphisms between 129 and C57BL/6J (blue). Cutler and colleagues [19] cataloged copy number variations in 41 inbred strains (green). Locke and colleagues [20] identified genome regions with high copy number variation calls in 351 different mouse strains and wild-caught mice (yellow). Morgan and colleagues [21] used a combination of wild and inbred mouse strains to define copy number variable regions (orange). Lilue and colleagues [16] used de novo assembly of 16 mouse strains (red). The gene families supported by multiple studies are named above. N/A indicates no protein coding genes in the region. SSDR, strain-specific diversity region.

https://doi.org/10.1371/journal.pgen.1008446.g001

thumbnail
Download:
Fig 2. A summary of SSDRs in 16 mouse strains.

(A) Proportion of sequence and coding genes in SSDRs for the classical and wild-derived inbred mouse strains. (B) Summary of annotated genes encoded in SSDRs. For gene (families) with known function, only 3.1% have complete sequences (introns and intergenic regions) for multiple mouse strains (green), 9.8% have coding regions from multiple mouse strains (yellow), and all others are based on a single mouse strain (red). (C) Top 10 PANTHER protein classes overrepresented in mouse SSDRs. X-axis indicates times underrepresentation or overrepresentation. Numbers after each protein class indicate corrected FDR(log10 value). CDS, coding sequence; FDR, false discovery rate; SSDR, strain-specific diversity region.

https://doi.org/10.1371/journal.pgen.1008446.g002

Immune-related genes in SSDRs

SSDRs are highly enriched for immunity and infection response related genes (Fig 2C). Examples include MHC, natural killer gene complex [22,23], T-cell receptors [24], and immunoglobulin variable regions [25], which play central roles in non-self-recognition and adaptive immunity. Other loci include oligoadenylate-synthetase 1 (Oas1) complex [26], AIM2-like receptors [27], and Schlafen gene family [28] for virus innate immune response; NOD-like receptors—Nlrp1—for anthrax lethal toxin resistance [29]; immunity-related GTPases (IRGs) for intracellular pathogen resistance [13]; and α and β defensins for immunomodulatory and antimicrobial function in intestinal crypts [30]. An interesting example is the Intelectin (Itln) members encoded on chromosome 1. Intelectin is known to be highly up-regulated in the immune response to parasitic infections, e.g., Trichinella spiralis. In the C57BL/6J reference genome, only one Itln allele can be found; however, BALB/cJ has two (Itln1 and Itlnb) [31], and strain 129S7 has up to 6 Itln alleles (Itln1 to Itln6) [32]. Similarly, the adjacent gene Natural Killer Cell Receptor 2B4 (Cd244) shows similar patterns of CNVs in recent de novo assemblies [16].

Recent de novo assemblies have also highlighted many loci with polymorphisms previously unreported in mice. Apolipoprotein L (APOL) members encoded on chromosome 15 show high levels of CNVs in classical and wild-derived mouse strains [16]. There are very few studies in murine ApoL members; however, their orthologues in humans (APOL1) are polymorphic [33]. Some alleles confer resistance to Trypanosoma brucei brucei in humans but at the same time lead to chronic kidney disease [33]. Skint gene family members, named after ‘skin thickness’ because they regulate epidermal γδ T cells, are associated with chronic wound healing deficiencies in humans [34,35]. A single SNP was reported in mouse strain FVBTac, which causes selective deficiency for epidermal Vγ5+Vδ1+ T cells [35]. However, the polymorphism of the Skint family appears to be much more complex than previously reported among mouse strains. Eosinophil-Associated RNases (Ears) encoded on chromosome 14 are orthologues of human eosinophil-derived neurotoxin (EDN) and eosinophil cationic protein (ECP), which are highly charged cytotoxic proteins released from activated eosinophil granules [34]. Mouse Ears can promote virus clearance [36] and play a role in the Schistosoma resistance [37]. Although evidence of positive selection has been found for Ears members in the reference genome, their diversity between inbred strains is poorly documented. At least three haplotypes can be found in classical inbred mouse strains (haplotype1: C57BL/6J, C57BL/6NJ, 129S1, AKR/J, BALB/cJ, A/J, CBA/J, DBA/2J, C3H/HeJ; haplotype2: NZO/HILtJ,LP/J; haplotype3: FVB/NJ, NOD/ShiLtJ), and four wild-derived strains all carry divergent sequences [16]. Signal-regulatory protein beta 1 members (SIRPB1) are cell surface glycoproteins expressed in leukocytes, which positively regulate neutrophil transepithelial migration [38]. A CNV in SIRPB1 has been reported in humans to be associated with autoimmune thyroid diseases [39] and impulsive-disinhibited personality [40]. In the GRCm38 reference genome, the Sirpb1 locus remains incomplete [16]. However, the de novo assembly of other inbred mouse strains, especially C57BL/6NJ, has partially improved the reference genome and confirmed significant conservation and high diversity across the strains compared to the C57BL/6J haplotype at the Sirpb1 locus [16]. Thus, the newly published draft genomes of multiple mouse strains will further facilitate the use of the house mouse for studying human disease.

Both mice and humans carry very large interferon inducible GTPases (GVIN). Their open reading frame is almost 8,000 base pair in length, encoded by a single colossal exon. They are highly expressed in lymph nodes and whole blood in humans [41] and are inducible by both type I and type II interferons (IFNs) in mice [42]. Although the function of GVIN members remains unknown, they are thought to play a role in pathogen immunity [42]. Alpha-1 antitrypsin (AAT) encoded by gene serine protease inhibitor A1 (SERPINA1) is the most abundant antiprotease in humans. It inhibits neutrophil elastase and regulates serine proteases during acute inflammatory responses, especially in the lungs where it protects the fragile alveolar tissues from proteolytic degradation [43]. Human SERPIN family members are highly polymorphic with 1/2,500 newborns in Western Europe carrying the PiZ or PiS allele that causes acute or chronic lung and liver disease [44]. The trade-off of these adverse alleles is still unclear; however, pathogens are reported to manipulate immunity regulation of host as evasion strategies, and SERPINA is a potential target [45]. Among mouse strains, both the SerpinA and SerpinB gene families are highly polymorphic, and mouse strains with different Serpin haplotype may confer a good model for human AAT diversity.

Many more immune-related genes or gene families are encoded in the SSDRs in mice, e.g., IFNs, guanylate-binding proteins (GBP), Ly6 members, orosomucoids (Orm), paired-Ig-like receptor (Pira/Pirb), interferon-induced proteins with tetratricopeptide repeats (IFIT), and CD200 receptors. A summary of these loci can be found in S1 Table.

Sensory and kin selection

Key to rodent survival is the ability to detect and avoid potentially harmful compounds by smell and taste. The polymorphisms of Tas2r members are believed to match the profiles of bitter chemicals that the mouse population encounter in their diets. Signatures of positive selection have been detected for the human bitter-taste receptor TAS2R16 [46]. The majority of bitter-taste receptors encoded on mouse chromosome 6 are lineage-specific [47]. Indeed, variations in aversion to chemical substances were observed in BXD mice [48] and between mouse strains C3HeB/FeJ and SWR/J [49]. Both phenotypes were mapped to mouse Tas2r loci on chromosome 6 [50].

Olfactory receptors (ORs) are the largest gene superfamily in house mice and most vertebrates [51]. There are 1,296 OR genes distributed in 27 clusters on the Celera mouse genome except chromosome 12 and Y [51]. As one of the most ancient animal senses, olfaction is important to recognise food, identify mates and offspring, and avoid predators or chemical dangers. Polymorphism in ORs in inbred mouse strains is well studied. Multiple OR members from strains 129S1/SvI, 129X1/SvI, 129S6/SvEvTac, A/J, AKR/J, BALB/c, C57BL/6, and DBA/2J were amplified from genomic DNA [51,52] and sequences available [53]. The de novo genome assemblies, especially of wild-derived inbred strains, have greatly boosted the abundance of novel OR genes. Taking strain CAST/EiJ as an example, 1,249 OR candidates have been annotated, 37 of which are not present in the reference mouse genome. In addition, multiple OR pseudogenes in GRCm38 are conserved with CAST/EiJ and vice versa [16]. Strain-specific polymorphisms can be found in 23 OR clusters among 16 mouse strains sequenced. Similar polymorphisms can be found in Taar7d, Taar7e, Taar8a, Taar8b and Taar8c in wild-derived mouse strains. These members are reported as ORs to recognise ethological odors [54]. The human genome encodes 950 OR genes with high diversity, comparable to the mouse genome [55].

The mouse genome contains many other lineage-specific gene family expansions compared to humans. Many of these genes are associated with reproduction, possibly caused by mating competition and kin selection [56]. One remarkable example is of vomeronasal receptors (VRs) that are mainly expressed in the vomeronasal organ and believed to detect pheromones for sexual recognition. Based on structural differences, VRs are classified into two superfamilies, Vmn1r and Vmn2r, and sum to more than 360 members encoded as clusters on multiple chromosomes in GRCm38 reference genome [57]. The dynamic evolution of VRs and the driving force behind it have been largely discussed in the last decades [58–62]. Wynn and colleagues [57] interrogated around 50% VR genes/alleles from 17 inbred mouse strains and found a significantly higher coding sequence variation with nonrandom distribution in the VRs, especially among three house mouse subspecies and between M. Musculus and M. Spretus. These results suggest that VRs may contribute to reproductive isolation between closely related subspecies [57].

As ligands for VRs [63], the major urinary proteins (Mups) are a set of 18–19 kDa communication proteins abundant in mouse urine and other secretions, including lacrimal, parotid, submaxillary, sublingual, preputial, and mammary glands [64,65]. Mups may either directly behave as pheromones or bind small molecule pheromones to stabilize them by a slow-release pattern [66,67]. In the house mouse, Mups are encoded by a gene-dense cluster on chromosome 4 with at least 19 Mup members per haplotype. Wild mice are reported to express complex ‘barcode’ patterns of Mups, which may provide gender, social dominance, and kinship information to other individuals, facilitating inbreeding avoidance and aiding pup identification. However, wild type individual variation on Mup locus is thought to have been lost during derivation of the classical laboratory strains [68]. It was proposed Mup alleles are also highly conserved between individual wild mice [69]. However, previous research based on PCR amplification could not assay novel haplotypes with highly diverse Mup members. De novo assemblies of 16 mouse strains have confirmed the sequence diversity in all four wild-derived strains [16].

Another important group of pheromone proteins are exocrine gland secreting peptides (ESPs). They may regulate mouse social behaviours via VR activation. Esp1 is reported to mediate Bruce effect in mice [70], and Esp22 secreted by juvenile mice may inhibit adult male mating behaviour[71]. ESPs are encoded by a gene cluster close to the class I MHC. In the GRCm38 reference genome, 38 Esp members are annotated, in which 14 appear to be pseudogenes [72,73]. Although most members in Esp family have high sequence diversity between mouse strains, their polymorphisms have not been widely reported. In the human genome, Mups and ESPs are not present, and all except five V1R genes are disrupted by deleterious mutations [74].

Behaviour and neuron development

Sensory receptors may affect the behaviour of the house mouse directly or indirectly, similar to VRs and ORs [75,76]. The modification of behaviour may also be achieved by regulating the development and connection of neuron cells. For example, protocadherin gamma (Pcdhg) members are encoded in a mouse SSDR on chromosome 18. Many Pcdhg members show high polymorphism among mouse strains. Pcdhga genes are found exclusively in vertebrates and predominantly expressed in the nervous system [77]. They may provide a synaptic address code for neuronal connectivity or a single cell barcode for self-recognition and self-avoidance, and their isoform diversity is necessary for postnatal development of neurons [77].

In humans, male and female specific brain function dimorphisms causing mental impairment have been linked to the X chromosome [78]. This is partially related to X-linked lymphocyte-regulated (Xlr) members [79]. Xlr3b and 4b are paternally imprinted in the cortex and other brain regions, which regulate the expression of other genes[79]. Xlr genes are encoded on rapidly evolving gene clusters. Among 16 mouse strains, very few SSDRs can be found on the X chromosome, but the two Xlr loci give strong signatures of CNVs and novel loci in the wild-derived strains.

One of the most remarkable SSDRs is on chromosome 12 (17–25 mega base pairs [mbp]). This 7 Mbp region encodes hippocalcin-like 1 (Hpcal1) and their homologues that belong to the neuronal calcium sensors. Humans have only single copy of HPCAL1, which is mainly expressed in retinal photoreceptors, neurons, and neuroendocrine cells [80]. Knockdown of HPCAL1 in neuroblastoma cells led to impaired neurite outgrowth and inhibited sympathetic neuronal differentiation [81]. In house mice, this gene has been duplicated into 50–100 copies. The mouse Hpcal1 complex is extremely repetitive containing several recent duplications of hundreds of kilobases. Current draft de novo assemblies do not accurately represent the Hpcal1 locus in any mouse strain, although it appears that in 12 classical strains at least eight different haplotypes can be observed, and four wild-derived strains contain a further four [16]. To date, the function of Hpcal1 homologs and the purpose of the rapid expansion of the loci remains unknown. Further candidate genes that potentially have functions in neuron development and regulation include Mas and related G Protein-Coupled Receptors (Mrgpra), Angiopoietins (Ang), and Neuronal apoptosis inhibitory proteins (Naip) [82–84].

Sexual reproduction and other biology processes

Sexual reproduction is a complex process from gamete recognition to maternal-fetal interaction. Many genes related to sperm-egg interaction show positive selection and polymorphism, which may reflect the evolutionary pressure from species recognition or inbreeding avoidance [85]. The a disintegrin and metalloprotease (Adam) gene family are important sperm surface proteins. Rapid evolution can be found within their sperm-egg adhesion domains [86]. Three Adam members are found in a SSDR, namely Adam20, Adam25, and Adam26a. The divergence can be only observed in M. Spretus, which indicates a potential role in hybridization avoidance [16]. Other sperm specific gene families, however, are polymorphic among classical laboratory mice. Sperm-associated glutamate (E)-rich protein (Speer) members are encoded in a gene-dense cluster on chromosome 5. At least three of them are expressed solely in the adult mouse testis. Speer homologs are not present in most other mammal species including humans [87], and their function remains unclear. Female specific genes can be also found in the SSDRs. Pregnancy-specific glycoproteins (Psg) are members of immunoglobulin superfamily. In humans, PSGs may be the most abundant trophoblastic proteins in maternal blood during pregnancy [88]. Human PSGs play an essential role in the regulation of maternal immunity, by protecting a fetus from immune responses in case of infection, inflammation, and trauma [89]. The polymorphism of Psg members in mice is possibly caused by a combination of immune tolerance and host–pathogen coevolution. Four haplotypes of the Psg complex can be found in classical inbred strains, and all wild-derived mouse strains have novel haplotypes.

Many other candidates in the list of mouse strain-specific diversity genes have various function or unknown function (see S1 Table). Variations in keratins and keratin associated proteins (Krtap) may affect the hair content characteristics of mouse individuals [90]. Polymorphisms of Hydroxysteroid sulfotransferase enzymes (Sult) may reflect challenge from chemical Metabolism. Variation of zinc finger proteins (Zfp) are thought to repress transposable elements in an evolutionary arms race [91].

Conclusions

Isogenic inbred mice have held a unique position as the key mammalian model in evolutionary, genetics, genomics, and biomedical research for over a century. Sequencing and functional studies have documented the extent of genetic polymorphism residing amongst the strains, both shared and unique to each strain. Genetic variation between mouse strains is not evenly distributed across the genome. In most regions, mice are >99.5% identical, but in SSDRs (around 0.5%–2.8% of the mouse genome), the difference is often higher than interspecies diversity between mouse and rat [13]. This scale of diversity cannot be easily represented using the reference genome with SNPs, indels, and SVs. SSDRs are overrepresented with genes associated with immunity, sensory, sexual reproduction, and behavioral phenotypes [16]. The selective pressures driving diversity and CNV includes host–pathogen coevolution (e.g., red queen hypothesis) [92], kin selection [93], mating preference [94], and even selective sweeps due to strong positive selection [95,96]. Many of these genes have direct orthologues in the human genome and are therefore important for understanding health and disease, drug development, and vaccine development. Multiple well-annotated reference genomes will allow researchers to use the appropriate strain for biological rather than historical reasons.

While this review has focused primarily on the limitations of our knowledge of diversity in protein coding regions of the mouse genome, there are other functional elements in which our knowledge is even more limited, e.g., long noncoding RNAs (ncRNAs); piRNAs; and transcription controlling elements such as promoters, enhancers, silencers, and insulators. Multiple reference quality chromosome sequences will provide the foundation for future mapping studies to interrogate these elements. The dramatic drop in second-generation sequencing costs has resulted in genome-wide catalogs of genetic variants for hundreds of mouse strains, but the process of producing a reference quality genome sequence that includes fully resolved novel haplotypes remains costly. Recent advances in third generation sequencing platforms, such as Pacific Biosciences and Oxford Nanopore, can produce mammalian genomes that are an order of magnitude more contiguous[97]. We expect that the representation of many SSDRs in mouse strains will be greatly improved by third generation sequencing platforms.

Human genome-wide association studies (GWAS) have discovered many loci associated with complex disease and traits. Knowledge from model organisms, combined with fine mapping techniques and functional studies, are used to identify causative genes and mechanisms. Mouse SSDR regions are enriched for genes with disease functions with known orthologs in the human genome. The completion of the mouse pan-genome that incorporates all known genetic variants and novel haplotypes will enable the functional characterization of many unresolved quantitative trait loci (QTLs) associated with human disease.

One interesting question is what the origins of these highly diverse haplotypes in the mouse genome are. To date, only a few of these loci have been studied in detail in both inbred and wild mice. Trachtulec and colleagues [98] constructed a haplotype map of the Hst1 region and H2 haplotypes for five mouse subspecies and found that trans-species SNPs were rare, concluding that the haplotypes are unlikely to have arisen by recombination during inbreeding. Lilue and colleagues[13] studied the polymorphic alleles of IRG proteins in inbred laboratory mice that have also been found in European wild mice, suggesting that these alleles arose prior to inbreeding, whilst other more ancient alleles are shared across mouse subspecies. The combination of multiple reference quality genomes for the primary mouse subspecies and availability of larger numbers of sequenced wild mice from ancestral populations will enable a comprehensive analysis of the origins of all SSDRs.

Supporting information

S1 Table. The gene families, publication identifiers, human orthologs, and mouse gene names for the SSDR regions in the mouse genome.

SSDR, strain-specific diversity region.

https://doi.org/10.1371/journal.pgen.1008446.s001

(XLSX)

S1 Data. Coordinates on GRCm38 for the SSDR regions per strain (BED format).

SSDR, strain-specific diversity region.

https://doi.org/10.1371/journal.pgen.1008446.s002

(GZ)

References

  1. 1. Beck JA, Lloyd S, Hafezparast M, Lennon-Pierce M, Eppig JT, Festing MF, et al. Genealogies of mouse inbred strains. Nat Genet. 2000;24: 23–25. pmid:10615122
  2. 2. Taft RA, Davisson M, Wiles MV. Know thy mouse. Trends Genet TIG. 2006;22: 649–653. pmid:17007958
  3. 3. Mouse Genome Sequencing Consortium, Waterston RH, Lindblad-Toh K, Birney E, Rogers J, Abril JF, et al. Initial sequencing and comparative analysis of the mouse genome. Nature. 2002;420: 520–562. pmid:12466850
  4. 4. Keane TM, Goodstadt L, Danecek P, White MA, Wong K, Yalcin B, et al. Mouse genomic variation and its effect on phenotypes and gene regulation. Nature. 2011;477: 289–294. pmid:21921910
  5. 5. Yalcin B, Wong K, Agam A, Goodson M, Keane TM, Gan X, et al. Sequence-based characterization of structural variation in the mouse genome. Nature. 2011;477: 326–329. pmid:21921916
  6. 6. van der Weyden L, Adams DJ, Bradley A. Tools for targeted manipulation of the mouse genome. Physiol Genomics. 2002;11: 133–164. pmid:12464689
  7. 7. Yang H, Bell TA, Churchill GA, Pardo-Manuel de Villena F. On the subspecific origin of the laboratory mouse. Nat Genet. 2007;39: 1100–1107. pmid:17660819
  8. 8. Yang H, Wang JR, Didion JP, Buus RJ, Bell TA, Welsh CE, et al. Subspecific origin and haplotype diversity in the laboratory mouse. Nat Genet. 2011;43: 648–655. pmid:21623374
  9. 9. Fischer Lindahl K. On naming H2 haplotypes: functional significance of MHC class Ib alleles. Immunogenetics. 1997;46: 53–62. pmid:9148789
  10. 10. Flaherty L, Elliott E, Tine JA, Walsh AC, Waters JB. Immunogenetics of the Q and TL regions of the mouse. Crit Rev Immunol. 1990;10: 131–175. pmid:2076187
  11. 11. Guénet JL, Bonhomme F. Wild mice: an ever-increasing contribution to a popular mammalian model. Trends Genet TIG. 2003;19: 24–31. pmid:12493245
  12. 12. Hassan MA, Olijnik A-A, Frickel E-M, Saeij JP. Clonal and atypical Toxoplasma strain differences in virulence vary with mouse sub-species. Int J Parasitol. 2019;49: 63–70. pmid:30471286
  13. 13. Lilue J, Müller UB, Steinfeldt T, Howard JC. Reciprocal virulence and resistance polymorphism in the relationship between Toxoplasma gondii and the house mouse. eLife. 2013;2: e01298. pmid:24175088
  14. 14. Song Y, Endepols S, Klemann N, Richter D, Matuschka F-R, Shih C-H, et al. Adaptive introgression of anticoagulant rodent poison resistance by hybridization between old world mice. Curr Biol CB. 2011;21: 1296–1301. pmid:21782438
  15. 15. Bagot S, Campino S, Penha-Gonçalves C, Pied S, Cazenave P-A, Holmberg D. Identification of two cerebral malaria resistance loci using an inbred wild-derived mouse strain. Proc Natl Acad Sci U S A. 2002;99: 9919–9923. pmid:12114535
  16. 16. Lilue J, Doran AG, Fiddes IT, Abrudan M, Armstrong J, Bennett R, et al. Sixteen diverse laboratory mouse reference genomes define strain-specific haplotypes and novel functional loci. Nat Genet. 2018;50: 1574–1583. pmid:30275530
  17. 17. Doran AG, Wong K, Flint J, Adams DJ, Hunter KW, Keane TM. Deep genome sequencing and variation analysis of 13 inbred mouse strains defines candidate phenotypic alleles, private variation and homozygous truncating mutations. Genome Biol. 2016;17: 167. pmid:27480531
  18. 18. Li J, Jiang T, Mao J-H, Balmain A, Peterson L, Harris C, et al. Genomic segmental polymorphisms in inbred mouse strains. Nat Genet. 2004;36: 952–954. pmid:15322544
  19. 19. Cutler G, Marshall LA, Chin N, Baribault H, Kassner PD. Significant gene content variation characterizes the genomes of inbred mouse strains. Genome Res. 2007;17: 1743–1754. pmid:17989247
  20. 20. Locke MEO, Milojevic M, Eitutis ST, Patel N, Wishart AE, Daley M, et al. Genomic copy number variation in Mus musculus. BMC Genomics. 2015;16: 497. pmid:26141061
  21. 21. Morgan AP, Holt JM, McMullan RC, Bell TA, Clayshulte AM-F, Didion JP, et al. The evolutionary fates of a large segmental duplication in mouse [Internet]. Genetics; 2016 Mar.
  22. 22. Carlyle JR, Mesci A, Fine JH, Chen P, Bélanger S, Tai L-H, et al. Evolution of the Ly49 and Nkrp1 recognition systems. Semin Immunol. 2008;20: 321–330. pmid:18595730
  23. 23. Brown MG, Scalzo AA. NK gene complex dynamics and selection for NK cell receptors. Semin Immunol. 2008;20: 361–368. pmid:18640056
  24. 24. Nobuhara H, Kuida K, Furutani M, Shiroishi T, Moriwaki K, Yanagi Y, et al. Polymorphism of T-cell receptor genes among laboratory and wild mice: diverse origins of laboratory mice. Immunogenetics. 1989;30: 405–413. pmid:2574156
  25. 25. Barstad P, Farnsworth V, Weigert M, Cohn M, Hood L. Mouse immunoglobulin heavy chains are coded by multiple germ line variable region genes. Proc Natl Acad Sci U S A. 1974;71: 4096–4100. pmid:4215076
  26. 26. Green R, Wilkins C, Thomas S, Sekine A, Hendrick DM, Voss K, et al. Oas1b-dependent Immune Transcriptional Profiles of West Nile Virus Infection in the Collaborative Cross. G3 Bethesda Md. 2017;7: 1665–1682. pmid:28592649
  27. 27. Nakaya Y, Lilue J, Stavrou S, Moran EA, Ross SR. AIM2-Like Receptors Positively and Negatively Regulate the Interferon Response Induced by Cytosolic DNA. mBio. 2017;8. pmid:28679751
  28. 28. Mavrommatis E, Fish EN, Platanias LC. The schlafen family of proteins and their regulation by interferons. J Interferon Cytokine Res Off J Int Soc Interferon Cytokine Res. 2013;33: 206–210. pmid:23570387
  29. 29. Sastalla I, Crown D, Masters SL, McKenzie A, Leppla SH, Moayeri M. Transcriptional analysis of the three Nlrp1 paralogs in mice. BMC Genomics. 2013;14: 188. pmid:23506131
  30. 30. Shanahan MT, Tanabe H, Ouellette AJ. Strain-specific polymorphisms in Paneth cell α-defensins of C57BL/6 mice and evidence of vestigial myeloid α-defensin pseudogenes. Infect Immun. 2011;79: 459–473. pmid:21041494
  31. 31. Pemberton AD, Knight PA, Gamble J, Colledge WH, Lee J-K, Pierce M, et al. Innate BALB/c enteric epithelial responses to Trichinella spiralis: inducible expression of a novel goblet cell lectin, intelectin-2, and its natural deletion in C57BL/10 mice. J Immunol Baltim Md 1950. 2004;173: 1894–1901. pmid:15265922
  32. 32. Lu ZH, di Domenico A, Wright SH, Knight PA, Whitelaw CBA, Pemberton AD. Strain-specific copy number variation in the intelectin locus on the 129 mouse chromosome 1. BMC Genomics. 2011;12: 110. pmid:21324158
  33. 33. Vanhamme L, Paturiaux-Hanocq F, Poelvoorde P, Nolan DP, Lins L, Van Den Abbeele J, et al. Apolipoprotein L-I is the trypanosome lytic factor of human serum. Nature. 2003;422: 83–87. pmid:12621437
  34. 34. Barbee SD, Woodward MJ, Turchinovich G, Mention J-J, Lewis JM, Boyden LM, et al. Skint-1 is a highly specific, unique selecting component for epidermal T cells. Proc Natl Acad Sci U S A. 2011;108: 3330–3335. pmid:21300860
  35. 35. Boyden LM, Lewis JM, Barbee SD, Bas A, Girardi M, Hayday AC, et al. Skint1, the prototype of a newly identified immunoglobulin superfamily gene cluster, positively selects epidermal gammadelta T cells. Nat Genet. 2008;40: 656–662. pmid:18408721
  36. 36. Percopo CM, Dyer KD, Ochkur SI, Luo JL, Fischer ER, Lee JJ, et al. Activated mouse eosinophils protect against lethal respiratory virus infection. Blood. 2014;123: 743–752. pmid:24297871
  37. 37. Nitto T, Dyer KD, Mejia RA, Byström J, Wynn TA, Rosenberg HF. Characterization of the divergent eosinophil ribonuclease, mEar 6, and its expression in response to Schistosoma mansoni infection in vivo. Genes Immun. 2004;5: 668–674. pmid:15526002
  38. 38. Liu Y, Soto I, Tong Q, Chin A, Bühring H-J, Wu T, et al. SIRPbeta1 is expressed as a disulfide-linked homodimer in leukocytes and positively regulates neutrophil transepithelial migration. J Biol Chem. 2005;280: 36132–36140. pmid:16081415
  39. 39. Jin X, Guan Y, Shen H, Pang Y, Liu L, Jia Q, et al. Copy Number Variation of Immune-Related Genes and Their Association with Iodine in Adults with Autoimmune Thyroid Diseases. Int J Endocrinol. 2018;2018: 1705478. pmid:29713342
  40. 40. Kajimoto N, Kirpekar SM, Wakade AR. An investigation of spontaneous potentials recorded from the smooth-muscle cells of the guinea-pig seminal vesicle. J Physiol. 1972;224: 105–119. pmid:5039969
  41. 41. Thierry-Mieg D, Thierry-Mieg J. AceView: a comprehensive cDNA-supported gene and transcripts annotation. Genome Biol. 2006;7 Suppl 1: S12.1–14. pmid:16925834
  42. 42. Klamp T, Boehm U, Schenk D, Pfeffer K, Howard JC. A giant GTPase, very large inducible GTPase-1, is inducible by IFNs. J Immunol Baltim Md 1950. 2003;171: 1255–1265. pmid:12874213
  43. 43. Bergin DA, Hurley K, McElvaney NG, Reeves EP. Alpha-1 antitrypsin: a potent anti-inflammatory and potential novel therapeutic agent. Arch Immunol Ther Exp (Warsz). 2012;60: 81–97. pmid:22349104
  44. 44. Fregonese L, Stolk J. Hereditary alpha-1-antitrypsin deficiency and its clinical consequences. Orphanet J Rare Dis. 2008;3: 16. pmid:18565211
  45. 45. Odendall C, Kagan JC. Activation and pathogenic manipulation of the sensors of the innate immune system. Microbes Infect. 2017;19: 229–237. pmid:28093320
  46. 46. Soranzo N, Bufe B, Sabeti PC, Wilson JF, Weale ME, Marguerie R, et al. Positive selection on a high-sensitivity allele of the human bitter-taste receptor TAS2R16. Curr Biol CB. 2005;15: 1257–1265. pmid:16051168
  47. 47. Lossow K, Hübner S, Roudnitzky N, Slack JP, Pollastro F, Behrens M, et al. Comprehensive Analysis of Mouse Bitter Taste Receptors Reveals Different Molecular Receptive Ranges for Orthologous Receptors in Mice and Humans. J Biol Chem. 2016;291: 15358–15377. pmid:27226572
  48. 48. Boughter JD, Raghow S, Nelson TM, Munger SD. Inbred mouse strains C57BL/6J and DBA/2J vary in sensitivity to a subset of bitter stimuli. BMC Genet. 2005;6: 36. pmid:15967025
  49. 49. Nelson TM, Munger SD, Boughter JD. Taste sensitivities to PROP and PTC vary independently in mice. Chem Senses. 2003;28: 695–704. pmid:14627538
  50. 50. Bachmanov AA, Bosak NP, Lin C, Matsumoto I, Ohmoto M, Reed DR, et al. Genetics of taste receptors. Curr Pharm Des. 2014;20: 2669–2683. pmid:23886383
  51. 51. Zhang X, Firestein S. The olfactory receptor gene superfamily of the mouse. Nat Neurosci. 2002;5: 124–133. pmid:11802173
  52. 52. Young JM, Friedman C, Williams EM, Ross JA, Tonnes-Priddy L, Trask BJ. Different evolutionary processes shaped the mouse and human olfactory receptor gene families. Hum Mol Genet. 2002;11: 535–546. pmid:11875048
  53. 53. Crasto C, Marenco L, Miller P, Shepherd G. Olfactory Receptor Database: a metadata-driven automated population from sources of gene and protein sequences. Nucleic Acids Res. 2002;30: 354–360. pmid:11752336
  54. 54. Ferrero DM, Wacker D, Roque MA, Baldwin MW, Stevens RC, Liberles SD. Agonists for 13 trace amine-associated receptors provide insight into the molecular basis of odor selectivity. ACS Chem Biol. 2012;7: 1184–1189. pmid:22545963
  55. 55. Young JM, Trask BJ. The sense of smell: genomics of vertebrate odorant receptors. Hum Mol Genet. 2002;11: 1153–1160. pmid:12015274
  56. 56. Silva L, Antunes A. Vomeronasal Receptors in Vertebrates and the Evolution of Pheromone Detection. Annu Rev Anim Biosci. 2017;5: 353–370. pmid:27912243
  57. 57. Wynn EH, Sánchez-Andrade G, Carss KJ, Logan DW. Genomic variation in the vomeronasal receptor gene repertoires of inbred mice. BMC Genomics. 2012;13: 415. pmid:22908939
  58. 58. Yoder AD, Larsen PA. The molecular evolutionary dynamics of the vomeronasal receptor (class 1) genes in primates: a gene family on the verge of a functional breakdown. Front Neuroanat. 2014;8: 153. pmid:25565978
  59. 59. Emes RD, Beatson SA, Ponting CP, Goodstadt L. Evolution and comparative genomics of odorant- and pheromone-associated genes in rodents. Genome Res. 2004;14: 591–602. pmid:15060000
  60. 60. Lane RP, Young J, Newman T, Trask BJ. Species specificity in rodent pheromone receptor repertoires. Genome Res. 2004;14: 603–608. pmid:15060001
  61. 61. Grus WE, Zhang J. Rapid turnover and species-specificity of vomeronasal pheromone receptor genes in mice and rats. Gene. 2004;340: 303–312. pmid:15475172
  62. 62. Park SH, Podlaha O, Grus WE, Zhang J. The microevolution of V1r vomeronasal receptor genes in mice. Genome Biol Evol. 2011;3: 401–412. pmid:21551350
  63. 63. Krieger J, Schmitt A, Löbel D, Gudermann T, Schultz G, Breer H, et al. Selective activation of G protein subtypes in the vomeronasal organ upon stimulation with urine-derived compounds. J Biol Chem. 1999;274: 4655–4662. pmid:9988702
  64. 64. Gubits RM, Lynch KR, Kulkarni AB, Dolan KP, Gresik EW, Hollander P, et al. Differential regulation of alpha 2u globulin gene expression in liver, lachrymal gland, and salivary gland. J Biol Chem. 1984;259: 12803–12809. pmid:6208189
  65. 65. Shahan K, Denaro M, Gilmartin M, Shi Y, Derman E. Expression of six mouse major urinary protein genes in the mammary, parotid, sublingual, submaxillary, and lachrymal glands and in the liver. Mol Cell Biol. 1987;7: 1947–1954. pmid:3600653
  66. 66. Hurst null, Robertson null, Tolladay null, Beynon null. Proteins in urine scent marks of male house mice extend the longevity of olfactory signals. Anim Behav. 1998;55: 1289–1297. pmid:9632512
  67. 67. Chamero P, Marton TF, Logan DW, Flanagan K, Cruz JR, Saghatelian A, et al. Identification of protein pheromones that promote aggressive behaviour. Nature. 2007;450: 899–902. pmid:18064011
  68. 68. Cheetham SA, Smith AL, Armstrong SD, Beynon RJ, Hurst JL. Limited variation in the major urinary proteins of laboratory mice. Physiol Behav. 2009;96: 253–261. pmid:18973768
  69. 69. Thoß M, Enk V, Yu H, Miller I, Luzynski KC, Balint B, et al. Diversity of major urinary proteins (MUPs) in wild house mice. Sci Rep. 2016;6: 38378. pmid:27922085
  70. 70. Hattori T, Osakada T, Masaoka T, Ooyama R, Horio N, Mogi K, et al. Exocrine Gland-Secreting Peptide 1 Is a Key Chemosensory Signal Responsible for the Bruce Effect in Mice. Curr Biol CB. 2017;27: 3197–3201.e3. pmid:29033330
  71. 71. Ferrero DM, Moeller LM, Osakada T, Horio N, Li Q, Roy DS, et al. A juvenile mouse pheromone inhibits sexual behaviour through the vomeronasal system. Nature. 2013;502: 368–371. pmid:24089208
  72. 72. Kimoto H, Sato K, Nodari F, Haga S, Holy TE, Touhara K. Sex- and strain-specific expression and vomeronasal activity of mouse ESP family peptides. Curr Biol CB. 2007;17: 1879–1884. pmid:17935991
  73. 73. Kimoto H, Haga S, Sato K, Touhara K. Sex-specific peptides from exocrine glands stimulate mouse vomeronasal sensory neurons. Nature. 2005;437: 898–901. pmid:16208374
  74. 74. Young JM, Massa HF, Hsu L, Trask BJ. Extreme variability among mammalian V1R gene families. Genome Res. 2010;20: 10–18. pmid:19952141
  75. 75. Ibarra-Soria X, Levitin MO, Logan DW. The genomic basis of vomeronasal-mediated behaviour. Mamm Genome Off J Int Mamm Genome Soc. 2014;25: 75–86. pmid:23884334
  76. 76. Glinka ME, Samuels BA, Diodato A, Teillon J, Feng Mei D, Shykind BM, et al. Olfactory deficits cause anxiety-like behaviors in mice. J Neurosci Off J Soc Neurosci. 2012;32: 6718–6725. pmid:22573694
  77. 77. Chen WV, Alvarez FJ, Lefebvre JL, Friedman B, Nwakeze C, Geiman E, et al. Functional significance of isoform diversification in the protocadherin gamma gene cluster. Neuron. 2012;75: 402–409. pmid:22884324
  78. 78. Zechner U, Wilda M, Kehrer-Sawatzki H, Vogel W, Fundele R, Hameister H. A high density of X-linked genes for general cognitive ability: a run-away process shaping human evolution? Trends Genet TIG. 2001;17: 697–701. pmid:11718922
  79. 79. Davies W, Isles A, Smith R, Karunadasa D, Burrmann D, Humby T, et al. Xlr3b is a new imprinted candidate for X-linked parent-of-origin effects on cognitive function in mice. Nat Genet. 2005;37: 625–629. pmid:15908950
  80. 80. Burgoyne RD. Neuronal calcium sensor proteins: generating diversity in neuronal Ca2+ signalling. Nat Rev Neurosci. 2007;8: 182–193. pmid:17311005
  81. 81. Wang W, Zhong Q, Teng L, Bhatnagar N, Sharma B, Zhang X, et al. Mutations that disrupt PHOXB interaction with the neuronal calcium sensor HPCAL1 impede cellular differentiation in neuroblastoma. Oncogene. 2014;33: 3316–3324. pmid:23873030
  82. 82. Geppetti P, Veldhuis NA, Lieu T, Bunnett NW. G Protein-Coupled Receptors: Dynamic Machines for Signaling Pain and Itch. Neuron. 2015;88: 635–649. pmid:26590341
  83. 83. Subramanian V, Crabtree B, Acharya KR. Human angiogenin is a neuroprotective factor and amyotrophic lateral sclerosis associated angiogenin variants affect neurite extension/pathfinding and survival of motor neurons. Hum Mol Genet. 2008;17: 130–149. pmid:17916583
  84. 84. Götz R, Karch C, Digby MR, Troppmair J, Rapp UR, Sendtner M. The neuronal apoptosis inhibitory protein suppresses neuronal differentiation and apoptosis in PC12 cells. Hum Mol Genet. 2000;9: 2479–2489. pmid:11030753
  85. 85. Firman RC, Gasparini C, Manier MK, Pizzari T. Postmating Female Control: 20 Years of Cryptic Female Choice. Trends Ecol Evol. 2017;32: 368–382. pmid:28318651
  86. 86. Civetta A. Positive selection within sperm-egg adhesion domains of fertilin: an ADAM gene with a potential role in fertilization. Mol Biol Evol. 2003;20: 21–29. pmid:12519902
  87. 87. Spiess A-N, Walther N, Müller N, Balvers M, Hansis C, Ivell R. SPEER—a new family of testis-specific genes from the mouse. Biol Reprod. 2003;68: 2044–2054. pmid:12606357
  88. 88. Moore T, Dveksler GS. Pregnancy-specific glycoproteins: complex gene families regulating maternal-fetal interactions. Int J Dev Biol. 2014;58: 273–280. pmid:25023693
  89. 89. Motrán CC, Díaz FL, Gruppi A, Slavin D, Chatton B, Bocco JL. Human pregnancy-specific glycoprotein 1a (PSG1a) induces alternative activation in human and mouse monocytes and suppresses the accessory cell-dependent T cell proliferation. J Leukoc Biol. 2002;72: 512–521. pmid:12223519
  90. 90. Wu D-D, Irwin DM, Zhang Y-P. Molecular evolution of the keratin associated protein gene family in mammals, role in the evolution of mammalian hair. BMC Evol Biol. 2008;8: 241. pmid:18721477
  91. 91. Imbeault M, Helleboid P-Y, Trono D. KRAB zinc-finger proteins contribute to the evolution of gene regulatory networks. Nature. 2017;543: 550–554. pmid:28273063
  92. 92. Morran LT, Schmidt OG, Gelarden IA, Parrish RC, Lively CM. Running with the Red Queen: Host-Parasite Coevolution Selects for Biparental Sex. Science. 2011;333: 216–218. pmid:21737739
  93. 93. Axelrod R, Hammond RA, Grafen A. Altruism via kin-selection strategies that rely on arbitrary tags with which they coevolve. Evol Int J Org Evol. 2004;58: 1833–1838.
  94. 94. Sherborne AL, Thom MD, Paterson S, Jury F, Ollier WER, Stockley P, et al. The genetic basis of inbreeding avoidance in house mice. Curr Biol CB. 2007;17: 2061–2066. pmid:17997307
  95. 95. Didion JP, Morgan AP, Clayshulte AM-F, Mcmullan RC, Yadgary L, Petkov PM, et al. A multi-megabase copy number gain causes maternal transmission ratio distortion on mouse chromosome 2. PLoS Genet. 2015;11: e1004850. pmid:25679959
  96. 96. Didion JP, Morgan AP, Yadgary L, Bell TA, McMullan RC, Ortiz de Solorzano L, et al. R2d2 Drives Selfish Sweeps in the House Mouse. Mol Biol Evol. 2016;33: 1381–1395. pmid:26882987
  97. 97. Gordon D, Huddleston J, Chaisson MJP, Hill CM, Kronenberg ZN, Munson KM, et al. Long-read sequence assembly of the gorilla genome. Science. 2016;352: aae0344. pmid:27034376
  98. 98. Trachtulec Z, Vlcek C, Mihola O, Gregorova S, Fotopulosova V, Forejt J. Fine Haplotype Structure of a Chromosome 17 Region in the Laboratory and Wild Mouse. Genetics. 2008;178: 1777–1784. pmid:18245833