RNAi screening in primary human hepatocytes of genes implicated in genome-wide association studies for roles in type 2 diabetes identifies roles for CAMK1D and CDKAL1, among others, in hepatic glucose regulation

CaMK1D signalling in AgRP neurons promotes ghrelin-mediated food intake

Hypothalamic AgRP/NPY neurons are key players in the control of feeding behaviour. Ghrelin, a major orexigenic hormone, activates AgRP/NPY neurons to stimulate food intake and adiposity. However, cell-autonomous ghrelin-dependent signalling mechanisms in AgRP/NPY neurons remain poorly defined. Here we show that calcium/calmodulin-dependent protein kinase ID (CaMK1D), a genetic hot spot in type 2 diabetes, is activated upon ghrelin stimulation and acts in AgRP/NPY neurons to mediate ghrelin-dependent food intake. Global Camk1d-knockout male mice are resistant to ghrelin, gain less body weight and are protected against high-fat-diet-induced obesity. Deletion of Camk1d in AgRP/NPY, but not in POMC, neurons is sufficient to recapitulate above phenotypes. In response to ghrelin, lack of CaMK1D attenuates phosphorylation of CREB and CREB-dependent expression of the orexigenic neuropeptides AgRP/NPY in fibre projections to the paraventricular nucleus (PVN). Hence, CaMK1D links ghrelin action to transcriptional control of orexigenic neuropeptide availability in AgRP neurons.

TRIM21 ameliorates hepatic glucose and lipid metabolic disorders in type 2 diabetes mellitus by ubiquitination of PEPCK1 and FASN

Hepatic glucose and lipid metabolism disorders promote the development and progression of type 2 diabetes mellitus (T2DM), yet the underlying mechanisms are not fully understood. Here, we identify tripartite motif-containing protein 21 (TRIM21), a class IV TRIM family member, as a pivotal regulator of hepatic metabolism in T2DM for the first time. Bioinformatic analysis suggests that TRIM21 expression is significantly reduced in T2DM patients. Intriguingly, in a mouse model of obese diabetes, TRIM21 expression is predominantly reduced in the liver rather than in other metabolic organs. It is further demonstrated that hepatic overexpression of TRIM21 significantly ameliorates glucose intolerance, insulin resistance, hepatic steatosis, and dyslipidemia in obese diabetic mice. In contrast, the knockdown of TRIM21 promotes glucose intolerance, insulin resistance, and triglyceride accumulation. Mechanistically, both phosphoenolpyruvate carboxykinase 1 (PEPCK1) and fatty acid synthase (FASN) are the hepatic targets of TRIM21. We revealed that TRIM21 promotes the degradation of PEPCK1 and FASN through a direct protein-protein interaction mediated K48-linked ubiquitination. Notably, overexpression of PEPCK1 and FASN essentially abolished the beneficial effects achieved by TRIM21 overexpression in obese diabetic mice. Overall, our data demonstrate that TRIM21 is a novel regulator of hepatic metabolic disorder, and suggest TRIM21 as a promising therapeutic target for T2DM.

Hepatocyte Nuclear Factor 4 alpha 2 Messenger RNA Reprograms Liver-Enriched Transcription Factors and Functional Proteins in End-Stage Cirrhotic Human Hepatocytes

