Genome-wide association mapping of blood cell traits in mice

Systems genetics approaches for understanding complex traits with relevance for human disease

Quantitative traits are often complex because of the contribution of many loci, with further complexity added by environmental factors. In medical research, systems genetics is a powerful approach for the study of complex traits, as it integrates intermediate phenotypes, such as RNA, protein, and metabolite levels, to understand molecular and physiological phenotypes linking discrete DNA sequence variation to complex clinical and physiological traits. The primary purpose of this review is to describe some of the resources and tools of systems genetics in humans and rodent models, so that researchers in many areas of biology and medicine can make use of the data.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Isometric skeletal muscle contractile properties in common strains of male laboratory mice

Assessing contractile function of skeletal muscle in murine models is a commonly employed laboratory technique that investigators utilize to measure the impact of genetic manipulations, drug efficacy, or other therapeutic interventions. Often overlooked is the potential for the strain of the mouse to influence the functional properties of the skeletal muscle. Thus, we sought to characterize commonly assessed isometric force measures in the hindlimb muscles across a variety of mouse strains. Using 6-8-week-old male mice, we measured isometric force, fatigue susceptibility, relaxation kinetics, muscle mass, myofiber cross-sectional area, and fiber type composition of the extensor digitorum longus (EDL) and soleus muscles in C57BL/6NJ, BALB/cJ, FVB/NJ, C57BL/6J, and C57BL/10 mice. The data demonstrate both unique differences and a number of similarities between both muscles in the various genetic backgrounds of mice. Soleus muscle specific force (i.e., force per unit size) exhibited higher variation across strains while specific force of the EDL muscle exhibited minimal variation. In contrast, absolute force differed only in a few mouse strains whereas analysis of muscle morphology revealed many distinctions when compared across all the groups. Collectively, the data suggest that the strain of the mouse can potentially influence the measured biological outcome and may possibly promote a synergistic effect with any genetic manipulation or therapeutic intervention. Thus, it is critical for the investigator to carefully consider the genetic background of the mouse used in the experimental design and precisely document the strain of mouse employed during publication.

Genetic Influence on Frequencies of Myeloid-Derived Cell Subpopulations in Mouse

Differences in frequencies of blood cell subpopulations were reported to influence the course of infections, atopic and autoimmune diseases, and cancer. We have discovered a unique mouse strain B10.O20 containing extremely high frequency of myeloid-derived cells (MDC) in spleen. B10.O20 carries 3.6% of genes of the strain O20 on the C57BL/10 genetic background. It contains much higher frequency of CD11b⁺Gr1⁺ cells in spleen than both its parents. B10.O20 carries O20-derived segments on chromosomes 1, 15, 17, and 18. Their linkage with frequencies of blood cell subpopulations in spleen was tested in F₂ hybrids between B10.O20 and C57BL/10. We found 3 novel loci controlling MDC frequencies: Mydc1, 2, and 3 on chromosomes 1, 15, and 17, respectively, and a locus controlling relative spleen weight (Rsw1) that co-localizes with Mydc3 and also influences proportion of white and red pulp in spleen. Mydc1 controls numbers of CD11b⁺Gr1⁺ cells. Interaction of Mydc2 and Mydc3 regulates frequency of CD11b⁺Gr1⁺ cells and neutrophils (Gr1⁺Siglec-F^- cells from CD11b⁺ cells). Interestingly, Mydc3/Rsw1 is orthologous with human segment 6q21 that was shown previously to determine counts of white blood cells. Bioinformatics analysis of genomic sequence of the chromosomal segments bearing these loci revealed polymorphisms between O20 and C57BL/10 that change RNA stability and genes’ functions, and we examined expression of relevant genes. This identified potential candidate genes Smap1, Vps52, Tnxb, and Rab44. Definition of genetic control of MDC can help to personalize therapy of diseases influenced by these cells.

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.

Plasma proteomic biomarker signature of age predicts health and life span

Older age is a strong shared risk factor for many chronic diseases, and there is increasing interest in identifying aging biomarkers. Here, a proteomic analysis of 1301 plasma proteins was conducted in 997 individuals between 21 and 102 years of age. We identified 651 proteins associated with age (506 over-represented, 145 underrepresented with age). Mediation analysis suggested a role for partial cis-epigenetic control of protein expression with age. Of the age-associated proteins, 33.5% and 45.3%, were associated with mortality and multimorbidity, respectively. There was enrichment of proteins associated with inflammation and extracellular matrix as well as senescence-associated secretory proteins. A 76-protein proteomic age signature predicted accumulation of chronic diseases and all-cause mortality. These data support the use of proteomic biomarkers to monitor aging trajectories and to identify individuals at higher risk of disease to be targeted for in depth diagnostic procedures and early interventions.

The Genetics of Variation of the Wave 1 Amplitude of the Mouse Auditory Brainstem Response

