👤 Simon C Ling

The vertebrate genome is spatially organized into topologically associating domains (TADs), primarily via cohesin-mediated loop extrusion which typically halts at convergent CTCF binding sites to establish domain boundaries. However, despite the essential roles of CTCF and cohesin in establishing TADs, a long-standing paradox persists: CTCF and cohesin binding sites dramatically outnumber observed TAD boundaries, suggesting the existence of undiscovered architectural factors. To identify such missing factors, we conducted high-resolution Show less

Activation of mitochondrial function and heat production in adipose tissue by the modification of dietary fat is a promising strategy against obesity. However, as an important source of lipids for ketogenic and daily diets, the function of fats extracted from different adipose tissue sites was largely unknown. In this study, we illustrated the function of fats extracted from adipose tissues with different "beigeing" properties in the ketogenic diet and identified lipid profiles of fats that facilitate energy expenditure. We found that the anti-obesity effect of ketogenic diets was potentiated by using "beigeing" fat [porcine subcutaneous adipose tissue (SAT)] as a major energy-providing ingredient. Through lipidomic analyses, phosphatidylserine (PS) was identified as a functional lipid activating thermogenesis in adipose tissue. Moreover, in vivo studies showed that PS induces adipose tissue thermogenesis and alleviates diet-induced obesity in mice. In vitro studies showed that PS promotes UCP1 expression and lipolysis of adipocytes. Mechanistically, PS promoted mitochondrial function in adipocytes via the ADCY3-cAMP-PKA-PGC1α pathway. In addition, PS-PGC1a binding may affect the stability of the PGC1α protein, which further augments PS-induced thermogenesis. These results demonstrated the efficacy of dietary SAT fats in diminishing lipid accumulation and the underlying molecular mechanism of PS in enhancing UCP1 expression and mitochondrial function. Thus, our findings suggest that as dietary fat, "beigeing" fat provides more beneficial lipids that contribute to the improvement of mitochondrial function, including PS, which may become a novel, nonpharmacological therapy to increase energy expenditure and counteract obesity and its related diseases. Show less

To explore the patterns of differentially expressed genes (DEGs) associated with different growth rates in rock carp (Procypris rabaudi), transcriptome sequencing was performed on the muscle, liver, and brain tissues of rock carp. Subsequently, bioinformatics analysis was conducted, and 2129, 1380, and 415 DEGs were identified in the muscle, liver, and brain tissues, respectively. GO enrichment and KEGG pathway analysis revealed that genes related to appetite regulation, protein degradation and digestion, lipid transport and metabolisms, and glycolysis/gluconeogenesis were upregulated in individuals with slower growth rates. Differential expression analysis identified 21 genes associated with feeding and metabolism across three tissues, including mc4r, npy, and npry in brain tissue; fatp, fabp, pparα, and apo in liver tissue; and prss, ctrl, and cela in muscle tissue. All these genes were upregulated in the slow-growing fish. Furthermore, weighted gene co-expression network analyses, including three modules (yellow, turquoise, and brown), significantly associated with growth. A network map that included these three modules enabled the identification of a series of hub genes, including rp13a, ube2o, h6pd, etc. These genes may be key candidate genes regulating the growth of rock carp. This study contributes to a deeper understanding of the growth control mechanism in rock carp and offers a scientific basis for efficient breeding and species improvement. Show less

It has been demonstrated that angiopoietin-like protein 4 (ANGPTL4) plays an important regulatory role in lipid metabolism and backfat deposition appears to vary in different pig breeds. However, the correlation between ANGPTL4 and backfat deposition have not been well characterized and the role of ANGPTL4 in regulating adipogenesis remains unclear. Therefore, this study aimed to investigate correlation between ANGPTL4 and backfat deposition and to explore the effects of ANGPTL4 on preadipocyte differentiation and the underlying mechanism. Our results showed that the backfat thickness and the ANGPTL4 gene expression of Laiwu pigs were significantly higher than those in DLY pigs and the ANGPTL4 gene expression was positively correlated with backfat thickness both in DLY pigs and Laiwu pigs. Moreover, an increase in ANGPTL4 expression and activation of autophagy were observed during the differentiation of stromal vascular fraction cells. In addition, knockdown of ANGPTL4 inhibited the differentiation of 3T3-L1 cells with decreased expression of LC3-II and ATG5 and increased expression of SQSTM1, suggesting the involvement of autophagy in ANGPTL4-mediated adipogenesis. In conclusion, these results suggested that ANGPTL4 is positively correlated with backfat deposition in pigs and knockdown of ANGPTL4 inhibits adipogenesis of preadipocyte via autophagy, providing new insights into the regulation of fat deposition and to improve the carcass quality and meat quality of porcine. Show less