The only definitive therapy for end-stage liver disease is whole-organ transplantation. The success of this intervention is severely limited by the complexity of the surgery, the cost of patient care, the need for long-term immunosuppression, and the shortage of donor organs. In rodents and humans, end-stage degeneration of hepatocyte function is associated with disruption of the liver-specific transcriptional network and a nearly complete loss of promoter P1-driven hepatocyte nuclear factor 4-alpha (P1-HNF4α) activity. Re-expression of HNF4α2, the predominant P1-HNF4α, reinstates the transcriptional network, normalizes the genes important for hepatocyte function, and reverses liver failure in rodents. In this study, we tested the effectiveness of supplementary expression of human HNF4α2 messenger RNA (mRNA) in primary human hepatocytes isolated from explanted livers of patients who underwent transplant for end-stage irreversibly decompensated liver failure (Child-Pugh B, C) resulting from alcohol-mediated cirrhosis and nonalcoholic steatohepatitis. Re-expression of HNF4α2 in decompensated cirrhotic human hepatocytes corrects the disrupted transcriptional network and normalizes the expression of genes important for hepatocyte function, improving liver-specific protein expression. End-stage liver disease in humans is associated with both loss of P1-HNF4α expression and failure of its localization to the nucleus. We found that while HNF4α2 re-expression increased the amount of P1-HNF4α protein in hepatocytes, it did not alter the ability of hepatocytes to localize P1-HNF4α to their nuclei. Conclusion: Re-expression of HNF4α2 mRNA in livers of patients with end-stage disease may be an effective therapy for terminal liver failure that would circumvent the need for organ transplantation. The efficacy of this strategy may be enhanced by discovering the cause for loss of nuclear P1-HNF4α localization in end-stage cirrhosis, a process not found in rodent studies.

Backbone resonance assignments of the catalytic and regulatory domains of Ca2+/calmodulin-dependent protein kinase 1D

The CaMK subfamily of Ser/Thr kinases are regulated by calmodulin interactions with their C-terminal regions. They are exemplified by Ca²⁺/calmodulin dependent protein kinase 1δ which is known as CaMK1D, CaMKIδ or CKLiK. CaMK1D mediates intracellular signalling downstream of Ca²⁺ influx and thereby exhibits amplifications of Ca²⁺signals and polymorphisms that have been implicated in breast cancer and diabetes. Here we report the backbone ¹H, ¹³C, ¹⁵N assignments of the 38 kDa human CaMK1D protein in its free state, including both the canonical bi-lobed kinase fold as well as the autoinhibitory and calmodulin binding domains.

Transcriptome-wide analysis of PGC-1α-binding RNAs identifies genes linked to glucagon metabolic action

The peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α) is a transcriptional coactivator that controls expression of metabolic/energetic genes, programming cellular responses to nutrient and environmental adaptations such as fasting, cold, or exercise. Unlike other coactivators, PGC-1α contains protein domains involved in RNA regulation such as serine/arginine (SR) and RNA recognition motifs (RRMs). However, the RNA targets of PGC-1α and how they pertain to metabolism are unknown. To address this, we performed enhanced ultraviolet (UV) cross-linking and immunoprecipitation followed by sequencing (eCLIP-seq) in primary hepatocytes induced with glucagon. A large fraction of RNAs bound to PGC-1α were intronic sequences of genes involved in transcriptional, signaling, or metabolic function linked to glucagon and fasting responses, but were not the canonical direct transcriptional PGC-1α targets such as OXPHOS or gluconeogenic genes. Among the top-scoring RNA sequences bound to PGC-1α were Foxo1, Camk1δ, Per1, Klf15, Pln4, Cluh, Trpc5, Gfra1, and Slc25a25 PGC-1α depletion decreased a fraction of these glucagon-induced messenger RNA (mRNA) transcript levels. Importantly, knockdown of several of these genes affected glucagon-dependent glucose production, a PGC-1α-regulated metabolic pathway. These studies show that PGC-1α binds to intronic RNA sequences, some of them controlling transcript levels associated with glucagon action.

Discovery of Highly Selective Inhibitors of Calmodulin-Dependent Kinases That Restore Insulin Sensitivity in the Diet-Induced Obesity in Vivo Mouse Model