This is the first genome-wide association study with the Hybrid Mouse Diversity Panel (HDMP) to define the genetic landscape of the variation in the suprathreshold wave 1 amplitude of the auditory brainstem response (ABR) both pre- and post-noise exposure. This measure is correlated with the density of the auditory neurons (AN) and/or the compliment of synaptic ribbons within the inner hair cells of the mouse cochlea. We analyzed suprathreshold ABR for 635 mice from 102 HMDP strains pre- and post-noise exposure (108 dB 10 kHz octave band noise exposure for 2 h) using auditory brainstem response (ABR) wave 1 suprathreshold amplitudes as part of a large survey (Myint et al., Hear Res 332:113-120, 2016). Genome-wide significance levels for pre- and post-exposure wave 1 amplitude across the HMDP were performed using FaST-LMM. Synaptic ribbon counts (Ctbp2 and mGluR2) were analyzed for the extreme strains within the HMDP. ABR wave 1 amplitude varied across all strains of the HMDP with differences ranging between 2.42 and 3.82-fold pre-exposure and between 2.43 and 7.5-fold post-exposure with several tone burst stimuli (4 kHz, 8 kHz, 12 kHz, 16 kHz, 24 kHz, and 32 kHz). Immunolabeling of paired synaptic ribbons and glutamate receptors of strains with the highest and lowest wave 1 values pre- and post-exposure revealed significant differences in functional synaptic ribbon counts. Genome-wide association analysis identified genome-wide significant threshold associations on chromosome 3 (24 kHz; JAX00105429; p < 1.12E-06) and chromosome 16 (16 kHz; JAX00424604; p < 9.02E-07) prior to noise exposure and significant associations on chromosomes 2 (32 kHz; JAX00497967; p < 3.68E-08) and 13 (8 kHz; JAX00049416; 1.07E-06) after noise exposure. In order to prioritize candidate genes, we generated cis-eQTLs from microarray profiling of RNA isolated from whole cochleae in 64 of the tested strains.This is the first report of a genome-wide association analysis, controlled for population structure, to explore the genetic landscape of suprathreshold wave 1 amplitude measurements of the mouse ABR. We have defined two genomic regions associated with wave 1 amplitude variation prior to noise exposure and an additional two associated with variation after noise exposure.

A GWAS approach identifies Dapp1 as a determinant of air pollution-induced airway hyperreactivity

Asthma is a chronic inflammatory disease of the airways with contributions from genes, environmental exposures, and their interactions. While genome-wide association studies (GWAS) in humans have identified ~200 susceptibility loci, the genetic factors that modulate risk of asthma through gene-environment (GxE) interactions remain poorly understood. Using the Hybrid Mouse Diversity Panel (HMDP), we sought to identify the genetic determinants of airway hyperreactivity (AHR) in response to diesel exhaust particles (DEP), a model traffic-related air pollutant. As measured by invasive plethysmography, AHR under control and DEP-exposed conditions varied 3-4-fold in over 100 inbred strains from the HMDP. A GWAS with linear mixed models mapped two loci significantly associated with lung resistance under control exposure to chromosomes 2 (p = 3.0x10^-6) and 19 (p = 5.6x10^-7). The chromosome 19 locus harbors Il33 and is syntenic to asthma association signals observed at the IL33 locus in humans. A GxE GWAS for post-DEP exposure lung resistance identified a significantly associated locus on chromosome 3 (p = 2.5x10^-6). Among the genes at this locus is Dapp1, an adaptor molecule expressed in immune-related and mucosal tissues, including the lung. Dapp1-deficient mice exhibited significantly lower AHR than control mice but only after DEP exposure, thus functionally validating Dapp1 as one of the genes underlying the GxE association at this locus. In summary, our results indicate that some of the genetic determinants for asthma-related phenotypes may be shared between mice and humans, as well as the existence of GxE interactions in mice that modulate lung function in response to air pollution exposures relevant to humans.

The genetic factors that modulate risk of asthma through gene-environment (GxE) interactions are poorly understood, due in large part to the inherent difficulties in carrying out such studies in humans. To address these challenges, we used the Hybrid Mouse Diversity Panel to elucidate the genetic architecture of asthma-related phenotypes in mice and identify loci that are associated with airway hyperreactivity (AHR) under control exposure conditions and in response to diesel exhaust particles (DEP), as a model traffic-related air pollutant. In the absence of exposure, we identified two loci on chromosomes 2 and 19 for AHR. The locus on chromosome 19 harbors Il33 and is syntenic to association signals observed for asthma at the IL33 locus in humans. In response to DEP exposure, we mapped AHR to a region on chromosome 3 and used a genetically modified mouse model to functionally demonstrate that Dapp1 is one of the genes underlying the GxE association at this locus. Collectively, our results support the concept that some of the genetic determinants for asthma-related phenotypes may be shared between mice and humans as well as the existence of GxE interactions in mice that modulate lung function in response to air pollution exposures relevant to humans.

Genetic Variation in the Magnitude and Longevity of the IgG Subclass Response to a Diphtheria-Tetanus-Acellular Pertussis (DTaP) Vaccine in Mice

