Genome-wide Association Study In Mice Research Articles

GWAS in Mice Maps Susceptibility to HIV-Associated Nephropathy to the Ssbp2 Locus.

To gain insight into the pathogenesis of collapsing glomerulopathy, a rare form of FSGS that often arises in the setting of viral infections, we performed a genome-wide association study (GWAS) among inbred mouse strains using a murine model of HIV-1 associated nephropathy (HIVAN). We first generated F1 hybrids between HIV-1 transgenic mice on the FVB/NJ background and 20 inbred laboratory strains. Analysis of histology, BUN, and urinary NGAL demonstrated marked phenotypic variation among the transgenic F1 hybrids, providing strong evidence for host genetic factors in the predisposition to nephropathy. A GWAS in 365 transgenic F1 hybrids generated from these 20 inbred strains was performed. We identified a genome-wide significant locus on chromosome 13-C3 and multiple additional suggestive loci. Crossannotation of the Chr. 13 locus, including single-cell transcriptomic analysis of wildtype and HIV-1 transgenic mouse kidneys, nominated Ssbp2 as the most likely candidate gene. Ssbp2 is highly expressed in podocytes, encodes a transcriptional cofactor that interacts with LDB1 and LMX1B, which are both previously implicated in FSGS. Consistent with these data, older Ssbp2 null mice spontaneously develop glomerulosclerosis, tubular casts, interstitial fibrosis, and inflammation, similar to the HIVAN mouse model. These findings demonstrate the utility of GWAS in mice to uncover host genetic factors for rare kidney traits and suggest Ssbp2 as susceptibility gene for HIVAN, potentially acting via the LDB1-LMX1B transcriptional network.

Noise Exposure and Distortion Product Otoacoustic Emission Suprathreshold Amplitudes: A Genome-Wide Association Study

Background: Although several candidate-gene association studies have been conducted to investigate noise-induced hearing loss (NIHL) in humans, most are underpowered, unreplicated, and account for only a fraction of the genetic risk. Mouse genome-wide association studies (GWASs) have revolutionized the field of genetics and have led to the discovery of hundreds of genes involved in complex traits. The hybrid mouse diversity panel (HMDP) is a collection of classic inbred and recombinant inbred strains whose genomes have been either genotyped at high resolution or sequenced. To further investigate the genetics of NIHL, we report the first GWAS based on distortion product otoacoustic emission (DPOAE) measurements and the HMDP. Methods: A total of 102 strains (n = 635) from the HMDP were evaluated based on DPOAE suprathreshold amplitudes before and after noise exposure. DPOAE amplitude variation was set at 60 and 70 dB SPL of the primary tones for each frequency separately (8, 11.3, 16, 22.6, and 32 kHz). These values provided an indirect assessment of outer hair cell integrity. Six-week-old mice were exposed for 2 h to 10 kHz octave-band noise at 108 dB SPL. To perform local expression quantitative trait locus (eQTL) analysis, gene expression microarray profiles were generated using cochlear RNA from 64 hybrid mouse strains (n = 3 arrays per strain). Results: Several new loci were identified and positional candidate-genes associated with NIHL were prioritized, especially after noise exposure (1 locus at baseline and 5 loci after exposure). A total of 35 candidate genes in these 6 loci were identified with at least 1 probe whose expression was regulated by a significant cis-eQTL in the cochlea. After careful analysis of the candidate genes based on cochlear gene expression, 2 candidate genes were prioritized: Eya1 (baseline) and Efr3a (post-exposure). Discussion and Conclusion: For the first time, an association analysis with correction for population structure was used to map several loci for hearing traits in inbred strains of mice based on DPOAE suprathreshold amplitudes before and after noise exposure. Our results identified a number of novel loci and candidate genes for susceptibility to NIHL, especially the Eya1 and Efr3a genes. Our findings validate the power of the HMDP for detecting NIHL susceptibility genes.

Lrp1 Regulation of Pulmonary Function. Follow-Up of Human GWAS in Mice.

