QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

doi:10.21203/rs.3.rs-1339100/v1

Download PDF

Research Article

QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

https://doi.org/10.21203/rs.3.rs-1339100/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background

The Berlin Fat Mouse Inbred line (BFMI) is a model for obesity and the metabolic syndrome. This study aimed to identify genetic variants associated with liver weight, liver triglycerides, and body weight using the obese BFMI sub-line BFMI861-S1. BFMI861-S1 mice are insulin resistant and store ectopic fat in the liver.

Methods

In generation 10, 58 males and 65 females of the advanced intercross line (AIL) BFMI861-S1xB6N were phenotyped under a standard diet over 20 weeks. QTL analysis was performed after genotyping with the MiniMUGA Genotyping Array. Whole-genome sequencing and gene expression data of the parental lines was used for the prioritization of positional candidate genes.

Results

Three QTL associated with liver weight, body weight, and subcutaneous adipose tissue (scAT) weight were identified. A highly significant QTL on chromosome (Chr) 1 (157–168 Mb) showed an association with liver weight. A QTL for body weight at 20 weeks was found on Chr 3 (34.1–40 Mb) overlapping with a QTL for scAT weight. In a multiple QTL mapping approach, an additional QTL affecting body weight at 16 weeks was identified on Chr 6 (9.5–26.1 Mb). Considering sequence variants and expression differences, Sec16b and Astn1 were prioritized as top positional candidate genes for the liver weight QTL on Chr 1; Met and Ica1 for the body weight QTL on Chr 6. Interestingly, all top candidate genes have previously been linked with metabolic traits.

Conclusions

This study shows once more the power of an advanced intercross line for fine mapping. QTL mapping combined with a detailed prioritization approach allowed us to identify additional and plausible candidate genes linked to metabolic traits in the BFMI861-S1xB6N AIL. By reidentifying known candidate genes in a different crossing population the causal link with specific traits is underlined and additional evidence is given for further investigations.

Obesity and its related pathologies such as insulin resistance, type 2 diabetes, and fatty liver are symptoms of an imbalanced energy homeostasis (1). A sedentary lifestyle as well as the (over)consumption of easily available energy-dense food contribute to this imbalance (1). However, genetic constitution sets the stage for the phenotypic characteristics. Genome-wide association studies (GWAS) in humans have revealed more than 300 single-nucleotide polymorphisms (SNPs) associated with obesity-related phenotypes (2). Nevertheless, identified loci and underlying genes explain only a minor proportion of the estimated heritability (2) (3).

Mice are ideal model organisms for studying genetic effects because the environmental conditions can be tightly controlled. Furthermore, different inbred mouse lines with distinct but well-defined genetic constitution are available, which can be used to improve our understanding of the genetic architecture of complex phenotypes like obesity. Crosses between inbred lines allow the generation of structured populations. These prerequisites make association studies between genetic loci and phenotypes feasible in relatively small populations with high statistical power(4). In particular, advanced intercross lines (AIL) allow high resolution QTL mapping by increasing recombination between any two loci (4) (5) (6) (7).

The Berlin Fat Mouse with its different inbred lines (BFMI) were generated initially by crossing various mice from different pet shops, subsequent selection for high body weight and high fat mass (8), and finally repeated inbreeding between pairs of full sibs to generate a mouse model for the investigation of body weight gain and body composition. The different BFMI sub-lines are genetically closely related. They are all obese, but show different features of the metabolic syndrome (9) (10).

Due to their unique genetic background in combination with the distinct obese phenotypes, the BFMI sublines allow the identification of diverse genetic contributors to the metabolic syndrome. In a cross between the obese line BFMI860-12 and C57BL/6N (B6N) as a lean strain, a major QTL for total fat mass was mapped on chromosome (Chr) 3 (jObes1) (11). The locus was further fine-mapped and characterized, which led to the identification of Bbs7 as causal gene (6), a gene that contributes to the Bardet-Biedl syndrome in humans. In an AIL between BFMI861-S1 and BFMI861-S2, two genetically very similar lines, we identified QTL and candidate genes responsible for differences in liver weight, liver triglycerides, gonadal adipose tissue weight, and body weight (7). The BFMI861-S1 line of this cross does not only carry the mutant jObes1 allele, it also shows the highest liver weight and liver triglyceride (TG) concentration among all BFMI sub-lines (9). Therefore, the BFMI861-S1 line is an interesting mouse model to study the genetic architecture of hepatic fat deposition in the context of obesity. In this study, we have generated the advanced intercross of BFMI861-S1 and B6N, where B6N is a lean counterpart to BFMI861-S1 to identify additional genes contributing to the specific phenotype of the BFMI861-S1 sub-line, in particular its fatty liver.

Mouse population

58 male and 65 female mice of the AIL BFMI861-S1xC57BL/6N (AIL BFMI861-S1xB6N) in generation 10 were genotyped and phenotyped. The AIL population was generated from an F2 population between an obese BFMI861-S1 male and lean B6N female. Beginning in generation F1, individuals were randomly mated to mice from the same generation using the program RandoMate (12).

Animal husbandry

All animal experiments were approved by the German Animal Welfare Authorities (approval no. G0099/16) and reported in accordance to ARRIVE guidelines. All methods were performed in accordance with the relevant guidelines and regulations.

Mice were maintained under conventional conditions and a 12:12 h light–dark cycle (lights on at 06:00) at a temperature of 22 ± 2°C. Animals had ad libitum access to food and water. Animals were fed a standard rodent diet containing 16.7 MJ/kg of metabolizable energy, 11% from fat, 26% from protein and 53% from carbohydrates (V1534-000, ssniff EF R/M; Ssniff Spezialdiäten GmbH, Soest, Germany).

