Searching for genomic region of high-fat diet-induced type 2 diabetes in mouse chromosome 2 by analysis of congenic strains

Analysis of a genetic region affecting mouse body weight

Genetic factors affect an individual's risk of developing obesity, but in most cases each genetic variant has a small effect. Discovery of genes that regulate obesity may provide clues about its underlying biological processes and point to new ways the disease can be treated. Preclinical animal models facilitate genetic discovery in obesity because environmental factors can be better controlled compared with the human population. We studied inbred mouse strains to identify novel genes affecting obesity and glucose metabolism. BTBR T+ Itpr3^tf/J (BTBR) mice are fatter and more glucose intolerant than C57BL/6J (B6) mice. Prior genetic studies of these strains identified an obesity locus on chromosome 2. Using congenic mice, we found that obesity was affected by a ∼316 kb region, with only two known genes, pyruvate dehydrogenase kinase 1 (Pdk1) and integrin α 6 (Itga6). Both genes had mutations affecting their amino acid sequence and reducing mRNA levels. Both genes have known functions that could modulate obesity, lipid metabolism, insulin secretion, and/or glucose homeostasis. We hypothesized that genetic variation in or near Pdk1 or Itga6 causing reduced Pdk1 and Itga6 expression would promote obesity and impaired glucose tolerance. We used knockout mice lacking Pdk1 or Itga6 fed an obesigenic diet to test this hypothesis. Under the conditions we studied, we were unable to detect an individual contribution of either Pdk1 or Itga6 to body weight. During our studies, with conditions outside our control, we were unable to reproduce some of our previous body weight data. However, we identified a previously unknown role for Pdk1 in cardiac cholesterol metabolism providing the basis for future investigations. The studies described in this paper highlight the importance and the challenge using physiological outcomes to study obesity genes in mice.

Evolution of the Membrane Transport Protein Domain

Membrane transport proteins are widely present in all living organisms, however, their function, transported substrate, and mechanism of action are unknown. Here we use diverse bioinformatics tools to investigate the evolution of MTPs, analyse domain organisation and loop topology, and study the comparative alignment of modelled 3D structures. Our results suggest a high level of conservancy between MTPs from different taxa on both amino acids and structural levels, which imply some degree of functional similarities. The presence of loop/s of different lengths in various positions suggests tax-on-specific adaptation to transported substrates, intracellular localisation, accessibility for post-translation modifications, and interaction with other proteins. The comparison of modelled structures proposes close relations and a common origin for MTP and Na/H exchanger. Further, a high level of amino acid similarity and identity between archaeal and bacterial MTPs and Na/H exchangers imply conservancy of ion transporting function at least for archaeal and bacterial MTPs.

Transgenerational effects of obesogenic diets in rodents: A meta-analysis

Obesity is a major health condition that affects millions worldwide. There is an increased interest in understanding the adverse outcomes associated with obesogenic diets. A multitude of studies have investigated the transgenerational impacts of maternal and parental obesogenic diets on subsequent generations of offspring, but results have largely been mixed. We conducted a systematic review and meta-analysis on rodent studies to elucidate how obesogenic diets impact the mean and variance of grand-offspring traits. Our study focused on transgenerational effects (i.e., F2 and F3 generations) in one-off and multigenerational exposure studies. From 33 included articles, we obtained 407 effect sizes representing pairwise comparisons of control and treatment grand-offspring groups pertaining to measures of body weight, adiposity, glucose, insulin, leptin, and triglycerides. We found evidence that male and female grand-offspring descended from grandparents exposed to an obesogenic diet displayed phenotypes consistent with metabolic syndrome, especially in cases where the obesogenic diet was continued across generations. Further, we found stronger evidence for the effects of grand-maternal than grand-paternal exposure on grand-offspring traits. A high-fat diet in one-off exposure studies did not seem to impact phenotypic variation, whereas in multigenerational exposure studies it reduced variation in several traits.

Identification and analysis of circulating long non-coding RNAs with high significance in diabetic cardiomyopathy