The type of IgG subclasses induced by vaccination is an important determinant of vaccine efficacy because the IgG subclasses vary in their biological function. The goal of this study was to determine the influence of the genetic background on the production and duration of vaccine-induced IgG subclasses. IgG1, IgG2b, and IgG3 titers against diphtheria toxoid (DT), pertussis toxin (PT), filamentous hemagglutinin (FHA), and pertactin (Prn) were measured in mice from 28 different inbred and wild-derived strains vaccinated with an aluminum hydroxide-adjuvanted DTaP vaccine. The titers and duration of vaccine-specific IgG subclass responses were different among mouse strains, indicating that genetic factors contribute to this variation. Statistical associations were used to identify potential mechanisms that contribute to antibody production and longevity. This analysis showed that the mechanisms guiding the magnitude of antibody production were antigen-dependent for IgG1 but antigen-independent for IgG2b and IgG3. However, the mechanisms driving the longevity of antibody titers were antigen-independent for IgG1, IgG2b, and IgG3. The ratio of IgG1 and IgG3 titers identified Th1 and Th2-prone mouse strains. TLR4-deficient C3H/HeJ mice had an enhanced IgG1 response compared with C3H/HeOuJ mice with intact TLR4. This work demonstrates that the genetic background contributes significantly to the magnitude and longevity of vaccine-induced IgG1, IgG2b, and IgG3 titers in mice.

Genomewide Association Study Identifies Cxcl Family Members as Partial Mediators of LPS-Induced Periodontitis

Periodontitis (PD) is characterized by bacterial infection and inflammation of tooth-supporting structures and can lead to tooth loss. PD affects ∼47% of the US population over age 30 years and has a heritability of about 50%. Although the host immunoinflammatory response and genetic background play a role, little is known of the underlying genetic factors. We examined natural genetic variation in lipopolysaccharide (LPS)-induced PD across a panel of inbred mouse strains, the hybrid mouse diversity panel (HMDP). We observed a strain-dependent sixfold difference in LPS-induced bone loss across the HMDP with a heritability of 53%. We performed a genomewide association study (GWAS) using FAST-LMM, which corrects for population structure, and identified loci significantly associated with PD. We examined candidate genes at a locus on chromosome 5, which suggested a relationship between LPS-induced bone loss and, together with expression data, identified Cxcl family members as associated with PD. We observed an increase in Cxcl10 protein, as well as immune cells and pro-inflammatory cytokines in C57BL/6J (high bone loss strain) but not in A/J (low bone loss strain) after LPS injections. Genetic deletion of CXCR3 (Cxcl9 and10 receptor) demonstrated a ∼50% reduction in bone loss and reduced osteoclasts after LPS injections. Furthermore, WT mice treated with AMG-487 (a CXCR3 antagonist) showed a ∼45% reduction in bone loss and decreased osteoclasts after LPS injections. We conclude that CXCR3 is a strong candidate for modulating the host response in individuals susceptible to PD. © 2018 American Society for Bone and Mineral Research.

Susceptibility of different mouse strains to peri-implantitis

BACKGROUND AND OBJECTIVE

Peri-implantitis (PI) is an inflammatory condition that affects the tissues surrounding dental implants. Although the pathogenesis of PI is not fully understood, evidence suggests that the etiology is multifactorial and may include a genetic component. The aim of this study was to investigate the role of genetics in the development of peri-implantitis.

MATERIAL AND METHODS

Four-week-old C57BL/6J, C3H/HeJ and A/J male mice had their left maxillary molars extracted. Implants were placed in the healed extraction sockets. Upon osseointegration, ligatures were placed around the implant head for 1 or 4 weeks to induce PI. Micro-computed tomography scanning was used to measure volumetric bone loss. Histological analyses were also performed to evaluate collagen organization and the presence of neutrophils and osteoclasts.

RESULTS

Radiographically, comparing the ligature-treated mice, C57BL/6J displayed the greatest amount of bone loss, followed by C3H/HeJ and A/J mice at 1 and 4 weeks. Histologically, at 1 week, C57BL/6J mice presented with the highest numbers of neutrophils and osteoclasts. At 4 weeks, C57BL/6J mice presented with the most active bone remodeling compared with the other two strains.

CONCLUSION

There were significant differences in the severity of peri-implantitis among the different mouse strains, suggesting that the genetic framework can affect implant survival and success. Future work is needed to dissect the genetic contribution to the development of peri-implantitis.

The Genetic Architecture of Coronary Artery Disease: Current Knowledge and Future Opportunities

PURPOSE OF REVIEW

We provide an overview of our current understanding of the genetic architecture of coronary artery disease (CAD) and discuss areas of research that provide excellent opportunities for further exploration.

RECENT FINDINGS

Large-scale studies in human populations, coupled with rapid advances in genetic technologies over the last decade, have clearly established the association of common genetic variation with risk of CAD. However, the effect sizes of the susceptibility alleles are for the most part modest and collectively explain only a small fraction of the overall heritability. By comparison, evidence that rare variants make a substantial contribution to risk of CAD has been somewhat disappointing thus far, suggesting that other biological mechanisms have yet to be discovered. Emerging data suggests that novel pathways involved in the development of CAD can be identified through complementary and integrative systems genetics strategies in mice or humans. There is also convincing evidence that gut bacteria play a previously unrecognized role in the development of CAD, particularly through metabolism of certain dietary nutrients that lead to proatherogenic metabolites in the circulation. A major effort is now underway to functionally understand the newly discovered genetic and biological associations for CAD, which could lead to the development of potentially novel therapeutic strategies. Other important areas of investigation for understanding the pathophysiology of CAD, including epistatic interactions between genes or with either sex and environmental factors, have not been studied on a broad scope and represent additional opportunities for future studies.