Polymorphisms in the region of the calmodulin-dependent kinase isoform D (CaMK1D) gene are associated with increased incidence of diabetes, with the most common polymorphism resulting in increased recognition by transcription factors and increased protein expression. While reducing CaMK1D expression has a potentially beneficial effect on glucose processing in human hepatocytes, there are no known selective inhibitors of CaMK1 kinases that can be used to validate or translate these findings. Here we describe the development of a series of potent, selective, and drug-like CaMK1 inhibitors that are able to provide significant free target cover in mouse models and are therefore useful as in vivo tool compounds. Our results show that a lead compound from this series improves insulin sensitivity and glucose control in the diet-induced obesity mouse model after both acute and chronic administration, providing the first in vivo validation of CaMK1D as a target for diabetes therapeutics.

An integrated investigation of oocyte developmental competence: expression of key genes in human cumulus cells, morphokinetics of early divisions, blastulation, and euploidy

PURPOSE

To investigate the association of cumulus cell (CC)-related expression of a selected cluster of key genes (PTGS2, CAMK1D, HAS2, STC1, and EFNB2) with embryo development to blastocyst.

METHODS

Exploratory study at a private clinic. Eighteen advanced maternal age patients were enrolled (37.3 ± 4.0 years). Seventy-five cumuli were collected, whose oocytes resulted in either developmental arrest (N = 33) or blastocyst formation (N = 42). The noninvasive CC gene expression was combined with time-lapse morphokinetic parameters and, for blastocysts, with qPCR-based aneuploidy testing on trophectoderm biopsies.

RESULTS

The detection rate was 100% for all transcripts, but STC1 (96%) and CAMK1D (89%). Among amplified assays, CC mean expression levels of CAMK1D, PTGS2, and HAS2 were higher from oocytes that developed to blastocyst. No difference in CC key gene expression was reported between euploid (N = 21) and aneuploid (N = 21) blastocysts. Some timings of early embryo development were faster in embryos developing to blastocyst (time of pronuclei appearance and fading, division to two- and four-cells, first and second cell cycles). However, the generalized linear models outlined increasing CAMK1D expression levels as the strongest parameter associated with oocytes' developmental potential from both a general (AUC = 0.78 among amplified samples) and an intrapatient perspectives (AUC = 0.9 among patients obtaining ≥ 2 zygotes from the cohort with different developmental outcomes).

CONCLUSIONS

CAMK1D level of expression in CCs associated with blastocyst development. If confirmed from larger studies in wider populations of patients, the investigation of CC key gene expression might suit IVF clinics not adopting blastocyst culture. Future investigations should clarify the role of CAMK1D in ovarian physiology and could provide novel insights on how oocytes gain competence during folliculogenesis.

Genetic Variants Associated With Obesity and Insulin Resistance in Hispanic Boys With Nonalcoholic Fatty Liver Disease

BACKGROUND AND OBJECTIVES

Nonalcoholic fatty liver disease (NAFLD) disproportionately affects Hispanic boys. Further, obesity and insulin resistance are major risk factors for NAFLD. No gene localization studies had been performed on children with biopsy-proven NAFLD. This study aims to identify genomic variants associated with increased adiposity and insulin resistance in a population of children with varying histologic severity of NAFLD.

METHODS

We conducted a genome-wide association scan (GWAS) including 624,297 single-nucleotide polymorphisms (SNPs) distributed among all 22 autosomal chromosomes in 234 Hispanic boys (up to 18 years of age) who were consecutively recruited in a prospective cohort study in the Nonalcoholic Steatohepatitis Clinical Research Network Studies. Traits were examined quantitatively using linear regression. SNPs with P value <10 and a minor allele frequency >5% were considered potentially significant.

RESULTS

Evaluated subjects had a median age of 12.0 years, body mass index (BMI) of 31.4, and hemoglobin A1C (Hgb A1C) of 5.3. The prevalence of NAFL, borderline NASH, and definite NASH were 23%, 53%, and 22%, respectively. The GWAS identified 10 SNPs that were associated with BMI z score, 6 within chromosome 2, and 1 within CAMK1D, which has a potential role in liver gluconeogenesis. In addition, the GWAS identified 9 novel variants associated with insulin resistance: HOMA-IR (6) and HbA1c (3).

CONCLUSIONS