Despite centuries of traditional use of silymarin for hepatoprotection, current randomized controlled trial (RCT) studies on the effectiveness of silymarin in managing metabolic dysfunction-associated steatotic liver disease (MASLD) are limited and inconclusive, particularly when it is administered alone. The low bioavailability of silymarin highlights the possible influence of gut microbiota on the effectiveness of silymarin; however, no human studies have investigated this aspect. To determine the potential efficacy of silymarin in improving MASLD indicators and to investigate the underlying mechanisms related to gut microbiota. In this 24-week randomized, double-blind, placebo-controlled trial, 83 patients with MASLD were randomized to either placebo (n = 41) or silymarin (103.2 mg/d, n = 42). At 0, 12, and 24 weeks, liver stiffness and hepatic steatosis were assessed using FibroScan, and blood samples were gathered for biochemical detection, while faecal samples were collected at 0 and 24 weeks for 16S rRNA sequencing. Silymarin supplementation significantly reduced liver stiffness (LSM, -0.21 ± 0.17 vs. 0.41 ± 0.17, P = 0.015) and serum levels of γ-glutamyl transpeptidase (GGT, -8.21 ± 3.01 vs. 1.23 ± 3.16, P = 0.042) and ApoB (-0.02 ± 0.03 vs. 0.07 ± 0.03, P = 0.023) but had no significant effect on the controlled attenuation parameter (CAP), other biochemical indicators (aminotransferases, total bilirubin, glucose and lipid parameters, hsCRP, SOD, and UA), physical measurements (DBP, SBP, BMI, WHR, BF%, and BMR), or APRI and FIB-4 indices. Gut microbiota analysis revealed increased species diversity and enrichment of Oscillospiraceae in the silymarin group. These findings suggest that silymarin supplementation could improve liver stiffness in MASLD patients, possibly by modulating the gut microbiota. The trial was registered at the Chinese Clinical Trial Registry (ChiCTR2200059043). Show less

Primary aldosteronism (PA), a significant cause of secondary hypertension affecting ∼10% of patients with severe hypertension, exacerbates cardiovascular, and cerebrovascular complications even after blood pressure control. PA is categorized into two main subtypes: unilateral aldosterone-producing adenomas (APA) and bilateral hyperaldosteronism (BHA), each requiring distinct treatment approaches. Accurate subtype classification is crucial for selecting the most effective treatment. The goal of this study was to develop novel blood-based proteomic biomarkers to differentiate between APA and BHA subtypes in patients with PA. Five subtyping differential protein biomarker candidates (APOC3, CD56, CHGA, KRT5, and AZGP1) were identified through targeted proteomic profiling of plasma. The subtyping efficiency of these biomarkers was assessed at both the tissue gene expression and blood protein expression levels. To explore the underlying biology of APA and BHA, significant differential pathways were investigated. The five-protein panel proved highly effective in distinguishing APA from BHA in both tissue and blood samples. By integrating these five protein biomarkers with aldosterone and renin, our blood-based predictive methods achieved remarkable receiver operating characteristic (ROC) area under the ROC curves of 0.986 (95% CI: 0.963-1.000) for differentiating essential hypertension from PA, and 0.922 (95% CI: 0.846-0.998) for subtyping APA versus BHA. These outcomes surpass the performance of the existing Kobayashi score subtyping system. Furthermore, the study validated differential pathways associated with the pathophysiology of PA, aligning with current scientific knowledge and opening new avenues for advancing PA care. The new blood-based biomarkers for PA subtyping hold the potential to significantly enhance clinical utility and advance the practice of PA care. Show less

Alzheimer's disease (AD) is a progressive neurodegenerative disease with the hallmark of aggregation of beta-amyloid (Aβ) into extracellular fibrillar deposition. Accumulating evidence suggests that soluble toxic Aβ oligomers exert diverse roles in neuronal cell death, oxidative stress, neuroinflammation, and the eventual pathogenesis of AD. Aβ is derived from the sequential cleavage of amyloid-β precursor protein (APP) by β-secretase (BACE1) and γ-secretase. The current effect of single targeting is not ideal for the treatment of AD. Therefore, developing multipotent agents with multiple properties, including anti-Aβ generation and anti-Aβ aggregation, is attracting more attention for AD treatment. Previous studies indicated that Quercetin was able to attenuate the effects of several pathogenetic factors in AD. Here, we showed that naturally synthesized Quercetin-3-O-glc-1-3-rham-1-6-glucoside (YCC31) could inhibit Aβ production by reducing β-secretase activity. Further investigations indicated that YCC31 could suppress toxic Aβ oligomer formation by directly binding to Aβ. Moreover, YCC31 could attenuate Aβ-mediated neuronal death, ROS and NO production, and pro-inflammatory cytokines release. Taken together, YCC31 targeting multiple pathogenetic factors deserves further investigation for drug development of AD. Show less

The coronavirus disease 2019 (COVID-19) pandemic has caused extensive disruption of public health worldwide. There were reports of COVID-19 patients having multiple complications. This study investigated COVID-19 from a genetic perspective. We conducted RNA sequencing (RNA-Seq) analysis of respiratory tract samples from 24 patients with COVID-19. Eight patients receiving mechanical ventilation or extracorporeal membrane oxygenation were regarded as severe cases; the remaining 16 patients were regarded as non-severe cases. After quality control, statistical analyses were performed by logistic regression and the Kolmogorov-Smirnov test to identify genes associated with disease severity. Six genes were associated with COVID-19 severity in both statistical tests, namely RNA sequencing analysis showed that severe acute respiratory syndrome coronavirus 2 infection is associated with the overexpression of genes involved in nervous system disorders. Show less

