Phosphoglucomutase Genetic Polymorphism and Body Mass

Characterization of gut microbial and metabolite alterations in faeces of Goto Kakizaki rats using metagenomic and untargeted metabolomic approach

BACKGROUND

In recent years, the incidence of type 2 diabetes (T2DM) has shown a rapid growth trend. Goto Kakizaki (GK) rats are a valuable model for the study of T2DM and share common glucose metabolism features with human T2DM patients. A series of studies have indicated that T2DM is associated with the gut microbiota composition and gut metabolites. We aimed to systematically characterize the faecal gut microbes and metabolites of GK rats and analyse the relationship between glucose and insulin resistance.

AIM

To evaluate the gut microbial and metabolite alterations in GK rat faeces based on metagenomics and untargeted metabolomics.

METHODS

Ten GK rats (model group) and Wistar rats (control group) were observed for 10 wk, and various glucose-related indexes, mainly including weight, fasting blood glucose (FBG) and insulin levels, homeostasis model assessment of insulin resistance (HOMA-IR) and homeostasis model assessment of β cell (HOMA-β) were assessed. The faecal gut microbiota was sequenced by metagenomics, and faecal metabolites were analysed by untargeted metabolomics. Multiple metabolic pathways were evaluated based on the differential metabolites identified, and the correlations between blood glucose and the gut microbiota and metabolites were analysed.

RESULTS

The model group displayed significant differences in weight, FBG and insulin levels, HOMA-IR and HOMA-β indexes (P < 0.05, P < 0.01) and a shift in the gut microbiota structure compared with the control group. The results demonstrated significantly decreased abundances of Prevotella sp. CAG:604 and Lactobacillus murinus (P < 0.05) and a significantly increased abundance of Allobaculum stercoricanis (P < 0.01) in the model group. A correlation analysis indicated that FBG and HOMA-IR were positively correlated with Allobaculum stercoricanis and negatively correlated with Lactobacillus murinus. An orthogonal partial least squares discriminant analysis suggested that the faecal metabolic profiles differed between the model and control groups. Fourteen potential metabolic biomarkers, including glycochenodeoxycholic acid, uric acid, 13(S)-hydroxyoctadecadienoic acid (HODE), N-acetylaspartate, β-sitostenone, sphinganine, 4-pyridoxic acid, and linoleic acid, were identified. Moreover, FBG and HOMA-IR were found to be positively correlated with glutathione, 13(S)-HODE, uric acid, 4-pyridoxic acid and allantoic acid and ne-gatively correlated with 3-α, 7-α, chenodeoxycholic acid glycine conjugate and 26-trihydroxy-5-β-cholestane (P < 0.05, P < 0.01). Allobaculum stercoricanis was positively correlated with linoleic acid and sphinganine (P < 0.01), and 2-methyl-3-hydroxy-5-formylpyridine-4-carboxylate was negatively associated with Prevotella sp. CAG:604 (P < 0.01). The metabolic pathways showing the largest differences were arginine biosynthesis; primary bile acid biosynthesis; purine metabolism; linoleic acid metabolism; alanine, aspartate and glutamate metabolism; and nitrogen metabolism.

CONCLUSION

Metagenomics and untargeted metabolomics indicated that disordered compositions of gut microbes and metabolites may be common defects in GK rats.

Birth weight and risk of type 2 diabetes: A dose-response meta-analysis of cohort studies

The association between birth weight and type 2 diabetes mellitus has been debated for several decades. The objective of this systematic review and meta-analysis was to quantitatively clarify the association between birth weight and risk of type 2 diabetes mellitus based on cohort studies. We searched PubMed, Web of Science, and Embase databases for cohort study articles on the association between birth weight and risk of type 2 diabetes mellitus published up to 1 March 2018. Random effects of generalized least square regression models were used to estimate relative risk (RR). Restricted cubic splines were conducted to model the dose-response relationship. We included 21 studies (19 articles) involving 1 041 879 individuals and 35 699 cases of type 2 diabetes mellitus, with follow-up ranged from 6 to 47 years. We identified significant decreasing trend for the highest versus lowest category of birth weight for the association with type 2 diabetes mellitus risk: The risk was reduced by 35% (RR, 0.65; 95% confidence interval [CI], 0.53-0.81) and by 12% (RR 0.88; 95% CI, 0.85-0.91) per 500-g increment in birth weight. Our results showed a dose-response relationship between birth weight and diabetes risk, which was nonlinear (P_nonlinearity  < 0.001) and L-shaped. With increasing birth weight (<5000 g), the risk of type 2 diabetes mellitus decreased substantially. The association between birth weight and type 2 diabetes mellitus was curvilinear and L-shaped.

