Carboxylesterase 1 gene duplication and mRNA expression in adipose tissue are linked to obesity and metabolic function

Prenylcysteine Oxidase 1 Is a Key Regulator of Adipogenesis

The process of adipogenesis involves the differentiation of preadipocytes into mature adipocytes. Excessive adipogenesis promotes obesity, a condition that increasingly threatens global health and contributes to the rapid rise of obesity-related diseases. We have recently shown that prenylcysteine oxidase 1 (PCYOX1) is a regulator of atherosclerosis-disease mechanisms, which acts through mechanisms not exclusively related to its pro-oxidant activity. To address the role of PCYOX1 in the adipogenic process, we extended our previous observations confirming that Pcyox1−/−/Apoe−/− mice fed a high-fat diet for 8 or 12 weeks showed significantly lower body weight, when compared to Pcyox1+/+/Apoe−/− mice, due to an evident reduction in visceral adipose content. We herein assessed the role of PCYOX1 in adipogenesis. Here, we found that PCYOX1 is expressed in adipose tissue, and, independently from its pro-oxidant enzymatic activity, is critical for adipogenesis. Pcyox1 gene silencing completely prevented the differentiation of 3T3-L1 preadipocytes, by acting as an upstream regulator of several key players, such as FABP4, PPARγ, C/EBPα. Proteomic analysis, performed by quantitative label-free mass spectrometry, further strengthened the role of PCYOX1 in adipogenesis by expanding the list of its downstream targets. Finally, the absence of Pcyox1 reduces the inflammatory markers in adipose tissue. These findings render PCYOX1 a novel adipogenic factor with possible pathophysiological or therapeutic potential.

Carboxylesterase 1 family knockout alters drug disposition and lipid metabolism

The mammalian carboxylesterase 1 (Ces1/CES1) family comprises several enzymes that hydrolyze many xenobiotic chemicals and endogenous lipids. To investigate the pharmacological and physiological roles of Ces1/CES1, we generated Ces1 cluster knockout (Ces1 ^-/- ) mice, and a hepatic human CES1 transgenic model in the Ces1 ^-/- background (TgCES1). Ces1 ^-/- mice displayed profoundly decreased conversion of the anticancer prodrug irinotecan to SN-38 in plasma and tissues. TgCES1 mice exhibited enhanced metabolism of irinotecan to SN-38 in liver and kidney. Ces1 and hCES1 activity increased irinotecan toxicity, likely by enhancing the formation of pharmacodynamically active SN-38. Ces1 ^-/- mice also showed markedly increased capecitabine plasma exposure, which was moderately decreased in TgCES1 mice. Ces1 ^-/- mice were overweight with increased adipose tissue, white adipose tissue inflammation (in males), a higher lipid load in brown adipose tissue, and impaired blood glucose tolerance (in males). These phenotypes were mostly reversed in TgCES1 mice. TgCES1 mice displayed increased triglyceride secretion from liver to plasma, together with higher triglyceride levels in the male liver. These results indicate that the carboxylesterase 1 family plays essential roles in drug and lipid metabolism and detoxification. Ces1 ^-/- and TgCES1 mice will provide excellent tools for further study of the in vivo functions of Ces1/CES1 enzymes.

Silencing long noncoding RNA-CES1P1 suppresses glomerular endothelial cell inflammation in diabetic nephropathy

Diabetic nephropathy (DN) has become the main cause of end-stage renal disease worldwide. Inflammation is associated with the occurrence and development of DN, and long noncoding RNAs (lncRNAs) are involved in the regulation of inflammatory processes. This study aims to determine the role and mechanism of lncRNA-CES1P1 in DN.C57BL/6 mice and human umbilical vein endothelial cells (HUVECs) were used for this experimental study. In vivo experimental intraperitoneal injection of streptozotocin (STZ) to construct a diabetes mellitus (DM) model in C57BL/6 mice caused increased expression of lncRNA-CES1P1, decreased expression of miR-214-3p in kidney tissue, and produced renal inflammation and proteinuria. Exogenous knockdown of lncRNA-CES1P1 expression decreased renal inflammatory infiltration. In vitro experiments using high glucose (HG) stimulation of HUVECs cell revealed increased expression of lncRNA-CES1P1, decreased expression of miR-214-3p, and increased expression of the inflammatory factors IL-17, IκB, NF-κB, and IL-6. Luciferase reporter assays showed direct targets of miR-214-3p interaction with lncRNA-CES1P1 and IL-17. These results suggest that hyperglycemia represses miR-214-3p by inducing lncRNA-CES1P1, which promotes the expression of the inflammatory factors IL-17, IκB, NF-κB and IL-6 ultimately leading to the development of DN. Interfering with lncRNA-CES1P1 can reduce hyperglycemia-induced DN.

