Searching for genetic factors of fatty liver in SMXA-5 mice by quantitative trait loci analysis under a high-fat diet

Discovery of genomic loci for liver health and steatosis reveals overlap with glutathione redox genetics

Non-alcoholic fatty liver disease (NAFLD) is the most common chronic liver condition in the United States, encompassing a wide spectrum of liver pathologies including steatosis, steatohepatitis, fibrosis, and cirrhosis. Despite its high prevalence, there are no medications currently approved by the Food and Drug Administration for the treatment of NAFLD. Recent work has suggested that NAFLD has a strong genetic component and identifying causative genes will improve our understanding of the molecular mechanisms contributing to NAFLD and yield targets for future therapeutic investigations. Oxidative stress is known to play an important role in NAFLD pathogenesis, yet the underlying mechanisms accounting for disturbances in redox status are not entirely understood. To better understand the relationship between the glutathione redox system and signs of NAFLD in a genetically-diverse population, we measured liver weight, serum biomarkers aspartate aminotransferase (AST) and alanine aminotransferase (ALT), and graded liver pathology in a large cohort of Diversity Outbred mice. We compared hepatic endpoints to those of the glutathione redox system previously measured in the livers and kidneys of the same mice, and we screened for statistical and genetic associations using the R/qtl2 software. We discovered several novel genetic loci associated with markers of liver health, including loci that were associated with both liver steatosis and glutathione redox status. Candidate genes within each locus point to possible new mechanisms underlying the complex relationship between NAFLD and the glutathione redox system, which could have translational implications for future studies targeting NAFLD pathology.

Neonatal lethality of mouse A/J-7SM consomic strain is caused by an insertion mutation in the Dchs1 gene

Homosomic mice of the A/J-7^SM consomic mouse strain that introduced the entire chromosome 7 (Chr 7) of SM/J into the A/J strain exhibited neonatal lethality. We tentatively maintained segregating inbred strains (A/J-7A^SM and A/J-7D^SM) in which the central portion of Chr 7 was heterozygous for the A/J and SM/J strains, and the centromeric and telomeric sides of Chr 7 were homozygous for the SM/J strain, instead of the A/J-7^SM strain. Based on the chromosomal constitution of Chr 7 in A/J-7A^SM and A/J-7D^SM mice, the causative gene for neonatal lethality in homosomic mice was suggested to be located within an approximately 1.620 Mb region between D7Mit125 (104.879 Mb) and D7Mit355 (106.499 Mb) on Chr 7. RT-PCR analysis revealed that homosomic mice lacked dachsous cadherin-related 1 (Dchs1), which is located within the D7Mit125 to D7Mit355 region and functions in the regulation of planar cell polarity. Screening for mutations in Dchs1 indicated that homosomic mice possessed an early transposable (ETn)-like sequence in intron 1 of Dchs1. Moreover, an allelism test between Dchs1 ETn-like-insertion alleles detected in homosomic mice and CRISPR/Cas9-induced Dchs1 deletion alleles revealed that Dchs1 is a causative gene for neonatal lethality in homosomic mice. Based on these results, we concluded that in the A/J-7^SM strain, ETn-like elements were inserted into intron 1 of SM/J-derived Dchs1 during strain development, which dramatically reduced Dchs1 expression, thus resulting in neonatal lethality in homosomic mice. Additionally, it was suggested that the timing of lethality in Dchs1 mutant mice is influenced by the genetic background.

Multi-Omic Approaches to Identify Genetic Factors in Metabolic Syndrome