Phenotyping

Animals were analyzed between the age of 4 (after weaning at 3 weeks) and 20 weeks with body weight being recorded weekly. To investigate glucose metabolism, an oral glucose tolerance test (GTT) (week 18) and an intraperitoneal insulin tolerance test (ITT) (week 20, 1 U insulin/kg body weight) were performed as previously described (13). The area under the curve (AUC) for blood glucose for GTT and ITT was calculated. At 20 weeks, the mice were anesthetized with isofluorane after a fasting period of 2 hours and sacrificed by decapitation (in the morning until 12 AM). Several tissues including gonadal adipose tissue (gonAT), subcutaneous inguinal adipose tissue (scAT), and liver were dissected, weighed, shock-frozen in liquid nitrogen and stored until further use at -80°C. Body length was measured (14) and the body mass index (BMI) was calculated using the DuBois Eq. (15). Liver TG were assessed as previously described (16). Plasma free fatty acids (FFA), cholesterol, and TG were measured as described in Schulz et al. 2011 (17). Plasma insulin and skeletal muscle fat % were measured as previously described (18).

Genotyping

DNA isolation was done by salt extraction (5 M NaCl with β-mercapthoethanol and proteinase K) and subsequent ethanol precipitation. Genotypes of all 123 mice were generated by Neogen GeneSeek (Lincoln, NE, USA) using the Mini Mouse Universal Genotyping Array (MiniMUGA; Illumina, San Diego, CA, USA). The MiniMUGA array contains probes targeting 10,171 known SNPs (markers) (19). Markers were removed when all genotypes were missing or when the marker was not segregating. In addition, to prevent spurious associations, we required that at least two of the genotype groups contained 10 observations each. In case one out of three genotype groups contained less than 10 individuals, these were set to N/A, but the marker was kept. After quality control 1,886 high quality markers were available and used for subsequent QTL analysis (supplementary Fig. 1). Genotypes are available in supplementary File 1.

QTL mapping

Linear models were used to investigate the influence of subfamily, litter size, and sex on each phenotype. Because the factors affected the phenotypes differently, different statistical models were used for mapping each phenotype, which included significant factors as fixed covariates accordingly.

To minimize the influence of population structure, the genomic inflation factor (λ) was computed. If the genomic inflation factor was above 1.05 (20), results were corrected using λ-correction. To account for multiple testing, significance thresholds were corrected using stringent Bonferroni correction. The number of independent SNPs was determined using the simpleM method (21). The threshold for significance was set using the number of independent SNPs (1,365) as the total number of tests performed. This resulted in a LOD score (as defined by -log10(pvalue)) after λ-correction above 5.1 to be deemed ‘genome-wide highly significant’ and above 4.4 to be ‘genome-wide significant’. QTL regions were defined by a 1.5 LOD drop from the top marker. Region start and end positions are defined by the first marker upstream and downstream, respectively, that have a drop of 1.5 from the LOD score of the top marker (22).

To discover additional QTLs for body weight, which might be hidden by the known strong effect QTL of the jObes1 locus, a variation of multiple QTL mapping (MQM) was used (23). The single QTL model was adjusted to compensate for the known effect of the jObes1 locus by including the top marker from the Chr 3 region (SNP gUNC5036315) as an additional cofactor into the model:

Body weight = sex + mother + gUNC5036315 + marker genotype + error

Gene expression analysis

Gene expression was measured in RNA isolated from liver of BFMI861-S1 male mice (n = 6) at 10 weeks. RNA was extracted as described in Hesse et al. 2018 (24). Gene expression was measured with the Clariom™ S Assay for mouse (Thermo Fisher Scientific) using the service of ATLAS Biolabs, Berlin, Germany. Gene expression data of male B6N mice (n = 5) measured with the Clariom™ S assay for mouse were downloaded from Gene Expression Omnibus (25). Probe intensities were log₂ transformed and quantile normalized. To test for expression differences between BFMI861-S1 and B6N mice, two-tailed t-tests were performed. False positives due to multiple testing were minimized using a Benjamini-Hochberg correction. For statistical analysis and for graphical presentation R: A Language and Environment for Statistical Computing (26) was used.

Whole genome sequencing

The BFMI861-S1 parental genome was paired-end sequenced using the Illumina HiSeq (Illumina) platform. Obtained DNA reads were trimmed using trimmomatic (27) after which trimmed reads were aligned to the mouse genome (MM10, GRCm38.p6) using the Burrows–Wheeler Aligner (BWA) software (28). Subsequently, SAM files were converted to BAM files, sorted, and indexed using Samtools (29). (Optical) Duplicate reads were removed using Picard tools v2.19.0 after which indel realignment and base recalibration was done using the GATK v4.1.0.0 (30), according to GATK best practices.

All sequence variants in BFMI861-S1 mice were called using BCFtools (29). Variants were further annotated using the Ensembl Variant Effect Predictor (VEP) (31).

Candidate genes prioritization

Prioritization of candidate genes in each QTL region was performed as described in Delpero et al. 2021 (7). In brief, genes in a QTL region containing sequence variants between the parental lines BFMI861-S1 and B6N were ranked according to the sum of scores for the functional annotation of coding and non-coding variants, gene expression data, and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Coding sequence variants leading to stop gain/stop loss codons and missense mutations located in functional protein domains were awarded a score of 3 points. A missense variant with either a deleterious or a tolerated SIFT (Sorting Intolerant From Tolerant) value obtained a score of 3 or 1, respectively. Non-coding variants were scored based on their location in potential functional sites. If a non-coding variant was located in the promoter or in a splice site, a score of 3 was awarded; if located in untranslated regions (UTRs), enhancers, or CTCF binding sites (involved in 3D structure of chromatin) the score was 1. Genes differentially expressed in the liver were scored with 2; genes annotated in relevant KEGG metabolic pathways with 1.