Pharmacogenomics variants are associated with BMI differences between individuals with bipolar and other psychiatric disorders

Aim: Regardless of the plethora of next-generation sequencing studies in the field of pharmacogenomics (PGx), the potential effect of covariate variables on PGx response within deeply phenotyped cohorts remains unexplored. Materials & methods: We explored with advanced statistical methods the potential influence of BMI, as a covariate variable, on PGx response in a Greek cohort with psychiatric disorders. Results: Nine PGx variants within UGT1A6, SLC22A4, GSTP1, CYP4B1, CES1, SLC29A3 and DPYD were associated with altered BMI in different psychiatric disorder groups. Carriers of rs2070959 (UGT1A6), rs199861210 (SLC29A3) and rs2297595 (DPYD) were also characterized by significant changes in the mean BMI, depending on the presence of psychiatric disorders. Conclusion: Specific PGx variants are significantly associated with BMI in a Greek cohort with psychiatric disorders.

Is platelet responsiveness to clopidogrel attenuated in overweight or obese patients and why? A reverse translational study in mice

BACKGROUND AND PURPOSE

Overweight or obese patients exhibit poorer platelet responses to clopidogrel. However, the mechanisms behind this phenotype remain to be elucidated. Here, we sought to discover whether and why obesity could affect the metabolic activation of and/or platelet response to clopidogrel in obese patients and high-fat diet-induced obese mice.

EXPERIMENTAL APPROACH

A post hoc stratified analysis of an observational clinical study was performed to investigate changes in residual platelet reactivity with increasing body weight in patients taking clopidogrel. Furthermore, high-fat diet-induced obese mice were used to reveal alterations in systemic exposure of clopidogrel thiol active metabolite H4, ADP-induced platelet activation and aggregation, the expression of genes involved in the metabolic activation of clopidogrel, count of circulating reticulated and mature platelets, and proliferation profiles of megakaryocytes in bone marrow. The relevant genes and potential signalling pathways were predicted and enriched according to the GEO datasets available from obese patients.

KEY RESULTS

Obese patients exhibited significantly attenuated antiplatelet effects of clopidogrel. In diet-induced obese mice, systemic exposure of clopidogrel active metabolite H4 was reduced but that of its hydrolytic metabolite was increased due to down-regulation of certain P450s but up-regulation of carboxylesterase-1 in the liver. Moreover, enhanced proliferation of megakaryocytes and elevated platelet count also contributed.

CONCLUSION AND IMPLICATIONS

Obesity attenuated metabolic activation of clopidogrel and increased counts of circulating reticulated and mature platelets, leading to impaired platelet responsiveness to the drug in mice, suggesting that clopidogrel dosage may need to be adjusted adequately in overweight or obese patients.

Physiologically-Based Pharmacokinetic-Pharmacodynamics Model Characterizing CYP2C19 Polymorphisms to Predict Clopidogrel Pharmacokinetics and Its Anti-Platelet Aggregation Effect Following Oral Administration to Coronary Artery Disease Patients With or Without Diabetes

Background and Objective: Clopidogrel (CLOP) is commonly used in coronary artery disease (CAD) patients with or without diabetes (DM), but these patients often suffer CLOP resistance, especially those with diabetes. This study was aimed to develop a physiologically-based pharmacokinetic-pharmacodynamic (PBPK-PD) model to describe the pharmacokinetics and pharmacodynamics of clopidogrel active metabolite (CLOP-AM) in CAD patients with or without DM.