A Suite of Tools for Biologists That Improve Accessibility and Visualization of Large Systems Genetics Datasets: Applications to the Hybrid Mouse Diversity Panel

In this chapter we address the recent explosion in large multilevel population studies such as the METSIM study in humans as well as large panels of animal models such as the Hybrid Mouse Diversity Panel or the BXD set of recombinant inbred strains. These studies have harnessed the increasing affordability of large-scale high-throughput profiling to gather massive quantities of data. These datasets, spread across different -omics levels (genome, transcriptome, etc.), different tissues (e.g. heart, plasma, bone) and different environmental factors (e.g. diet, drugs) each individually have led to a number of novel findings relevant to a variety of complex diseases and other phenotypes. The analysis of these results, however, is often limited to individuals with a comprehensive understanding of database languages such as SQL. In this chapter, we describe the development of a GUI-based database analysis suite, using the Hybrid Mouse Diversity Panel as an example to lay out a series of methods for visualization and integration of large systems genetics datasets. The database is based on the Shiny suite of tools in R, and is transferrable to other SQL-based datasets.

Genome-Wide Association Mapping for Female Infertility in Inbred Mice

The genetic factors underlying female infertility in humans are only partially understood. Here, we performed a genome-wide association study of female infertility in 25 inbred mouse strains by using publicly available SNP data. As a result, a total of four SNPs were identified after chromosome-wise multiple test correction. The first SNP rs29972765 is located in a gene desert on chromosome 18, about 72 kb upstream of Skor2 (SKI family transcriptional corepressor 2). The second SNP rs30415957 resides in the intron of Plce1 (phospholipase C epsilon 1). The remaining two SNPs (rs30768258 and rs31216810) are close to each other on chromosome 19, in the vicinity of Sorbs1 (sorbin and SH3 domain containing 1). Using quantitative RT-PCR, we found that Sorbs1 is highly expressed in the mouse uterus during embryo implantation. Knockdown of Sorbs1 by siRNA attenuates the induction of differentiation marker gene Prl8a2 (decidual prolactin-related protein) in an in vitro model of decidualization using mouse endometrial stromal cells, suggesting that Sorbs1 may be a potential candidate gene for female infertility in mice. Our results may represent an opportunity to further understand female infertility in humans.

Genome-Wide Association Analysis Identifies Dcc as an Essential Factor in the Innervation of the Peripheral Vestibular System in Inbred Mice

This study aimed to investigate the genetic causes of vestibular dysfunction. We used vestibular sensory-evoked potentials (VsEPs) to characterize the vestibular function of 35 inbred mouse strains selected from the Hybrid Mouse Diversity Panel and demonstrated strain-dependent phenotypic variation in vestibular function. Using these phenotypic data, we performed the first genome-wide association study controlling for population structure that has revealed two highly suggestive loci, one of which lies within a haplotype block containing five genes (Stard6, 4930503L19Rik, Poli, Mbd2, Dcc) on Chr. 18 (peak SNP rs29632020), one gene, deleted in colorectal carcinoma (Dcc) has a well-established role in nervous system development. An in-depth analysis of Dcc-deficient mice demonstrated elevation in mean VsEP threshold for Dcc (+/-) mice (-11.86 dB) compared to wild-type (-9.68 dB) littermates. Synaptic ribbon studies revealed Dcc (-/-) (P0) and Dcc (+/-) (6-week-old) mice showed lower density of the presynaptic marker (CtBP2) as compared to wild-type controls. Vestibular ganglion cell counts of Dcc (-/-) (P0) was lower than controls. Whole-mount preparations showed abnormal innervation of the utricle, saccule, and crista ampullaris at E14.5, E16.5, and E18.5. Postnatal studies were limited by the perinatal lethality in Dcc (-/-) mice. Expression analyses using in situ hybridization and immunohistochemistry showed Dcc expression in the mouse vestibular ganglion (E15.5), and utricle and crista ampullaris (6-week-old), respectively. In summary, we report the first GWAS for vestibular functional variation in inbred mice and provide evidence for the role of Dcc in the normal innervation of the peripheral vestibular system.

The Hybrid Mouse Diversity Panel: a resource for systems genetics analyses of metabolic and cardiovascular traits

The Hybrid Mouse Diversity Panel (HMDP) is a collection of approximately 100 well-characterized inbred strains of mice that can be used to analyze the genetic and environmental factors underlying complex traits. While not nearly as powerful for mapping genetic loci contributing to the traits as human genome-wide association studies, it has some important advantages. First, environmental factors can be controlled. Second, relevant tissues are accessible for global molecular phenotyping. Finally, because inbred strains are renewable, results from separate studies can be integrated. Thus far, the HMDP has been studied for traits relevant to obesity, diabetes, atherosclerosis, osteoporosis, heart failure, immune regulation, fatty liver disease, and host-gut microbiota interactions. High-throughput technologies have been used to examine the genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes of the mice under various environmental conditions. All of the published data are available and can be readily used to formulate hypotheses about genes, pathways and interactions.

Genetic Architecture of Atherosclerosis in Mice: A Systems Genetics Analysis of Common Inbred Strains