Human genome-wide association studies (GWASs) have identified more than 270 loci associated with pulmonary function; however, follow-up studies to determine causal genes at these loci are few. SNPs in low-density lipoprotein receptor-related protein 1 (LRP1) are associated with human pulmonary function in GWASs. Using murine models, we investigated the effect of genetic disruption of the Lrp1 gene in smooth muscle cells on pulmonary function in naive animals and after exposure to bacterial LPS or house dust mite extract. Disruption of Lrp1 in smooth muscle cells leads to an increase in tissue resistance, elastance, and tissue elastance at baseline. Furthermore, disruption of Lrp1 in smooth muscle increases airway responsiveness as measured by increased total lung resistance and airway resistance after methacholine. Immune cell counts in BAL fluid were increased in animals with Lrp1 disruption. The difference in airway responsiveness by genotype observed in naive animals was not observed after LPS or house dust mite extract exposure. To further explore the mechanisms contributing to changes in pulmonary function, we identified several ligands dysregulated with Lrp1 disruption in smooth muscle cells. These data suggest that dysregulation of LRP1 in smooth muscle cells affects baseline pulmonary function and airway responsiveness and helps establish LRP1 as the causal gene at this GWAS locus.

A genome-wide association study in mice reveals a role for Rhbdf2 in skeletal homeostasis

Low bone mass and an increased risk of fracture are predictors of osteoporosis. Individuals who share the same bone-mineral density (BMD) vary in their fracture risk, suggesting that microstructural architecture is an important determinant of skeletal strength. Here, we utilized the rich diversity of the Collaborative Cross mice to identify putative causal genes that contribute to the risk of fractures. Using microcomputed tomography, we examined key structural features that pertain to bone quality in the femoral cortical and trabecular compartments of male and female mice. We estimated the broad-sense heritability to be 50–60% for all examined traits, and we identified five quantitative trait loci (QTL) significantly associated with six traits. We refined each QTL by combining information inferred from the ancestry of the mice, ranging from RNA-Seq data and published literature to shortlist candidate genes. We found strong evidence for new candidate genes, particularly Rhbdf2, whose close association with the trabecular bone volume fraction and number was strongly suggested by our analyses. We confirmed our findings with mRNA expression assays of Rhbdf2 in extreme-phenotype mice, and by phenotyping bones of Rhbdf2 knockout mice. Our results indicate that Rhbdf2 plays a decisive role in bone mass accrual and microarchitecture.

A genome-wide association study in mice reveals a role for Rhbdf2 in skeletal homeostasis

Low bone mass and an increased risk of fracture are predictors of osteoporosis. Individuals who share the same bone-mineral density (BMD) vary in their fracture risk, suggesting that microstructural architecture is an important determinant of skeletal strength. Here, we utilized the rich diversity of the Collaborative Cross mice to identify putative causal genes that contribute to the risk of fractures. Using microcomputed tomography, we examined key structural features that pertain to bone quality in the femoral cortical and trabecular compartments of male and female mice. We estimated the broad-sense heritability to be 50–60% for all examined traits, and we identified five quantitative trait loci (QTL) significantly associated with six traits. We refined each QTL by combining information inferred from the ancestry of the mice, ranging from RNA-Seq data and published literature to shortlist candidate genes. We found strong evidence for new candidate genes, particularly Rhbdf2, whose close association with the trabecular bone volume fraction and number was strongly suggested by our analyses. We confirmed our findings with mRNA expression assays of Rhbdf2 in extreme-phenotype mice, and by phenotyping bones of Rhbdf2 knockout mice. Our results indicate that Rhbdf2 plays a decisive role in bone mass accrual and microarchitecture.

Mouse genome-wide association studies and systems genetics uncover the genetic architecture associated with hepatic pharmacokinetic and pharmacodynamic properties of a constrained ethyl antisense oligonucleotide targeting Malat1.

Antisense oligonucleotides (ASOs) have demonstrated variation of efficacy in patient populations. This has prompted our investigation into the contribution of genetic architecture to ASO pharmacokinetics (PK) and pharmacodynamics (PD). Genome wide association (GWA) and transcriptomic analysis in a hybrid mouse diversity panel (HMDP) were used to identify and validate novel genes involved in the uptake and efficacy of a single dose of a Malat1 constrained ethyl (cEt) modified ASO. The GWA of the HMDP identified two significant associations on chromosomes 4 and 10 with hepatic Malat1 ASO concentrations. Stabilin 2 (Stab2) and vesicle associated membrane protein 3 (Vamp3) were identified by cis-eQTL analysis. HMDP strains with lower Stab2 expression and Stab2 KO mice displayed significantly lower PK than strains with higher Stab2 expression and the wild type (WT) animals respectively, confirming the role of Stab2 in regulating hepatic Malat1 ASO uptake. GWA examining ASO efficacy uncovered three loci associated with Malat1 potency: Small Subunit Processome Component (Utp11l) on chromosome 4, Rho associated coiled-coil containing protein kinase 2 (Rock2) and Aci-reductone dioxygenase (Adi1) on chromosome 12. Our results demonstrate the utility of mouse GWAS using the HMDP in detecting genes capable of impacting the uptake of ASOs, and identifies genes critical for the activity of ASOs in vivo.