Methods: The PBPK-PD model was first established and validated in healthy subjects and then in CAD patients with or without DM. The influences of CYP2C19, CYP2C9, CYP3A4, carboxylesterase 1 (CES1), gastrointestinal transit rates (K_t,i) and platelets response to CLOP-AM (k_irre) on predicted pharmacokinetics and pharmacodynamics were investigated, followed with their individual and integrated effects on CLOP-AM pharmacokinetics due to changes in DM status.

Results: Most predictions fell within 0.5–2.0 folds of observations, indicating successful predictions. Sensitivity analysis showed that contributions of interested factors to pharmacodynamics were CES1> k_irre> K_t,i> CYP2C19 > CYP3A4> CYP2C9. Mimicked analysis showed that the decreased exposure of CLOP-AM by DM was mainly attributed to increased CES1 activity, followed by decreased CYP2C19 activity.

Conclusion: The pharmacokinetics and pharmacodynamics of CLOP-AM were successfully predicted using the developed PBPK-PD model. Clopidogrel resistance by DM was the integrated effects of altered K_t,i, CYP2C19, CYP3A4, CES1 and k_irre.

Carboxylesterase 1 and Precision Pharmacotherapy: Pharmacogenetics and Nongenetic Regulators

Carboxylesterase (CES) 1 is the most abundant drug-metabolizing enzyme in human livers, comprising approximately 1% of the entire liver proteome. CES1 is responsible for 80%-95% of total hydrolytic activity in the liver and plays a crucial role in the metabolism of a wide range of drugs (especially ester-prodrugs), pesticides, environmental pollutants, and endogenous compounds. Expression and activity of CES1 vary markedly among individuals, which is a major contributing factor to interindividual variability in the pharmacokinetics (PK) and pharmacodynamics (PD) of drugs metabolized by CES1. Both genetic and nongenetic factors contribute to CES1 variability. Here, we discuss genetic polymorphisms, including single-nucleotide polymorphisms (SNPs), and copy number variants and nongenetic contributors, such as developmental status, genders, and drug-drug interactions, that could influence CES1 functionality and the PK and PD of CES1 substrates. Currently, the loss-of-function SNP G143E (rs71647871) is the only clinically significant CES1 variant identified to date, and alcohol is the only potent CES1 inhibitor that could alter the therapeutic outcomes of CES1 substrate medications. However, G143E and alcohol can only explain a small portion of the interindividual variability in the CES1 function. A better understanding of the regulation of CES1 expression and activity and identification of biomarkers for CES1 function in vivo could lead to the development of a precision pharmacotherapy strategy to improve the efficacy and safety of many CES1 substrate drugs. SIGNIFICANCE STATEMENT: The clinical relevance of CES1 has been well demonstrated in various clinical trials. Genetic and nongenetic regulators can affect CES1 expression and activity, resulting in the alteration of the metabolism and clinical outcome of CES1 substrate drugs, such as methylphenidate and clopidogrel. Predicting the hepatic CES1 function can provide clinical guidance to optimize pharmacotherapy of numerous medications metabolized by CES1.

Adipocyte abundances of CES1, CRYAB, ENO1 and GANAB are modified in-vitro by glucose restriction and are associated with cellular remodelling during weight regain

Long-term weight loss maintenance is a problem of overweight and obesity. Changes of gene expression during weight loss (WL) by calorie restriction (CR) are linked to the risk of weight regain (WR). However, detailed information on genes/proteins involved in the mechanism is still lacking. Therefore, we developed an in-vitro model system for glucose restriction (GR) and refeeding (RF) to uncover proteome differences between GR with RF vs normal feeding, of which we explored the relation with WR after WL. Human Simpson-Golabi-Behmel Syndrome cells were subjected to changing levels of glucose to mimic the condition of CR and RF. Proteome profiling was performed by liquid chromatography tandem mass spectrometry. This in-vitro model revealed 44 proteins differentially expressed after GR and RF versus feeding including proteins of the focal adhesions. Four proteins showed a persistent up- or down-regulation: liver carboxylesterase (CES1), mitochondrial superoxide dismutase [Mn] (SOD2), alpha-crystallin B-chain (CRYAB), alpha-enolase (ENO1). In-vivo weight loss-induced RNA expression changes linked CES1, CRYAB and ENO1 to WR. Moreover, of these 44 proteins, CES1 and glucosidase II alpha subunit (GANAB) during follow up correlated with WR. Correlation clustering of in-vivo protein expression data indicated an interaction of these proteins with structural components of the focal adhesions and cytoplasmic filaments in the adipocytes.