Phenotypic variation and correlation analysis in the AIL

Animals of the AIL population showed high standard deviation for the collected phenotypes which were expected and are needed for QTL analysis. In detail, body weight was on average 21.9 ± 7.27 g at the end of the experiment (week 20). GonAT weight, scAT, and liver weight were on average 1.53 ± 1.0 g, 0.65 ± 0.32 g, and 1.92 ± 0.55 g, respectively. The areas under the curve for blood glucose during GTT and ITT were 21,002 ± 12,415 and 7,308 ± 3,195, respectively. Liver triglycerides and plasma triglycerides were on average 112 ± 68 µg TG / µg protein and 896 ± 506 ug/ml, respectively. Additional plasma parameters such as plasma cholesterol, plasma FFA, and plasma insulin were on average 44 ± 9 mg/dl, 0.23 ± 0.06 mmol/l, and 7 ± 13 ng/mL, respectively. Skeletal muscle fat % also showed high standard deviation in the AIL population (mean = 19 ± 5%) (supplementary Table 1).

In order to assess the relationship among the phenotypes measured in all AIL mice, Spearman correlation was computed between all the collected phenotypes. Most of the phenotypes (scAT weight, gonAT weight, liver weight, body weight, GTT AUC, and ITT AUC, plasma cholesterol, plasma insulin, BMI, body length, and skeletal muscle fat %) were positively correlated among each other (Table 1).

No significant correlation was detected between gonAT weight and liver weight, gonAT weight and plasma TG, and scAT weight and plasma TG. Liver TG showed a significant positive correlation with body weight at 20 weeks (r = 0.34), gonAT weight (r = 0.32), scAT weight (r = 0.44), and skeletal muscle fat % (r = 0.45). In addition, plasma FFA did not show any correlation with the collected phenotypes.

QTL mapping

For QTL analysis, different statistical models were used for mapping each phenotype (Table 2). The results revealed genome-wide significant loci on three different chromosomes (1, 3, and 6) associated with one or more of the investigated phenotypes (Table 3). Additionally, a suggestive QTL associated with liver TG was found on Chr 8.

In detail, the significant QTL on Chr 1 (157,132,066–168,495,457) with a LOD score of 4.96 was associated with liver weight (Fig. 1a). This region contains 89 annotated protein-coding positional candidate genes. The most significant SNP for liver weight in this region was “gUNC2036998” (Chr1:158,663,689). Interestingly, this SNP showed only two genotype classes (homozygous BFMI861-S1 and heterozygous). The liver of homozygous mice carrying the BFMI861-S1 allele was 17% heavier compared to the liver of heterozygous (Het) mice (mean BFMI861-S1 = 1.81 ± 0.25 g, mean Het = 1.55 ± 0.42 g) (Fig. 1b).

The highly significant region for body weight on Chr 3 (34,066,622 − 40,043,158) corresponded with the jObes1 locus that was identified in BFMI mice before (6). This QTL effect in the AIL BFMI861-S1xB6N persisted at all time points starting from week 9 until week 20 (Fig. 2a). The most significant association (LOD = 8.89) was body weight at week 14 with the top marker gUNC5036315 (Chr3:35,986,311) (Fig. 2b). This marker was 604 kbp away from the Bbs7 gene that had been identified recently as causal gene for obesity in BFMI mice (6). At the top marker locus, 14 weeks-old mice homozygous for the BFMI861-S1 allele were 10.09 g heavier than homozygous B6N counterparts. The same region affected also scAT weight (LOD = 5.8), and BMI (LOD = 4.94) with homozygous BFMI861-S1 mice carrying 98% more scAT than B6N homozygous mice (mean BFMI861-S1 = 0.95 ± 0.24 g, mean B6N = 0.48 ± 0.21 g) and 83% compared to heterozygous mice (mean Het = 0.52 ± 0.27 g). In addition, homozygous BFMI861-S1 mice showed 14% BMI increase compared to B6N homozygous mice (mean BFMI861-S1 = 4.20 ± 0.20 kg/m2, mean B6N = 3.67 ± 0.31 kg/m2) and 12% compared to heterozygous mice (mean Het = 3.76 ± 0.33 kg/m2). This region contains 30 annotated protein-coding genes.

When correcting for the top marker of the jObes1 locus on Chr 3 (gUNC5036315), an additional region associated with body weight at 16 weeks was detected on Chr 6 (0–17,553,096). This region contains 59 protein-coding genes. The most significant SNP of this region was gUNC10595065 (3,919,413; LOD = 5.41) (Fig. 1c). Heterozygous mice showed 9% increase in body weight compared to homozygous B6N mice (mean Het = 45.14 ± 2.91 g, mean B6N = 41.47 ± 3.24 g) and 4.5% increase compared to homozygous BFMI861-S1 mice (mean BFMI861-S1 = 43.33 ± 2.7 g) (Fig. 1d).

A suggestive QTL for liver TG was identified on Chr 8 (86,158,420 − 106,738,488). Due to the suggestive significance, the region is large containing 179 protein-coding genes. The top marker in this region was “S1H083826428” (Chr8:95,660,710; LOD = 3.93). On average, homozygous mice carrying the BFMI861-S1 allele at this marker showed 90% increase amounts of liver TG compared to homozygous B6N mice (mean BFMI861-S1 = 192 ± 54 µg TG / µg protein, mean B6N = 101 ± 53 µg TG / µg protein) and 83% increase compared to heterozygous mice (mean Het = 105 ± 61 µg TG / µg protein).

Candidate gene prioritization