This study of Hispanic boys with biopsy-proven NAFLD with increased risk for the metabolic syndrome revealed novel genetic variants that are associated with obesity and insulin resistance.

High-Content Screening Approaches That Minimize Confounding Factors in RNAi, CRISPR, and Small Molecule Screening

Screening arrayed libraries of reagents, particularly small molecules began as a vehicle for drug discovery, but the in last few years it has become a cornerstone of biological investigation, joining RNAi and CRISPR as methods for elucidating functional relationships that could not be anticipated, and illustrating the mechanisms behind basic and disease biology, and therapeutic resistance. However, these approaches share some common challenges, especially with respect to specificity or selectivity of the reagents as they are scaled to large protein families or the genome. High-content screening (HCS) has emerged as an important complement to screening, mostly the result of a wide array of specific molecular events, such as protein kinase and transcription factor activation, morphological changes associated with stem cell differentiation or the epithelial-mesenchymal transition of cancer cells. Beyond the range of cellular events that can be screened by HCS, image-based screening introduces new processes for differentiating between specific and nonspecific effects on cells. This chapter introduces these complexities and discusses strategies available in image-based screening that can mitigate the challenges they can bring to screening.

Contribution of independent and pleiotropic genetic effects in the metabolic syndrome in a hypertensive rat

Hypertension is a major risk factor for cardiovascular disease, Type 2 diabetes, and end organ failure, and is often found concomitant with disorders characteristic of the Metabolic Syndrome (MetS), including obesity, dyslipidemia, and insulin resistance. While the associated features often occur together, the pathway(s) or mechanism(s) linking hypertension in MetS are not well understood. Previous work determined that genetic variation on rat chromosome 17 (RNO17) contributes to several MetS-defining traits (including hypertension, obesity, and dyslipidemia) in the Lyon Hypertensive (LH) rat, a genetically determined MetS model. We hypothesized that at least some of the traits on RNO17 are controlled by a single gene with pleiotropic effects. To address this hypothesis, consomic and congenic strains were developed, whereby a defined fragment of RNO17 from the LH rat was substituted with the control Lyon Normotensive (LN) rat, and MetS phenotypes were measured in the resultant progeny. Compared to LH rats, LH-17^LN consomic rats have significantly reduced body weight, blood pressure, and lipid profiles. A congenic strain (LH-17^LNc), with a substituted fragment at the distal end of RNO17 (17q12.3; 74–97 Mb; rn4 assembly), showed differences from the LH rat in blood pressure and serum total cholesterol and triglycerides. Interestingly, there was no difference in body weight between the LH-17^LNc and the parental LH rat. These data indicate that blood pressure and serum lipids are regulated by a gene(s) in the distal congenic interval, and could be due to pleiotropy. The data also indicate that body weight is not determined by the same gene(s) at this locus. Interestingly, only two small haplotypes spanning a total of approximately 0.5 Mb differ between the LH and LN genomes in the congenic interval. Genes in these haplotypes are strong candidate genes for causing dyslipidemia in the LH rat. Overall, MetS, even in a simplified genetic model such as the LH-17^LN rat, is likely due to both independent and pleiotropic gene effects.

Maternal Obesity in Pregnancy Developmentally Programs Adipose Tissue Inflammation in Young, Lean Male Mice Offspring

Obesity during pregnancy has a long-term effect on the health of the offspring including risk of developing the metabolic syndrome. Using a mouse model of maternal diet-induced obesity, we employed a genome-wide approach to investigate the microRNA (miRNA) and miRNA transcription profile in adipose tissue to understand mechanisms through which this occurs. Male offspring of diet-induced obese mothers, fed a control diet from weaning, showed no differences in body weight or adiposity at 8 weeks of age. However, offspring from the obese dams had up-regulated cytokine (Tnfα; P < .05) and chemokine (Ccl2 and Ccl7; P < .05) signaling in their adipose tissue. This was accompanied by reduced expression of miR-706, which we showed can directly regulate translation of the inflammatory proteins IL-33 (41% up-regulated; P < .05) and calcium/calmodulin-dependent protein kinase 1D (30% up-regulated; P < .01). We conclude that exposure to obesity during development primes an inflammatory environment in adipose tissue that is independent of offspring adiposity. Programming of adipose tissue miRNAs that regulate expression of inflammatory signaling molecules may be a contributing mechanism.