Paroxysmal nocturnal haemoglobinuria (PNH) is a disorder resulting from erythrocyte membrane deficiencies caused by PIG-A gene mutations. While current treatments alleviate symptoms, they fail to address the underlying cause of the disease-the pathogenic PNH clones. In this study, we found that the expression of carbamoyl phosphate synthetase 1 (CPS1) was downregulated in PNH clones, and the level of CPS1 was negatively correlated with the proportion of PNH clones. Using PIG-A knockout K562 (K562 KO) cells, we demonstrated that CPS1 knockdown increased cell proliferation and altered cell metabolism, suggesting that CPS1 participates in PNH clonal proliferation through metabolic reprogramming. Furthermore, we observed an increase in the expression levels of the histone demethylase JMJD1C in PNH clones, and JMJD1C expression was negatively correlated with CPS1 expression. Knocking down JMJD1C in K562 KO cells upregulated CPS1 and H3K36me3 expression, decreased cell proliferation and increased cell apoptosis. Chromatin immunoprecipitation analysis further demonstrated that H3K36me3 regulated CPS1 expression. Finally, we demonstrated that histone demethylase inhibitor JIB-04 can suppressed K562 KO cell proliferation and reduced the proportion of PNH clones in PNH mice. In conclusion, aberrant regulation of the JMJD1C-H3K36me3-CPS1 axis contributes to PNH clonal proliferation. Targeting JMJD1C with a specific inhibitor unveils a potential strategy for treating PNH patients. Show less

It is not known why severe cystic fibrosis (CF) liver disease (CFLD) with portal hypertension occurs in only ~7% of people with CF. We aimed to identify genetic modifiers for severe CFLD to improve understanding of disease mechanisms. Whole-genome sequencing was available in 4082 people with CF with pancreatic insufficiency (n = 516 with severe CFLD; n = 3566 without CFLD). We tested ~15.9 million single nucleotide polymorphisms (SNPs) for association with severe CFLD versus no-CFLD, using pre-modulator clinical phenotypes including (1) genetic variant ( SERPINA1 ; Z allele) previously associated with severe CFLD; (2) candidate SNPs (n = 205) associated with non-CF liver diseases; (3) genome-wide association study of common/rare SNPs; (4) transcriptome-wide association; and (5) gene-level and pathway analyses. The Z allele was significantly associated with severe CFLD ( p = 1.1 × 10 -4 ). No significant candidate SNPs were identified. A genome-wide association study identified genome-wide significant SNPs in 2 loci and 2 suggestive loci. These 4 loci contained genes [significant, PKD1 ( p = 8.05 × 10 -10 ) and FNBP1 ( p = 4.74 × 10 -9 ); suggestive, DUSP6 ( p = 1.51 × 10 -7 ) and ANKUB1 ( p = 4.69 × 10 -7 )] relevant to severe CFLD pathophysiology. The transcriptome-wide association identified 3 genes [ CXCR1 ( p = 1.01 × 10 -6 ) , AAMP ( p = 1.07 × 10 -6 ), and TRBV24 ( p = 1.23 × 10 -5 )] involved in hepatic inflammation and innate immunity. Gene-ranked analyses identified pathways enriched in genes linked to multiple liver pathologies. These results identify loci/genes associated with severe CFLD that point to disease mechanisms involving hepatic fibrosis, inflammation, innate immune function, vascular pathology, intracellular signaling, actin cytoskeleton and tight junction integrity and mechanisms of hepatic steatosis and insulin resistance. These discoveries will facilitate mechanistic studies and the development of therapeutics for severe CFLD. Show less

Dendrobium nobile Lindl. (DNL) is a traditional Chinese ethnobotanical herb. Dendrobine (DNE) has been designated as a quality indicator for DNL in the Chinese Pharmacopoeia. DNE exhibits various pharmacological activities, including the reduction of blood lipids, regulation of blood sugar levels, as well as anti-inflammatory and antioxidant properties. The objective of this study is to explore the impact of DNE on lipid degeneration in nonalcoholic fatty liver disease (NAFLD) liver cells and elucidate its specific mechanism. The findings aim to offer theoretical support for the development of drugs related to DNL. We utilized male C57BL/6J mice, aged 6 weeks old, to establish a NAFLD model. This model allowed us to assess the impact of DNE on liver pathology and lipid levels in NAFLD mice. We investigated the mechanism of DNE's regulation of lipid metabolism through RNA-seq analysis. Furthermore, a NAFLD model was established using HepG2 cells to further evaluate the impact of DNE on the pathological changes of NAFLD liver cells. The potential mechanism of DNE's improvement was rapidly elucidated using HT-qPCR technology. These results were subsequently validated using mouse liver samples. Following the in vitro activation or inhibition of PPARα function, we observed changes in DNE's ability to ameliorate pathological changes in NAFLD hepatocytes. This mechanism was further verified through RT-qPCR and Western blot analysis. DNE demonstrated a capacity to enhance serum TC, TG, and liver TG levels in mice, concurrently mitigating liver lipid degeneration. RNA-seq analysis unveiled that DNE primarily modulates the expression of genes related to metabolic pathways in mouse liver. Utilizing HT-qPCR technology, it was observed that DNE markedly regulates the expression of genes associated with the PPAR signaling pathway in liver cells. Consistency was observed in the in vivo data, where DNE significantly up-regulated the expression of PPARα mRNA and its protein level in mouse liver. Additionally, the expression of fatty acid metabolism-related genes (ACOX1, CPT2, HMGCS2, LPL), regulated by PPARα, was significantly elevated following DNE treatment. In vitro experiments further demonstrated that DNE notably ameliorated lipid deposition, peroxidation, and inflammation levels in NAFLD hepatocytes, particularly when administered in conjunction with fenofibrate. Notably, the PPARα inhibitor GW6471 attenuated these effects of DNE. In summary, DNE exerts its influence on the expression of genes associated with downstream fat metabolism by regulating PPARα. This regulatory mechanism enhances liver lipid metabolism, mitigates lipid degeneration in hepatocytes, and ultimately ameliorates the pathological changes in NAFLD hepatocytes. Show less