Common forms of atherosclerosis involve multiple genetic and environmental factors. While human genome-wide association studies have identified numerous loci contributing to coronary artery disease and its risk factors, these studies are unable to control environmental factors or examine detailed molecular traits in relevant tissues. We now report a study of natural variations contributing to atherosclerosis and related traits in over 100 inbred strains of mice from the Hybrid Mouse Diversity Panel (HMDP). The mice were made hyperlipidemic by transgenic expression of human apolipoprotein E-Leiden (APOE-Leiden) and human cholesteryl ester transfer protein (CETP). The mice were examined for lesion size and morphology as well as plasma lipid, insulin and glucose levels, and blood cell profiles. A subset of mice was studied for plasma levels of metabolites and cytokines. We also measured global transcript levels in aorta and liver. Finally, the uptake of acetylated LDL by macrophages from HMDP mice was quantitatively examined. Loci contributing to the traits were mapped using association analysis, and relationships among traits were examined using correlation and statistical modeling. A number of conclusions emerged. First, relationships among atherosclerosis and the risk factors in mice resemble those found in humans. Second, a number of trait-loci were identified, including some overlapping with previous human and mouse studies. Third, gene expression data enabled enrichment analysis of pathways contributing to atherosclerosis and prioritization of candidate genes at associated loci in both mice and humans. Fourth, the data provided a number of mechanistic inferences; for example, we detected no association between macrophage uptake of acetylated LDL and atherosclerosis. Fifth, broad sense heritability for atherosclerosis was much larger than narrow sense heritability, indicating an important role for gene-by-gene interactions. Sixth, stepwise linear regression showed that the combined variations in plasma metabolites, including LDL/VLDL-cholesterol, trimethylamine N-oxide (TMAO), arginine, glucose and insulin, account for approximately 30 to 40% of the variation in atherosclerotic lesion area. Overall, our data provide a rich resource for studies of complex interactions underlying atherosclerosis.

While recent genetic association studies in human populations have succeeded in identifying genetic loci that contribute to coronary artery disease (CAD) and related phenotypes, these loci explain only a small fraction of the genetic variation in CAD and associated traits. Here, we present a complementary approach using association analysis of atherosclerotic traits among inbred strains of mice. A strength of this approach is that it enables in-depth phenotypic characterization including gene expression and metabolic profiling across a variety of tissues, and integration of these molecular phenotypes with coronary artery disease itself. A striking finding was the large fraction of atherosclerosis that was explained by genetic interactions. Association analysis allowed us to identify genetic loci for atherosclerotic lesion area as well as transcript, cytokine and metabolite levels, and relationships among the traits were examined by correlation and network modeling. The plasma metabolites associated with atherosclerosis in mice, namely, LDL/VLDL-cholesterol, TMAO, arginine, glucose and insulin, overlapped with those observed in humans and accounted for approximately 30 to 40% of the observed variation in atherosclerotic lesion area. In summary, our data provide a detailed overview of the genetic architecture of atherosclerosis in mice and a rich resource for studies of the complex genetic and metabolic interactions that underlie the disease.

Genetic Mapping of Novel Loci Affecting Canine Blood Phenotypes

Since the publication of the dog genome and the construction of high-quality genome-wide SNP arrays, thousands of dogs have been genotyped for disease studies. For many of these dogs, additional clinical phenotypes are available, such as hematological and clinical chemistry results collected during routine veterinary care. Little is known about the genetic basis of variation in blood phenotypes, but this variation may play an important role in the etiology and progression of many diseases. From a cohort of dogs that had been previously genotyped on a semi-custom Illumina CanineHD array for various genome-wide association studies (GWAS) at Cornell University Hospital for Animals, we chose 353 clinically healthy, adult dogs for our analysis of clinical pathologic test results (14 hematological tests and 25 clinical chemistry tests). After correcting for age, body weight and sex, genetic associations were identified for amylase, segmented neutrophils, urea nitrogen, glucose, and mean corpuscular hemoglobin. Additionally, a strong genetic association (P = 8.1×10⁻¹³) was evident between a region of canine chromosome 13 (CFA13) and alanine aminotransferase (ALT), explaining 23% of the variation in ALT levels. This region of CFA13 encompasses the GPT gene that encodes the transferase. Dogs homozygous for the derived allele exhibit lower ALT activity, making increased ALT activity a less useful marker of hepatic injury in these individuals. Overall, these associations provide a roadmap for identifying causal variants that could improve interpretation of clinical blood tests and understanding of genetic risk factors associated with diseases such as canine diabetes and anemia, and demonstrate the utility of holistic phenotyping of dogs genotyped for disease mapping studies.

High-Density Genotypes of Inbred Mouse Strains: Improved Power and Precision of Association Mapping

Human genome-wide association studies have identified thousands of loci associated with disease phenotypes. Genome-wide association studies also have become feasible using rodent models and these have some important advantages over human studies, including controlled environment, access to tissues for molecular profiling, reproducible genotypes, and a wide array of techniques for experimental validation. Association mapping with common mouse inbred strains generally requires 100 or more strains to achieve sufficient power and mapping resolution; in contrast, sample sizes for human studies typically are one or more orders of magnitude greater than this. To enable well-powered studies in mice, we have generated high-density genotypes for ∼175 inbred strains of mice using the Mouse Diversity Array. These new data increase marker density by 1.9-fold, have reduced missing data rates, and provide more accurate identification of heterozygous regions compared with previous genotype data. We report the discovery of new loci from previously reported association mapping studies using the new genotype data. The data are freely available for download, and Web-based tools provide easy access for association mapping and viewing of the underlying intensity data for individual loci.