Genome-Wide Association Mapping of the Antibody Response to Diphtheria, Tetanus and Acellular Pertussis Vaccine in Mice.

The objective of the current study was to investigate the genetics of antibody responses to an acellular pertussis vaccine by a genome-wide association study in mice. Female mice of 28 inbred strains received this vaccine at 6, 8, and 12 weeks of age. The antibody titer and avidity of immunoglobulin (Ig) G specific for diphtheria toxin, pertussis toxin, filamentous hemagglutinin and pertactin were measured at 14 and 24 weeks of age. The magnitude, longevity and avidity of IgG differed significantly among mouse strains. There was significant correlation between antigen-specific IgGs for longevity but not for magnitude and avidity. Association mapping and analysis with PolyPhen software identified 6 genetic markers associated with longevity for all 4 antigens, although the expression levels of these genes did not correlate with longevity phenotype. This study provides novel insights into the genetic basis and potential candidate genes for differences in the IgG responses to vaccination.

Genome-wide association study of behavioral, physiological and gene expression traits in outbred CFW mice.

Although mice are the most widely used model organism, genetic studies have suffered from limited mapping resolution due to extensive linkage disequilibrium (LD) that is characteristic of crosses among inbred strains. Carworth Farms White (CFW) mice are a commercially available outbred mouse population that exhibit rapid LD decay compared to other available mouse populations. We performed a genome-wide association study (GWAS) of behavioral, physiological and gene expression phenotypes using 1,200 male CFW mice. We used genotyping-by-sequencing (GBS) to obtain genotypes at 92,734 single nucleotide polymorphisms (SNPs). We also measured gene expression using RNA-Sequencing in three brain regions. Our study identified numerous behavioral, physiological and expression quantitative trait loci (QTLs). We integrated the behavioral QTL and eQTL results to implicate specific genes, including Azi2 in sensitivity to methamphetamine and Zmynd11 in anxiety-like behavior. The combination of CFW mice, GBS and RNA-Sequencing constitutes a powerful approach to GWAS in mice.

Genome-wide association study in mice identifies loci affecting liver-related phenotypes including Sel1l influencing serum bile acids.

Using publicly available data from inbred mouse strains, we conducted a genome-wide association study to identify loci that accounted for liver-related phenotypes between C57BL/6J and A/J mice fed a Paigen diet. We confirmed genome-wide significant associations for hepatic cholesterol (chromosome 10A2) and serum total bile acid concentration (chromosome 12E) and identified a new locus for liver inflammation (chromosome 7C). Analysis of consomic mice confirmed that chromosome 12 A/J alleles accounted for the variance in serum total bile acid concentrations and had pleiotropic effects on liver mass, serum cholesterol, and serum alanine aminotransferase activity. Using an affected-only haplotype analysis among strains, we refined the chromosome 12E signal to a 1.95 Mb linkage disequilibrium block containing only one gene, sel-1 suppressor of lin-12-like (Sel1l). RNA sequencing and immunoblotting demonstrated that the risk allele locally conferred reduced expression of SEL1L in liver and distantly down-regulated pathways associated with hepatocyte nuclear factor 1 homeobox A (Hnf1a) and hepatocyte nuclear factor 4A (Hnf4a), known modifiers of bile acid transporters and metabolic traits. Consistent with these data, knockdown of SEL1L in HepG2 cells resulted in reduced HNF1A and HNF4A and increased bile acids in culture media; it further captured multiple molecular signatures observed in consomic mouse livers with reduced SEL1L. Finally, dogs harboring a SEL1L mutation and Sel1l(+/-) mice fed a Paigen diet had significantly increased serum total bile acid concentrations, providing independent confirmation linking SEL1L to bile acid metabolism. Genetic analyses of inbred mouse strains identified loci affecting different liver-related traits and implicated Sel1l as a significant determinant of serum bile acid concentration. (Hepatology 2016;63:1943-1956).