Human gastrointestinal tumours have been shown to contain massive numbers of tumour infiltrating regulatory T cells (Tregs), the presence of which are closely related to tumour immunity. This study was designed to develop new Treg-related prognostic biomarkers to monitor the prognosis of patients with gastric cancer (GC). Treg-related prognostic genes were screened from Treg-related differentially expressed genes in GC patients by using Cox regression analysis, based on which a prognostic model was constructed. Then, combined with RiskScore, survival curve, survival status assessment and ROC analysis, these genes were used to verify the accuracy of the model, whose independent prognostic ability was also evaluated. Six Treg-related prognostic genes (CHRDL1, APOC3, NPTX1, TREML4, MCEMP1, GH2) in GC were identified, and a 6-gene Treg-related prognostic model was constructed. Survival analysis revealed that patients had a higher survival rate in the low-risk group. Combining clinicopathological features, we performed univariate and multivariate regression analyses, with results establishing that the RiskScore was an independent prognostic factor. Predicted 1-, 3- and 5-year survival rates of GC patients had a good fit with the actual survival rates according to nomogram results. In addition patients in the low-risk group had higher tumour mutational burden (TMB) values. Gene Set Enrichment Analysis (GSEA) demonstrated that genes in the high-risk group were significantly enriched in pathways related to immune inflammation, tumour proliferation and migration. In general, we constructed a 6-gene Treg-associated GC prognostic model with good prediction accuracy, where RiskScore could act as an independent prognostic factor. This model is expected to provide a reference for clinicians to estimate the prognosis of GC patients. Show less

Abdominal fat deposition depends on both the proliferation of preadipocytes and their maturation into adipocytes, which is a well-orchestrated multistep process involving many regulatory molecules. Circular RNAs (circRNAs) have emergingly been implicated in mammalian adipogenesis. However, circRNA-mediated regulation in chicken adipogenesis remains unclear. Our previous circRNA sequencing data identified a differentially expressed novel circRNA, 8:27,886,180|27,889,657, during the adipogenic differentiation of chicken abdominal preadipocytes. This study aimed to investigate the regulatory role of circDOCK7 in the proliferation and adipogenic differentiation of chicken abdominal preadipocytes, and explore its molecular mechanisms of competing endogenous RNA underlying chicken adipogenesis. Our results showed that 8:27,886,180|27,889,657 is an exonic circRNA derived from the head-to-tail splicing of exons 19-22 of the dedicator of cytokinesis 7 (DOCK7) gene, abbreviated as circDOCK7. CircDOCK7 is mainly distributed in the cytoplasm of chicken abdominal preadipocytes and is stable because of its RNase R resistance and longer half-life. CircDOCK7 is significantly upregulated in the abdominal fat tissues of fat chickens compared to lean chickens, and its expression gradually increases during the proliferation and adipogenic differentiation of chicken abdominal preadipocytes. Functionally, the gain- and loss-of-function experiments showed that circDOCK7 promoted proliferation, G0/G1- to S-phase progression, and glucose uptake capacity of chicken abdominal preadipocytes, in parallel with adipogenic differentiation characterized by remarkably increased intracellular lipid droplet accumulation and triglyceride and acetyl coenzyme A content in differentiated chicken abdominal preadipocytes. Mechanistically, a pull-down assay and a dual-luciferase reporter assay confirmed that circDOCK7 interacted with gga-miR-301b-3p, which was identified as an inhibitor of chicken abdominal adipogenesis. Moreover, the ACSL1 gene was demonstrated to be a direct target of gga-miR-301b-3p. Chicken ACSL1 protein is localized in the endoplasmic reticulum and mitochondria of chicken abdominal preadipocytes and acts as an adipogenesis accelerator. Rescue experiments showed that circDOCK7 could counteract the inhibitory effects of gga-miR-301b-3p on ACSL1 mRNA abundance as well as the proliferation and adipogenic differentiation of chicken abdominal preadipocytes. CircDOCK7 serves as a miRNA sponge that directly sequesters gga-miR-301b-3p away from the ACSL1 gene, thus augmenting adipogenesis in chickens. These findings may elucidate a new regulatory mechanism underlying abdominal fat deposition in chickens. Show less

Dysregulation of RNA-binding proteins (RBP) is one of the characteristics of cancer. Investigating the biological functions and molecular mechanisms of abnormal RBPs can help uncover new cancer biomarkers and treatment strategies. To identify oncogenic RBPs in triple-negative breast cancer (TNBC), we employed an in vivo CRISPR screen and a TNBC progression model, which revealed small nuclear ribonucleoprotein polypeptide C (SNRPC), a subunit of the U1 small nuclear ribonucleoprotein particle (U1 snRNP), as a key modulator of TNBC progression. SNRPC was frequently upregulated, which corresponded to poor prognosis in patients with TNBC. SNRPC ablation significantly impaired the proliferation, migration, and invasion of TNBC cells in vitro and in vivo. In addition, SNRPC was essential for the stability of U1 snRNP and contributed to the RNA Pol II-controlled transcriptional program. Knockdown of SNRPC decreased RNA Pol II enrichment on a subset of oncogenes (TNFAIP2, E2F2, and CDK4) and reduced their expression levels. Furthermore, SNRPC deletion was confirmed to inhibit TNBC progression partially through regulation of the TNFAIP2-Rac1-β-catenin signaling pathway. Taken together, this data suggests that SNRPC plays an oncogenic role in TNBC, is a marker of poor prognosis, and may be a valuable therapeutic target for patients with intractable TNBC. A functional CRISPR screen identifies SNRPC as an RNA-binding protein that promotes the aggressiveness of breast cancer by facilitating Pol II-controlled transcription of oncogenes. Show less

