Two Novel Candidate Genes for Insulin Secretion Identified by Comparative Genomics of Multiple Backcross Mouse Populations

Identifying target ion channel-related genes to construct a diagnosis model for insulinoma

Background: Insulinoma is the most common functional pancreatic neuroendocrine tumor (PNET) with abnormal insulin hypersecretion. The etiopathogenesis of insulinoma remains indefinable. Based on multiple bioinformatics methods and machine learning algorithms, this study proposed exploring the molecular mechanism from ion channel-related genes to establish a genetic diagnosis model for insulinoma. Methods: The mRNA expression profile dataset of GSE73338 was applied to the analysis, which contains 17 insulinoma samples, 63 nonfunctional PNET (NFPNET) samples, and four normal islet samples. Differently expressed ion channel-related genes (DEICRGs) enrichment analyses were performed. We utilized the protein-protein interaction (PPI) analysis and machine learning of LASSO and support vector machine-recursive feature elimination (SVM-RFE) to identify the target genes. Based on these target genes, a nomogram diagnostic model was constructed and verified by a receiver operating characteristic (ROC) curve. Moreover, immune infiltration analysis, single-gene gene set enrichment analysis (GSEA), and gene set variation analysis (GSVA) were executed. Finally, a drug-gene interaction network was constructed. Results: We identified 29 DEICRGs, and enrichment analyses indicated they were primarily enriched in ion transport, cellular ion homeostasis, pancreatic secretion, and lysosome. Moreover, the PPI network and machine learning recognized three target genes (MCOLN1, ATP6V0E1, and ATP4A). Based on these target genes, we constructed an efficiently predictable diagnosis model for identifying insulinomas with a nomogram and validated it with the ROC curve (AUC = 0.801, 95% CI 0.674-0.898). Then, single-gene GSEA analysis revealed that these target genes had a significantly positive correlation with insulin secretion and lysosome. In contrast, the TGF-beta signaling pathway was negatively associated with them. Furthermore, statistically significant discrepancies in immune infiltration were revealed. Conclusion: We identified three ion channel-related genes and constructed an efficiently predictable diagnosis model to offer a novel approach for diagnosing insulinoma.

E96V Mutation in the Kdelr3 Gene Is Associated with Type 2 Diabetes Susceptibility in Obese NZO Mice

Type 2 diabetes (T2D) represents a multifactorial metabolic disease with a strong genetic predisposition. Despite elaborate efforts in identifying the genetic variants determining individual susceptibility towards T2D, the majority of genetic factors driving disease development remain poorly understood. With the aim to identify novel T2D risk genes we previously generated an N2 outcross population using the two inbred mouse strains New Zealand obese (NZO) and C3HeB/FeJ (C3H). A linkage study performed in this population led to the identification of the novel T2D-associated quantitative trait locus (QTL) Nbg15 (NZO blood glucose on chromosome 15, Logarithm of odds (LOD) 6.6). In this study we used a combined approach of positional cloning, gene expression analyses and in silico predictions of DNA polymorphism on gene/protein function to dissect the genetic variants linking Nbg15 to the development of T2D. Moreover, we have generated congenic strains that associated the distal sublocus of Nbg15 to mechanisms altering pancreatic beta cell function. In this sublocus, Cbx6, Fam135b and Kdelr3 were nominated as potential causative genes associated with the Nbg15 driven effects. Moreover, a putative mutation in the Kdelr3 gene from NZO was identified, negatively influencing adaptive responses associated with pancreatic beta cell death and induction of endoplasmic reticulum stress. Importantly, knockdown of Kdelr3 in cultured Min6 beta cells altered insulin granules maturation and pro-insulin levels, pointing towards a crucial role of this gene in islets function and T2D susceptibility.

Accelerated aging in cyclophilin B deficient mice downstream of p21‐Cip1 /Waf1

Loss of bone mass and strength is a common problem of advanced age in humans. Defective bone is also a primary finding in osteogenesis imperfecta (OI), a genetic condition most commonly caused by autosomal dominant mutations in the type I collagen genes. Although altered collagen has been proposed to correlate with cellular processes that underlie aging, the causal relationships between them in vivo have not yet been completely explored. Whether aging plays a promoting role in OI development or whether OI contributes to aging, also remains unknown. The PpiB gene encodes cyclophilin B (CypB), a prolyl isomerase residing in the endoplasmic reticulum required for normal assembly of collagen. Germline deletion or mutations of CypB in mice or humans cause autosomal recessive OI (type IX). Here, we show that mice lacking CypB develop early onset of aging‐associated phenotypes, including kyphosis, fat reduction and weight loss, as well as abnormal teeth, skin, and muscle. Elevated senescence‐associated beta‐galactosidase (SA‐β‐Gal) activity was observed in fat tissues and in bone marrow–derived multipotent stromal cells. Protein levels of the cyclin‐dependent kinase (cdk)‐inhibitor p21‐Cip1/Waf1, a well known senescence marker, were significantly elevated in CypB‐deficient primary cells and mouse tissues. Importantly, loss of p21 in CypB knockout mice attenuated SA‐β‐Gal activity and delayed the development of kyphosis. In addition, less adipose tissue depot and higher SA‐β‐Gal activity were observed in a second OI model, Cola2^oim mutant mice. A potential upregulation of p21 was also revealed in a limited number of these mice. These findings suggest that some of the features in OI patients may be mediated in part through activation of the p21‐dependent pathway, one of which is closely associated with senescence and aging. This study provides new mechanistic insight into relationships between OI and aging and raises the possibility of using senolytics drugs to treat OI in the future. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Comparative genomic analyses of multiple backcross mouse populations suggest SGCG as a novel potential obesity-modifier gene