Metabolic syndrome (MetS) is a highly heritable disease and a major public health burden worldwide. MetS diagnosis criteria are met by the simultaneous presence of any three of the following: high triglycerides, low HDL/high LDL cholesterol, insulin resistance, hypertension, and central obesity. These diseases act synergistically in people suffering from MetS and dramatically increase risk of morbidity and mortality due to stroke and cardiovascular disease, as well as certain cancers. Each of these component features is itself a complex disease, as is MetS. As a genetically complex disease, genetic risk factors for MetS are numerous, but not very powerful individually, often requiring specific environmental stressors for the disease to manifest. When taken together, all sequence variants that contribute to MetS disease risk explain only a fraction of the heritable variance, suggesting additional, novel loci have yet to be discovered. In this article, we will give a brief overview on the genetic concepts needed to interpret genome-wide association studies (GWAS) and quantitative trait locus (QTL) data, summarize the state of the field of MetS physiological genomics, and to introduce tools and resources that can be used by the physiologist to integrate genomics into their own research on MetS and any of its component features. There is a wealth of phenotypic and molecular data in animal models and humans that can be leveraged as outlined in this article. Integrating these multi-omic QTL data for complex diseases such as MetS provides a means to unravel the pathways and mechanisms leading to complex disease and promise for novel treatments. © 2022 American Physiological Society. Compr Physiol 12:1-40, 2022.

Ablation of Iah1, a candidate gene for diet-induced fatty liver, does not affect liver lipid accumulation in mice

Nonalcoholic fatty liver disease (NAFLD) is a pathological condition caused by excess triglyceride deposition in the liver. The SMXA-5 severe fatty liver mouse model has been established from the SM/J and A/J strains. To explore the genetic factors involved in fatty liver development in SMXA-5 mice, we had previously performed quantitative trait locus (QTL) analysis, using (SM/J×SMXA-5)F2 intercross mice, and identified Fl1sa on chromosome 12 (centromere-53.06 Mb) as a significant QTL for fatty liver. Furthermore, isoamyl acetate-hydrolyzing esterase 1 homolog (Iah1) was selected as the most likely candidate gene for Fl1sa. Iah1 gene expression in fatty liver-resistant A/J-12^SM mice was significantly higher than in fatty liver-susceptible A/J mice. These data indicated that the Iah1 gene might be associated with fatty liver development. However, the function of murine Iah1 remains unknown. Therefore, in this study, we created Iah1 knockout (KO) mice with two different backgrounds [C57BL/6N (B6) and A/J-12^SM (A12)] to investigate the relationship between Iah1 and liver lipid accumulation. Liver triglyceride accumulation in Iah1-KO mice of B6 or A12 background did not differ from their respective Iah1-wild type mice under a high-fat diet. These results indicated that loss of Iah1 did not contribute to fatty liver. On the other hands, adipose tissue dysfunction causes lipid accumulation in ectopic tissues (liver, skeletal muscle, and pancreas). To investigate the effect of Iah1 deficiency on white adipose tissue, we performed DNA microarray analysis of epididymal fat in Iah1-KO mice of A12 background. This result showed that Iah1 deficiency might decrease adipokines Sfrp4 and Metrnl gene expression in epididymal fat. This study demonstrated that Iah1 deficiency did not cause liver lipid accumulation and that Iah1 was not a suitable candidate gene for Fl1sa.

Genetic dissection of the fatty liver QTL Fl1sa by using congenic mice and identification of candidate genes in the liver and epididymal fat

BACKGROUND

Nonalcoholic fatty liver disease (NAFLD) is a multifactorial disease caused by interactions between environmental and genetic factors. The SMXA-5 mouse is a high-fat diet-induced fatty liver model established from SM/J and A/J strains. We have previously identified Fl1sa, a quantitative trait locus (QTL) for fatty liver on chromosome 12 (centromere-53.06 Mb) of SMXA-5 mice. However, the chromosomal region containing Fl1sa was too broad. The aim of this study was to narrow the Fl1sa region by genetic dissection using novel congenic mice and to identify candidate genes within the narrowed Fl1sa region.

RESULTS