Multiple genome-wide association studies (GWAS) have identified SNPs in the 8q24 locus near TRIB1 that are significantly associated with plasma lipids and other markers of cardiometabolic health, and prior studies have revealed the roles of hepatic and myeloid Trib1 in plasma lipid regulation and atherosclerosis. The same 8q24 SNPs are additionally associated with plasma adiponectin levels in humans, implicating TRIB1 in adipocyte biology. Here, we hypothesize that TRIB1 in adipose tissue regulates plasma adiponectin, lipids, and metabolic health. We investigate the metabolic phenotype of adipocyte-specific Trib1 knockout mice (Trib1_ASKO) fed on chow and high-fat diet (HFD). Through secretomics of adipose tissue explants and RNA-seq of adipocytes and livers from these mice, we further investigate the mechanism of TRIB1 in adipose tissue. Trib1_ASKO mice have an improved metabolic phenotype with increased plasma adiponectin levels, improved glucose tolerance, and decreased plasma lipids. Trib1_ASKO adipocytes have increased adiponectin production and secretion independent of the known TRIB1 function of regulating proteasomal degradation. RNA-seq analysis of adipocytes and livers from Trib1_ASKO mice indicates that alterations in adipocyte function underlie the observed plasma lipid changes. Adipose tissue explant secretomics further reveals that Trib1_ASKO adipose tissue has decreased ANGPTL4 production, and we demonstrate an accompanying increase in the lipoprotein lipase (LPL) activity that likely underlies the triglyceride phenotype. This study shows that adipocyte Trib1 regulates multiple aspects of metabolic health, confirming previously observed genetic associations in humans and shedding light on the further mechanisms by which TRIB1 regulates plasma lipids and metabolic health. Show less

Hepatocellular carcinoma (HCC) is a globally prevailing cancer with a low 5-year survival rate. Little is known about its intricate gene expression profile. Single-cell RNA sequencing is an indispensable tool to explore the genetic characteristics of HCC at a more detailed level. In this study, we profiled the gene expression of single cells from human HCC tumor and para-tumor tissues using the Smart-seq 2 sequencing method. Based on differentially expressed genes, we identified heterogeneous subclones in HCC tissues, including five HCC and two hepatocyte subclones. We then carried out hub-gene co-network and functional annotations analysis followed pseudo-time analysis with regulated transcriptional factor co-networks to determine HCC cellular trajectory. We found that MLX interacting protein like (MLXIPL) was commonly upregulated in the single cells and tissues and associated with a poor survival rate in HCC. Mechanistically, MLXIPL activation is crucial for promoting cell proliferation and inhibits cell apoptosis by accelerating cell glycolysis. Taken together, our work identifies the heterogeneity of HCC subclones, and suggests MLXIPL might be a promising therapeutic target for HCC. Show less

PFKFB3 regulates glycolysis in tumor cells, might function as an oncogene, and is associated with cancer metastasis. However, its role in gastric cancer (GC) remains largely unknown. PFKFB3 expression was assessed by immunohistochemistry (IHC) in GC tissues and paired paracancerous histological normal tissues (PCHNTs). The associations of PFKFB3 expression with clinical features and HIF-1α, Ki-67, E-cadherin, Snail, and Vimentin expression levels were assessed. A series of in vivo and in vitro experiments were performed to investigate the effects of PFKFB3 on the growth, migration, and invasion of GC cells. We found that PFKFB3 expression was significantly higher in GC tissues compared with PCHNTs (P = 0.000). PFKFB3 expression was positively correlated with tumor size (P = 0.000), differentiation (P = 0.025), venous invasion (P = 0.084), nerve invasion (P = 0.014), lymphatic invasion (P = 0.000), local invasion (P = 0.000), invasive depth (P = 0.000), nodal metastasis (P = 0.000), tumor-node-metastasis stage (P = 0.000), and patient survival (P = 0.000). Notably, PFKFB3 upregulation was highly correlated with increased epithelial-mesenchymal transition (EMT) in GC samples. PFKFB3 overexpression positively modulated cell proliferation, migration, and EMT in GC cells in vitro, with concomitant activation of NF-κB signaling. Administration of an NF-κB inhibitor attenuated PFKFB3-induced EMT in GC cells. PFKFB3 overexpression promoted tumor development and EMT in nude mice, which were attenuated by PFK-15, a PFKFB3 inhibitor. PFKFB3 could potentiate malignancy in GC cells through NF-κB pathway-mediated EMT, suggesting PFKFB3 represents a potential target for GC therapy. Show less

