REALGAR: a web app of integrated respiratory omics data.

MOTIVATION: In the post genome-wide association study (GWAS) era, omics techniques have characterized information beyond genomic variants to include cell and tissue type-specific gene transcription, transcription factor binding sites, expression quantitative trait loci (eQTL) and many other biological layers. Analysis of omics data and its integration has in turn improved the functional interpretation of disease-associated genetic variants. Over 170 000 transcriptomic and epigenomic datasets corresponding to studies of various cell and tissue types under specific disease, treatment and exposure conditions are available in the Gene Expression Omnibus resource. Although these datasets are valuable to guide the design of experimental validation studies to understand the function of disease-associated genetic loci, in their raw form, they are not helpful to experimental researchers who lack adequate computational resources or experience analyzing omics data. We sought to create an integrated re-source of tissue-specific results from omics studies that is guided by disease-specific knowledge to facilitate the design of experiments that can provide biologically meaningful insights into genetic associations. RESULTS: We designed the Reducing Associations by Linking Genes and omics Results web app to provide multi-layered omics information based on results from GWAS, transcriptomic, epigenomic and eQTL studies for gene-centric analysis and visualization. With a focus on asthma datasets, the integrated omics results it contains facilitate the formulation of hypotheses related to airways disease-associated genes and can be addressed with experimental validation studies. AVAILABILITY AND IMPLEMENTATION: The REALGAR web app is available at: http://realgar.org/. The source code is available at: https://github.com/HimesGroup/realgar. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multiomics analysis identifies BIRC3 as a novel glucocorticoid response-associated gene.

BACKGROUND: Inhaled corticosteroid (ICS) response among patients with asthma is influenced by genetics, but biologically actionable insights based on associations have not been found. Various glucocorticoid response omics data sets are available to interrogate their biological effects. OBJECTIVE: We sought to identify functionally relevant ICS-response genetic associations by integrating complementary multiomics data sets. METHODS: Variants with P values less than 10-4 from a previous ICS-response genome-wide association study were reranked on the basis of integrative scores determined from (1) glucocorticoid receptor- and (2) RNA polymerase II-binding regions inferred from ChIP-Seq data for 3 airway cell types, (3) glucocorticoid response element motifs, (4) differentially expressed genes in response to glucocorticoid exposure according to 20 transcriptomic data sets, and (5) expression quantitative trait loci from GTEx. Candidate variants were tested for association with ICS response and asthma in 6 independent studies. RESULTS: Four variants had significant (q value < 0.05) multiomics integrative scores. These variants were in a locus consisting of 52 variants in high linkage disequilibrium (r2 ≥ 0.8) near glucocorticoid receptor-binding sites by the gene BIRC3. Variants were also BIRC3 expression quantitative trait loci in lung, and 2 were within/near putative glucocorticoid response element motifs. BIRC3 had increased RNA polymerase II occupancy and gene expression, with glucocorticoid exposure in 2 ChIP-Seq and 13 transcriptomic data sets. Some BIRC3 variants in the 52-variant locus were associated (P < .05) with ICS response in 3 independent studies and others with asthma in 1 study. CONCLUSIONS: BIRC3 should be prioritized for further functional studies of ICS response.

Genetic Variation and Sensory Perception of a Pediatric Formulation of Ibuprofen: Can a Medicine Taste Too Good for Some?

There is wide variation in how individuals perceive the chemosensory attributes of liquid formulations of ibuprofen, encompassing both adults and children. To understand personal variation in the taste and chemesthesis properties of this medicine, and how to measure it, our first scientific strategy centered on utilizing trained adult panelists, due to the complex and time-consuming psychophysical tasks needed at this initial stage. We conducted a double-blind cohort study in which panelists underwent whole-genome-wide genotyping and psychophysically evaluated an over-the-counter pediatric medicine containing ibuprofen. Associations between sensory phenotypes and genetic variation near/within irritant and taste receptor genes were determined. Panelists who experienced the urge to cough or throat sensations found the medicine less palatable and sweet, and more irritating. Perceptions varied with genetic ancestry; panelists of African genetic ancestry had fewer chemesthetic sensations, rating the medicine sweeter, less irritating, and more palatable than did those of European genetic ancestry. We discovered a novel association between TRPA1 rs11988795 and tingling sensations, independent of ancestry. We also determined for the first time that just tasting the medicine allowed predictions of perceptions after swallowing, simplifying future psychophysical studies on diverse populations of different age groups needed to understand genetic, cultural-dietary, and epigenetic factors that influence individual perceptions of palatability and, in turn, adherence and the risk of accidental ingestion.