Diabetic cardiomyopathy (DCM) lacks diagnostic biomarkers. Circulating long non-coding RNAs (lncRNAs) can serve as valuable diagnostic biomarkers in cardiovascular disease. To seek potential lncRNAs as a diagnostic biomarker for DCM, we investigated the genome-wide expression profiling of circulating lncRNAs and mRNAs in type 2 diabetic db/db mice with and without DCM and performed bioinformatic analyses of the deregulated lncRNA-mRNA co-expression network. Db/db mice had obesity and hyperglycemia with normal cardiac function at 6 weeks of age (diabetes without DCM) but with an impaired cardiac function at 20 weeks of age (DCM) on an isolated Langendorff apparatus. Compared with the age-matched controls, 152 circulating lncRNAs, 127 mRNAs and 3355 lncRNAs, 2580 mRNAs were deregulated in db/db mice without and with DCM, respectively. The lncRNA-mRNA co-expression network analysis showed that five deregulated lncRNAs, XLOC015617, AK035192, Gm10435, TCR-α chain, and MouselincRNA0135, have the maximum connections with differentially expressed mRNAs. Bioinformatic analysis revealed that these five lncRNAs were highly associated with the development and motion of myofilaments, regulation of inflammatory and immune responses, and apoptosis. This finding was validated by the ultrastructural examination of myocardial samples from the db/db mice with DCM using electron microscopy and changes in the expression of myocardial tumor necrosis factor-α and phosphorylated p38 mitogen-activated protein kinase in db/db mice with DCM. These results indicate that XLOC015617, AK035192, Gm10435, TCR-α chain, and MouselincRNA0135 are crucial circulating lncRNAs in the pathogenesis of DCM. These five circulating lncRNAs may have high potential as a diagnostic biomarker for DCM.

Pierre CL, Macias‐Velasco JF, Hughes JW, Kunzmann M, Schmidt H, Wayhart JP, Lawson HA. Spontaneous restoration of functional β-cell mass in obese SM/J mice

Maintenance of functional β-cell mass is critical to preventing diabetes, but the physiological mechanisms that cause β-cell populations to thrive or fail in the context of obesity are unknown. High fat-fed SM/J mice spontaneously transition from hyperglycemic-obese to normoglycemic-obese with age, providing a unique opportunity to study β-cell adaptation. Here, we characterize insulin homeostasis, islet morphology, and β-cell function during SM/J's diabetic remission. As they resolve hyperglycemia, obese SM/J mice dramatically increase circulating and pancreatic insulin levels while improving insulin sensitivity. Immunostaining of pancreatic sections reveals that obese SM/J mice selectively increase β-cell mass but not α-cell mass. Obese SM/J mice do not show elevated β-cell mitotic index, but rather elevated α-cell mitotic index. Functional assessment of isolated islets reveals that obese SM/J mice increase glucose-stimulated insulin secretion, decrease basal insulin secretion, and increase islet insulin content. These results establish that β-cell mass expansion and improved β-cell function underlie the resolution of hyperglycemia, indicating that obese SM/J mice are a valuable tool for exploring how functional β-cell mass can be recovered in the context of obesity.

Diethylnitrosamine enhances hepatic tumorigenic pathways in mice fed with high fat diet (Hfd)

Obesity has been implicated in the genesis of metabolic syndromes including insulin resistance and Type 2 Diabetes Mellitus (T2DM). Given the association between T2DM and the risk of hepatocellular carcinoma (HCC), our specific goal was to determine whether the liver of HFD-induced T2DM mice is more sensitive to the carcinogen diethylnitrosamine (DEN), due to a modification of the molecular pathways implicated in the early stages of HCC pathogenesis. C57BL/6 male mice (five-week-old) were divided into 4 groups: C, C + DEN, HFD and HFD + DEN. Mice were euthanized twenty-five weeks after DEN-injection. Livers of HDF-fed mice showed a higher proliferative index than Control groups. In line with this, HFD groups showed an increase of nuclear β-catenin, and interestingly, DEN treatment led to a slight increase in the expression of this protein in HFD group. Based on these results, and to confirm this effect, we analyzed β-catenin target genes, finding that DEN treatment in HFD group led to a significant increase of Vegf, c-myc, c-jun and cyclin D1 expression levels. According to our results, the expression of TCF4 showed to be significantly increased in HFD + DEN vs. HFD. In this regard, the β-catenin/TCF4 complex enhanced its association with pSmads 2/3, as we observed an increase of nuclear Smads expression in HFD + DEN, suggesting a possible role of TGF-β1/Smads signaling pathway in this phenomenon. Our results show that the liver of HFD fed model that resembles early T2DM pathology in mice, is more sensitive to DEN, by inducing both Wnt/β-catenin and TGF β1/Smads tumorigenic pathways.

Burly1 is a mouse QTL for lean body mass that maps to a 0.8-Mb region of chromosome 2