Genome-wide investigation of schizophrenia associated plasma Ndel1 enzyme activity

Ndel1 is a DISC1-interacting oligopeptidase that cleaves in vitro neuropeptides as neurotensin and bradykinin, and which has been associated with both neuronal migration and neurite outgrowth. We previously reported that plasma Ndel1 enzyme activity is lower in patients with schizophrenia (SCZ) compared to healthy controls (HCs). To our knowledge, no previous study has investigated the genetic factors associated with the plasma Ndel1 enzyme activity. In the current analyses, samples from 83 SCZ patients and 92 control subjects that were assayed for plasma Ndel1 enzyme activity were genotyped on Illumina Omni Express arrays. A genetic relationship matrix using genome-wide information was then used for ancestry correction, and association statistics were calculated genome-wide. Ndel1 enzyme activity was significantly lower in patients with SCZ (t=4.9; p<0.001) and was found to be associated with CAMK1D, MAGI2, CCDC25, and GABGR3, at a level of suggestive significance (p<10(-6)), independent of the clinical status. Then, we performed a model to investigate the observed differences for case/control measures. 2 SNPs at region 1p22.2 reached the p<10(-7) level. ZFPM2 and MAD1L1 were the only two genes with more than one hit at 10(-6) order of p value. Therefore, Ndel1 enzyme activity is a complex trait influenced by many different genetic variants that may contribute to SCZ physiopathology.

Genome-wide association studies in the Japanese population identify seven novel loci for type 2 diabetes

Genome-wide association studies (GWAS) have identified more than 80 susceptibility loci for type 2 diabetes (T2D), but most of its heritability still remains to be elucidated. In this study, we conducted a meta-analysis of GWAS for T2D in the Japanese population. Combined data from discovery and subsequent validation analyses (23,399 T2D cases and 31,722 controls) identify 7 new loci with genome-wide significance (P<5 × 10(-8)), rs1116357 near CCDC85A, rs147538848 in FAM60A, rs1575972 near DMRTA1, rs9309245 near ASB3, rs67156297 near ATP8B2, rs7107784 near MIR4686 and rs67839313 near INAFM2. Of these, the association of 4 loci with T2D is replicated in multi-ethnic populations other than Japanese (up to 65,936 T2Ds and 158,030 controls, P<0.007). These results indicate that expansion of single ethnic GWAS is still useful to identify novel susceptibility loci to complex traits not only for ethnicity-specific loci but also for common loci across different ethnicities.

miRNAs: early prognostic biomarkers for Type 2 diabetes mellitus?

Type 2 diabetes mellitus (T2DM) has reached epidemic proportions and is associated with peripheral insulin resistance. The currently used therapies aim to delay progression of T2DM. Their efficacy could drastically be improved if implemented at earlier stages. Classical diagnostic markers (blood glucose and HbA1C) are generally detected once metabolic imbalance has already set in. Therefore, development of biomarkers for early diagnosis would help identify individuals at risk for developing T2DM. Along with genetic predisposition, epigenetics also plays a major role in T2DM development. In this review, we discuss the potential role of early diagnostic markers such as circulating miRNAs, studies done so far and challenges to be considered while taking into account the novel role of miRNAs as prognostic biomarkers.

Association of PCK1 with Body Mass Index and Other Metabolic Features in Patients With Psychotropic Treatments