We established two congenic strains, R2 and R3, from parental A/J-12^SM and A/J strains. R2 and R3 strains have genomic intervals of centromere-29.20 Mb and 29.20-46.75 Mb of chromosome 12 derived from SM/J, respectively. Liver triglyceride content in R2 and R3 mice was significantly lower than that in A/J mice fed with a high-fat diet for 7 weeks. This result suggests that at least one of the genes responsible for fatty liver exists within the two chromosomal regions centromere-29.20 Mb (R2) and 29.20-46.75 Mb (R3). We found that liver triglyceride accumulation is inversely correlated with epididymal fat weight among the parental and congenic strains. Therefore, the ectopic fat accumulation in the liver may be due to organ-organ interactions between the liver and epididymal fat. To identify candidate genes in Fl1sa, we performed a DNA microarray analysis using the liver and epididymal fat in A/J and A/J-12^SM mice fed with a high-fat diet for 7 weeks. In epididymal fat, mRNA levels of Zfp125 (in R2) and Nrcam (in R3) were significantly different in A/J-12^SM mice from those in A/J mice. In the liver, mRNA levels of Iah1 (in R2) and Rrm2 (in R2) were significantly different in A/J-12^SM mice from those in A/J mice.

CONCLUSIONS

In this study, using congenic mice analysis, we narrowed the chromosomal region containing Fl1sa to two regions of mouse chromosome 12. We then identified 4 candidate genes in Fl1sa: Iah1 and Rrm2 from the liver and Zfp125 and Nrcam from epididymal fat.

Detection of differentially expressed candidate genes for a fatty liver QTL on mouse chromosome 12

Background

The SMXA-5 mouse is an animal model of high-fat diet-induced fatty liver. The major QTL for fatty liver, Fl1sa on chromosome 12, was identified in a SM/J × SMXA-5 intercross. The SMXA-5 genome consists of the SM/J and A/J genomes, and the A/J allele of Fl1sa is a fatty liver-susceptibility allele. The existence of the responsible genes for fatty liver within Fl1sa was confirmed in A/J-12^SM consomic mice. The aim of this study was to identify candidate genes for Fl1sa, and to investigate whether the identified genes affect the lipid metabolism.

Results

A/J-12^SM mice showed a significantly lower liver triglyceride content compared to A/J mice when fed the high-fat diet for 7 weeks. We detected differences in the accumulation of liver lipids in response to the high-fat diet between A/J and A/J-12^SM consomic mice. To identify candidate genes for Fl1sa, we performed DNA microarray analysis using the livers of A/J-12^SM and A/J mice fed the high-fat diet. The mRNA levels of three genes (Iah1, Rrm2, Prkd1) in the chromosomal region of Fl1sa were significantly different between the strains. Iah1 mRNA levels in the liver, kidney, and lung were significantly higher in A/J-12^SM mice than in A/J mice. The hepatic Iah1 mRNA level in A/J-12^SM mice was 3.2-fold higher than that in A/J mice. To examine the effect of Iah1 on hepatic lipid metabolism, we constructed a stable cell line expressing the mouse Iah1 protein in mouse hepatoma Hepa1-6 cells. Overexpression of Iah1 in Hepa1-6 cells suppressed the mRNA levels of Cd36 and Dgat2, which play important roles in triglyceride synthesis and lipid metabolism.

Conclusions

These results demonstrated that Fl1sa on the proximal region of chromosome 12 affected fatty liver in mice on a high-fat diet. Iah1 (isoamyl acetate-hydrolyzing esterase 1 homolog) was identified as one of the candidate genes for Fl1sa. This study revealed that the mouse Iah1 gene regulated the expression of genes related to lipid metabolism in the liver.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0385-2) contains supplementary material, which is available to authorized users.

The genetic architecture of NAFLD among inbred strains of mice

To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis.

Prenatal polycyclic aromatic hydrocarbon, adiposity, peroxisome proliferator-activated receptor (PPAR) γ methylation in offspring, grand-offspring mice

Rationale

Greater levels of prenatal exposure to polycyclic aromatic hydrocarbon (PAH) have been associated with childhood obesity in epidemiological studies. However, the underlying mechanisms are unclear.

Objectives

We hypothesized that prenatal PAH over-exposure during gestation would lead to weight gain and increased fat mass in offspring and grand-offspring mice. Further, we hypothesized that altered adipose gene expression and DNA methylation in genes important to adipocyte differentiation would be affected.

Materials and Methods