Copy number variation profiling in pharmacogenes using panel-based exome resequencing and correlation to human liver expression

Structural variants including copy number variations (CNV) have gained widespread attention, especially in pharmacogenomics but for several genes functional relevance and clinical evidence are still lacking. Detection of CNVs in next-generation sequencing data is challenging but offers widespread applications. We developed a cohort-based CNV detection workflow to extract CNVs from read counts of targeted NGS of 340 genes involved in absorption, distribution, metabolism and excretion (ADME) of drugs. We applied our method to 150 human liver tissue samples and correlated identified CNVs to mRNA expression levels. In total, we identified 445 deletions (73%) and 167 duplications (27%) in 36 pharmacogenes including all well-known CNVs of CYPs, GSTs, SULTs, UGTs, numerous described rare CNVs of CYP2E1, SLC16A3 or UGT2B15 as well as novel observations, e.g., for SLC22A12, SLC22A17 and GPS2 (G Protein Pathway Suppressor 2). We were able to fine-map complex CNVs of CYP2A6 and CYP2D6 with exon resolution. Correlation analysis confirmed known expression patterns for common CNVs and suggested an influence on expression variability for some rare CNVs. Our straightforward CNV detection workflow can be easily applied to any NGS coverage data and helped to analyze CNVs in an ADME-NGS panel of 340 pharmacogenes to improve genotype-phenotype correlations.

Sex differences in gene expression related to antipsychotic induced weight gain

Antipsychotics are crucial for the treatment of schizophrenia and contribute to weight gain in psychosis, particularly during early phases. Antipsychotic Induced Weight Gain (AIWG) might contribute to reduce the quality of life, drug compliance and to increase mortality. To characterize sex differences of gene expression related to AIWG, we sequenced total mRNA from blood samples of schizophrenia patients, before and after 3 months of antipsychotic-treatment. We analyzed schizophrenia patients according to their sex (38 males and 39 females) and their BMI increase after medication, characterizing the differential gene expression before and after medication. Individuals in each group were categorized in patients who gain weight and those whose do not gain weight. The “weight gain” groups included patients with an increase of body mass index (BMI) > 1.0 points (27 males and 23 females with a median BMI increase of 2.68 and 2.32 respectively). The “no weight gain” groups included patients with a change of BMI between < 1.0 and > -1.0 points (11 males and 16 females with a median BMI increase of 0.21 and 0.16 respectively). The males had 331 genes with significant differential expression in the weight gain group and 24 genes in the no weight gain group. The females had 119 genes with significant differential expression in the weight gain group and 75 genes in the no weight gain group. Both weight gain groups were significantly enriched with “obesity” genes (Fisher; p = 1.1E-09 and p = 0.0001 respectively), according to the Gene Reference into Function (GeneRIF) database.In conclusion, we characterized genes with differential expression associated to AIWG that are specific to males, to females and common to both sexes. These genes are good candidates to depict the biological processes involved in AIWG and provide additional evidence of the genetic links between weight gain and the immune system.

Breed Differences in Pig Liver Esterase (PLE) between Tongcheng (Chinese Local Breed) and Large White Pigs

Human carboxylesterases has been proven to be age and race-related and a sound basis of clinical medication. PLE involve in signal transduction and highly catalyze hydrolysis. Therefore, the expression level of PLE most probably exist age and breed difference and lead to significant differences of pharmacology and physiology. Four age groups of Tongcheng (TC) and Large White (LW) pigs were selected to explore PLE breed and age differences, and it was found that PLE mRNA was most abundant in liver in both breeds. In liver, PLE levels and hydrolytic activities increased with age, and PLE levels (except for 3 month) and the hydrolytic activities were higher in LW than in TC across all age groups. Abundance of PLE isoenzymes was obvious different between breeds and among age groups. The most abundant PLE isoenzyme in LW and TC pigs was PLE-A1 (all age groups) and PLE-B9 (three early age groups) or PLE-G3 (adult groups), respectively. 103 new PLE isoenzymes were found, and 55 high-frequency PLE isoenzymes were accordingly classified into seven categories (A-G). The results of this research provide a necessary basis not only for clinical medication of pigs but also for pig breeding purposes.