Weight gain is a major health problem among psychiatric populations. It implicates several receptors and hormones involved in energy balance and metabolism. Phosphoenolpyruvate carboxykinase 1 is a rate-controlling enzyme involved in gluconeogenesis, glyceroneogenesis and cataplerosis and has been related to obesity and diabetes phenotypes in animals and humans. The aim of this study was to investigate the association of phosphoenolpyruvate carboxykinase 1 polymorphisms with metabolic traits in psychiatric patients treated with psychotropic drugs inducing weight gain and in general population samples. One polymorphism (rs11552145G > A) significantly associated with body mass index in the psychiatric discovery sample (n = 478) was replicated in 2 other psychiatric samples (n1 = 168, n2 = 188), with AA-genotype carriers having lower body mass index as compared to G-allele carriers. Stronger associations were found among women younger than 45 years carrying AA-genotype as compared to G-allele carriers (-2.25 kg/m, n = 151, P = 0.009) and in the discovery sample (-2.20 kg/m, n = 423, P = 0.0004). In the discovery sample for which metabolic parameters were available, AA-genotype showed lower waist circumference (-6.86 cm, P = 0.008) and triglycerides levels (-5.58 mg/100 mL, P < 0.002) when compared to G-allele carriers. Finally, waist-to-hip ratio was associated with rs6070157 (proxy of rs11552145, r = 0.99) in a population-based sample (N = 123,865, P = 0.022). Our results suggest an association of rs11552145G > A polymorphism with metabolic-related traits, especially in psychiatric populations and in women younger than 45 years.

Identification of a regulatory variant that binds FOXA1 and FOXA2 at the CDC123/CAMK1D type 2 diabetes GWAS locus

Many of the type 2 diabetes loci identified through genome-wide association studies localize to non-protein-coding intronic and intergenic regions and likely contain variants that regulate gene transcription. The CDC123/CAMK1D type 2 diabetes association signal on chromosome 10 spans an intergenic region between CDC123 and CAMK1D and also overlaps the CDC123 3′UTR. To gain insight into the molecular mechanisms underlying the association signal, we used open chromatin, histone modifications and transcription factor ChIP-seq data sets from type 2 diabetes-relevant cell types to identify SNPs overlapping predicted regulatory regions. Two regions containing type 2 diabetes-associated variants were tested for enhancer activity using luciferase reporter assays. One SNP, rs11257655, displayed allelic differences in transcriptional enhancer activity in 832/13 and MIN6 insulinoma cells as well as in human HepG2 hepatocellular carcinoma cells. The rs11257655 risk allele T showed greater transcriptional activity than the non-risk allele C in all cell types tested. Using electromobility shift and supershift assays we demonstrated that the rs11257655 risk allele showed allele-specific binding to FOXA1 and FOXA2. We validated FOXA1 and FOXA2 enrichment at the rs11257655 risk allele using allele-specific ChIP in human islets. These results suggest that rs11257655 affects transcriptional activity through altered binding of a protein complex that includes FOXA1 and FOXA2, providing a potential molecular mechanism at this GWAS locus.

GWAS have identified more than 1200 loci contributing to risk of disease, including more than 70 loci associated with type 2 diabetes. With a majority of associated variants localized to non-coding regions of the genome, focus has moved to identifying the functional variants explaining the association signals. One mechanism by which variants may act is to affect activity of enhancer elements regulating target gene expression. In this study, we take advantage of recent advances in genome-wide annotation of human regulatory elements to prioritize candidate functional variants at the CDC123/CAMK1D locus. We identify two T2D-associated variants that overlap predicted regulatory enhancer elements. We demonstrate that one variant, rs11257655, shows allele-specific transcriptional enhancer activity in mammalian cell lines relevant to type 2 diabetes. We also show differential protein-DNA binding suggesting that the rs11257655 type 2 diabetes- risk allele increased transcriptional activity through binding a protein complex that includes FOXA1 and FOXA2. This study demonstrates that genome-wide maps of regulatory elements are a useful resource to guide identification of variants differentially affecting transcriptional activity and provides insight into molecular mechanisms underlying a T2D susceptibility locus.