Within the confidence intervals of the four QTL (including the QTL on Chr 8 suggestively associated with liver TG) 357 protein coding positional candidate genes were located. 152 genes were polymorphic between BFMI861-S1 and B6N in protein-coding and/or regulatory regions; 29 on Chr 1, 38 on Chr 3, 22 on Chr 6, and 63 on Chr 8.

In order to identify the most likely candidate genes for each QTL, the 152 polymorphic positional candidate genes were scored according to the decision tree (supplementary Table 2). After applying the prioritization criteria, two genes (Astn1 and Sec16b) located in the region on Chr 1 associated with liver weight ranked with the highest score of 12 and 10, respectively (Table 4). Astn1 and Sec16b carry deleterious missense variants according to the variant effect predictor. In addition, both Sec16b and Astn1 show variants in the promoter region and 5-prime and 3-prime UTRs. However, despite variants in regulatory regions Sec16b and Astn1 did not show gene expression differences in the liver.

Three genes on Chr 3 (Frem2, Bbs7, and Noct), with a score of 13, 12, and 12, respectively ranked as top candidate genes. Among the candidate genes located in the region on Chr 3 (jObes1) associated with body weight from week 9 to 20, Bbs7, Noct, and Frem2 all carried missense variants in domains and regulatory region variants. In addition, Bbs7 and Noct were both downregulated (p = 0.01346 and p = 0.00876, respectively) in liver of BFMI861-S1 mice compare to B6N, while Frem2 did not show differences in the expression.

In the region on Chr 6 associated with body weight two genes (Met and Ica1) with scores of 10 and 9, respectively ranked as top candidate genes. Met carries a tolerated missense variant in the IPT (Ig-like, plexins, transcription factors) domain and variants in regulatory regions such as enhancers and untranslated regions. Ica1 showed variants only in regulatory regions (promoter, CTCF binds, enhancers, and untranslated regions). Met and Ica1 were both downregulated in the liver of BFMI861-S1 mice (p = 8.2E-06 and 0.00200, respectively).

In the region affecting liver TG on Chr 8, Fto and Lpcat2 were identified as the top candidate genes with scores of 16 and 13, respectively. Fto carried a stop gain variant and additional variants in different regulatory regions (promoter, CTCF binds, and enhancer) and was downregulated in the liver of BFMI861-S1 mice (p = 0.01092). Lpcat2 instead carried one deleterious missense variant and regulatory region variants (promoter and enhancer) but did not show expression differences in the liver between BFMI861-S1 and B6N mice.

Genome wide association studies (GWAS) on obesity-associated phenotypes identified loci that account in sum only for a small percentage of the total variance of the examined population (2). Therefore, studies are needed that better allow the identification of genetic effects than most populations do. In order to unravel the genetics behind obesity and hepatic fat deposition, we investigated an AIL population generated from a cross between the obese BFMI861-S1 mouse line and the lean reference strain B6N.

The AIL population used in this study has the advantage of having a high resolution for QTL mapping. Because the examined AIL accumulated recombination over 10 generations, the physical length of the QTL regions is relatively short and as such the number of positional candidate genes is low (4). In our population, the number of positional candidate genes could be further reduced by removing regions in the genome that are identical between BFMI861-S1 and B6N. Excluding non-polymorphic genes reduced the number of protein coding candidate genes from 357 to 152. In addition, the application of the decision tree led to the prioritization of the most likely candidate genes among the 152 polymorphic genes.

In the region on Chr 1 associated with liver weight, the top candidate genes are Sec16b and Astn1. Sec16b is required for secretory cargo traffic from the endoplasmic reticulum to the Golgi apparatus (32). Previously, the gene has been linked to increased fat storage in both mice and humans. A human GWAS associated Sec16b with differences in body composition (33). In mice, dysfunctional Sec16b was associated with increased body weight (34). The BFMI861-S1 mice of our study carry a deleterious missense variant of Sec16b leading to an impaired protein variant and, in addition, SNPs in the promoter region and 3-prime and 5-prime UTRs. Since Sec16b was not differentially expressed between the parental lines BFMI861-S1 and B6N, we hypothesize that the detected deleterious missense variant is responsible for the dysfunction of the encoded protein. Therefore, we consider Sec16b as a very strong candidate responsible for the increased liver weight in BFMI861-S1 mice. Astn1, the second prioritized candidate gene in the Chr 1 region, is a neuronal adhesion molecule required for the migration of young postmitotic neuroblasts along glial fibers (35). The gene has not been associated with obesity related phenotypes yet. In BFMI861-S1 mice, Astn1 carries one deleterious missense variant in the fibronectin type 3 domain. This domain is known to be responsible for interactions with other extracellular matrix (ECM) or cell surface proteins (36). Therefore, the identified deleterious mutation could reduce the interaction ability of the encoded protein. Additional SNPs between the BFMI861-S1 and B6N in the promoter, 5-prime UTR, enhancers, and CTCF binding sites of Astn1 could be responsible for the observed downregulation of the gene in the liver of BFMI861-S1 mice versus B6N. This finding confirms that the selected prioritization approach is useful to identify so far unknown candidate genes for the phenotype under investigation which should be considered for follow-up studies.

In the QTL interval for body weight on Chr 3, the prioritization approach identified Frem2, Bbs7, and Noct as the top candidate genes. Bbs7 has been identify to be the causal gene for elevated fat mass and obesity in all BFMI lines before (6). Among diverse sequence variants between BFMI and B6N in the Bbs7 gene, it has been clarified that a large intronic deletion is partially responsible for the high fat content in BFMI mice (37). This is additional evidence for the prioritization approach being suitable for the correct identification of positional candidate genes in a defined confidence interval.