Isoflavones enhance pharmacokinetic exposure of active lovastatin acid via the upregulation of carboxylesterase in high-fat diet mice after oral administration of Xuezhikang capsules

Xuezhikang capsule (XZK) is a traditional Chinese medicine that contains lovastatin (Lv) for hyperlipidemia treatment, although it has fewer side effects than Lv. However, the pharmacokinetic mechanisms contributing to its distinct efficacy and low side effects are unclear. Mice were fed a high-fat diet (HFD) for 6 weeks to induce hyperlipidemia. We first conducted the pharmacokinetic studies in HFD mice following oral administration of Lv (10 mg/kg, i.g.) and found that HFD remarkably decreased the active form of Lv (the lovastatin acid, LvA) exposure in the circulation system, especially in the targeting organ liver, with a declined conversion from Lv to LvA, whereas the Lv (responsible for myotoxicity) exposure in muscle markedly increased. Then we compared the pharmacokinetic profiles of Lv in HFD mice after the oral administration of XZK (1200 mg/kg, i.g.) or an equivalent dose of Lv (10 mg/kg, i.g.). A higher exposure of LvA and lower exposure of Lv were observed after XZK administration, suggesting a pharmacokinetic interaction of some ingredients in XZK. Further studies revealed that HFD promoted the inflammation and inhibited carboxylesterase (CES) activities in the intestine and the liver, thus contributing to the lower transformation of Lv into LvA. In contrast, XZK inhibited the inflammation and upregulated CES in the intestine and the liver. Finally, we evaluated the effects of monacolins and phytosterols, the fractional extracts of isoflavones, on inflammatory LS174T or HepG2 cells, which showed that isoflavones inhibited inflammation, upregulated CES, and markedly enhanced the conversion of Lv into LvA. For the first time, we provide evidence that isoflavones and Lv in XZK act in concert to enhance the efficacy and reduce the side effects of Lv.

Carboxylesterase 1 genes: systematic review and evaluation of existing genotyping procedures

The carboxylesterase 1 gene (CES1) encodes a hydrolase that metabolizes commonly used drugs. The CES1-related pseudogene, carboxylesterase 1 pseudogene 1 (CES1P1), has been implicated in gene exchange with CES1 and in the formation of hybrid genes including the carboxylesterase 1A2 gene (CES1A2). Hence, the CES1 region is complex. Using in silico PCR and alignment, we assessed the specificity of PCR-assisted procedures for genotyping CES1, CES1A2 and CES1P1 in studies identified in PubMed. We identified 33 such studies and excluded those that were not the first to use a procedure or lacked sequence information. After this 17 studies remained. Ten of these used haplotype-specific amplification, restriction enzyme treatment or amplicon sequencing, and included five that were predicted to lack specificity. All procedures for genotyping of single nucleotide polymorphisms in eight studies lacked specificity. One of these studies also used amplicon sequencing, thus being present in the group above. Some primers and their intended targets were mismatched. We provide experimental evidence that one of the procedures lacked specificity. Additionally, a complex pattern of segmental duplications in the CES1 region was revealed. In conclusion, many procedures for CES1, CES1A2 and CES1P1 genotyping appear to lack specificity. Knowledge about the segmental duplications may improve the typing of these genes.

Reappraisal of the genetic diversity and pharmacogenetic assessment of CES1

The CES1 gene encodes a hydrolase that metabolizes important drugs. Variants generated by exchange of segments with CES1P1 complicate genotyping of CES1. Using a highly specific procedure we examined DNA samples from 200 Caucasians and identified 46 single nucleotide variants (SNVs) in CES1 and 21 SNVs in CES1A2, a hybrid composed of CES1 and CES1P1. Several of these SNVs were novel. The frequencies of SNVs with a potential functional impact were below 0.02 suggesting limited pharmacogenetic potential for CES1 genotyping. In silico PCR revealed that the majority of the primer pairs for amplification of CES1 or CES1A2 in three previous studies lacked specificity, which partially explains a limited overlap with our findings.

Gene copy number variation analysis reveals dosage-insensitive expression of CYP2E1

