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.

GenomeMUSter mouse genetic variation service enables multi-trait, multi-population data integration and analyses.

Hundreds of inbred laboratory mouse strains and intercross populations have been used to functionalize 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 the Collaborative Cross, expanded BXD and inbred wild-derived strains add to set of complex disease mouse models, genetic mapping resources and sensitized backgrounds against which to evaluate engineered mutations. The genome sequences of many inbred strains, along with dense genotypes from others could allow integrated analysis of trait - variant associations across populations, but these analyses are not feasible due to the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense data resource by harmonizing multiple variant datasets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extensible to other model organism species. The result is a web- and programmatically-accessible data service called GenomeMUSter ( https://muster.jax.org ), comprising allelic data covering 657 strains at 106.8M segregating sites. Interoperation with phenotype databases, analytic tools and other resources enable a wealth of applications including multi-trait, multi-population meta-analysis. We demonstrate this in a cross-species comparison of the meta-analysis of Type 2 Diabetes and of substance use disorders, resulting in the more specific characterization of the role of human variant effects in light of mouse phenotype data. 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.

Keratinocyte-specific deletion of SHARPIN induces atopic dermatitis-like inflammation in mice.

Spontaneous mutations in the SHANK-associated RH domain interacting protein (Sharpin) resulted in a severe autoinflammatory type of chronic proliferative dermatitis, inflammation in other organs, and lymphoid organ defects. To determine whether cell-type restricted loss of Sharpin causes similar lesions, a conditional null mutant was created. Ubiquitously expressing cre-recombinase recapitulated the phenotype seen in spontaneous mutant mice. Limiting expression to keratinocytes (using a Krt14-cre) induced a chronic eosinophilic dermatitis, but no inflammation in other organs or lymphoid organ defects. The dermatitis was associated with a markedly increased concentration of serum IgE and IL18. Crosses with S100a4-cre resulted in milder skin lesions and moderate to severe arthritis. This conditional null mutant will enable more detailed studies on the role of SHARPIN in regulating NFkB and inflammation, while the Krt14-Sharpin-/- provides a new model to study atopic dermatitis.

2-deoxy D-glucose treatment does not elicit a hair growth response in alopecia areata.

2-deoxy D-glucose (2DG) was tested for efficacy in treating alopecia areata using the C3H/HeJ skin graft model. 2DG has proven to be efficacious in treatment of various mouse models of autoimmunity with minimal serious side effects noted. This agent has been shown to normalize abnormally activated T-cell populations while also preventing cell surface expression of NKG2D; key factors defining alopecia areata disease progression. Daily oral ingestion of 2DG via drinking water to mice with patchy or diffuse alopecia areata for 16 weeks failed to prevent expansion of alopecia or cause regrowth of hair in treated mice. Histologically, there were no differences between treated and control groups. These results indicate that, while 2DG is effective for some autoimmune diseases, it was not efficacious for the cell-mediated autoimmune mouse disease, alopecia areata.

Gain of function p.E138A alteration in Card14 leads to psoriasiform skin inflammation and implicates genetic modifiers in disease severity.

Psoriasis (PS) is a common inflammatory and incurable skin disease affecting 2-3% of the human population. Although genome-wide association studies implicate more than 60 loci, the full complement of genetic factors leading to disease is not known. Rare, highly penetrant, gain-of-function, dominantly acting mutations within the human caspase recruitment domain family, member 14 (CARD14) gene lead to the development of PS and psoriatic arthritis (PSA) (a familial p.G117S and de-novo p.E138A alteration). These residues are conserved in mouse and orthologous Knock-In (KI) mutations within Card14 were created. The Card14tm.1.1Sun allele (G117S) resulted in no clinically or histologically evident phenotype of the skin or joints in young adult or old mice. However, mice carrying the Card14tm2.1Sun mutant allele (E138A) were runted and developed thick, white, scaly skin soon after birth, dying within two weeks or less. The skin hyperplasia and inflammation was remarkable similarity to human PS at the clinical, histological, and transcriptomic levels. For example, the skin was markedly acanthotic and exhibited orthokeratotic hyperkeratosis with minimal inflammation and no pustules and transcripts affecting critical pathways of epidermal differentiation and components of the IL17 axis (IL23, IL17A, IL17C, TNF and IL22) were altered. Similar changes were seen in a set of orthologous microRNAs previously associated with PS suggesting conservation across species. Crossing the Card14tm2.1Sun/WT mice to C57BL/6NJ, FVB/NJ, CBA/J, C3H/HeJ, and 129S1/SvImJ generated progeny with epidermal acanthosis and marked orthokeratotic hyperkeratosis regardless of the hybrid strain. Of these hybrid lines, only the FVB;B6N(129S4) mice survived to 250 days of age or older and has led to recombinant inbred lines homozygous for Card14E138A that are fecund and have scaly skin disease. This implicates that modifiers of PS severity exist in mice, as in the familial forms of the disease in humans.

Quantitative Trait Locus and Integrative Genomics Revealed Candidate Modifier Genes for Ectopic Mineralization in Mouse Models of Pseudoxanthoma Elasticum.

Pseudoxanthoma elasticum, a prototype of heritable multisystem ectopic mineralization disorders, is caused by mutations in the ABCC6 gene encoding a putative efflux transporter, ABCC6. The phenotypic spectrum of pseudoxanthoma elasticum varies, and the correlation between genotype and phenotype has not been established. To identify genetic modifiers, we performed quantitative trait locus analysis in inbred mouse strains that carry the same hypomorphic allele in Abcc6 yet with highly variable ectopic mineralization phenotypes of pseudoxanthoma elasticum. Abcc6 was confirmed as a major determinant for ectopic mineralization in multiple tissues. Integrative analysis using functional genomics tools that included GeneWeaver, String, and Mouse Genome Informatics identified a total of nine additional candidate modifier genes that could influence the organ-specific ectopic mineralization phenotypes. Integration of the candidate genes into the existing ectopic mineralization gene network expands the current knowledge on the complexity of the network that, as a whole, governs ectopic mineralization in soft connective tissues.

Hyaline Arteriolosclerosis in 30 Strains of Aged Inbred Mice.

During a screen for vascular phenotypes in aged laboratory mice, a unique discrete phenotype of hyaline arteriolosclerosis of the intertubular arteries and arterioles of the testes was identified in several inbred strains. Lesions were limited to the testes and did not occur as part of any renal, systemic, or pulmonary arteriopathy or vasculitis phenotype. There was no evidence of systemic or pulmonary hypertension, and lesions did not occur in ovaries of females. Frequency was highest in males of the SM/J (27/30, 90%) and WSB/EiJ (19/26, 73%) strains, aged 383 to 847 days. Lesions were sporadically present in males from several other inbred strains at a much lower (<20%) frequency. The risk of testicular hyaline arteriolosclerosis is at least partially underpinned by a genetic predisposition that is not associated with other vascular lesions (including vasculitis), separating out the etiology of this form and site of arteriolosclerosis from other related conditions that often co-occur in other strains of mice and in humans. Because of their genetic uniformity and controlled dietary and environmental conditions, mice are an excellent model to dissect the pathogenesis of human disease conditions. In this study, a discrete genetically driven phenotype of testicular hyaline arteriolosclerosis in aging mice was identified. These observations open the possibility of identifying the underlying genetic variant(s) associated with the predisposition and therefore allowing future interrogation of the pathogenesis of this condition.

Quantitative evaluation of ontology design patterns for combining pathology and anatomy ontologies.

Data are increasingly annotated with multiple ontologies to capture rich information about the features of the subject under investigation. Analysis may be performed over each ontology separately, but recently there has been a move to combine multiple ontologies to provide more powerful analytical possibilities. However, it is often not clear how to combine ontologies or how to assess or evaluate the potential design patterns available. Here we use a large and well-characterized dataset of anatomic pathology descriptions from a major study of aging mice. We show how different design patterns based on the MPATH and MA ontologies provide orthogonal axes of analysis, and perform differently in over-representation and semantic similarity applications. We discuss how such a data-driven approach might be used generally to generate and evaluate ontology design patterns.

Nail abnormalities identified in an ageing study of 30 inbred mouse strains.

In a large-scale ageing study, 30 inbred mouse strains were systematically screened for histologic evidence of lesions in all organ systems. Ten strains were diagnosed with similar nail abnormalities. The highest frequency was noted in NON/ShiLtJ mice. Lesions identified fell into two main categories: acute to chronic penetration of the third phalangeal bone through the hyponychium with associated inflammation and bone remodelling or metaplasia of the nail matrix and nail bed associated with severe orthokeratotic hyperkeratosis replacing the nail plate. Penetration of the distal phalanx through the hyponychium appeared to be the initiating feature resulting in nail abnormalities. The accompanying acute to subacute inflammatory response was associated with osteolysis of the distal phalanx. Evaluation of young NON/ShiLtJ mice revealed that these lesions were not often found, or affected only one digit. The only other nail unit abnormality identified was sporadic subungual epidermoid inclusion cysts which closely resembled similar lesions in human patients. These abnormalities, being age-related developments, may have contributed to weight loss due to impacts upon feeding and should be a consideration for future research due to the potential to interact with other experimental factors in ageing studies using the affected strains of mice.

Systematic screening for skin, hair, and nail abnormalities in a large-scale knockout mouse program.

The International Knockout Mouse Consortium was formed in 2007 to inactivate ("knockout") all protein-coding genes in the mouse genome in embryonic stem cells. Production and characterization of these mice, now underway, has generated and phenotyped 3,100 strains with knockout alleles. Skin and adnexa diseases are best defined at the gross clinical level and by histopathology. Representative retired breeders had skin collected from the back, abdomen, eyelids, muzzle, ears, tail, and lower limbs including the nails. To date, 169 novel mutant lines were reviewed and of these, only one was found to have a relatively minor sebaceous gland abnormality associated with follicular dystrophy. The B6N(Cg)-Far2tm2b(KOMP)Wtsi/2J strain, had lesions affecting sebaceous glands with what appeared to be a secondary follicular dystrophy. A second line, B6N(Cg)-Ppp1r9btm1.1(KOMP)Vlcg/J, had follicular dystrophy limited to many but not all mystacial vibrissae in heterozygous but not homozygous mutant mice, suggesting that this was a nonspecific background lesion. We discuss potential reasons for the low frequency of skin and adnexal phenotypes in mice from this project in comparison to those seen in human Mendelian diseases, and suggest alternative approaches to identification of human disease-relevant models.

Mouse genome-wide association study identifies polymorphisms on chromosomes 4, 11, and 15 for age-related cardiac fibrosis.

Dystrophic cardiac calcinosis (DCC), also called epicardial and myocardial fibrosis and mineralization, has been detected in mice of a number of laboratory inbred strains, most commonly C3H/HeJ and DBA/2J. In previous mouse breeding studies between these DCC susceptible and the DCC-resistant strain C57BL/6J, 4 genetic loci harboring genes involved in DCC inheritance were identified and subsequently termed Dyscalc loci 1 through 4. Here, we report susceptibility to cardiac fibrosis, a sub-phenotype of DCC, at 12 and 20 months of age and close to natural death in a survey of 28 inbred mouse strains. Eight strains showed cardiac fibrosis with highest frequency and severity in the moribund mice. Using genotype and phenotype information of the 28 investigated strains, we performed genome-wide association studies (GWAS) and identified the most significant associations on chromosome (Chr) 15 at 72 million base pairs (Mb) (P < 10(-13)) and Chr 4 at 122 Mb (P < 10(-11)) and 134 Mb (P < 10(-7)). At the Chr 15 locus, Col22a1 and Kcnk9 were identified. Both have been reported to be morphologically and functionally important in the heart muscle. The strongest Chr 4 associations were located approximately 6 Mb away from the Dyscalc 2 quantitative trait locus peak within the boundaries of the Extl1 gene and in close proximity to the Trim63 and Cap1 genes. In addition, a single-nucleotide polymorphism association was found on chromosome 11. This study provides evidence for more than the previously reported 4 genetic loci determining cardiac fibrosis and DCC. The study also highlights the power of GWAS in the mouse for dissecting complex genetic traits.

Dermal lymphatic dilation in a mouse model of alopecia areata.

Mouse models of various types of inflammatory skin disease are often accompanied by increased dermal angiogenesis. The C3H/HeJ inbred strain spontaneously develops alopecia areata (AA), a cell mediated autoimmune disorder that can be controllably expanded using full thickness skin grafts to young unaffected mice. This provides a reproducible and progressive model for AA in which the vascularization of the skin can be examined. Mice receiving skin grafts from AA or normal mice were evaluated at 5, 10, 15, and 20 weeks after engraftment. Lymphatics are often overlooked as they are small slit-like structures above the hair follicle that resemble artifact-like separation of collagen bundles with some fixatives. Lymphatics are easily detected using lymphatic vessel endothelial hyaluronan receptor 1 (LYVE1) by immunohistochemistry to label their endothelial cells. Using LYVE1, there were no changes in distribution or numbers of lymphatics although they were more prominent (dilated) in the mice with AA. Lyve1 transcripts were not significantly upregulated except at 10 weeks after skin grafting when clinical signs of AA first become apparent. Other genes involved with vascular growth and dilation or movement of immune cells were dysregulated, mostly upregulated. These findings emphasize aspects of AA not commonly considered and provide potential targets for therapeutic intervention.

Genetic determinants of fibro-osseous lesions in aged inbred mice.

Fibro-osseous lesions in mice are progressive aging changes in which the bone marrow is replaced to various degrees by fibrovascular stroma and bony trabeculae in a wide variety of bones. The frequency and severity varied greatly among 28 different inbred mouse stains, predominantly affecting females, ranging from 0% for 10 strains to 100% for KK/HlJ and NZW/LacJ female mice. Few lesions were observed in male mice and for 23 of the strains, no lesions were observed in males for any of the cohorts. There were no significant correlations between strain-specific severities of fibro-osseous lesions and ovarian (r=0.11; P=0.57) or endometrial (r=0.03; P=0.89) cyst formation frequency or abnormalities in parathyroid glands. Frequency of fibro-osseous lesions was most strongly associated (P<10(-6)) with genome variations on chromosome (Chr) 8 at 90.6 and 90.8Mb (rs33108071, rs33500669; P=5.0·10(-10), 1.3·10(-6)), Chr 15 at 23.6 and 23.8Mb (rs32087871, rs45770368; P=7.3·10(-7), 2.7·10(-6)), and Chr 19 at 33.2, 33.4, and 33.6Mb (rs311004232, rs30524929, rs30448815; P=2.8·10(-6), 2.8·10(-6), 2.8·10(-6)) in genome-wide association studies (GWAS). The relatively large number of candidate genes identified in the GWAS analyses suggests that this may be an extremely complex polygenic disease. These results indicate that fibro-osseous lesions are surprisingly common in many inbred strains of laboratory mice as they age. While this presents little problem in most studies that utilize young animals, it may complicate aging studies, particularly those focused on bone.

Crisp1 and alopecia areata in C3H/HeJ mice.

Alopecia areata (AA), a cell mediated autoimmune disease, is the second most common form of hair loss in humans. While the autoimmune disease is responsible for the underlying pathogenesis, the alopecia phenotype is ultimately due to hair shaft fragility and breakage associated with structural deficits. Quantitative trait genetic analyses using the C3H/HeJ mouse AA model identified cysteine-rich secretory protein 1 (Crisp1), a hair shaft structural protein, as a candidate gene within the major AA locus. Crisp1 transcripts in the skin at various times during disease development were barely detectable. In situ hybridization identified Crisp1 expression within the medulla of hair shafts from clinically normal strains of mice but not C3H/HeJ mice with AA. Follow-up work with 5-day-old C3H/HeJ mice with normal hair also had essentially no expression of Crisp1. Other non-inflammatory based follicular dystrophy mouse models with similar hair shaft abnormalities also have little or no Crisp1 expression. Shotgun proteomics, used to determine strain difference in hair proteins, confirmed that there was very little CRISP1 within normal C3H/HeJ mouse hair in comparison to 11 other strains. However, mutant mice with hair medulla defects also had undetectable levels of CRISP1 in their hair. Crisp1 null mice had normal skin, hair follicles, and hair shafts indicating that the lack of the CRISP1 protein does not translate directly into defects in the hair shaft or hair follicle. These results suggest that CRISP1 may be an important structural component of mouse hair and that its strain-specific dysregulation may indicate a predisposition to hair shaft disease such as AA.

Alopecia areata: updates from the mouse perspective.

Alopecia areata (AA) is a cell-mediated autoimmune disease that targets actively growing hair follicles in mammals, including humans and mice. Development of the C3H/HeJ spontaneous mouse model AA nearly 20 years ago provided a much needed tool to test the hypotheses and ultimately serve as a preclinical model for drug testing. Discoveries in both human AA patients and the mouse model supported each other and lead to discoveries on the incredibly complex genetic basis of this disease. The discovery that A/J, MRL/MpJ, SJL/J, and SWR/J strains also develop AA now allows genome-wide association mapping studies to expand the list of genes underlying this disease. Potential new targets for unraveling the pathogenesis of AA include the role of retinoic acid metabolism in the severity of disease and hair shaft proteins that may be either the inciting antigen or ultimate target of the immune reaction leading to breakage of the shaft causing clinical alopecia. Comparing these model systems with human and mouse clinical disease, for both discovery and validation of the discoveries, continues to resolve the complex questions surrounding AA.

The mouse pathology ontology, MPATH; structure and applications.

BACKGROUND: The capture and use of disease-related anatomic pathology data for both model organism phenotyping and human clinical practice requires a relatively simple nomenclature and coding system that can be integrated into data collection platforms (such as computerized medical record-keeping systems) to enable the pathologist to rapidly screen and accurately record observations. The MPATH ontology was originally constructed in 2,000 by a committee of pathologists for the annotation of rodent histopathology images, but is now widely used for coding and analysis of disease and phenotype data for rodents, humans and zebrafish. CONSTRUCTION AND CONTENT: MPATH is divided into two main branches describing pathological processes and structures based on traditional histopathological principles. It does not aim to include definitive diagnoses, which would generally be regarded as disease concepts. It contains 888 core pathology terms in an almost exclusively is_a hierarchy nine layers deep. Currently, 86% of the terms have textual definitions and contain relationships as well as logical axioms to other ontologies such the Gene Ontology. APPLICATION AND UTILITY: MPATH was originally devised for the annotation of histopathological images from mice but is now being used much more widely in the recording of diagnostic and phenotypic data from both mice and humans, and in the construction of logical definitions for phenotype and disease ontologies. We discuss the use of MPATH to generate cross-products with qualifiers derived from a subset of the Phenotype and Trait Ontology (PATO) and its application to large-scale high-throughput phenotyping studies. MPATH provides a largely species-agnostic ontology for the descriptions of anatomic pathology, which can be applied to most amniotes and is now finding extensive use in species other than mice. It enables investigators to interrogate large datasets at a variety of depths, use semantic analysis to identify the relations between diseases in different species and integrate pathology data with other data types, such as pharmacogenomics.

The mouse as a model for understanding chronic diseases of aging: the histopathologic basis of aging in inbred mice.

Inbred mice provide a unique tool to study aging populations because of the genetic homogeneity within an inbred strain, their short life span, and the tools for analysis which are available. A large-scale longitudinal and cross-sectional aging study was conducted on 30 inbred strains to determine, using histopathology, the type and diversity of diseases mice develop as they age. These data provide tools that when linked with modern in silico genetic mapping tools, can begin to unravel the complex genetics of many of the common chronic diseases associated with aging in humans and other mammals. In addition, novel disease models were discovered in some strains, such as rhabdomyosarcoma in old A/J mice, to diseases affecting many but not all strains including pseudoxanthoma elasticum, pulmonary adenoma, alopecia areata, and many others. This extensive data set is now available online and provides a useful tool to help better understand strain-specific background diseases that can complicate interpretation of genetically engineered mice and other manipulatable mouse studies that utilize these strains.

Recombinant human hepatitis B vaccine initiating alopecia areata: testing the hypothesis using the C3H/HeJ mouse model.

Untoward effects of human vaccines suggest that recombinant hepatitis B vaccine may induce alopecia areata (AA) in some patients. Similar untoward immunological effects may also account for AA-like diseases in domestic species. In this study, the C3H/HeJ spontaneous adult onset AA mouse model was used to test the role, if any, of recombinant hepatitis B vaccine on the initiation or activation of AA. Initial experiments demonstrated no effect on induction of AA in young adult female C3H/HeJ mice (P = 0.5689). By contrast, older females, those at the age when AA first begins to appear in this strain, had a significant increase (P = 0.0264) in the time of onset of AA, suggesting that the vaccine may initiate disease in mice predisposed to AA. However, larger vaccine trials, which included diphtheria and tetanus toxoids as additional controls, did not support these initial result findings and suggest that AA associated with vaccination may be within the normal background levels of the given population.

Integrating mouse anatomy and pathology ontologies into a phenotyping database: tools for data capture and training.

The Mouse Disease Information System (MoDIS) is a data capture system for pathology data from laboratory mice designed to support phenotyping studies. The system integrates the mouse anatomy (MA) and mouse pathology (MPATH) ontologies into a Microsoft Access database facilitating the coding of organ, tissue, and disease process to recognized semantic standards. Grading of disease severity provides scores for all lesions that can then be used for quantitative trait locus (QTL) analyses and haplotype association gene mapping. Direct linkage to the Pathbase online database provides reference definitions for disease terms and access to photomicrographic images of similar diagnoses in other mutant mice. MoDIS is an open source and freely available program (http://research.jax.org/faculty/sundberg/index.html). This provides a valuable tool for setting up a mouse pathology phenotyping program.