Pregnant dams were exposed to a nebulized PAH mixture versus negative control aerosol 5 days a week, for 3 weeks. Body weight was recorded from postnatal day (PND) 21 through PND60. Body composition, adipose cell size, gene expression of peroxisome proliferator-activated receptor (PPAR) γ, CCAAT/enhancer-binding proteins (C/EBP) α, cyclooxygenase (Cox)-2, fatty acid synthase (FAS) and adiponectin, and DNA methylation of PPAR γ, were assayed in both the offspring and grand-offspring adipose tissue.

Findings

Offspring of dams exposed to greater PAH during gestation had increased weight, fat mass, as well as higher gene expression of PPAR γ, C/EBP α, Cox2, FAS and adiponectin and lower DNA methylation of PPAR γ. Similar differences in phenotype and DNA methylation extended through the grand-offspring mice.

Conclusions

Greater prenatal PAH exposure was associated with increased weight, fat mass, adipose gene expression and epigenetic changes in progeny.

ATR-FTIR spectroscopy reveals genomic loci regulating the tissue response in high fat diet fed BXD recombinant inbred mouse strains

Background

Obesity-associated organ-specific pathological states can be ensued from the dysregulation of the functions of the adipose tissues, liver and muscle. However, the influence of genetic differences underlying gross-compositional differences in these tissues is largely unknown. In the present study, the analytical method of ATR-FTIR spectroscopy has been combined with a genetic approach to identify genetic differences responsible for phenotypic alterations in adipose, liver and muscle tissues.

Results

Mice from 29 BXD recombinant inbred mouse strains were put on high fat diet and gross-compositional changes in adipose, liver and muscle tissues were measured by ATR-FTIR spectroscopy. The analysis of genotype-phenotype correlations revealed significant quantitative trait loci (QTL) on chromosome 12 for the content of fat and collagen, collagen integrity, and the lipid to protein ratio in adipose tissue and on chromosome 17 for lipid to protein ratio in liver. Using gene expression and sequence information, we suggest Rsad2 (viperin) and Colec11 (collectin-11) on chromosome 12 as potential quantitative trait candidate genes. Rsad2 may act as a modulator of lipid droplet contents and lipid biosynthesis; Colec11 might play a role in apoptopic cell clearance and maintenance of adipose tissue. An increased level of Rsad2 transcripts in adipose tissue of DBA/2J compared to C57BL/6J mice suggests a cis-acting genetic variant leading to differential gene activation.

Conclusion

The results demonstrate that the analytical method of ATR-FTIR spectroscopy effectively contributed to decompose the macromolecular composition of tissues that accumulate fat and to link this information with genetic determinants. The candidate genes in the QTL regions may contribute to obesity-related diseases in humans, in particular if the results can be verified in a bigger BXD cohort.

Biglycan deletion alters adiponectin expression in murine adipose tissue and 3T3-L1 adipocytes

Obesity promotes increased secretion of a number of inflammatory factors from adipose tissue. These factors include cytokines and very lately, extracellular matrix components (ECM). Biglycan, a small leucine rich proteoglycan ECM protein, is up-regulated in obesity and has recently been recognized as a pro-inflammatory molecule. However, it is unknown whether biglycan contributes to adipose tissue dysfunction. In the present study, we characterized biglycan expression in various adipose depots in wild-type mice fed a low fat diet (LFD) or obesity-inducing high fat diet (HFD). High fat feeding induced biglycan mRNA expression in multiple adipose depots. Adiponectin is an adipokine with anti-inflammatory and insulin sensitizing effects. Due to the importance of adiponectin, we examined the effect of biglycan on adiponectin expression. Comparison of adiponectin expression in biglycan knockout (bgn(-/0)) and wild-type (bgn(+/0)) reveals higher adiponectin mRNA and protein in epididymal white adipose tissue in bgn(-/0) mice, as well higher serum concentration of adiponectin, and lower serum insulin concentration. On the contrary, knockdown of biglycan in 3T3-L1 adipocytes led to decreased expression and secretion of adiponectin. Furthermore, treatment of 3T3-L1 adipocytes with conditioned medium from biglycan treated macrophages resulted in an increase in adiponectin mRNA expression. These data suggest a link between biglycan and adiponectin expression. However, the difference in the pattern of regulation between in vivo and in vitro settings reveals the complexity of this relationship.