Transcriptomic analysis between Normal and high-intake feeding geese provides insight into adipose deposition and susceptibility to fatty liver in migratory birds

BACKGROUND

Dysregulation of adipogenesis causes metabolic diseases, like obesity and fatty liver. Migratory birds such as geese have a high tolerance of massive energy intake and exhibit little pathological development. Domesticated goose breeds, derivatives of the wild greyleg goose (Anser anser) or swan goose (Anser cygnoides), have high tolerance of energy intake resembling their ancestor species. Thus, goose is potentially a model species to study mechanisms associated with adipogenesis.

RESULTS

Phenotypically, goose liver exhibited higher fat accumulation than adipose tissues during fattening (liver increased by 3.35 fold than 1.65 fold in adipose), showing a priority of fat accumulation in liver. We found the number of differentially expressed genes in liver (13.97%) was nearly twice the number of that in adipose (6.60%). These differentially expressed genes in liver function in several important lipid metabolism pathways, immune response, regulation of cancer, while in adipose, terms closely related to protein binding, gluconeogenesis were enriched. Typically, genes like MDH2 and SCD, which have key roles in glycolysis and fatty acids metabolism, had higher fold change in liver than in adipose tissues. Three hundred two differentially expressed long noncoding RNAs involved in regulation of metabolism in liver were also identified. For example, lncRNA XLOC_292762, which was 5.7 kb downstream of FERMT2, a gene involved phosphatidylinositol-3,4,5-trisphosphate binding, was significantly down-regulated after the high-intake feeding period. Further investigation of documented obesity-related orthologous genes in goose suggested that understanding the evolutionary split from mammals in adipogenesis will make goose fatty liver a better resource for future research.

CONCLUSIONS

Our research reveals that goose uses liver as the major tissue to regulate a distinct lipid synthesis and degradation flux and the dynamic expression network analyses showed numerous layers of positive responses to both massive energy intake and possible pathological development. Our results offer insights into goose adipogenesis and provide a new perspective for research in human metabolic dysregulation.

Proteomics-Based Identification of the Molecular Signatures of Liver Tissues from Aged Rats following Eight Weeks of Medium-Intensity Exercise

Physical activity has emerged as a powerful intervention that promotes healthy aging by maintaining the functional capacity of critical organ systems. Here, by combining functional and proteomics analyses, we examined how hepatic phenotypes might respond to exercise treatment in aged rats. 16 male aged (20 months old) SD rats were divided into exercise and parallel control groups at random; the exercise group had 8 weeks of treadmill training with medium intensity. Whole protein samples of the liver were extracted from both groups and separated by two-dimensional gel electrophoresis. Alternatively objective protein spots with >2-fold difference in expression were selected for enzymological extraction and MS/MS identification. Results show increased activity of the manganese superoxide dismutase and elevated glutathione levels in the livers of exercise-treated animals, but malondialdehyde contents obviously decreased in the liver of the exercise group. Proteomics-based identification of differentially expressed proteins provided an integrated view of the metabolic adaptations occurring in the liver proteome during exercise, which significantly altered the expression of several proteins involved in key liver metabolic pathways including mitochondrial sulfur, glycolysis, methionine, and protein metabolism. These findings indicate that exercise may be beneficial to aged rats through modulation of hepatic protein expression profiles.

The evolutionary fate of alternatively spliced homologous exons after gene duplication

Alternative splicing and gene duplication are the two main processes responsible for expanding protein functional diversity. Although gene duplication can generate new genes and alternative splicing can introduce variation through alternative gene products, the interplay between the two processes is complex and poorly understood. Here, we have carried out a study of the evolution of alternatively spliced exons after gene duplication to better understand the interaction between the two processes. We created a manually curated set of 97 human genes with mutually exclusively spliced homologous exons and analyzed the evolution of these exons across five distantly related vertebrates (lamprey, spotted gar, zebrafish, fugu, and coelacanth). Most of these exons had an ancient origin (more than 400 Ma). We found examples supporting two extreme evolutionary models for the behaviour of homologous axons after gene duplication. We observed 11 events in which gene duplication was accompanied by splice isoform separation, that is, each paralog specifically conserved just one distinct ancestral homologous exon. At other extreme, we identified genes in which the homologous exons were always conserved within paralogs, suggesting that the alternative splicing event cannot easily be separated from the function in these genes. That many homologous exons fall in between these two extremes highlights the diversity of biological systems and suggests that the subtle balance between alternative splicing and gene duplication is adjusted to the specific cellular context of each gene.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

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