Gene copy number variants (CNVs) of CYP2E1 have been described but not functionally characterized. Here we investigated effects of CNVs on hepatic and lymphoblastoid CYP2E1 expression. Using available single-nuleotide polymorphism microarray data and quantitative PCR, CYP2E1 gene duplication and deletion carriers were identified. CYP2E1 mRNA, protein and enzyme activity (chlorzoxazone-6-hydroxylation) phenotypes of CYP2E1 were not associated with gene copy number. Analysis of gene expression in lymphoblastoid cell lines in relation to CNV confirmed this finding in an extrahepatic tissue and for other ethnicities. Further analyses identified a linked haplotype cluster with possible influence on gene expression. In summary, our data suggest a homeostatic, gene dosage-insensitive regulation of CYP2E1 expression by unknown gene dosage compensation mechanisms. This is in striking contrast to well-known structural variations of CYP2A6 and CYP2D6 that have a strong impact on expression and activity. These findings are important in the context of pharmacogenetic prediction.

Genome-wide expression in visceral adipose tissue from obese prepubertal children

Characterization of the genes expressed in adipose tissue (AT) is key to understanding the pathogenesis of obesity and to developing treatments for this condition. Our objective was to compare the gene expression in visceral AT (VAT) between obese and normal-weight prepubertal children. A total of fifteen obese and sixteen normal-weight children undergoing abdominal elective surgery were selected. RNA was extracted from VAT biopsies. Microarray experiments were independently performed for each sample (six obese and five normal-weight samples). Validation by quantitative PCR (qPCR) was performed on an additional 10 obese and 10 normal-weight VAT samples. Of 1276 differentially expressed genes (p < 0.05), 245 were more than two-fold higher in obese children than in normal-weight children. As validated by qPCR, expression was upregulated in genes involved in lipid and amino acid metabolism (CES1, NPRR3 and BHMT2), oxidative stress and extracellular matrix regulation (TNMD and NQO1), adipogenesis (CRYAB and AFF1) and inflammation (ANXA1); by contrast, only CALCRL gene expression was confirmed to be downregulated. In conclusion, this study in prepubertal children demonstrates the up- and down-regulation of genes that encode molecules that were previously proposed to influence the pathogenesis of adulthood obesity, as well as previously unreported dysregulated genes that may be candidate genes in the aetiology of obesity.

Association between copy number variation losses and alcohol dependence across African American and European American ethnic groups

BACKGROUND

Copy number variations (CNVs) are structural genetic mutations consisting of segmental gains or losses in DNA sequence. Although CNVs contribute substantially to genomic variation, few genetic and imaging studies report association of CNVs with alcohol dependence (AD). Our purpose is to find evidence of this association across ethnic populations and genders. This work is the first AD-CNV study across ethnic groups and the first to include the African American (AA) population.

METHODS

This study considers 2 CNV data sets, one for discovery (2,345 samples) and the other for validation (239 samples), both including subjects with AD and healthy controls of European and African ancestry. Our analysis assesses the association between AD and CNV losses across ethnic groups and gender by examining the effect of overall losses across the whole genome, collective losses within individual cytogenetic bands, and specific losses in CNV regions.

RESULTS

Results from the discovery data set showed an association between CNV losses within 16q12.2 and AD diagnosis (p = 4.53 × 10(-3) ). An overlapping CNV region from the validation data set exhibited the same direction of effect with respect to AD (p = 0.051). This CNV region affects the genes CES1p1 and CES1, which are members of the carboxylesterase (CES) family. The enzyme encoded by CES1 is a major liver enzyme that typically catalyzes the decomposition of ester into alcohol and carboxylic acid and is involved in drug or xenobiotics, fatty acid, and cholesterol metabolisms. In addition, the most significantly associated CNV region was located at 9p21.2 (p = 1.9 × 10(-3) ) in our discovery data set. Although not observed in the validation data set, probably due to small sample size, this result might hold potential connection to AD given its connection with neuronal death. In contrast, we did not find any association between AD and the overall total losses or the collective losses within individual cytogenetic bands.

CONCLUSIONS

Overall, our study provides evidence that the specific CNVs at 16q12.2 contribute to the development of alcoholism in AA and European American populations.