To nominate novel disease genes for obesity and type 2 diabetes (T2D), we recently generated two mouse backcross populations of the T2D-susceptible New Zealand Obese (NZO/HI) mouse strain and two genetically different, lean and T2D-resistant strains, 129P2/OlaHsd and C3HeB/FeJ. Comparative linkage analysis of our two female backcross populations identified seven novel body fat-associated quantitative trait loci (QTL). Only the locus Nbw14 (NZO body weight on chromosome 14) showed linkage to obesity-related traits in both backcross populations, indicating that the causal gene variant is likely specific for the NZO strain as NZO allele carriers in both crosses displayed elevated body weight and fat mass. To identify candidate genes for Nbw14, we used a combined approach of gene expression and haplotype analysis to filter for NZO-specific gene variants in gonadal white adipose tissue (gWAT), defined as the main QTL-target tissue. Only two genes, Arl11 and Sgcg, fulfilled our candidate criteria. In addition, expression QTL analysis revealed cis-signals for both genes within the Nbw14 locus. Moreover, retroviral overexpression of Sgcg in 3 T3-L1 adipocytes resulted in increased insulin-stimulated glucose uptake. In humans, mRNA levels of SGCG correlated with BMI and body fat mass exclusively in diabetic subjects, suggesting that SGCG may present a novel marker for metabolically unhealthy obesity. In conclusion, our comparative-cross analysis could substantially improve the mapping resolution of the obesity locus Nbw14. Future studies will shine light on the mechanism by which Sgcg may protect from the development of obesity.

Nudix hydrolase NUDT19 regulates mitochondrial function and ATP production in murine hepatocytes

Changes in intracellular CoA levels are known to contribute to the development of non-alcoholic fatty liver disease (NAFLD) in type 2 diabetes (T2D) in human and rodents. However, the underlying genetic basis is still poorly understood. Due to their diverse susceptibility towards metabolic diseases, mouse inbred strains have been proven to serve as powerful tools for the identification of novel genetic factors that underlie the pathophysiology of NAFLD and diabetes. Transcriptome analysis of mouse liver samples revealed the nucleoside diphosphate linked moiety X-type motif Nudt19 as novel candidate gene responsible for NAFLD and T2D development. Knockdown (KD) of Nudt19 increased mitochondrial and glycolytic ATP production rates in Hepa 1-6 cells by 41% and 10%, respectively. The enforced utilization of glutamine or fatty acids as energy substrate reduced uncoupled respiration by 41% and 47%, respectively, in non-target (NT) siRNA transfected cells. This reduction was prevented upon Nudt19 KD. Furthermore, incubation with palmitate or oleate respectively increased mitochondrial ATP production by 31% and 20%, and uncoupled respiration by 23% and 30% in Nudt19 KD cells, but not in NT cells. The enhanced fatty acid oxidation in Nudt19 KD cells was accompanied by a 1.3-fold increased abundance of Pdk4. This study is the first to describe Nudt19 as regulator of hepatic lipid metabolism and potential mediator of NAFLD and T2D development.

Identification of Novel Genes Involved in Hyperglycemia in Mice

Current attempts to prevent and manage type 2 diabetes have been moderately effective, and a better understanding of the molecular roots of this complex disease is important to develop more successful and precise treatment options. Recently, we initiated the collective diabetes cross, where four mouse inbred strains differing in their diabetes susceptibility were crossed with the obese and diabetes-prone NZO strain and identified the quantitative trait loci (QTL) Nidd13/NZO, a genomic region on chromosome 13 that correlates with hyperglycemia in NZO allele carriers compared to B6 controls. Subsequent analysis of the critical region, harboring 644 genes, included expression studies in pancreatic islets of congenic Nidd13/NZO mice, integration of single-cell data from parental NZO and B6 islets as well as haplotype analysis. Finally, of the five genes (Acot12, S100z, Ankrd55, Rnf180, and Iqgap2) within the polymorphic haplotype block that are differently expressed in islets of B6 compared to NZO mice, we identified the calcium-binding protein S100z gene to affect islet cell proliferation as well as apoptosis when overexpressed in MIN6 cells. In summary, we define S100z as the most striking gene to be causal for the diabetes QTL Nidd13/NZO by affecting β-cell proliferation and apoptosis. Thus, S100z is an entirely novel diabetes gene regulating islet cell function.

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.