The Genetic Landscape of Hematopoietic Stem Cell Frequency in Mice

Prior 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.

•

Genetic variation across mouse strains influences hematopoietic stem cell frequency

•

This variation can be exploited for genome-wide association studies

•

Hopx is a regulator of hematopoietic stem/progenitor cell function

•

This approach can be used to identify genetic determinants of other stem cell systems

Contrasting genetic architectures in different mouse reference populations used for studying complex traits

Quantitative trait loci (QTLs) are being used to study genetic networks, protein functions, and systems properties that underlie phenotypic variation and disease risk in humans, model organisms, agricultural species, and natural populations. The challenges are many, beginning with the seemingly simple tasks of mapping QTLs and identifying their underlying genetic determinants. Various specialized resources have been developed to study complex traits in many model organisms. In the mouse, remarkably different pictures of genetic architectures are emerging. Chromosome Substitution Strains (CSSs) reveal many QTLs, large phenotypic effects, pervasive epistasis, and readily identified genetic variants. In contrast, other resources as well as genome-wide association studies (GWAS) in humans and other species reveal genetic architectures dominated with a relatively modest number of QTLs that have small individual and combined phenotypic effects. These contrasting architectures are the result of intrinsic differences in the study designs underlying different resources. The CSSs examine context-dependent phenotypic effects independently among individual genotypes, whereas with GWAS and other mouse resources, the average effect of each QTL is assessed among many individuals with heterogeneous genetic backgrounds. We argue that variation of genetic architectures among individuals is as important as population averages. Each of these important resources has particular merits and specific applications for these individual and population perspectives. Collectively, these resources together with high-throughput genotyping, sequencing and genetic engineering technologies, and information repositories highlight the power of the mouse for genetic, functional, and systems studies of complex traits and disease models.

Genome-wide association study identifies nox3 as a critical gene for susceptibility to noise-induced hearing loss

In the United States, roughly 10% of the population is exposed daily to hazardous levels of noise in the workplace. Twin studies estimate heritability for noise-induced hearing loss (NIHL) of approximately 36%, and strain specific variation in sensitivity has been demonstrated in mice. Based upon the difficulties inherent to the study of NIHL in humans, we have turned to the study of this complex trait in mice. We exposed 5 week-old mice from the Hybrid Mouse Diversity Panel (HMDP) to a 10 kHz octave band noise at 108 dB for 2 hours and assessed the permanent threshold shift 2 weeks post exposure using frequency specific stimuli. These data were then used in a genome-wide association study (GWAS) using the Efficient Mixed Model Analysis (EMMA) to control for population structure. In this manuscript we describe our GWAS, with an emphasis on a significant peak for susceptibility to NIHL on chromosome 17 within a haplotype block containing NADPH oxidase-3 (Nox3). Our peak was detected after an 8 kHz tone burst stimulus. Nox3 mutants and heterozygotes were then tested to validate our GWAS. The mutants and heterozygotes demonstrated a greater susceptibility to NIHL specifically at 8 kHz both on measures of distortion product otoacoustic emissions (DPOAE) and on auditory brainstem response (ABR). We demonstrate that this sensitivity resides within the synaptic ribbons of the cochlea in the mutant animals specifically at 8 kHz. Our work is the first GWAS for NIHL in mice and elucidates the power of our approach to identify tonotopic genetic susceptibility to NIHL.

Noise-induced hearing loss (NIHL) is the most common work-related disease in the world and the second cause of hearing loss. Although several candidate gene association studies for NIHL in humans have been conducted, each are underpowered, un-replicated, and account for only a fraction of the genetic risk. Buoyed by the prospects and successes of human association studies, several groups have proposed mouse genome-wide association studies. The environment can be carefully controlled, facilitating the study of complex traits like NIHL. In this manuscript, we describe, for the first time, an association analysis with correction for population structure for the mapping of several loci for susceptibility to NIHL in inbred strains of mice. We identify Nox3 as the associated gene for susceptibility to NIHL that the genetic susceptibility is frequency specific and that it occurs at the level of the cochlear synaptic ribbon.

Epigenomes: the missing heritability in human cardiovascular disease?

Cardiovascular disease is a tremendous burden on human health and results from malfunction of various networks of biological molecules in the context of environmental stress. Despite strong evidence of heritability, many common forms of heart disease (heart failure in particular) have not yielded to genome-wide association studies to identify causative mutations acting via the disruption of individual molecules. Increasing evidence suggests, however, that genetic variation in noncoding regions is strongly linked to disease susceptibility. We hypothesize that epigenomic variation may engender different chromatin environments in the absence of (or in parallel with) changes in protein or mRNA sequence and abundance. In this manner, distinct-genetically encoded-chromatin environments can exhibit distinct responses to environmental stresses that cause heart failure, explaining a significant portion of the altered susceptibility that is observed in human disease.