Inhibition of ABCC1 Decreases cAMP Egress and Promotes Human Airway Smooth Muscle Cell Relaxation.

In most living cells, the second-messenger roles for adenosine 3',5'-cyclic monophosphate (cAMP) are short-lived, confined to the intracellular space, and tightly controlled by the binary switch-like actions of Gαs (stimulatory G protein)-activated adenylyl cyclase (cAMP production) and cAMP-specific PDE (cAMP breakdown). Here, by using human airway smooth muscle (HASM) cells in culture as a model, we report that activation of the cell-surface ß2AR (ß2-adrenoceptor), a Gs-coupled GPCR (G protein-coupled receptor), evokes cAMP egress to the extracellular space. Increased extracellular cAMP levels ([cAMP]e) are long-lived in culture and are induced by receptor-dependent and receptor-independent mechanisms in such a way as to define a universal response class of increased intracellular cAMP levels ([cAMP]i). We find that HASM cells express multiple ATP-binding cassette (ABC) membrane transporters, with ABCC1 (ABC subfamily member C 1) being the most highly enriched transcript mapped to MRPs (multidrug resistance-associated proteins). We show that pharmacological inhibition or downregulation of ABCC1 with siRNA markedly reduces ß2AR-evoked cAMP release from HASM cells. Furthermore, inhibition of ABCC1 activity or expression decreases basal tone and increases ß-agonist-induced HASM cellular relaxation. These findings identify a previously unrecognized role for ABCC1 in the homeostatic regulation of [cAMP]i in HASM that may be conserved traits of the Gs-GPCRs (Gs-coupled family of GPCRs). Hence, the general features of this activation mechanism may uncover new disease-modifying targets in the treatment of airflow obstruction in asthma. Surprisingly, we find that serum cAMP levels are elevated in a small cohort of patients with asthma as compared with control subjects, which warrants further investigation.

Epinephrine evokes shortening of human airway smooth muscle cells following ß2 adrenergic receptor desensitization.

Epinephrine (EPI), an endogenous catecholamine involved in the body's fight-or-flight responses to stress, activates α1-adrenergic receptors (α1ARs) expressed on various organs to evoke a wide range of physiological functions, including vasoconstriction. In the smooth muscle of human bronchi, however, the functional role of EPI on α1ARs remains controversial. Classically, evidence suggests that EPI promotes bronchodilation by stimulating ß2-adrenergic receptors (ß2ARs). Conventionally, the selective ß2AR agonism of EPI was thought to be, in part, due to a predominance of ß2ARs and/or a sparse, or lack of α1AR activity in human airway smooth muscle (HASM) cells. Surprisingly, we find that HASM cells express a high abundance of ADRA1B (the α1AR subtype B) and identify a spontaneous "switch-like" activation of α1ARs that evokes intracellular calcium, myosin light chain phosphorylation, and HASM cell shortening. The switch-like responses, and related EPI-induced biochemical and mechanical signals, emerged upon pharmacological inhibition of ß2ARs and/or under experimental conditions that induce ß2AR tachyphylaxis. EPI-induced procontractile effects were abrogated by an α1AR antagonist, doxazosin mesylate (DM). These data collectively uncover a previously unrecognized feed-forward mechanism driving bronchospasm via two distinct classes of G protein-coupled receptors (GPCRs) and provide a basis for reexamining α1AR inhibition for the management of stress/exercise-induced asthma and/or ß2-agonist insensitivity in patients with difficult-to-control, disease subtypes.

CEBPD modulates the airway smooth muscle transcriptomic response to glucocorticoids.