Compromised catalysis and potential folding defects in in vitro studies of missense mutants associated with hereditary phosphoglucomutase 1 deficiency

Recent studies have identified phosphoglucomutase 1 (PGM1) deficiency as an inherited metabolic disorder in humans. Affected patients show multiple disease phenotypes, including dilated cardiomyopathy, exercise intolerance, and hepatopathy, reflecting the central role of the enzyme in glucose metabolism. We present here the first in vitro biochemical characterization of 13 missense mutations involved in PGM1 deficiency. The biochemical phenotypes of the PGM1 mutants cluster into two groups: those with compromised catalysis and those with possible folding defects. Relative to the recombinant wild-type enzyme, certain missense mutants show greatly decreased expression of soluble protein and/or increased aggregation. In contrast, other missense variants are well behaved in solution, but show dramatic reductions in enzyme activity, with kcat/Km often <1.5% of wild-type. Modest changes in protein conformation and flexibility are also apparent in some of the catalytically impaired variants. In the case of the G291R mutant, severely compromised activity is linked to the inability of a key active site serine to be phosphorylated, a prerequisite for catalysis. Our results complement previous in vivo studies, which suggest that both protein misfolding and catalytic impairment may play a role in PGM1 deficiency.

Genome-wide linkage and copy number variation analysis reveals 710 kb duplication on chromosome 1p31.3 responsible for autosomal dominant omphalocele

BACKGROUND

Omphalocele is a congenital birth defect characterised by the presence of internal organs located outside of the ventral abdominal wall. The purpose of this study was to identify the underlying genetic mechanisms of a large autosomal dominant Caucasian family with omphalocele.

METHODS AND FINDINGS

A genetic linkage study was conducted in a large family with an autosomal dominant transmission of an omphalocele using a genome-wide single nucleotide polymorphism (SNP) array. The analysis revealed significant evidence of linkage (non-parametric NPL = 6.93, p=0.0001; parametric logarithm of odds (LOD) = 2.70 under a fully penetrant dominant model) at chromosome band 1p31.3. Haplotype analysis narrowed the locus to a 2.74 Mb region between markers rs2886770 (63014807 bp) and rs1343981 (65757349 bp). Molecular characterisation of this interval using array comparative genomic hybridisation followed by quantitative microsphere hybridisation analysis revealed a 710 kb duplication located at 63.5-64.2 Mb. All affected individuals who had an omphalocele and shared the haplotype were positive for this duplicated region, while the duplication was absent from all normal individuals of this family. Multipoint linkage analysis using the duplication as a marker yielded a maximum LOD score of 3.2 at 1p31.3 under a dominant model. The 710 kb duplication at 1p31.3 band contains seven known genes including FOXD3, ALG6, ITGB3BP, KIAA1799, DLEU2L, PGM1, and the proximal portion of ROR1. Importantly, this duplication is absent from the database of genomic variants.

CONCLUSIONS

The present study suggests that development of an omphalocele in this family is controlled by overexpression of one or more genes in the duplicated region. To the authors' knowledge, this is the first reported association of an inherited omphalocele condition with a chromosomal rearrangement.

Multivariate linkage scan for metabolic syndrome traits in families with type 2 diabetes

The purpose of this study was to evaluate evidence for linkage to interrelated quantitative features of the metabolic syndrome (MetS). Data on eight quantitative MetS traits (body weight, waist circumference, systolic and diastolic blood pressure, high-density lipoprotein (HDL) cholesterol, triglycerides (TGs), and fasting glucose and insulin measurements) and a 10 cM genome scan were available for 78 white families (n = 532 subjects). These data were used to conduct multipoint, multivariate linkage analyses, including tests for coincident linkage and complete pleiotropy. The strongest evidence for linkage from the bivariate analyses was observed on chromosome 1 (1p22.2) (HDL-TG; univariate lod score equivalent (lod(eq) = 3.99)) with stronger results from the trivariate analysis at the same location (HDL-TG-Insulin; lod(eq) = 4.32). Seven additional susceptibility regions (lod(eq) scores >1.9) were observed (1p36, 1q23, 2q21.2, 8q23.3, 14q23.2, 14q32.11, and 20p11.21). The results from this study indicate that several correlated traits of the MetS are influenced by the same gene(s) that account for some of the clustering of the MetS features.