Alternative exon splicing and differential expression in pancreatic islets reveals candidate genes and pathways implicated in early diabetes development

Type 2 diabetes (T2D) has a strong genetic component. Most of the gene variants driving the pathogenesis of T2D seem to target pancreatic β-cell function. To identify novel gene variants acting at early stage of the disease, we analyzed whole transcriptome data to identify differential expression (DE) and alternative exon splicing (AS) transcripts in pancreatic islets collected from two metabolically diverse mouse strains at 6 weeks of age after three weeks of high-fat-diet intervention. Our analysis revealed 1218 DE and 436 AS genes in islets from NZO/Hl vs C3HeB/FeJ. Whereas some of the revealed genes present well-established markers for β-cell failure, such as Cd36 or Aldh1a3, we identified numerous DE/AS genes that have not been described in context with β-cell function before. The gene Lgals2, previously associated with human T2D development, was DE as well as AS and localizes in a quantitative trait locus (QTL) for blood glucose on Chr.15 that we reported recently in our N₂(NZOxC3H) population. In addition, pathway enrichment analysis of DE and AS genes showed an overlap of only half of the revealed pathways, indicating that DE and AS in large parts influence different pathways in T2D development. PPARG and adipogenesis pathways, two well-established metabolic pathways, were overrepresented for both DE and AS genes, probably as an adaptive mechanism to cope for increased cellular stress. Our results provide guidance for the identification of novel T2D candidate genes and demonstrate the presence of numerous AS transcripts possibly involved in islet function and maintenance of glucose homeostasis.

Immunity-related GTPase induces lipophagy to prevent excess hepatic lipid accumulation

BACKGROUND & AIMS

Currently, only a few genetic variants explain the heritability of fatty liver disease. Quantitative trait loci (QTL) analysis of mouse strains has identified the susceptibility locus Ltg/NZO (liver triglycerides from New Zealand obese [NZO] alleles) on chromosome 18 as associating with increased hepatic triglycerides. Herein, we aimed to identify genomic variants responsible for this association.

METHODS

Recombinant congenic mice carrying 5.3 Mbp of Ltg/NZO were fed a high-fat diet and characterized for liver fat. Bioinformatic analysis, mRNA profiles and electrophoretic mobility shift assays were performed to identify genes responsible for the Ltg/NZO phenotype. Candidate genes were manipulated in vivo by injecting specific microRNAs into C57BL/6 mice. Pulldown coupled with mass spectrometry-based proteomics and immunoprecipitation were performed to identify interaction partners of IFGGA2.

RESULTS

Through positional cloning, we identified 2 immunity-related GTPases (Ifgga2, Ifgga4) that prevent hepatic lipid storage. Expression of both murine genes and the human orthologue IRGM was significantly lower in fatty livers. Accordingly, liver-specific suppression of either Ifgga2 or Ifgga4 led to a 3-4-fold greater increase in hepatic fat content. In the liver of low-fat diet-fed mice, IFGGA2 localized to endosomes/lysosomes, while on a high-fat diet it associated with lipid droplets. Pulldown experiments and proteomics identified the lipase ATGL as a binding partner of IFGGA2 which was confirmed by co-immunoprecipitation. Both proteins partially co-localized with the autophagic marker LC3B. Ifgga2 suppression in hepatocytes reduced the amount of LC3B-II, whereas overexpression of Ifgga2 increased the association of LC3B with lipid droplets and decreased triglyceride storage.

CONCLUSION

IFGGA2 interacts with ATGL and protects against hepatic steatosis, most likely by enhancing the binding of LC3B to lipid droplets.

LAY SUMMARY

The genetic basis of non-alcoholic fatty liver disease remains incompletely defined. Herein, we identified members of the immunity-related GTPase family in mice and humans that act as regulators of hepatic fat accumulation, with links to autophagy. Overexpression of the gene Ifgga2 was shown to reduce hepatic lipid storage and could be a therapeutic target for the treatment of fatty liver disease.

Identification of Novel Potential Type 2 Diabetes Genes Mediating β-Cell Loss and Hyperglycemia Using Positional Cloning