Only if we accounted for the strong effect of jObes1 by including the top marker from the Chr 3 region (gUNC5036315) as a cofactor into the model we detected another QTL for body weight on Chr 6. This QTL affecting body weight at 16 weeks partially overlaps with a previously identified QTL for body weight at 10 weeks that was identified in the F2 population BFMI860-12xB6N, which used BFMI860-12 instead of BFMI861-S1 from the current study (11). In the confidence interval of the QTL identified in our cross, Met and Ica1 are the highest scored top candidate genes. Met encodes a receptor tyrosine kinase involved in the transduction of signals from the extracellular matrix to the cytoplasm by binding to hepatocyte growth factor/HGF ligand (38). In BFMI861-S1 mice Met carries a tolerated missense variant in the Ig-like, plexins, transcription factors domain which is involved in the control of cell dissociation, motility, and invasion of extracellular matrices (39). In addition, Met carries SNPs in regulatory regions (enhancers and CTCF binds) which could contribute to the observed downregulation in the liver of BFMI861-S1 mice compare to B6N. According to literature, Met is known to be involved in pancreatic-cell death and diabetes (34) (40). Therefore, we consider this gene to be an interesting candidate for the increased body weight of heterozygous mice compared to homozygous BFMI861-S1 and homozygous B6N mice. Ica1 which is the second prioritized gene of the Chr 6 region encodes for Islet cell autoantigen 1 and plays a role in neurotransmitter secretion (41). Ica1 carries several SNPs in the promoter and 5-prime UTR which might cause the downregulation in liver of BFMI861-S1 mice. Ica1 is also known to be associated with type 1 diabetes mellitus in non-obese diabetic mice (42) as well as with glucose homeostasis (34). Therefore, Ica1 could cause the higher body weight that we observed in heterozygous mice compared to homozygous BFMI861-S1 and homozygous B6N mice.

In the confidence interval of the suggestive liver TG QTL on Chr 8, Fto and Lpcat2 are the most likely candidate genes responsible for hepatic fat accumulation in BFMI861-S1 mice. The same genomic region has previously been associated with liver TG in the Collaborative Cross (43). In this QTL, we suggested Fto to be responsible for the increased amount of TG in the liver of BFMI861-S1 mice. Fto has been associated with the body mass index in humans and therefore this gene has been examined extensively (44) (45). Fto codes for an RNA demethylase. In mice, both Fto knockout and Fto missense variants are responsible for fat accumulation and hypertriglyceridemia (46). The BFMI861-S1 allele carriers of our study carry a stop gain variant at FTO amino acid position 314, which leads to a premature stop codon, and thereby to a shortened protein. Furthermore, BFMI861-S1 mice contain additional missense and promoter variants. The BFMI861-S1 Fto haplotype likely contributes as a whole to the observed downregulation in the liver of BFMI861-S1 mice. These findings led us to consider Fto as the main contributor to the hepatic fat accumulation in BFMI861-S1 mice. Lpcat2, the other prioritized candidate gene is an acetyltransferase (47). In BFMI861-S1 mice Lpcat2 carries a missense variant at amino acid position 59 located in the transmembrane helical domain which could affect the function of the protein. In addition, Lpcat2 carries variants in regulatory regions but is not differentially expressed in liver of BFMI861-S1 mice compare to B6N. According to literature, Lpcat proteins are associated with polyunsaturated fatty acid accumulation (48). Moreover, knockdown of both Lpcat1 and Lpcat2 leads to an increase in lipid droplets size (49). The occurrence of several SNPs in functional regions of Lpcat2 in BFMI861-S1 mice and its known function in fat accumulation let us to consider this gene as another potential contributor to liver TG accumulation.

In conclusion, our approach identified strong candidate genes that are likely involved in the development of obesity and fatty liver disease in our Berlin Fat Mouse model. However, although we have prioritized candidate genes using the available information, we cannot completely rule out that one of the other polymorphic genes was wrongly discarded, or that non-protein coding regions might be causal. Nevertheless, the natural mutations found in this study in the Berlin Fat Mouse Inbred line BFMI861-S1 contribute to our understanding of which genes impact obesity and hepatic fat storage. These findings help to clarify and support the role of known candidates. The examined mouse model and the applied gene prioritization approach allow to unravel the effects of the identified QTL regions and to link genes with observed phenotypes. Additional studies on the candidate genes should be performed to discover by which molecular mechanisms they contribute to the development of obesity and liver associated diseases in mouse models but also in humans.

Data availability

DNA sequencing data were deposited at the NCBI Sequence Read Archive (SRA) under BioProject ID: PRJNA717237 and is available at: https://www.ncbi.nlm.nih.gov/bioproject/717237.

BFMI861-S1 gene expression data measured with Clariom™ S Assay for mouse were deposited at the ArrayExpress Archive (accession E-MTAB-11512).

Acknowledgement

DH was funded by the Deutsche Forschungsgemeinschaft (DFG, HE8165/1-1). The project was supported by the DDG (Deutsche Diabetes Gesellschaft). We thank Marion Bütow, Ulf Kiessling, and Ines Walter for assistance in animal care and Maximilian Sprechert for technical assistance.

Authors contribution

MD: Data collection, data analysis, and manuscript draft. DA: Method discussion, data analysis, and manuscript revision. AF: Data collection. GAB: Data discussion, supervision, and manuscript feedback. DH: Experimental design, supervision, method discussion, and manuscript draft, feedback, and revision.