Lack of association between PCK1 polymorphisms and obesity, physical activity, and fitness in European Youth Heart Study (EYHS)

Phosphoenolpyruvate carboxykinase-1 (PCK1) is the rate-limiting enzyme in the hepatic gluconeogenic pathway. Studies have shown that overexpression of Pck1 in mice results in obesity-related traits and higher levels of physical activity (PA). Therefore, our aims were to investigate whether common genetic variation in the PCK1 gene influences obesity-related traits, PA, and fitness, and to examine whether PA and fitness attenuate the influence of the PCK1 polymorphisms on obesity in children. Analyses were undertaken on data from Danish and Estonian children (958 boys and 1,104 girls) from the European Youth Heart Study (EYHS), a school-based, cross-sectional study of children (mean ± s.d. age: 9.6 ± 0.4 years) and adolescents (15.5 ± 0.5 years). We genotyped eight polymorphisms that captured the common genetic variations in the PCK1 gene. The association between the PCK1 polymorphisms and BMI, waist circumference (WC), sum of four skinfolds, PA, and fitness was tested using an additive model adjusted for age, age-group, gender, maturity, and country. Interactions were tested by including interaction terms in the model. None of the polymorphisms were significantly associated with BMI, WC, sum of four skinfolds, PA, and fitness, and also with the risk of being overweight or obese (P > 0.05). The interactions between the polymorphisms and age-group, gender, PA, and fitness were not statistically significant. This is the first study to comprehensively examine the association of PCK1 polymorphisms with obesity, PA, and fitness. Despite strong evidence from animal studies, our study in the EYHS cohort failed to identify an association of PCK1 polymorphisms with obesity, PA, and fitness.

Detecting purely epistatic multi-locus interactions by an omnibus permutation test on ensembles of two-locus analyses

BACKGROUND

Purely epistatic multi-locus interactions cannot generally be detected via single-locus analysis in case-control studies of complex diseases. Recently, many two-locus and multi-locus analysis techniques have been shown to be promising for the epistasis detection. However, exhaustive multi-locus analysis requires prohibitively large computational efforts when problems involve large-scale or genome-wide data. Furthermore, there is no explicit proof that a combination of multiple two-locus analyses can lead to the correct identification of multi-locus interactions.

RESULTS

The proposed 2LOmb algorithm performs an omnibus permutation test on ensembles of two-locus analyses. The algorithm consists of four main steps: two-locus analysis, a permutation test, global p-value determination and a progressive search for the best ensemble. 2LOmb is benchmarked against an exhaustive two-locus analysis technique, a set association approach, a correlation-based feature selection (CFS) technique and a tuned ReliefF (TuRF) technique. The simulation results indicate that 2LOmb produces a low false-positive error. Moreover, 2LOmb has the best performance in terms of an ability to identify all causative single nucleotide polymorphisms (SNPs) and a low number of output SNPs in purely epistatic two-, three- and four-locus interaction problems. The interaction models constructed from the 2LOmb outputs via a multifactor dimensionality reduction (MDR) method are also included for the confirmation of epistasis detection. 2LOmb is subsequently applied to a type 2 diabetes mellitus (T2D) data set, which is obtained as a part of the UK genome-wide genetic epidemiology study by the Wellcome Trust Case Control Consortium (WTCCC). After primarily screening for SNPs that locate within or near 372 candidate genes and exhibit no marginal single-locus effects, the T2D data set is reduced to 7,065 SNPs from 370 genes. The 2LOmb search in the reduced T2D data reveals that four intronic SNPs in PGM1 (phosphoglucomutase 1), two intronic SNPs in LMX1A (LIM homeobox transcription factor 1, alpha), two intronic SNPs in PARK2 (Parkinson disease (autosomal recessive, juvenile) 2, parkin) and three intronic SNPs in GYS2 (glycogen synthase 2 (liver)) are associated with the disease. The 2LOmb result suggests that there is no interaction between each pair of the identified genes that can be described by purely epistatic two-locus interaction models. Moreover, there are no interactions between these four genes that can be described by purely epistatic multi-locus interaction models with marginal two-locus effects. The findings provide an alternative explanation for the aetiology of T2D in a UK population.

CONCLUSION

An omnibus permutation test on ensembles of two-locus analyses can detect purely epistatic multi-locus interactions with marginal two-locus effects. The study also reveals that SNPs from large-scale or genome-wide case-control data which are discarded after single-locus analysis detects no association can still be useful for genetic epidemiology studies.