To fine map a mouse QTL for lean body mass (Burly1), we used information from intercross, backcross, consomic, and congenic mice derived from the C57BL/6ByJ (host) and 129P3/J (donor) strains. The results from these mapping populations were concordant and showed that Burly1 is located between 151.9 and 152.7 Mb (rs33197365 to rs3700604) on mouse chromosome 2. The congenic region harboring Burly1 contains 26 protein-coding genes, 11 noncoding RNA elements (e.g., lncRNA), and 4 pseudogenes, with 1949 predicted functional variants. Of the protein-coding genes, 7 have missense variants, including genes that may contribute to lean body weight, such as Angpt41, Slc52c3, and Rem1. Lean body mass was increased by the B6-derived variant relative to the 129-derived allele. Burly1 influenced lean body weight at all ages but not food intake or locomotor activity. However, congenic mice with the B6 allele produced more heat per kilogram of lean body weight than did controls, pointing to a genotype effect on lean mass metabolism. These results show the value of integrating information from several mapping populations to refine the map location of body composition QTLs and to identify a short list of candidate genes.

Expression profile analysis of long non-coding RNAs involved in the metformin-inhibited gluconeogenesis of primary mouse hepatocytes

Long non-coding RNAs (lncRNAs) have been demonstrated to regulate metabolic tissue development and function, including adipogenesis, hepatic lipid metabolism, islet function and energy balance. However, the role of lncRNAs in gluconeogenesis remains completely unknown. Metformin reduces glucose output mainly via the inhibition of gluconeogenesis. In the present study, the lncRNA expression profile of primary mouse hepatocytes exposed to cyclic adenosine monophosphate (cAMP), a gluconeogenic stimulus, with or without metformin was analyzed by microarray. Among the 22,016 lncRNAs that were identified, 456 were upregulated and 409 were downregulated by cAMP (fold-change ≥2.0). Furthermore, the cAMP-induced upregulation of 189 lncRNAs and downregulation of 167 lncRNAs was attenuated by metformin. The expression levels of eight lncRNAs were validated by reverse transcription-quantitative polymerase chain reaction, and the results were consistent with those of the microarray analysis. Among them, two lncRNAs NR_027710 and ENSMUST00000138573, were identified to have an association with two protein coding genes, namely peroxisome proliferator-activated receptor-γ coactivator-1α, a critical transcriptional coactivator in gluconeogenesis, and G protein-coupled receptor 155, respectively. The two protein coding genes exhibited similar expression patterns to their associated lncRNAs. The findings of the present study suggest that lncRNAs are potentially involved in the regulation of gluconeogenesis.

A Systems Approach Implicates a Brain Mitochondrial Oxidative Homeostasis Co-expression Network in Genetic Vulnerability to Alcohol Withdrawal

Genetic factors significantly affect vulnerability to alcohol dependence (alcoholism). We previously identified quantitative trait loci on distal mouse chromosome 1 with large effects on predisposition to alcohol physiological dependence and associated withdrawal following both chronic and acute alcohol exposure in mice (Alcdp1 and Alcw1, respectively). We fine-mapped these loci to a 1.1–1.7 Mb interval syntenic with human 1q23.2-23.3. Alcw1/Alcdp1 interval genes show remarkable genetic variation among mice derived from the C57BL/6J and DBA/2J strains, the two most widely studied genetic animal models for alcohol-related traits. Here, we report the creation of a novel recombinant Alcw1/Alcdp1 congenic model (R2) in which the Alcw1/Alcdp1 interval from a donor C57BL/6J strain is introgressed onto a uniform, inbred DBA/2J genetic background. As expected, R2 mice demonstrate significantly less severe alcohol withdrawal compared to wild-type littermates. Additionally, comparing R2 and background strain animals, as well as reciprocal congenic (R8) and appropriate background strain animals, we assessed Alcw1/Alcdp1 dependent brain gene expression using microarray and quantitative PCR analyses. To our knowledge this includes the first Weighted Gene Co-expression Network Analysis using reciprocal congenic models. Importantly, this allows detection of co-expression patterns limited to one or common to both genetic backgrounds with high or low predisposition to alcohol withdrawal severity. The gene expression patterns (modules) in common contain genes related to oxidative phosphorylation, building upon human and animal model studies that implicate involvement of oxidative phosphorylation in alcohol use disorders (AUDs). Finally, we demonstrate that administration of N-acetylcysteine, an FDA-approved antioxidant, significantly reduces symptoms of alcohol withdrawal (convulsions) in mice, thus validating a phenotypic role for this network. Taken together, these studies support the importance of mitochondrial oxidative homeostasis in alcohol withdrawal and identify this network as a valuable therapeutic target in human AUDs.