Blüher M. Obesity: global epidemiology and pathogenesis. Nat Rev Endocrinol. 2019 May 27;15(5):288–98. Available from: http://www.nature.com/articles/s41574-019-0176-8
Goodarzi MO. Genetics of obesity: what genetic association studies have taught us about the biology of obesity and its complications. Lancet Diabetes Endocrinol. 2018 Mar;6(3):223–36. Available from: https://linkinghub.elsevier.com/retrieve/pii/S2213858717302000
Albuquerque D, Nóbrega C, Manco L, Padez C. The contribution of genetics and environment to obesity. Br Med Bull. 2017 Sep 1;123(1):159–73. Available from: https://academic.oup.com/bmb/article/123/1/159/3930933
Darvasi A, Soller M. Advanced intercross lines, an experimental population for fine genetic mapping. Genetics. 1995 Nov;141(3):1199–207. Available from: http://www.ncbi.nlm.nih.gov/pubmed/8582624
Rockman M V., Kruglyak L. Breeding Designs for Recombinant Inbred Advanced Intercross Lines. Genetics. 2008 Jun;179(2):1069–78. Available from: http://www.genetics.org/lookup/doi/10.1534/genetics.107.083873
Arends D, Heise S, Kärst S, Trost J, Brockmann GA. Fine mapping a major obesity locus (jObes1) using a Berlin Fat Mouse × B6N advanced intercross population. Int J Obes. 2016 Nov;40(11):1784–8.
Delpero M, Arends D, Sprechert M, Krause F, Kluth O, Schürmann A, et al. Identification of four novel QTL linked to the metabolic syndrome in the Berlin Fat Mouse. Int J Obes. 2021 Oct 23; Available from: https://www.nature.com/articles/s41366-021-00991-3
Hantschel C, Wagener A, Neuschl C, Teupser D, Brockmann GA. Features of the metabolic syndrome in the Berlin Fat Mouse as a model for human obesity. Obes Facts. 2011;4(4):270–7. Available from: http://www.ncbi.nlm.nih.gov/pubmed/21921649
Heise S, Trost J, Arends D, Wirth E, Schäfer N, Köhrle J, et al. High Variability of Insulin Sensitivity in Closely Related Obese Mouse Inbred Strains. Exp Clin Endocrinol Diabetes. 2016 Jul;124(09):519–28.
Wagener A, Schmitt AO, Aksu S, Schlote W, Neuschl C, Brockmann GA. Genetic, sex, and diet effects on body weight and obesity in the Berlin Fat Mouse Inbred lines. Physiol Genomics. 2006 Dec;27(3):264–70. Available from: https://www.physiology.org/doi/10.1152/physiolgenomics.00225.2005
Neuschl C, Hantschel C, Wagener A, Schmitt AO, Illig T, Brockmann GA. A unique genetic defect on chromosome 3 is responsible for juvenile obesity in the Berlin Fat Mouse. Int J Obes. 2010 Dec;34(12):1706–14.
Schmitt AO, Bortfeldt R, Neuschl C, Brockmann GA. RandoMate: a program for the generation of random mating schemes for small laboratory animals. Mamm Genome. 2009 May;20(5):321–5. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19444508
Hesse D, Jaschke A, Kanzleiter T, Witte N, Augustin R, Hommel A, et al. GTPase ARFRP1 Is Essential for Normal Hepatic Glycogen Storage and Insulin-Like Growth Factor 1 Secretion. Mol Cell Biol. 2012 Nov 1;32(21):4363–74. Available from: http://mcb.asm.org/cgi/doi/10.1128/MCB.00522-12
Hesse D, Dunn M, Heldmaier G, Klingenspor M, Rozman J. Behavioural mechanisms affecting energy regulation in mice prone or resistant to diet- induced obesity. Physiol Behav. 2010 Mar;99(3):370–80. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0031938409003874
Gargiulo S, Gramanzini M, Megna R, Greco A, Albanese S, Manfredi C, et al. Evaluation of Growth Patterns and Body Composition in C57Bl/6J Mice Using Dual Energy X-Ray Absorptiometry. Biomed Res Int. 2014;2014:1–11. Available from: http://www.hindawi.com/journals/bmri/2014/253067/
Hesse D, Radloff K, Jaschke A, Lagerpusch M, Chung B, Tailleux A, et al. Hepatic trans-Golgi action coordinated by the GTPase ARFRP1 is crucial for lipoprotein lipidation and assembly. J Lipid Res. 2014 Jan;55(1):41–52. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24186947
Schulz N, Himmelbauer H, Rath M, van Weeghel M, Houten S, Kulik W, et al. Role of Medium- and Short-Chain L-3-Hydroxyacyl-CoA Dehydrogenase in the Regulation of Body Weight and Thermogenesis. Endocrinology. 2011 Dec 1;152(12):4641–51. Available from: https://academic.oup.com/endo/article/152/12/4641/2457348
Heise S, Trost J, Arends D, Wirth EK, Schäfer N, Köhrle J, et al. High Variability of Insulin Sensitivity in Closely Related Obese Mouse Inbred Strains. Exp Clin Endocrinol Diabetes. 2016 Oct;124(9):519–28. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27437914
Sigmon JS, Blanchard MW, Baric RS, Bell TA, Brennan J, Brockmann GA, et al. Content and Performance of the MiniMUGA Genotyping Array: A New Tool To Improve Rigor and Reproducibility in Mouse Research. Genetics. 2020 Dec 1;216(4):905–30. Available from: https://academic.oup.com/genetics/article/216/4/905/6065869
Devlin B, Roeder K. Genomic control for association studies. Biometrics. 1999 Dec;55(4):997–1004. Available from: http://www.ncbi.nlm.nih.gov/pubmed/11315092
Gao X, Becker LC, Becker DM, Starmer JD, Province MA. Avoiding the high Bonferroni penalty in genome-wide association studies. Genet Epidemiol. 2010;34(1):100–5.
Dupuis J, Siegmund D. Statistical methods for mapping quantitative trait loci from a dense set of markers. Genetics. 1999 Jan;151(1):373–86. Available from: http://www.ncbi.nlm.nih.gov/pubmed/9872974
Arends D, Prins P, Jansen RC, Broman KW. R/qtl: high-throughput multiple QTL mapping: Fig. 1. Bioinformatics. 2010 Dec 1;26(23):2990–2. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btq565
Hesse D, Trost J, Schäfer N, Schwerbel K, Hoeflich A, Schürmann A, et al. Effect of adipocyte-derived IGF-I on adipose tissue mass and glucose metabolism in the Berlin Fat Mouse. Growth Factors. 2018 Mar 4;36(1–2):78–88. Available from: https://www.tandfonline.com/doi/full/10.1080/08977194.2018.1497621
Clough E, Barrett T. The Gene Expression Omnibus Database. In 2016. p. 93–110. Available from: http://link.springer.com/10.1007/978-1-4939-3578-9_5
R Core Team (2018). R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014 Aug 1;30(15):2114–20. Available from: https://academic.oup.com/bioinformatics/article-lookup/doi/10.1093/bioinformatics/btu170
Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009 Jul;25(14):1754–60.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009 Aug 15;25(16):2078–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/19505943
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010 Sep;20(9):1297–303. Available from: http://www.ncbi.nlm.nih.gov/pubmed/20644199
McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GRS, Thormann A, et al. The Ensembl Variant Effect Predictor. Genome Biol. 2016 Dec 6;17(1):122. Available from: http://genomebiology.biomedcentral.com/articles/10.1186/s13059-016-0974-4
Bhattacharyya D, Glick BS. Two Mammalian Sec16 Homologues Have Nonredundant Functions in Endoplasmic Reticulum (ER) Export and Transitional ER Organization. Barr F, editor. Mol Biol Cell. 2007 Mar;18(3):839–49. Available from: https://www.molbiolcell.org/doi/10.1091/mbc.e06-08-0707
Sahibdeen V, Crowther NJ, Soodyall H, Hendry LM, Munthali RJ, Hazelhurst S, et al. Genetic variants in SEC16B are associated with body composition in black South Africans. Nutr Diabetes. 2018 Dec 19;8(1):43. Available from: http://www.nature.com/articles/s41387-018-0050-0
Bult CJ, Blake JA, Smith CL, Kadin JA, Richardson JE, Anagnostopoulos A, et al. Mouse Genome Database (MGD) 2019. Nucleic Acids Res. 2019 Jan 8;47(D1):D801–6. Available from: https://academic.oup.com/nar/article/47/D1/D801/5165331
Zheng C, Heintz N, Hatten ME. CNS Gene Encoding Astrotactin, Which Supports Neuronal Migration Along Glial Fibers. Science (80-). 1996 Apr 19;272(5260):417–9. Available from: https://www.science.org/doi/10.1126/science.272.5260.417
Kim H-S, McKnite A, Xie Y, Christian JL. Fibronectin type III and intracellular domains of Toll-like receptor 4 interactor with leucine-rich repeats (Tril) are required for developmental signaling. Bronner M, editor. Mol Biol Cell. 2018 Mar;29(5):523–31. Available from: https://www.molbiolcell.org/doi/10.1091/mbc.E17-07-0446
Florian Krause, Kourosh Mohebian, Manuel Delpero, Deike Hesse, Ralf Kühn, Danny Arends, Dr. Gudrun A. Brockmann (in press). A deletion containing a CTCF-element in intron 8 of the Bbs7 gene is partially responsible for juvenile obesity in the Berlin Fat Mouse. Mamm Genome. 2021;
Bladt F, Riethmacher D, Isenmann S, Aguzzi A, Birchmeier C. Essential role for the c-met receptor in the migration of myogenic precursor cells into the limb bud. Nature. 1995 Aug;376(6543):768–71. Available from: http://www.nature.com/articles/376768a0
Collesi C, Santoro MM, Gaudino G, Comoglio PM. A splicing variant of the RON transcript induces constitutive tyrosine kinase activity and an invasive phenotype. Mol Cell Biol. 1996 Oct;16(10):5518–26. Available from: https://journals.asm.org/doi/10.1128/MCB.16.10.5518
Mellado-Gil J, Rosa TC, Demirci C, Gonzalez-Pertusa JA, Velazquez-Garcia S, Ernst S, et al. Disruption of Hepatocyte Growth Factor/c-Met Signaling Enhances Pancreatic -Cell Death and Accelerates the Onset of Diabetes. Diabetes. 2011 Feb 1;60(2):525–36. Available from: http://diabetes.diabetesjournals.org/cgi/doi/10.2337/db09-1305
Pilon M, Peng X-R, Spence AM, Plasterk RHA, Dosch H-M. The Diabetes Autoantigen ICA69 and Its Caenorhabditis elegans Homologue, ric-19, Are Conserved Regulators of Neuroendocrine Secretion. Kimble J, editor. Mol Biol Cell. 2000 Oct;11(10):3277–88. Available from: https://www.molbiolcell.org/doi/10.1091/mbc.11.10.3277
Bonner SM, Pietropaolo SL, Fan Y, Chang Y, Sethupathy P, Morran MP, et al. Sequence Variation in Promoter of Ica1 Gene, Which Encodes Protein Implicated in Type 1 Diabetes, Causes Transcription Factor Autoimmune Regulator (AIRE) to Increase Its Binding and Down-regulate Expression. J Biol Chem. 2012 May;287(21):17882–93. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0021925820503086
Abu-Toamih Atamni HJ, Botzman M, Mott R, Gat-Viks I, Iraqi FA. Mapping novel genetic loci associated with female liver weight variations using Collaborative Cross mice. Anim Model Exp Med. 2018 Sep;1(3):212–20. Available from: http://doi.wiley.com/10.1002/ame2.12036
Lan N, Lu Y, Zhang Y, Pu S, Xi H, Nie X, et al. FTO – A Common Genetic Basis for Obesity and Cancer. Front Genet. 2020 Nov 16;11. Available from: https://www.frontiersin.org/articles/10.3389/fgene.2020.559138/full
Zhao X, Yang Y, Sun B-F, Zhao Y-L, Yang Y-G. FTO and Obesity: Mechanisms of Association. Curr Diab Rep. 2014 May 14;14(5):486. Available from: http://link.springer.com/10.1007/s11892-014-0486-0
Church C, Lee S, Bagg EAL, McTaggart JS, Deacon R, Gerken T, et al. A Mouse Model for the Metabolic Effects of the Human Fat Mass and Obesity Associated FTO Gene. Barsh GS, editor. PLoS Genet. 2009 Aug 14;5(8):e1000599. Available from: https://dx.plos.org/10.1371/journal.pgen.1000599
Shindou H, Hishikawa D, Nakanishi H, Harayama T, Ishii S, Taguchi R, et al. A Single Enzyme Catalyzes Both Platelet-activating Factor Production and Membrane Biogenesis of Inflammatory Cells. J Biol Chem. 2007 Mar;282(9):6532–9. Available from: https://linkinghub.elsevier.com/retrieve/pii/S0021925820522631
Harayama T, Eto M, Shindou H, Kita Y, Otsubo E, Hishikawa D, et al. Lysophospholipid Acyltransferases Mediate Phosphatidylcholine Diversification to Achieve the Physical Properties Required In Vivo. Cell Metab. 2014 Aug;20(2):295–305. Available from: https://linkinghub.elsevier.com/retrieve/pii/S1550413114002277
Moessinger C, Klizaite K, Steinhagen A, Philippou-Massier J, Shevchenko A, Hoch M, et al. Two different pathways of phosphatidylcholine synthesis, the Kennedy Pathway and the Lands Cycle, differentially regulate cellular triacylglycerol storage. BMC Cell Biol. 2014 Dec 10;15(1):43. Available from: https://bmcmolcellbiol.biomedcentral.com/articles/10.1186/s12860-014-0043-3