Although the immune checkpoint role of programmed death ligand 1 (PD-L1) has been established and targeted in cancer immunotherapy, the tumor-intrinsic role of PD-L1 is less appreciated in tumor biology and therapeutics development, partly because of the incomplete mechanistic understanding. Here we demonstrate a potentially novel mechanism by which PD-L1 promotes the epithelial-mesenchymal transition (EMT) in triple-negative breast cancer (TNBC) cells by suppressing the destruction of the EMT transcription factor Snail. PD-L1 directly binds to and inhibits the tyrosine phosphatase PTP1B, thus preserving p38-MAPK activity that phosphorylates and inhibits glycogen synthase kinase 3β (GSK3β). Via this mechanism, PD-L1 prevents the GSK3β-mediated phosphorylation, ubiquitination, and degradation of Snail and consequently promotes the EMT and metastatic potential of TNBC. Significantly, PD-L1 antibodies that confine the tumor-intrinsic PD-L1/Snail pathway restricted TNBC progression in immunodeficient mice. More importantly, targeting both tumor-intrinsic and tumor-extrinsic functions of PD-L1 showed strong synergistic tumor suppression effect in an immunocompetent TNBC mouse model. Our findings support that PD-L1 intrinsically facilitates TNBC progression by promoting the EMT, and this potentially novel PD-L1 signaling pathway could be targeted for better clinical management of PD-L1-overexpressing TNBCs. Show less

SNAI1 is an epithelial-mesenchymal transition (EMT) inducer, which has been indicated to play a role in the progression of cancers. We aimed to evaluate the expression and prognostic roles of SNAI1 in head and neck carcinoma (HNC). The study involved two major phases. In the in silico phase, the SNAI1 expression and its association with clinical features as well as its prognostic values were assessed; then, the target genes of SNAI1 were predicted and the relationship between SNAI1 expression and immune cell infiltration was evaluated. In the validation phase, a cohort of a tissue microarray (47 cases) and a cohort of HNC patients (68 cases) were enrolled. SNAI1 was detected by using an immunochemistry assay. The in silico analysis showed that overexpression of SNAI1 in HNC tissues may be correlated with metastatic lymph node numbers and may predict poor outcomes. Six genes, including CREB3L1, MITF, KLF9, RARA, KLF7, and ETV1, were predicted to be the target genes of SNAI1. The expression of SNAI1 was negatively correlated with tumor purity of HNC, while it was positively correlated with the infiltration of diverse immune cells, such as B cells and macrophages. In the validation phase, the relationships of SNAI1 expression with lymph node metastasis and poor prognosis were verified. Overexpression of SNAI1 might promote lymph node metastasis through complex molecular mechanisms and act as a prognostic indicator in HNC. SNAI1 expression may have a correlation with immune cell infiltrates. Future studies are needed to address these points. Show less

Diabetes, dyslipidemia and hypertension are important metabolic diseases that impose a great burden on many populations worldwide. However, certain population strata have reduced prevalence for all three diseases, but the underlying mechanisms are poorly understood. We sought to identify the phenotypic, genomic and metabolomic characteristics of the low-prevalence population to gain insights into possible innate non-susceptibility against metabolic diseases. We performed k-means cluster analysis of 16,792 subjects using anthropometric and clinical biochemistry data collected by the Taiwan Biobank. Nuclear magnetic resonance spectra-based metabolome analysis was carried out for 217 subjects with normal body mass index, good exercise habits and healthy lifestyles. We found that the gene APOA5 was significantly associated with reduced prevalence of disease, and lesser associations included the genes HIF1A, LIMA1, LPL, MLXIPL, and TRPC4. Blood plasma of subjects belonging to the low disease prevalence cluster exhibited lowered levels of the GlycA inflammation marker, very low-density lipoprotein and low-density lipoprotein cholesterol, triglycerides, valine and leucine compared to controls. Literature mining revealed that these genes and metabolites are biochemically linked, with the linkage between lipoprotein metabolism and inflammation being particularly prominent. The combination of phenomic, genomic and metabolomic analysis may also be applied towards the study of metabolic disease prevalence in other populations. Show less

Clinical diagnosis of SLE is currently challenging due to its heterogeneity. Many autoantibodies are associated with SLE and are considered potential diagnostic markers, but systematic screening and validation of such autoantibodies is lacking. This study aimed to systematically discover new autoantibodies that may be good biomarkers for use in SLE diagnosis. Sera from 15 SLE patients and 5 healthy volunteers were analysed using human proteome microarrays to identify candidate SLE-related autoantibodies. The results were validated by screening of sera from 107 SLE patients, 94 healthy volunteers and 60 disease controls using focussed arrays comprised of autoantigens corresponding to the identified candidate antibodies. Logistic regression was used to derive and validate autoantibody panels that can discriminate SLE disease. Extensive ELISA screening of sera from 294 SLE patients and 461 controls was performed to validate one of the newly discovered autoantibodies. A total of 31, 11 and 18 autoantibodies were identified to be expressed at significantly higher levels in the SLE group than in the healthy volunteers, disease controls and healthy volunteers plus disease control groups, respectively, with 25, 7 and 13 of these differentially expressed autoantibodies being previously unreported. Diagnostic panels comprising anti-RPLP2, anti-SNRPC and anti-PARP1, and anti-RPLP2, anti-PARP1, anti-MAK16 and anti- RPL7A were selected. Performance of the newly discovered anti-MAK16 autoantibody was confirmed by ELISA. Some associations were seen with clinical characteristics of SLE patients, such as disease activity with the level of anti-PARP1 and rash with the level of anti-RPLP2, anti-MAK16 and anti- RPL7A. The combined autoantibody panels identified here show promise for the diagnosis of SLE and for differential diagnosis of other major rheumatic immune diseases. Show less