BACKGROUND: CCAAT/Enhancer Binding Protein D (CEBPD), a pleiotropic glucocorticoid-responsive transcription factor, modulates inflammatory responses. Of relevance to asthma, expression of CEBPD in airway smooth muscle (ASM) increases with glucocorticoid exposure. We sought to characterize CEBPD-mediated transcriptomic responses to glucocorticoid exposure in ASM by measuring changes observed after knockdown of CEBPD and its impact on asthma-related ASM function. METHODS: Primary ASM cells derived from four donors were transfected with CEBPD or non-targeting (NT) siRNA and exposed to vehicle control, budesonide (100 nM, 18 h), TNFα (10 ng/ml, 18 h), or both budesonide and TNFα. Subsequently, RNA-Seq was used to measure gene expression levels, and pairwise differential expression results were obtained for exposures versus vehicle and knockdown versus control conditions. Weighted gene co-expression analysis was performed to identify groups of genes with similar expression patterns across the various experimental conditions (i.e., CEBPD knockdown status, exposures). RESULTS: CEBPD knockdown altered expression of 3037 genes under at least one exposure (q-value < 0.05). Co-expression analysis identified sets of 197, 152 and 290 genes that were correlated with CEBPD knockdown status, TNFα exposure status, and both, respectively. JAK-STAT signaling pathway genes, including IL6R and SOCS3, were among those influenced by both TNFα and CEBPD knockdown. Immunoblot assays revealed that budesonide-induced IL-6R protein expression and augmented IL-6-induced STAT3 phosphorylation levels were attenuated by CEBPD knockdown in ASM. CONCLUSIONS: CEBPD modulates glucocorticoid responses in ASM, in part via modulation of IL-6 receptor signaling.

Gene Expression Profiles in Children With Suspected Sepsis.

STUDY OBJECTIVE: Sepsis recognition is a clinical challenge in children. We aim to determine whether peripheral blood gene expression profiles are associated with pathogen type and sepsis severity in children with suspected sepsis. METHODS: This was a prospective pilot observational study in a tertiary pediatric emergency department with a convenience sample of children enrolled. Participants were older than 56 days and younger than 18 years, had suspected sepsis, and had not received broad-spectrum antibiotics in the previous 4 hours. Primary outcome was source pathogen, defined as confirmed bacterial source from sterile body fluid or confirmed viral source. Secondary outcome was sepsis severity, defined as maximum therapy required for shock reversal in the first 3 hospital days. We drew peripheral blood for ribonucleic acid isolation at the sepsis protocol activation, obtained gene expression measures with the GeneChip Human Gene 2.0 ST Array, and conducted differential expression analysis. RESULTS: We collected ribonucleic acid samples from a convenience sample of 122 children with suspected sepsis and 12 healthy controls. We compared the 66 children (54%) with confirmed bacterial or viral infection and found 558 differentially expressed genes, many related to interferon signaling or viral immunity. We did not find statistically significant gene expression differences in patients according to sepsis severity. CONCLUSION: The study demonstrates feasibility of evaluating gene expression profiling data in children evaluated for sepsis in the pediatric emergency department setting. Our results suggest that gene expression profiling may facilitate identification of source pathogen in children with suspected sepsis, which could ultimately lead to improved tailoring of sepsis treatment and antimicrobial stewardship.

Airway Smooth Muscle-Specific Transcriptomic Signatures of Glucocorticoid Exposure.

Glucocorticoids, commonly used asthma controller medications, decrease symptoms in most patients, but some remain symptomatic despite high-dose treatment. The physiological basis underlying the glucocorticoid response, especially in asthma patients with severe, refractory disease, is not fully understood. We sought to identify differences between the transcriptomic response of airway smooth muscle (ASM) cells derived from donors with fatal asthma and donors without asthma to glucocorticoid exposure and to compare ASM-specific changes with those observed in other cell types. In cells derived from nine donors with fatal asthma and eight donors without asthma, RNA sequencing was used to measure ASM transcriptome changes after exposure to budesonide (100 nM 24 h) or control vehicle (DMSO). Differential expression results were obtained for this dataset, as well as 13 publicly available glucocorticoid-response transcriptomic datasets corresponding to seven cell types. Specific genes were differentially expressed in response to glucocorticoid exposure (7,835 and 6,957 in ASM cells derived from donors with fatal asthma and donors without asthma, respectively; adjusted P value < 0.05). Transcriptomic changes in response to glucocorticoid exposure were similar in ASM derived from donors with fatal asthma and donors without asthma, with enriched ontological pathways that included cytokine- and chemokine-related categories. A comparison of glucocorticoid-induced changes in the nonasthma ASM transcriptome with those observed in six other cell types showed that ASM has a distinct glucocorticoid-response signature that is also present in ASM cells from donors with fatal asthma.