Establishment of consomic strains derived from A/J and SM/J mice for genetic analysis of complex traits

Consomic strains, in which one chromosome is derived from a donor strain and the other chromosomes are derived from the recipient strain, provide a powerful tool for the dissection of complex genetic traits. In this study we established ten consomic strains (A-2(SM), A-6(SM), A-11(SM), A-12(SM), A-13(SM), A-15(SM), A-17(SM), A-18(SM), A-19(SM), A-Y(SM)) using the SM/J strain as the donor and the A/J strain as the recipient; these are the parental strains of a set of SMXA recombinant inbred (RI) strains that we had developed previously. We analyzed body weights and blood lipid levels in the consomic and parental strains. The mean values for each trait showed a continuous range of variation in the consomic strains suggesting that they are controlled by multiple genes. We previously identified suggestive QTLs for body weight on chromosome 6 in SMXA RI strains and (SM/J × A/J)F(2) mice. The observation that the A-6(SM) consomic strain had a significantly lower mean body weight than the A/J strain supports the presence of this QTL on chromosome 6. Similarly, the higher blood triglyceride level in the A-11(SM) strain shows the existence of a previously mapped QTL on chromosome 11, and the A-12(SM) strain provides evidence of a QTL for blood total cholesterol level on chromosome 12. These consomic strains, along with the previously developed set of SMXA RI strains from A/J and SM/J mice, offer an invaluable and powerful resource for the analysis of complex genetic traits in mice.

Chromosome substitution strains: gene discovery, functional analysis, and systems studies

Laboratory mice are valuable in biomedical research in part because of the extraordinary diversity of genetic resources that are available for studies of complex genetic traits and as models for human biology and disease. Chromosome substitution strains (CSSs) are important in this resource portfolio because of their demonstrated use for gene discovery, genetic and epigenetic studies, functional characterizations, and systems analysis. CSSs are made by replacing a single chromosome in a host strain with the corresponding chromosome from a donor strain. A complete CSS panel involves a total of 22 engineered inbred strains, one for each of the 19 autosomes, one each for the X and Y chromosomes, and one for mitochondria. A genome survey simply involves comparing each phenotype for each of the CSSs with the phenotypes of the host strain. The CSS panels that are available for laboratory mice have been used to dissect a remarkable variety of phenotypes and to characterize an impressive array of disease models. These surveys have revealed considerable phenotypic diversity even among closely related progenitor strains, evidence for strong epistasis and for heritable epigenetic changes. Perhaps most importantly, and presumably because of their unique genetic constitution, CSSs, and congenic strains derived from them, the genetic variants underlying quantitative trait loci (QTLs) are readily identified and functionally characterized. Together these studies show that CSSs are important resource for laboratory mice.

Quantitative trait loci affecting liver fat content in mice

Nonalcoholic fatty liver disease, a condition in which excess fat accumulates in the liver, is strongly associated with the metabolic syndrome, including obesity and other related conditions. This disease has the potential to progress from steatosis to steatohepatitis, fibrosis, and cirrhosis. The recent increase in the prevalence of the metabolic syndrome is largely driven by changes in diet and activity levels. Individual variation in the response to this obesogenic environment, however, is attributable in part to genetic variation between individuals, but very few mammalian genetic loci have been identified with effects on fat accumulation in the liver. To study the genetic basis for variation in liver fat content in response to dietary fat, liver fat proportion was determined using quantitative magnetic resonance imaging in 478 mice from 16 LG/J X SM/J recombinant inbred strains fed either a high-fat (42% kcal from fat) or low-fat (15% kcal from fat) diet. An analysis of variance confirmed that there is a genetic basis for variation in liver fat content within the population with significant effects of sex and diet. Three quantitative trail loci that contribute to liver fat content also were mapped.

Genetic analysis of abdominal fat distribution in SM/J and A/J mice