Impaired insulin secretion from pancreatic islets is a hallmark of type 2 diabetes (T2D). Altered chromatin structure may contribute to the disease. We therefore studied the impact of T2D on open chromatin in human pancreatic islets. We used assay for transposase-accessible chromatin using sequencing (ATAC-seq) to profile open chromatin in islets from T2D and non-diabetic donors. We identified 57,105 and 53,284 ATAC-seq peaks representing open chromatin regions in islets of non-diabetic and diabetic donors, respectively. The majority of ATAC-seq peaks mapped near transcription start sites. Additionally, peaks were enriched in enhancer regions and in regions where islet-specific transcription factors (TFs), e.g. FOXA2, MAFB, NKX2.2, NKX6.1 and PDX1, bind. Islet ATAC-seq peaks overlap with 13 SNPs associated with T2D (e.g. rs7903146, rs2237897, rs757209, rs11708067 and rs878521 near TCF7L2, KCNQ1, HNF1B, ADCY5 and GCK, respectively) and with additional 67 SNPs in LD with known T2D SNPs (e.g. SNPs annotated to GIPR, KCNJ11, GLIS3, IGF2BP2, FTO and PPARG). There was enrichment of open chromatin regions near highly expressed genes in human islets. Moreover, 1,078 open chromatin peaks, annotated to 898 genes, differed in prevalence between diabetic and non-diabetic islet donors. Some of these peaks are annotated to candidate genes for T2D and islet dysfunction (e.g. HHEX, HMGA2, GLIS3, MTNR1B and PARK2) and some overlap with SNPs associated with T2D (e.g. rs3821943 near WFS1 and rs508419 near ANK1). Enhancer regions and motifs specific to key TFs including BACH2, FOXO1, FOXA2, NEUROD1, MAFA and PDX1 were enriched in differential islet ATAC-seq peaks of T2D versus non-diabetic donors. Our study provides new understanding into how T2D alters the chromatin landscape, and thereby accessibility for TFs and gene expression, in human pancreatic islets. Show less

Many single nucleotide polymorphisms (SNPs) have been reported to be associated with migraine susceptibility. However, evidences for their associations with migraine endophenotypes or subtypes are scarce. We aimed to investigate the associations of pre-identified migraine susceptibility loci in Taiwanese with migraine endophenotypes or subtypes, including chronic migraine and allodynia. The associations of six SNPs identified from our previous study, including TRPM8 rs10166942, LRP1 rs1172113, DLG2 rs655484, GFRA1 rs3781545, UPP2 rs7565931, and GPR39 rs10803531, and migraine endophenotypes, including chronic migraine and allodynia were tested. Significant associations in the discovery cohort were validated in the replication cohort. The adjusted odds ratios (aOR) were calculated after controlling for confounders. In total, 1904 patients (mean age 37.5 ± 12.2 years old, female ratio: 77.7%) including 1077 in the discovery cohort and 827 in the replication cohort were recruited. Of them, 584 (30.7%) had chronic migraine. Of the 6 investigated SNPs, TRPM8 rs10166942 T allele-carrying patients were more likely to have chronic migraine than non-T allele carriers in both discovery and replication cohorts and combined samples (33.7% vs. 25.8%, p = 0.004, aOR = 1.62). In addition, T allele carriers reported more allodynic symptoms than non-T allele carriers (3.5 ± 3.7 vs. 2.6 ± 2.8, p < 0.001). However, allodynia severity did not differ between episodic and chronic migraine patients. No further correlations between genetic variants and endophenotypes were noted for the other SNPs. TRPM8 may contribute to the pathogenesis of chronic migraine. However, our study did not support allodynia as a link between them. The underlying mechanisms deserve further investigations. Show less