A functional splice variant associated with decreased asthma risk abolishes the ability of gasdermin B to induce epithelial cell pyroptosis.

BACKGROUND: Genetic variants in the chromosomal region 17q21 are consistently associated with asthma. However, mechanistic studies have not yet linked any of the associated variants to a function that could influence asthma, and as a result, the identity of the asthma gene(s) remains elusive. OBJECTIVES: We sought to identify and characterize functional variants in the 17q21 locus. METHODS: We used the Exome Aggregation Consortium browser to identify coding (amino acid-changing) variants in the 17q21 locus. We obtained asthma association measures for these variants in both the Genetic Epidemiology Research in Adult Health and Aging (GERA) cohort (16,274 cases and 38,269 matched controls) and the EVE Consortium study (5,303 asthma cases and 12,560 individuals). Gene expression and protein localization were determined by quantitative RT-PCR and fluorescence immunostaining, respectively. Molecular and cellular studies were performed to determine the functional effects of coding variants. RESULTS: Two coding variants (rs2305480 and rs11078928) of the gasdermin B (GSDMB) gene in the 17q21 locus were associated with lower asthma risk in both GERA (odds ratio, 0.92; P = 1.01 × 10-6) and EVE (odds ratio, 0.85; joint PEVE = 1.31 × 10-13). In GERA, rs11078928 had a minor allele frequency (MAF) of 0.45 in unaffected (nonasthmatic) controls and 0.43 in asthma cases. For European Americans in EVE, the MAF of rs2305480 was 0.45 for controls and 0.39 for cases; for all EVE subjects, the MAF was 0.32 for controls and 0.27 for cases. GSDMB is highly expressed in differentiated airway epithelial cells, including the ciliated cells. We found that, when the GSDMB protein is cleaved by inflammatory caspase-1 to release its N-terminal fragment, potent pyroptotic cell death is induced. The splice variant rs11078928 deletes the entire exon 6, which encodes 13 amino acids in the critical N-terminus, and abolishes the pyroptotic activity of the GSDMB protein. CONCLUSIONS: Our study identified a functional asthma variant in the GSDMB gene of the 17q21 locus and implicates GSDMB-mediated epithelial cell pyroptosis in pathogenesis.

Multiregion sequencing reveals the intratumor heterogeneity of driver mutations in TP53-driven non-small cell lung cancer.

Intratumor heterogeneity (ITH) in non-small cell lung cancer (NSCLC) may account for resistance after a period of targeted therapies because drugs destroy only a portion of tumor cells. The recognition of ITH helps identify high-risk patients to make effective treatment decisions. However, ITH studies are confounded by interpatient heterogeneity in NSCLC and a large amount of passenger mutations. To address these issues, we recruited NSCLC patients carrying TP53 mutations and selected driver mutations within recurrently mutated genes in NSCLC. A total of 12-paired normal-tumor tissues were subjected to whole-genome/whole-exome sequencing. From these, 367 non-silent mutations were selected as driver mutations and deeply sequenced in 61 intratumoral microdissections. We identified a universal prevalence of heterogeneity in all 12 tumors, indicating branched evolution. Although TP53 mutations were observed in single biopsy of all 12 tumors, most tumors consist of both TP53 mutated and non-mutated cells in separate regions within the same tumor. This suggests the late molecular timing of the acquisition of TP53 mutations; therefore, the detection of TP53 mutations in a single biopsy may simply not reflect the early malignant potential. In addition, we identified regions of loss of heterozygosity surrounding TP53 and CDKN2A mutations in tumor 711, which also exhibited heterogeneity in different regional samples. Because the ITH of driver mutations likely has clinical consequences, further efforts are needed to limit the impact of ITH and to improve therapeutic efficiency, which will benefit NSCLC patients receiving targeted treatments.

Rare-variant extensions of the transmission disequilibrium test: application to autism exome sequence data.