Each abdominal fat depot, such as mesenteric or epididymal, differently contributes to the development of insulin resistance. The aim of this study was to identify the genetic regions that contribute to fat accumulation in epididymal/mesenteric fat and to examine whether or not the genetic regions that affect glucose metabolism and body fat distribution are coincident. We previously mapped a major quantitative trait locus (QTL) (T2dm2sa) for impaired glucose tolerance on chromosome 2 and revealed that SM.A-T2dm2sa congenic mice showed not only glucose tolerance but also fat accumulation. In the present study, to identify the loci/genes that control the accumulation of abdominal fat, we performed QTL analyses of epididymal/mesenteric fat weight by using (A/J x SM.A-T2dm2sa)F2 mice in which the effect of T2dm2sa was excluded. As a result, two highly significant QTLs for mesenteric fat, as well as three significant QTLs for epididymal/mesenteric fat, were mapped on the different chromosomal regions. This suggests that the fat accumulations in individual fat depots are controlled by distinct genomic regions. Our comparison of these QTLs for abdominal fat distribution with those for glucose metabolism revealed that the major genetic factors affecting body fat distribution do not coincide with genetic factors affecting glucose metabolism in (A/J x SM.A-T2dm2sa)F2.

Body and major organ weights of A/J-Chr 11(SM) consomic mice

The body and major organ weights of A/J-Chr 11(SM) consomic mice were compared to those of the progenitor strains, A/J and SM/J. The weights of the body and organs, except for those of the brain and uterus, were significantly greater in A/J mice than in SM/J mice. However, those of consomic mice were highly variable. For example, the average body weight of consomic mice was significantly greater than that of SM/J mice and nearly equal to that of A/J mice. Chromosome 11 of SM/J mice induced various significant changes of the organ weights of A/J mice, especially those of the brain, lung, kidney, adrenal, and ovary, demonstrating the importance of this chromosome in determining the phenotypes.

A Cmv2 QTL on chromosome X affects MCMV resistance in New Zealand male mice

NK cell-mediated resistance to viruses is subject to genetic control in humans and mice. Here we used classical and quantitative genetic strategies to examine NK-mediated murine cytomegalovirus (MCMV) control in genealogically related New Zealand white (NZW) and black (NZB) mice. NZW mice display NK cell-dependent MCMV resistance while NZB NK cells fail to limit viral replication after infection. Unlike Ly49H(+) NK resistance in C57BL/6 mice, NZW NK-mediated MCMV control was Ly49H-independent. Instead, MCMV resistance in NZW (Cmv2) involves multiple genetic factors. To establish the genetic basis of Cmv2 resistance, we further characterized a major chromosome X-linked resistance locus (DXMit216) responsible for innate MCMV control in NZW x NZB crosses. We found that the DXMit216 locus affects early MCMV control in New Zealand F(2) crosses and demonstrate that the NZB-derived DXMit216 allele enhances viral resistance in F(2) males. The evolutionary conservation of the DXMit216 region in mice and humans suggests that a Cmv2-related mechanism may affect human antiviral responses.

The metabolic syndrome

The "metabolic syndrome" (MetS) is a clustering of components that reflect overnutrition, sedentary lifestyles, and resultant excess adiposity. The MetS includes the clustering of abdominal obesity, insulin resistance, dyslipidemia, and elevated blood pressure and is associated with other comorbidities including the prothrombotic state, proinflammatory state, nonalcoholic fatty liver disease, and reproductive disorders. Because the MetS is a cluster of different conditions, and not a single disease, the development of multiple concurrent definitions has resulted. The prevalence of the MetS is increasing to epidemic proportions not only in the United States and the remainder of the urbanized world but also in developing nations. Most studies show that the MetS is associated with an approximate doubling of cardiovascular disease risk and a 5-fold increased risk for incident type 2 diabetes mellitus. Although it is unclear whether there is a unifying pathophysiological mechanism resulting in the MetS, abdominal adiposity and insulin resistance appear to be central to the MetS and its individual components. Lifestyle modification and weight loss should, therefore, be at the core of treating or preventing the MetS and its components. In addition, there is a general consensus that other cardiac risk factors should be aggressively managed in individuals with the MetS. Finally, in 2008 the MetS is an evolving concept that continues to be data driven and evidence based with revisions forthcoming.