A systems genetics study of swine illustrates mechanisms underlying human phenotypic traits

Background

The pig, which shares greater similarities with human than with mouse, is important for agriculture and for studying human diseases. However, similarities in the genetic architecture and molecular regulations underlying phenotypic variations in humans and swine have not been systematically assessed.

Results

We systematically surveyed ~500 F2 pigs genetically and phenotypically. By comparing candidates for anemia traits identified in swine genome-wide SNP association and human genome-wide association studies (GWAS), we showed that both sets of candidates are related to the biological process “cellular lipid metabolism” in liver. Human height is a complex heritable trait; by integrating genome-wide SNP data and human adipose Bayesian causal network, which closely represents bone transcriptional regulations, we identified PLAG1 as a causal gene for limb bone length. This finding is consistent with GWAS findings for human height and supports the common genetic architecture between swine and humans. By leveraging a human protein-protein interaction network, we identified two putative candidate causal genes TGFB3 and DAB2IP and the known regulators MESP1 and MESP2 as responsible for the variation in rib number and identified the potential underlying molecular mechanisms. In mice, knockout of Tgfb3 and Tgfb2 together decreases rib number.

Conclusion

Our findings show that integrative network analyses reveal causal regulators underlying the genetic association of complex traits in swine and that these causal regulators have similar effects in humans. Thus, swine are a potentially good animal model for studying some complex human traits that are not under intense selection.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1240-y) contains supplementary material, which is available to authorized users.

In-silico QTL mapping of postpubertal mammary ductal development in the mouse uncovers potential human breast cancer risk loci

Genetic background plays a dominant role in mammary gland development and breast cancer (BrCa). Despite this, the role of genetics is only partially understood. This study used strain-dependent variation in an inbred mouse mapping panel, to identify quantitative trait loci (QTL) underlying structural variation in mammary ductal development, and determined if these QTL correlated with genomic intervals conferring BrCa susceptibility in humans. For about half of the traits, developmental variation among the complete set of strains in this study was greater (P < 0.05) than that of previously studied strains, or strains in current common use for mammary gland biology. Correlations were also detected with previously reported variation in mammary tumor latency and metastasis. In-silico genome-wide association identified 20 mammary development QTL (Mdq). Of these, five were syntenic with previously reported human BrCa loci. The most significant (P = 1 × 10(-11)) association of the study was on MMU6 and contained the genes Plxna4, Plxna4os1, and Chchd3. On MMU5, a QTL was detected (P = 8 × 10(-7)) that was syntenic to a human BrCa locus on h12q24.5 containing the genes Tbx3 and Tbx5. Intersection of linked SNP (r(2) > 0.8) with genomic and epigenomic features, and intersection of candidate genes with gene expression and survival data from human BrCa highlighted several for further study. These results support the conclusion that mammary tumorigenesis and normal ductal development are influenced by common genetic factors and that further studies of genetically diverse mice can improve our understanding of BrCa in humans.

Quantitative trait Loci for resistance to the congenital nephropathy in tensin 2-deficient mice

The ICR-derived glomerulonephritis (ICGN) mouse is a chronic kidney disease (CKD) model that is characterized histologically by glomerulosclerosis, vascular sclerosis and tubulointerstitial fibrosis, and clinically by proteinuria and anemia, which are common symptoms and pathological changes associated with a variety of kidney diseases. Previously, we performed a quantitative trait locus (QTL) analysis to identify the causative genes for proteinuria in ICGN mice, and found a deletion mutation of the tensin 2 gene (Tns2^nph, MGI no: 2447990). Interestingly, the congenic strain carrying the Tns2^nph mutation on a C57BL/6J (B6) genetic background exhibited milder phenotypes than did ICGN mice, indicating the presence of several modifier genes controlling the disease phenotype. In this study, to identify the modifier/resistant loci for CKD progression in Tns2-deficient mice, we performed QTL analysis using backcross progenies from susceptible ICGN and resistant B6 mice. We identified a significant locus on chromosome (Chr) 2 (LOD = 5.36; 31 cM) and two suggestive loci on Chrs 10 (LOD = 2.27; 64 cM) and X (LOD = 2.65; 67 cM) with linkage to the severity of tubulointerstitial injury. One significant locus on Chr 13 (LOD = 3.49; approximately 14 cM) and one suggestive locus on Chr 2 (LOD = 2.41; 51 cM) were identified as QTLs for the severity of glomerulosclerosis. Suggestive locus in BUN was also detected in the same Chr 2 region (LOD = 2.34; 51 cM). A locus on Chr 2 (36 cM) was significantly linked with HGB (LOD = 4.47) and HCT (LOD = 3.58). Four novel epistatic loci controlling either HGB or tubulointerstitial injury were discovered. Further genetic analysis should lead to identification of CKD modifier gene(s), aiding early diagnosis and providing novel approaches to the discovery of drugs for the treatment and possible prevention of kidney disease.

Genome-wide association study for age-related hearing loss (AHL) in the mouse: a meta-analysis