Many population-based rare-variant (RV) association tests, which aggregate variants across a region, have been developed to analyze sequence data. A drawback of analyzing population-based data is that it is difficult to adequately control for population substructure and admixture, and spurious associations can occur. For RVs, this problem can be substantial, because the spectrum of rare variation can differ greatly between populations. A solution is to analyze parent-child trio data, by using the transmission disequilibrium test (TDT), which is robust to population substructure and admixture. We extended the TDT to test for RV associations using four commonly used methods. We demonstrate that for all RV-TDT methods, using proper analysis strategies, type I error is well-controlled even when there are high levels of population substructure or admixture. For trio data, unlike for population-based data, RV allele-counting association methods will lead to inflated type I errors. However type I errors can be properly controlled by obtaining p values empirically through haplotype permutation. The power of the RV-TDT methods was evaluated and compared to the analysis of case-control data with a number of genetic and disease models. The RV-TDT was also used to analyze exome data from 199 Simons Simplex Collection autism trios and an association was observed with variants in ABCA7. Given the problem of adequately controlling for population substructure and admixture in RV association studies and the growing number of sequence-based trio studies, the RV-TDT is extremely beneficial to elucidate the involvement of RVs in the etiology of complex traits.

Using omics approaches to understand pulmonary diseases.

Omics approaches are high-throughput unbiased technologies that provide snapshots of various aspects of biological systems and include: 1) genomics, the measure of DNA variation; 2) transcriptomics, the measure of RNA expression; 3) epigenomics, the measure of DNA alterations not involving sequence variation that influence RNA expression; 4) proteomics, the measure of protein expression or its chemical modifications; and 5) metabolomics, the measure of metabolite levels. Our understanding of pulmonary diseases has increased as a result of applying these omics approaches to characterize patients, uncover mechanisms underlying drug responsiveness, and identify effects of environmental exposures and interventions. As more tissue- and cell-specific omics data is analyzed and integrated for diverse patients under various conditions, there will be increased identification of key mechanisms that underlie pulmonary biological processes, disease endotypes, and novel therapeutics that are efficacious in select individuals. We provide a synopsis of how omics approaches have advanced our understanding of asthma, chronic obstructive pulmonary disease (COPD), acute respiratory distress syndrome (ARDS), idiopathic pulmonary fibrosis (IPF), and pulmonary arterial hypertension (PAH), and we highlight ongoing work that will facilitate pulmonary disease precision medicine.

DNA methylation profiling in the thalamus and hippocampus of postnatal malnourished mice, including effects related to long-term potentiation.

BACKGROUND: DNA methylation has been viewed as the most highly characterized epigenetic mark for genome regulation and development. Postnatal brains appear to exhibit stimulus-induced methylation changes because of factors such as environment, lifestyle, and diet (nutrition). The purpose of this study was to examine how extensively the brain DNA methylome is regulated by nutrition in early life. RESULTS: By quantifying the total amount of 5-methylcytosine (5mC) in the thalamus and the hippocampus of postnatal malnourished mice and normal mice, we found the two regions showed differences in global DNA methylation status. The methylation level in the thalamus was much higher than that in the hippocampus. Then, we used a next-generation sequencing (NGS)-based method (MSCC) to detect the whole genome methylation of the two regions in malnourished mice and normal mice. Notably, we found that in the thalamus, 500 discriminable variations existed and that approximately 60% were related to neuronal development or psychiatric diseases. Pathway analyses of the corresponding genes highlighted changes for 9 genes related to long-term potentiation (5.3-fold enrichment, P = 0.033). CONCLUSIONS: Our findings may help to indicate the genome-wide DNA methylation status of different brain regions and the effects of malnutrition on brain DNA methylation. The results also indicate that postnatal malnutrition may increase the risk of psychiatric disorders.

GPC5 rs2352028 polymorphism and risk of lung cancer in Han Chinese.

rs2352028 in GPC5 has been reported to be associated with the risk of lung cancer in never-smokers. We performed a replication study in 1,045 lung cancer patients and 1,094 controls of Han Chinese origin. We found no association between rs2352028 and lung cancer/adenocarcinoma in never-smokers, but a p value of .04 (under the recessive model) was obtained between this SNP and overall lung cancer/adenocarcinoma. Our data and a recent meta-analysis suspected the possibility of rs2352028 being a risk variant of lung cancer risk in never-smokers. Our findings suggested that rs2352028 might confer a slight risk to lung cancer/adenocarcinoma.