Identification of DEP domain-containing proteins by a machine learning method and experimental analysis of their expression in human HCC tissues

The Dishevelled/EGL-10/Pleckstrin (DEP) domain-containing (DEPDC) proteins have seven members. However, whether this superfamily can be distinguished from other proteins based only on the amino acid sequences, remains unknown. Here, we describe a computational method to segregate DEPDCs and non-DEPDCs. First, we examined the Pfam numbers of the known DEPDCs and used the longest sequences for each Pfam to construct a phylogenetic tree. Subsequently, we extracted 188-dimensional (188D) and 20D features of DEPDCs and non-DEPDCs and classified them with random forest classifier. We also mined the motifs of human DEPDCs to find the related domains. Finally, we designed experimental verification methods of human DEPDC expression at the mRNA level in hepatocellular carcinoma (HCC) and adjacent normal tissues. The phylogenetic analysis showed that the DEPDCs superfamily can be divided into three clusters. Moreover, the 188D and 20D features can both be used to effectively distinguish the two protein types. Motif analysis revealed that the DEP and RhoGAP domain was common in human DEPDCs, human HCC and the adjacent tissues that widely expressed DEPDCs. However, their regulation was not identical. In conclusion, we successfully constructed a binary classifier for DEPDCs and experimentally verified their expression in human HCC tissues.

Selection signature analysis reveals genes associated with tail type in Chinese indigenous sheep

Fat-tailed sheep have commercial value because consumers prefer high-protein and low-fat food and producers care about feed conversion rate. However, fat-tailed sheep still have some scientific significance, as the fat tail is commonly regarded as a characteristic of environmental adaptability. Finding the candidate genes associated with fat tail formation is essential for breeding and conservation. To identify these candidate genes, we applied FST and hapFLK approaches in fat- and thin-tailed sheep with available 50K SNP genotype data. These two methods found 6.24 Mb of overlapped regions and 43 genes that may associated with fat tail development. Gene annotation showed that HOXA11, BMP2, PPP1CC, SP3, SP9, WDR92, PROKR1 and ETAA1 may play important roles in fat tail formation. These findings provide insight into tail fat development and a guide for molecular breeding and conservation.

A GWA study reveals genetic loci for body conformation traits in Chinese Laiwu pigs and its implications for human BMI

Pigs share numerous physiological and phenotypic similarities with human and thus have been considered as a good model in nonrodent mammals for the study of genetic basis of human obesity. Researches on candidate genes for obesity traits have successfully identified some common genes between humans and pigs. However, few studies have assessed how many similarities exist between the genetic architecture of obesity in pigs and humans by large-scale comparative genomics. Here, we performed a genome-wide association study (GWAS) using the porcine 60 K SNP Beadchip for BMI and other four conformation traits at three different ages in a Chinese Laiwu pig population, which shows a large variability in fat deposition. In total, 35 SNPs were found to be significant at Bonferroni-corrected 5 % chromosome-wise level (P = 2.13 × 10^-5) and 88 SNPs had suggestive (P < 10^-4) association with the conformation traits. Some SNPs showed age-dependent association. Intriguingly, out of 32 regions associated with BMI in pigs, 18 were homologous with the loci for BMI in humans. Furthermore, five closest genes to GWAS peaks including HIF1AN, SMYD3, COX10, SLMAP, and GBE1 have been already associated with BMI in humans, which makes them very promising candidates for these QTLs. The result of GO analysis provided strong support to the fact that mitochondria and synapse play important roles in obesity susceptibility, which is consistent with previous findings on human obesity, and it also implicated new gene sets related to chromatin modification and Ig-like C2-type 5 domain. Therefore, these results not only provide new insights into the genetic architecture of BMI in pigs but also highlight that humans and pigs share the significant overlap of obesity-related genes.

An Overview of Murine High Fat Diet as a Model for Type 2 Diabetes Mellitus

Type 2 diabetes mellitus (T2DM) is a worldwide epidemic, which by all predictions will only increase. To help in combating the devastating array of phenotypes associated with T2DM a highly reproducible and human disease-similar mouse model is required for researchers. The current options are genetic manipulations to cause T2DM symptoms or diet induced obesity and T2DM symptoms. These methods to model human T2DM have their benefits and their detractions. As far as modeling the majority of T2DM cases, HFD establishes the proper etiological, pathological, and treatment options. A limitation of HFD is that it requires months of feeding to achieve the full spectrum of T2DM symptoms and no standard protocol has been established. This paper will attempt to rectify the last limitation and argue for a standard group of HFD protocols and standard analysis procedures.