Type 2 diabetes (T2D) is a complex metabolic disease regulated by an interaction of genetic predisposition and environmental factors. To understand the genetic contribution in the development of diabetes, mice varying in their disease susceptibility were crossed with the obese and diabetes-prone New Zealand obese (NZO) mouse. Subsequent whole-genome sequence scans revealed one major quantitative trait loci (QTL), Nidd/DBA on chromosome 4, linked to elevated blood glucose and reduced plasma insulin and low levels of pancreatic insulin. Phenotypical characterization of congenic mice carrying 13.6 Mbp of the critical fragment of DBA mice displayed severe hyperglycemia and impaired glucose clearance at week 10, decreased glucose response in week 13, and loss of β-cells and pancreatic insulin in week 16. To identify the responsible gene variant(s), further congenic mice were generated and phenotyped, which resulted in a fragment of 3.3 Mbp that was sufficient to induce hyperglycemia. By combining transcriptome analysis and haplotype mapping, the number of putative responsible variant(s) was narrowed from initial 284 to 18 genes, including gene models and non-coding RNAs. Consideration of haplotype blocks reduced the number of candidate genes to four (Kti12, Osbpl9, Ttc39a, and Calr4) as potential T2D candidates as they display a differential expression in pancreatic islets and/or sequence variation. In conclusion, the integration of comparative analysis of multiple inbred populations such as haplotype mapping, transcriptomics, and sequence data substantially improved the mapping resolution of the diabetes QTL Nidd/DBA. Future studies are necessary to understand the exact role of the different candidates in β-cell function and their contribution in maintaining glycemic control.

H+/K+ATPase Expression in the Larynx of Laryngopharyngeal Reflux and Laryngeal Cancer Patients

OBJECTIVES

The gastric H+/K+ ATPase proton pump has previously been shown to be expressed in the human larynx, however its contribution to laryngopharyngeal reflux (LPR) signs, symptoms and associated diseases such as laryngeal cancer is unknown. Proton pump expression in the larynx of patients with LPR and laryngeal cancer was investigated herein. A human hypopharyngeal cell line expressing the proton pump was generated to investigate its effects.

STUDY DESIGN

In-vitro translational.

METHODS

Laryngeal biopsies were obtained from three LPR and eight LSCC patients. ATP4A, ATP4B and HRPT1 were assayed via qPCR. Human hypopharyngeal FaDu cell lines stably expressing proton pump were created using lentiviral transduction and examined via transmission electron microscopy and qPCR for genes associated with inflammation or laryngeal cancer.

RESULTS

Expression of ATP4A and ATP4B was detected in 3/3 LPR, 4/8 LSCC-tumor and 3/8 LSCC-adjacent specimens. Expression of ATP4A and ATP4B in FaDu elicited mitochondrial damage and expression of IL1B, PTGS2, and TNFA (P < .0001); expression of ATP4B alone did not.

CONCLUSIONS

Gastric proton pump subunits are expressed in the larynx of LPR and LSCC patients. Mitochondrial damage and changes in gene expression observed in cells expressing the full proton pump, absent in those expressing a single subunit, suggest that acid secretion by functional proton pumps expressed in upper airway mucosa may elicit local cell and molecular changes associated with inflammation and cancer.

LEVEL OF EVIDENCE

NA Laryngoscope, 131:130-135, 2021.

Distinct alterations of gut morphology and microbiota characterize accelerated diabetes onset in nonobese diabetic mice

The rising prevalence of type 1 diabetes (T1D) over the past decades has been linked to lifestyle changes, but the underlying mechanisms are largely unknown. Recent findings point to gut-associated mechanisms in the control of T1D pathogenesis. In nonobese diabetic (NOD) mice, a model of T1D, diabetes development accelerates after deletion of the Toll-like receptor 4 (TLR4). We hypothesized that altered intestinal functions contribute to metabolic alterations, which favor accelerated diabetes development in TLR4-deficient (TLR4^-/-) NOD mice. In 70-90-day-old normoglycemic (prediabetic) female NOD TLR4^+/+ and NOD TLR4^-/- mice, gut morphology and microbiome composition were analyzed. Parameters of lipid metabolism, glucose homeostasis, and mitochondrial respiratory activity were measured in vivo and ex vivo Compared with NOD TLR4^+/+ mice, NOD TLR4^-/- animals showed lower muscle mass of the small intestine, higher abundance of Bacteroidetes, and lower Firmicutes in the large intestine, along with lower levels of circulating short-chain fatty acids (SCFA). These changes are associated with higher body weight, hyperlipidemia, and severe insulin and glucose intolerance, all occurring before the onset of diabetes. These mice also exhibited insulin resistance-related abnormalities of energy metabolism, such as lower total respiratory exchange rates and higher hepatic oxidative capacity. Distinct alterations of gut morphology and microbiota composition associated with reduction of circulating SCFA may contribute to metabolic disorders promoting the progression of insulin-deficient diabetes/T1D development.