A Common variant near the melanocortin 4 receptor is associated with low-density lipoprotein cholesterol and total cholesterol in the Chinese Han population.

A number of recent genome-wide association studies in European populations have reported that variant rs17782313 is significantly associated with obesity and body mass index (BMI). The purpose of the present study was to evaluate the association of rs17782313 with obesity and BMI in the Chinese Han population. We also sought to extend previous studies by determining whether this SNP is associated with plasma lipid levels in the Chinese Han population. Rs17782313 was genotyped in two independent Chinese Han cohorts (Cohort1: n = 2533; Cohort2: n = 2105). In our study, rs17782313 did not show significant association with either obesity or BMI in the Chinese Han population, but showed evidence for association with LDL-C (P ~ 0.003) and TC (P ~ 0.001). Our findings indicate that the variant rs17782313 near MC4R is likely to have an impact on plasma lipid levels of LDL-C and TC in the Chinese Han population.

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness.

Variability in response to short-acting ß2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.

Insights into glucocorticoid responses derived from omics studies.

Glucocorticoid drugs are commonly used in the treatment of several conditions, including autoimmune diseases, asthma and cancer. Despite their widespread use and knowledge of biological pathways via which they act, much remains to be learned about the cell type-specific mechanisms of glucocorticoid action and the reasons why patients respond differently to them. In recent years, human and in vitro studies have addressed these questions with genomics, transcriptomics and other omics approaches. Here, we summarize key insights derived from omics studies of glucocorticoid response, and we identify existing knowledge gaps related to mechanisms of glucocorticoid action that future studies can address.

Variations in/nearby genes coding for JAZF1, TSPAN8/LGR5 and HHEX-IDE and risk of type 2 diabetes in Han Chinese.

Several genetic loci (JAZF1, CDC123/CAMK1D, TSPAN8/LGR5, ADAMTS9, VEGFA and HHEX-IDE) were identified to be significantly related to the risk of type 2 diabetes and quantitative metabolic traits in European populations. Here, we aimed to evaluate the impacts of these novel loci on type 2 diabetes risk in a population-based case-control study of Han Chinese (1912 cases and 2041 controls). We genotyped 13 single-nucleotide polymorphisms (SNPs) in/near these genes and examined the differences in allele/genotype frequency between cases and controls. We found that both IDE rs11187007 and HHEX rs1111875 were associated with type 2 diabetes risk (for both variants: odds ratio (OR)=1.15, 95% confidence interval (CI) 1.04-1.28, P=0.009). In a meta-analysis where we pooled our data with the three previous studies conducted in East Asians, we found that the variants of JAZF1 rs864745 (1.09 (1.03-1.16); P=3.49 × 10(-3)) and TSPAN8/LGR5 rs7961581 (1.11(1.05-1.17); P=1.89 × 10(-4)) were significantly associated with type 2 diabetes risk. In addition, the meta-analysis (7207 cases and 8260 controls) also showed that HHEX rs1111875 did have effects on type 2 diabetes in Chinese population (OR=1.15(1.10-1.21); P=1.93 × 10(-8)). This large population-based study and meta-analysis further confirmed the modest effects of the JAZF1, TSPAN8/LGR5 and HHEX-IDE loci on type 2 diabetes in Chinese and other East Asians.

Chaperone-mediated autophagy regulates the pluripotency of embryonic stem cells.

Embryonic stem cells can propagate indefinitely in a pluripotent state, able to differentiate into all types of specialized cells when restored to the embryo. What sustains their pluripotency during propagation remains unclear. Here, we show that core pluripotency factors OCT4 and SOX2 suppress chaperone-mediated autophagy (CMA), a selective form of autophagy, until the initiation of differentiation. Low CMA activity promotes embryonic stem cell self-renewal, whereas its up-regulation enhances differentiation. CMA degrades isocitrate dehydrogenases IDH1 and IDH2 and reduces levels of intracellular α-ketoglutarate, an obligatory cofactor for various histone and DNA demethylases involved in pluripotency. These findings suggest that CMA mediates the effect of core pluripotency factors on metabolism, shaping the epigenetic landscape of stem cells and governing the balance between self-renewal and differentiation.