Non-alcoholic fatty liver disease has been associated with increased mRNA expression of FADS2 in the liver and estimated activity of delta-6 desaturase in serum, encoded by the FADS2 gene. Since DNA methylation in the FADS1/2/3 gene cluster has been previously linked with genetic variants and desaturase activities, we now aimed to discover factors regulating DNA methylation of the CpG sites annotated to FADS1/2 genes. DNA methylation levels in the CpG sites annotated to FADS2 and FADS1 were analyzed from liver samples of 95 obese participants of the Kuopio Obesity Surgery Study (34 men and 61 women, age 49.5 ± 7.7 years, BMI 43.0 ± 5.7 kg/m We found a negative correlation between DNA methylation levels of cg06781209 and cg07999042 and hepatic FADS2 mRNA expression (both p < 0.05), and with estimated delta-6 desaturase activity based on both liver and serum fatty acids (all p < 0.05). Interestingly, the methylation level of cg07999042 (p = 0.001) but not of cg06781209 (p = 0.874) was associated with FADS2 variant rs174616. Genetic variants of FADS2 may contribute to the pathogenesis of non-alcoholic fatty liver disease by modifying DNA methylation. Show less

Impaired adipose tissue insulin signalling is a critical feature of insulin resistance. Here we identify a pathway linking the lipolytic enzyme hormone-sensitive lipase (HSL) to insulin action via the glucose-responsive transcription factor ChREBP and its target, the fatty acid elongase ELOVL6. Genetic inhibition of HSL in human adipocytes and mouse adipose tissue results in enhanced insulin sensitivity and induction of ELOVL6. ELOVL6 promotes an increase in phospholipid oleic acid, which modifies plasma membrane fluidity and enhances insulin signalling. HSL deficiency-mediated effects are suppressed by gene silencing of ChREBP and ELOVL6. Mechanistically, physical interaction between HSL, independent of lipase activity, and the isoform activated by glucose metabolism ChREBPα impairs ChREBPα translocation into the nucleus and induction of ChREBPβ, the isoform with high transcriptional activity that is strongly associated with whole-body insulin sensitivity. Targeting the HSL-ChREBP interaction may allow therapeutic strategies for the restoration of insulin sensitivity. Show less

Aberrant expression of RBM6 has been implicated in the development of human malignancies. However, the bio-function of RBM6 in laryngocarcinoma is still almost blank. Here we identified that RBM6 was downregulated in laryngocarcinoma tissues, as well as laryngocarcinoma cell lines. Notably, the expression level of RBM6 was lower in laryngocarcinoma patients at stage3/4 than that in laryngocarcinoma patients at stage1/2. Upregulation of RBM6 suppressed the proliferation of TU212 and Hep-2 cells, as shown by decreased cell viability and Ki67 level. In parallel, overexpression of RBM6 inhibited invasion and promoted apoptosis of TU212 and Hep-2 cells, as evidenced by downregulation of MMP-2 and MMP-9 protein expression and upregulation of cleaved caspase-3 protein expression. In vivo, RBM6 overexpression repressed the laryngocarcinoma tumor growth. EGFR mRNA level was higher in the laryngocarcinoma tissues than that in the adjacent normal tissues. Moreover, upregulation of RBM6 reduced the expression of EGFR, ERK and p-ERK in vitro and in vivo. Our data suggest that RBM6 as a tumor suppressor represses the growth and progression in laryngocarcinoma. Show less

Ankylosing spondylitis (AS) is a debilitating autoimmune disease affecting tens of millions of people in the world. The genetics of AS is unclear. Analysis of rare AS pedigrees might facilitate our understanding of AS pathogenesis. We used genome-wide linkage analysis and whole-exome sequencing in combination with variant co-segregation verification and haplotype analysis to study an AS pedigree and a sporadic AS patient. We identified a missense variant in the ankyrin repeat and death domain containing 1B gene ANKDD1B from a Han Chinese pedigree with dominantly inherited AS. This variant (p.L87V) co-segregates with all male patients of the pedigree. In females, the penetrance of the symptoms is incomplete with one identified patient out of 5 carriers, consistent with the reduced frequency of AS in females of the general population. We further identified a distinct missense variant affecting a conserved amino acid (p.R102L) of ANKDD1B in a male from 30 sporadic early onset AS patients. Both variants are absent in 500 normal controls. We determined the haplotypes of four major known AS risk loci, including HLA-B*27, 2p15, ERAP1 and IL23R, and found that only HLA-B*27 is strongly associated with patients in our cohort. Together these results suggest that ANKDD1B variants might be associated with AS and genetic analyses of more AS patients are warranted to verify this association. Show less

Metabolic syndrome (MetS), a cluster of metabolic disturbances that increase the risk for cardiovascular disease and diabetes, was because of genetic susceptibility and environmental risk factors. To identify the genetic variants associated with MetS and metabolic components, we conducted a genome-wide association study followed by replications in totally 12,720 participants from the north, north-eastern and eastern China. In combined analyses, independent of the top known signal at rs651821 on APOA5, we newly identified a secondary triglyceride-associated signal at rs180326 on BUD13 (P Show less

Little is known about the extent to which interactions between genetics and epigenetics may affect the risk of complex metabolic diseases and/or their intermediary phenotypes. We performed a genome-wide DNA methylation quantitative trait locus (mQTL) analysis in human adipose tissue of 119 men, where 592,794 single nucleotide polymorphisms (SNPs) were related to DNA methylation of 477,891 CpG sites, covering 99% of RefSeq genes. SNPs in significant mQTLs were further related to gene expression in adipose tissue and obesity related traits. We found 101,911 SNP-CpG pairs (mQTLs) in cis and 5,342 SNP-CpG pairs in trans showing significant associations between genotype and DNA methylation in adipose tissue after correction for multiple testing, where cis is defined as distance less than 500 kb between a SNP and CpG site. These mQTLs include reported obesity, lipid and type 2 diabetes loci, e.g. ADCY3/POMC, APOA5, CETP, FADS2, GCKR, SORT1 and LEPR. Significant mQTLs were overrepresented in intergenic regions meanwhile underrepresented in promoter regions and CpG islands. We further identified 635 SNPs in significant cis-mQTLs associated with expression of 86 genes in adipose tissue including CHRNA5, G6PC2, GPX7, RPL27A, THNSL2 and ZFP57. SNPs in significant mQTLs were also associated with body mass index (BMI), lipid traits and glucose and insulin levels in our study cohort and public available consortia data. Importantly, the Causal Inference Test (CIT) demonstrates how genetic variants mediate their effects on metabolic traits (e.g. BMI, cholesterol, high-density lipoprotein (HDL), hemoglobin A1c (HbA1c) and homeostatic model assessment of insulin resistance (HOMA-IR)) via altered DNA methylation in human adipose tissue. This study identifies genome-wide interactions between genetic and epigenetic variation in both cis and trans positions influencing gene expression in adipose tissue and in vivo (dys)metabolic traits associated with the development of obesity and diabetes. Show less