The Genetic Landscape of Hematopoietic Stem Cell Frequency in Mice

SummaryPrior efforts to identify regulators of hematopoietic stem cell physiology have relied mainly on candidate gene approaches with genetically modified mice. Here we used a genome-wide association study (GWAS) strategy with the hybrid mouse diversity panel to identify the genetic determinants of hematopoietic stem/progenitor cell (HSPC) frequency. Among 108 strains, we observed ∼120- to 300-fold variation in three HSPC populations. A GWAS analysis identified several loci that were significantly associated with HSPC frequency, including a locus on chromosome 5 harboring the homeodomain-only protein gene (Hopx). Hopx previously had been implicated in cardiac development but was not known to influence HSPC biology. Analysis of the HSPC pool in Hopx−/− mice demonstrated significantly reduced cell frequencies and impaired engraftment in competitive repopulation assays, thus providing functional validation of this positional candidate gene. These results demonstrate the power of GWAS in mice to identify genetic determinants of the hematopoietic system.

Comparative Genome-Wide Association Studies in Mice and Humans for Trimethylamine N-Oxide, a Proatherogenic Metabolite of Choline and L-Carnitine

Comparative Genome-Wide Association Studies in Mice and Humans for Trimethylamine N-Oxide, a Proatherogenic Metabolite of Choline and L-Carnitine

Genome-wide association studies in mice

Genome-wide association studies (GWASs) have transformed the field of human genetics and have led to the discovery of hundreds of genes that are implicated in human disease. The technological advances that drove this revolution are now poised to transform genetic studies in model organisms, including mice. However, the design of GWASs in mouse strains is fundamentally different from the design of human GWASs, creating new challenges and opportunities. This Review gives an overview of the novel study designs for mouse GWASs, which dramatically improve both the statistical power and resolution compared to classical gene-mapping approaches.

Non-MHC Risk Alleles in Rheumatoid Arthritis and in the Syntenic Chromosome Regions of Corresponding Animal Models

Rheumatoid arthritis (RA) is a polygenic autoimmune disease primarily affecting the synovial joints. Numerous animal models show similarities to RA in humans; some of them not only mimic the clinical phenotypes but also demonstrate the involvement of homologous genomic regions in RA. This paper compares corresponding non-MHC genomic regions identified in rodent and human genome-wide association studies (GWAS). To date, over 30 non-MHC RA-associated loci have been identified in humans, and over 100 arthritis-associated loci have been identified in rodent models of RA. The genomic regions associated with the disease are designated by the name(s) of the gene having the most frequent and consistent RA-associated SNPs or a function suggesting their involvement in inflammatory or autoimmune processes. Animal studies on rats and mice preferentially have used single sequence length polymorphism (SSLP) markers to identify disease-associated qualitative and quantitative trait loci (QTLs) in the genome of F2 hybrids of arthritis-susceptible and arthritis-resistant rodent strains. Mouse GWAS appear to be far ahead of rat studies, and significantly more mouse QTLs correspond to human RA risk alleles.

Dark Matter: Are Mice the Solution to Missing Heritability?

Genome-wide association studies (GWAS) in humans have identified hundreds of single nucleotide polymorphisms associated with complex traits, yet for most traits studied, the sum total of all these identified variants fail to explain a significant portion of the heritable variation. Reasons for this “missing heritability” are thought to include the existence of rare causative variants not captured by current genotyping arrays, structural variants that go undetected by existing technology, insufficient power to identify multi-gene interactions, small sample sizes, and the influence of environmental and epigenetic effects. As genotyping technologies have evolved it has become inexpensive and relatively straightforward to perform GWAS in mice. Mice offer a powerful tool for elucidating the genetic architecture of behavioral and physiological traits, and are complementary to human studies. Unlike F2 crosses of inbred strains, advanced intercross lines, heterogeneous stocks, outbred, and wild-caught mice have more rapid breakdown of linkage disequilibrium which allow for increasingly high resolution mapping. Because some of these populations are created using a small number of founder chromosomes they are not expected to harbor rare alleles. We discuss the differences between these mouse populations and examine their potential to overcome some of the pitfalls that have plagued human GWAS studies.