Table 1 to 4 is only available as a download in the Supplemental Files section.

No competing interests reported.

Table1.xlsx
Spearman correlation coefficients between the collected phenotypes in the AIL population. Bold indicates significant correlation after multiple testing correction (p < 9.10E-04). Abbreviations: gonAT, gonadal adipose tissue; scAT, subcutaneous adipose tissue; TG, triglycerides; ITT, insulin tolerance test; GTT, glucose tolerance test; AUC, area under the curve; FFA, free fatty acids; BMI, body mass index; SMuscle fat %, skeletal muscle fat percentage.
Table2.xlsx
values for effects of covariates on each phenotype.Abbreviations: gonAT, gonadal adipose tissue; scAT, subcutaneous adipose tissue; Gluc, blood glucose concentration; ITT, insulin tolerance test; GTT, glucose tolerance test; AUC, area under the curve; FFA, free fatty acids; TG, triglycerides; BMI, body mass index; SMuscle fat %, skeletal muscle fat percentage.
Table3.xlsx
QTL identified for different phenotypes in the AIL BFMI861-S1xB6N population.Abbreviations: ScAT, subcutaneous adipose tissue; BMI, body mass index; QTL, quantitative phenotypes locus; Chr, chromosome number; Start, Top, and Stop, position of the start of the QTL confidence interval, position of the SNP with the highest LOD score, and position of the end of the QTL confidence interval in base pairs, respectively. Positions are given according to the Mouse Genome Version MM10, GRCm38.p6. SNP, single-nucleotide polymorphism. The conﬁdence interval gives the 1.5 LOD drop region of the top SNP position. A LOD score above 5.1 was deemed to be ‘genome-wide highly significant’ and above 4.4 was deemed ‘genome-wide significant’. BH, Bonferroni correction; LOD, logarithm (base 10) of odds; Var %, percentage of total variance in AIL explained. The Mean columns show the phenotypic mean adjusted for significant covariates of homozygous BFMI861-S1, heterozygous, and homozygous B6N animals, respectively.
Table4.xlsx
Top candidate genes after applying the prioritization criteria. Bold indicates differentially expressed in liver after Benjamini-Hochberg correction.Abbreviations: scAT, subcutaneous adipose tissue; TG, triglycerides; BMI, body mass index; Chr, chromosome; FC, fold change.
supplementaryFigure1.png
supplementaryTable1.xlsx
supplementaryTable2.xlsx

Download PDF

Editorial decision: Major revision
09 May, 2022
Reviews received at journal
07 May, 2022
Reviewers agreed at journal
24 Apr, 2022
Reviewers invited by journal
03 Apr, 2022
Editor assigned by journal
03 Apr, 2022
Editor invited by journal
09 Mar, 2022
Submission checks completed at journal
09 Mar, 2022
First submitted to journal
08 Feb, 2022

You are reading this latest preprint version

QTL-mapping in the obese Berlin Fat Mouse identifies additional candidate genes for obesity and fatty liver disease

Status:

Version 1

Abstract

Background

Methods

Results

Conclusions

Figures

Introduction

Material And Methods

Mouse population

Animal husbandry

Phenotyping

Genotyping

QTL mapping

Gene expression analysis

Whole genome sequencing

Candidate genes prioritization

Results

Phenotypic variation and correlation analysis in the AIL

QTL mapping

Candidate gene prioritization

Discussion And Conclusion

Declarations

References

Tables

Additional Declarations

Supplementary Files

Status:

Version 1