Age-related hearing loss (AHL) is characterized by a symmetric sensorineural hearing loss primarily in high frequencies and individuals have different levels of susceptibility to AHL. Heritability studies have shown that the sources of this variance are both genetic and environmental, with approximately half of the variance attributable to hereditary factors as reported by Huag and Tang (Eur Arch Otorhinolaryngol 267(8):1179-1191, 2010). Only a limited number of large-scale association studies for AHL have been undertaken in humans, to date. An alternate and complementary approach to these human studies is through the use of mouse models. Advantages of mouse models include that the environment can be more carefully controlled, measurements can be replicated in genetically identical animals, and the proportion of the variability explained by genetic variation is increased. Complex traits in mouse strains have been shown to have higher heritability and genetic loci often have stronger effects on the trait compared to humans. Motivated by these advantages, we have performed the first genome-wide association study of its kind in the mouse by combining several data sets in a meta-analysis to identify loci associated with age-related hearing loss. We identified five genome-wide significant loci (<10(-6)). One of these loci confirmed a previously identified locus (ahl8) on distal chromosome 11 and greatly narrowed the candidate region. Specifically, the most significant associated SNP is located 450 kb upstream of Fscn2. These data confirm the utility of this approach and provide new high-resolution mapping information about variation within the mouse genome associated with hearing loss.

FVB/NJ mice demonstrate a youthful sensitivity to noise-induced hearing loss and provide a useful genetic model for the study of neural hearing loss

The hybrid mouse diversity panel (HMDP), a panel of 100 strains, has been employed in genome wide association studies (GWAS) to study complex traits in mice. Hearing is a complex trait and the CBA/CaJ mouse strain is a widely used model for age-related hearing loss (ARHI) and noise induced hearing loss (NIHL). The CBA/CaJ strain's youthful sensitivity to noise and limited age-related loss led us to attempt to identify additional strains segregating a similar phenotype for our panel. FVB/NJ is part of the HMDP and has been previously described as having a similar ARHI phenotype to CBA/CaJ. For these reasons, we have studied the FVB/NJ mouse for ARHI and NIHL phenotypes in hopes of incorporating its phenotype into HMDP studies. We demonstrate that FVB/NJ exhibits ARHI at an earlier age than CBA/CaJ and young FVB/NJ mice are vulnerable to NIHL up until 10 to 12 weeks. This suggests that FVB/NJ may be used as an additional genetic model for neural forms of progressive hearing loss and for the study of youthful sensitivity to noise.

Systems genetics approaches to understand complex traits

Systems genetics is an approach to understand the flow of biological information that underlies complex traits. It uses a range of experimental and statistical methods to quantitate and integrate intermediate phenotypes, such as transcript, protein or metabolite levels, in populations that vary for traits of interest. Systems genetics studies have provided the first global view of the molecular architecture of complex traits and are useful for the identification of genes, pathways and networks that underlie common human diseases. Given the urgent need to understand how the thousands of loci that have been identified in genome-wide association studies contribute to disease susceptibility, systems genetics is likely to become an increasingly important approach to understanding both biology and disease.

The systems genetics resource: a web application to mine global data for complex disease traits

The Systems Genetics Resource (SGR) (http://systems.genetics.ucla.edu) is a new open-access web application and database that contains genotypes and clinical and intermediate phenotypes from both human and mouse studies. The mouse data include studies using crosses between specific inbred strains and studies using the Hybrid Mouse Diversity Panel. SGR is designed to assist researchers studying genes and pathways contributing to complex disease traits, including obesity, diabetes, atherosclerosis, heart failure, osteoporosis, and lipoprotein metabolism. Over the next few years, we hope to add data relevant to deafness, addiction, hepatic steatosis, toxin responses, and vascular injury. The intermediate phenotypes include expression array data for a variety of tissues and cultured cells, metabolite levels, and protein levels. Pre-computed tables of genetic loci controlling intermediate and clinical phenotypes, as well as phenotype correlations, are accessed via a user-friendly web interface. The web site includes detailed protocols for all of the studies. Data from published studies are freely available; unpublished studies have restricted access during their embargo period.

Continuous glucose monitoring in infants of very low birth weight

OBJECTIVES

To evaluate the feasibility and efficacy of a continuous glucose monitoring system (CGMS) in a population of infants of very low birth weight (VLBW).

STUDY DESIGN

Infants weighing <or=1,500 g and of <or=32 weeks of gestation were recruited within 24 h of delivery. A subcutaneous sensor connected to a CGMS was inserted and maintained for 7 days or until dysfunction. Therapeutic management followed the usual standard protocols.

RESULTS

38 patients (21 male) were included over 17 months. Their mean gestational age was 27.5 +/- 2.0 weeks and their mean birth weight was 958.3 +/- 205.5 g. Their perinatal histories and complications during admission were unremarkable for extremely premature babies. Continuous monitoring lasted an average of 7.84 +/- 1.99 days per patient. Hyperglycaemia was detected in 22 (57.90%) patients and it lasted a mean of 20.33 +/- 30.13 h, while 14 (36.8%) presented with hypoglycaemia for a mean of 2.45 +/- 2.3 h.

CONCLUSIONS

The CGMS gave a safe and useful estimate of glucose levels in VLBW infants, revealing abnormal glucose levels at a much higher rate than expected by usual sampling. However, it was not able to provide real-time glucose concentration data. CGMS may be very useful in providing information on the role of hyper- and hypoglycaemia on short- and long-term outcomes in VLBW infants.