Novel Lipid Long Intervening Noncoding RNA, Oligodendrocyte Maturation-Associated Long Intergenic Noncoding RNA, Regulates the Liver Steatosis Gene Stearoyl-Coenzyme A Desaturase As an Enhancer RNA

Single nucleus RNA-sequencing integrated into risk variant colocalization discovers 17 cell-type-specific abdominal obesity genes for metabolic dysfunction-associated steatotic liver disease

BACKGROUND

Abdominal obesity increases the risk for non-alcoholic fatty liver disease (NAFLD), now known as metabolic dysfunction-associated steatotic liver disease (MASLD).

METHODS

To elucidate the directional cell-type level biological mechanisms underlying the association between abdominal obesity and MASLD, we integrated adipose and liver single nucleus RNA-sequencing and bulk cis-expression quantitative trait locus (eQTL) data with the UK Biobank genome-wide association study (GWAS) data using colocalization. Then we used colocalized cis-eQTL variants as instrumental variables in Mendelian randomization (MR) analyses, followed by functional validation experiments on the target genes of the cis-eQTL variants.

FINDINGS

We identified 17 colocalized abdominal obesity GWAS variants, regulating 17 adipose cell-type marker genes. Incorporating these 17 variants into MR discovers a putative tissue-of-origin, cell-type-aware causal effect of abdominal obesity on MASLD consistently with multiple MR methods without significant evidence for pleiotropy or heterogeneity. Single cell data confirm the adipocyte-enriched mean expression of the 17 genes. Our cellular experiments across human adipogenesis identify risk variant -specific epigenetic and transcriptional mechanisms. Knocking down two of the 17 genes, PPP2R5A and SH3PXD2B, shows a marked decrease in adipocyte lipidation and significantly alters adipocyte function and adipogenesis regulator genes, including DGAT2, LPL, ADIPOQ, PPARG, and SREBF1. Furthermore, the 17 genes capture a characteristic MASLD expression signature in subcutaneous adipose tissue.

INTERPRETATION

Overall, we discover a significant cell-type level effect of abdominal obesity on MASLD and trace its biological effect to adipogenesis.

FUNDING

NIH grants R01HG010505, R01DK132775, and R01HL170604; the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation program (Grant No. 802825), Academy of Finland (Grants Nos. 333021), the Finnish Foundation for Cardiovascular Research the Sigrid Jusélius Foundation and the Jane and Aatos Erkko Foundation; American Association for the Study of Liver Diseases (AASLD) Advanced Transplant Hepatology award and NIH/NIDDK (P30DK41301) Pilot and Feasibility award; NIH/NIEHS F32 award (F32ES034668); Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2021), the Academy of Finland grant (Contract no. 138006); Academy of Finland (Grant Nos 335443, 314383, 272376 and 266286), Sigrid Jusélius Foundation, Finnish Medical Foundation, Finnish Diabetes Research Foundation, Novo Nordisk Foundation (#NNF20OC0060547, NNF17OC0027232, NNF10OC1013354) and Government Research Funds to Helsinki University Hospital; Orion Research Foundation, Maud Kuistila Foundation, Finish Medical Foundation, and University of Helsinki.

Hepatokine ITIH3 protects against hepatic steatosis by downregulating mitochondrial bioenergetics and de novo lipogenesis

Recent studies demonstrate that liver secretory proteins, also known as hepatokines, regulate normal development, obesity, and simple steatosis to non-alcoholic steatohepatitis (NASH) progression. Using a panel of ∼100 diverse inbred strains of mice and a cohort of bariatric surgery patients, we found that one such hepatokine, inter-trypsin inhibitor heavy chain 3 (ITIH3), was progressively lower in severe non-alcoholic fatty liver disease (NAFLD) disease states highlighting an inverse relationship between Itih3/ITIH3 expression and NAFLD severity. Follow-up animal and cell culture models demonstrated that hepatic ITIH3 overexpression lowered liver triglyceride and lipid droplet accumulation, respectively. Conversely, ITIH3 knockdown in mice increased the liver triglyceride in two independent NAFLD models. Mechanistically, ITIH3 reduced mitochondrial respiration and this, in turn, reduced liver triglycerides, via downregulated de novo lipogenesis. This was accompanied by increased STAT1 signaling and Stat3 expression, both of which are known to protect against NAFLD/NASH. Our findings indicate hepatokine ITIH3 as a potential biomarker and/or treatment for NAFLD.

Key genes involved in nonalcoholic steatohepatitis improvement after bariatric surgery

Background

Nonalcoholic steatohepatitis (NASH) is the advanced stage of nonalcoholic fatty liver disease (NAFLD), one of the most prevalent chronic liver diseases. The effectiveness of bariatric surgery in treating NASH and preventing or even reversing liver fibrosis has been demonstrated in numerous clinical studies, but the underlying mechanisms and crucial variables remain unknown.

Methods

Using the GSE135251 dataset, we examined the gene expression levels of NASH and healthy livers. Then, the differentially expressed genes (DEGs) of patients with NASH, at baseline and one year after bariatric surgery, were identified in GSE83452. We overlapped the hub genes performed by protein-protein interaction (PPI) networks and DEGs with different expression trends in both datasets to obtain key genes. Genomic enrichment analysis (GSEA) and genomic variation analysis (GSVA) were performed to search for signaling pathways of key genes. Meanwhile, key molecules that regulate the key genes are found through the construction of the ceRNA network. NASH mice were induced by a high-fat diet (HFD) and underwent sleeve gastrectomy (SG). We then cross-linked the DEGs in clinical and animal samples using quantitative polymerase chain reaction (qPCR) and validated the key genes.

Results

Seven key genes (FASN, SCD, CD68, HMGCS1, SQLE, CXCL10, IGF1) with different expression trends in GSE135251 and GSE83452 were obtained with the top 30 hub genes selected by PPI. The expression of seven key genes in mice after SG was validated by qPCR. Combined with the qPCR results from NASH mice, the four genes FASN, SCD, HMGCS1, and CXCL10 are consistent with the biological analysis. The GSEA results showed that the 'cholesterol homeostasis' pathway was enriched in the FASN, SCD, HMGCS1, and SQLE high-expression groups. The high-expression groups of CD68 and CXCL10 were extremely enriched in inflammation-related pathways. The construction of the ceRNA network obtained microRNAs and ceRNAs that can regulate seven key genes expression.

Conclusion

In summary, this study contributes to our understanding of the mechanisms by which bariatric surgery improves NASH, and to the development of potential biomarkers for the treatment of NASH.

Breast Cancer Chemoresistance: Insights into the Regulatory Role of lncRNA

Long noncoding RNAs (lncRNAs) are a subclass of noncoding RNAs composed of more than 200 nucleotides without the ability to encode functional proteins. Given their involvement in critical cellular processes such as gene expression regulation, transcription, and translation, lncRNAs play a significant role in organism homeostasis. Breast cancer (BC) is the second most common cancer worldwide and evidence has shown a relationship between aberrant lncRNA expression and BC development. One of the main obstacles in BC control is multidrug chemoresistance, which is associated with the deregulation of multiple mechanisms such as efflux transporter activity, mitochondrial metabolism reprogramming, and epigenetic regulation as well as apoptosis and autophagy. Studies have shown the involvement of a large number of lncRNAs in the regulation of such pathways. However, the underlying mechanism is not clearly elucidated. In this review, we present the principal mechanisms associated with BC chemoresistance that can be directly or indirectly regulated by lncRNA, highlighting the importance of lncRNA in controlling BC chemoresistance. Understanding these mechanisms in deep detail may interest the clinical outcome of BC patients and could be used as therapeutic targets to overcome BC therapy resistance.

PTHrP Regulates Fatty Acid Metabolism via Novel lncRNA in Breast Cancer Initiation and Progression Models

Parathyroid hormone-related peptide (PTHrP) is the primary cause of malignancy-associated hypercalcemia (MAH). We previously showed that PTHrP ablation, in the MMTV-PyMT murine model of breast cancer (BC) progression, can dramatically prolong tumor latency, slow tumor growth, and prevent metastatic spread. However, the signaling mechanisms using lineage tracing have not yet been carefully analyzed. Here, we generated Pthrpflox/flox; Cre+ mT/mG mice (KO) and Pthrpwt/wt; Cre+ mT/mG tumor mice (WT) to examine the signaling pathways under the control of PTHrP from the early to late stages of tumorigenesis. GFP+ mammary epithelial cells were further enriched for subsequent RNA sequencing (RNAseq) analyses. We observed significant upregulation of cell cycle signaling and fatty acid metabolism in PTHrP WT tumors, which are linked to tumor initiation and progression. Next, we observed that the expression levels of a novel lncRNA, GM50337, along with stearoyl-Coenzyme A desaturase 1 (Scd1) are significantly upregulated in PTHrP WT but not in KO tumors. We further validated a potential human orthologue lncRNA, OLMALINC, together with SCD1 that can be regulated via PTHrP in human BC cell lines. In conclusion, these novel findings could be used to develop targeted strategies for the treatment of BC and its metastatic complications.

Cross-tissue omics analysis discovers ten adipose genes encoding secreted proteins in obesity-related non-alcoholic fatty liver disease

BACKGROUND

Non-alcoholic fatty liver disease (NAFLD) is a fast-growing, underdiagnosed, epidemic. We hypothesise that obesity-related inflammation compromises adipose tissue functions, preventing efficient fat storage, and thus driving ectopic fat accumulation into the liver.

METHODS

To identify adipose-based mechanisms and potential serum biomarker candidates (SBCs) for NAFLD, we utilise dual-tissue RNA-sequencing (RNA-seq) data in adipose tissue and liver, paired with histology-based NAFLD diagnosis, from the same individuals in a cohort of obese individuals. We first scan for genes that are differentially expressed (DE) for NAFLD in obese individuals' subcutaneous adipose tissue but not in their liver; encode proteins secreted to serum; and show preferential adipose expression. Then the identified genes are filtered to key adipose-origin NAFLD genes by best subset analysis, knockdown experiments during human preadipocyte differentiation, recombinant protein treatment experiments in human liver HepG2 cells, and genetic analysis.

FINDINGS

We discover a set of genes, including 10 SBCs, that may modulate NAFLD pathogenesis by impacting adipose tissue function. Based on best subset analysis, we further follow-up on two SBCs CCDC80 and SOD3 by knockdown in human preadipocytes and subsequent differentiation experiments, which show that they modulate crucial adipogenesis genes, LPL, SREBPF1, and LEP. We also show that treatment of the liver HepG2 cells with the CCDC80 and SOD3 recombinant proteins impacts genes related to steatosis and lipid processing, including PPARA, NFE2L2, and RNF128. Finally, utilizing the adipose NAFLD DE gene cis-regulatory variants associated with serum triglycerides (TGs) in extensive genome-wide association studies (GWASs), we demonstrate a unidirectional effect of serum TGs on NAFLD with Mendelian Randomization (MR) analysis. We also demonstrate that a single SNP regulating one of the SBC genes, rs2845885, produces a significant MR result by itself. This supports the conclusion that genetically regulated adipose expression of the NAFLD DE genes may contribute to NAFLD through changes in serum TG levels.

INTERPRETATION

Our results from the dual-tissue transcriptomics screening improve the understanding of obesity-related NAFLD by providing a targeted set of 10 adipose tissue-active genes as new serum biomarker candidates for the currently grossly underdiagnosed fatty liver disease.

FUNDING

The work was supported by NIH grants R01HG010505 and R01DK132775. The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The KOBS study (J. P.) was supported by the Finnish Diabetes Research Foundation, Kuopio University Hospital Project grant (EVO/VTR grants 2005-2019), and the Academy of Finland grant (Contract no. 138006). This study was funded by the European Research Council under the European Union's Horizon 2020 research and innovation program (Grant No. 802825 to M. U. K.). K. H. P. was funded by the Academy of Finland (grant numbers 272376, 266286, 314383, and 335443), the Finnish Medical Foundation, Gyllenberg Foundation, Novo Nordisk Foundation (grant numbers NNF10OC1013354, NNF17OC0027232, and NNF20OC0060547), Finnish Diabetes Research Foundation, Finnish Foundation for Cardiovascular Research, University of Helsinki, and Helsinki University Hospital and Government Research Funds. I. S. was funded by the Instrumentarium Science Foundation. Personal grants to U. T. A. were received from the Matti and Vappu Maukonen Foundation, Ella och Georg Ehrnrooths Stiftelse and the Finnish Foundation for Cardiovascular Research.

CPT1B, a metabolic molecule, is also an independent risk factor in CN-AML

BACKGROUND

Fatty acid oxidation has been considered as an important energy source for tumorigenesis and development. Several studies have investigated the role of CPT1A, a kind of fatty acid oxidation rate-limiting enzyme, in AML. However, prognostic value and regulatory network of another subtype, CPT1B in AML remains elusive. This study aims to clarify the independent prognostic role of CPT1B in CN-AML based on clinical data and molecular level data (mRNA, miRNA and lncRNA).

OBJECTIVE

The aim of this study is to investigate the prognostic value of CPT1B in AML patients.

METHODS

First, we analyzed the CPT1B expression in AML cohort via the online database "GEPIA". Subsequently, miRNA-mRNA and ceRNA networks were constructed to help predict the role of CPT1B in AML. Several molecules which showed the prognostic value and metabolic function of CPT1B were identified. Finally, the expression of CPT1B in our own cohort of 324 CN-AML patients was analyzed to clarify the results.

RESULTS

It was found that CPT1B was markedly higher in AML patients compared to normal people and this upregulation was associated with the poor clinical outcome. Several molecules revealed the possible regulatory mechanism of CPT1B in AML.

CONCLUSION

CPT1B is a potential prognostic factor and a therapeutic target for AML treatment.

Cholesterol-Related lncRNAs as Response Predictors of Atorvastatin Treatment in Chilean Hypercholesterolemic Patients: A Pilot Study

Statins are currently the treatment of choice for hypercholesterolemia. However, wide interindividual variability has been observed in the response to treatment. Recent studies have reported the role of lncRNAs in the metabolism of lipids; nevertheless, there are few studies to date that show their role in the response to treatment with statins. Thus, the aim of this study was to assess the levels of expression of three lncRNAs (RP1-13D10.2; MANTIS; lncHR1) associated with genes involved in cholesterol homeostasis in leukocyte cells of hypercholesterolemic patients after treatment with atorvastatin and compare them with levels in subjects with normal cholesterol levels. A secondary aim was to assess the levels of expression in monocytic THP-1 cells differentiated to macrophages. The study included 20 subjects with normal cholesterol (NC) levels and 20 individuals with hypercholesterolemia (HC). The HC patients were treated with atorvastatin (20 mg/day/4 weeks). THP-1 cells were differentiated to macrophages with PMA and treated with different doses of atorvastatin for 24 h. Expression of lncRNAs was determined by RT-qPCR. The lncRNAs RP1-13D10.2 (p < 0.0001), MANTIS (p = 0.0013) and lncHR1 (p < 0.0001) presented increased expression in HC subjects compared with NC subjects. Furthermore, atorvastatin had a negative regulatory effect on the expression of lncHR1 (p < 0.0001) in HC subjects after treatment. In vitro, all the lncRNAs showed significant differences in expression after atorvastatin treatment. Our findings show that the lncRNAs tested present differential expression in HC patients and play a role in the variability reported in the response to atorvastatin treatment. Further research is needed to clarify the biological impact of these lncRNAs on cholesterol homeostasis and treatment with statins.

Enhancer RNA (eRNA) in Human Diseases

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

Downregulation of hepatic lipopolysaccharide binding protein improves lipogenesis-induced liver lipid accumulation

Circulating lipopolysaccharide-binding protein (LBP) is increased in individuals with liver steatosis. We aimed to evaluate the possible impact of liver LBP downregulation using lipid nanoparticle-containing chemically modified LBP small interfering RNA (siRNA) (LNP-Lbp UNA-siRNA) on the development of fatty liver. Weekly LNP-Lbp UNA-siRNA was administered to mice fed a standard chow diet, a high-fat and high-sucrose diet, and a methionine- and choline-deficient diet (MCD). In mice fed a high-fat and high-sucrose diet, which displayed induced liver lipogenesis, LBP downregulation led to reduced liver lipid accumulation, lipogenesis (mainly stearoyl-coenzyme A desaturase 1 [Scd1]) and lipid peroxidation-associated oxidative stress markers. LNP-Lbp UNA-siRNA also resulted in significantly decreased blood glucose levels during an insulin tolerance test. In mice fed a standard chow diet or an MCD, in which liver lipogenesis was not induced or was inhibited (especially Scd1 mRNA), liver LBP downregulation did not impact on liver steatosis. The link between hepatocyte LBP and lipogenesis was further confirmed in palmitate-treated Hepa1-6 cells, in primary human hepatocytes, and in subjects with morbid obesity. Altogether, these data indicate that siRNA against liver Lbp mRNA constitutes a potential target therapy for obesity-associated fatty liver through the modulation of hepatic Scd1.

Bioinformatics analysis of lncRNAs in the occurrence and development of osteosarcoma

BACKGROUND

Osteosarcoma (OS) is a disease with high mortality in children and adolescents, and metastasis is one of its important clinical features. However, the molecular mechanism of OS occurrence is not completely clear. Thus, we screened potential biomarkers of OS and analyze their prognostic value.

METHODS

The Cancer Genome Atlas (TCGA) datasets were used to analyze the differential lncRNAs in patients with OS of different immune score and the lncRNAs expressed by immune cells. Cox regression was used to develop the prognosis prediction model and specify the prognosis outcomes. Risk-proportional regression model was constructed, and the samples were divided into high and low groups based on the risk scores for the survival analysis. The areas under the receiver operating characteristic (ROC) curve were calculated and the risk-score model was verified. Finally, using 4 gene sets (comprising chemokines, immune checkpoint blockades, immune activity-related genes, and immune cells), and 4 analysis tools (CIBERSORT, TIMER, XCELL and MCP) to evaluated tumor immune infiltration.

RESULTS

Twenty-nine long non-coding ribonucleic acids (lncRNAs) were obtained from the intersection of the screened lncRNAs. Caspase recruitment domain-containing protein 8-antisense RNA 1 (CARD8-AS1), lncRNA five prime to Xist (FTX), KAT8 regulatory NSL complex unit 1-antisense RNA 1 (KANSL1-AS1), Neuroplastin Intronic Transcript 1 (NPTN-IT1), oligodendrocyte maturation-associated long intervening non-coding RNA (OLMALINC) and RPARP Antisense RNA 1 (RPARP-AS1) were found to be correlated with survival. Univariate and multivariate regression analysis showed risk score [HR (hazard ratio) 3.5, P value 0.0043; HR 3.7, P value 0.0033] and metastasis (HR 4.7, P value 6.60E-05; HR 4.8, P value 8.36E-05) were the key factors of patients with OS. The areas under curves (AUCs) of the 1-, 3-, and 5-year ROC curves of the prognostic model were 0.715, 0.729, and 0.771. The low-risk patients tended to have a high abundance of immune cells.

CONCLUSIONS

This study showed that a risk score based on 6 lncRNAs has potential value in the prognosis of OS, and patients with low-risk scores have high immune cell infiltration and good prognosis. This study may enrich understandings of underlying mechanisms related to the occurrence and development of OS.

Multi-omics to characterize the functional relationships of R-loops with epigenetic modifications, RNAPII transcription and gene expression

Abnormal accumulation of R-loops results in replication stress, genome instability, chromatin alterations and gene silencing. Little research has been done to characterize functional relationships among R-loops, histone marks, RNA polymerase II (RNAPII) transcription and gene regulation. We built extremely randomized trees (ETs) models to predict the genome-wide R-loops using RNAPII and multiple histone modifications chromatin immunoprecipitation (ChIP)-seq, DNase-seq, Global Run-On sequencing (GRO-seq) and R-loop profiling data. We compared the performance of ET models to multiple machine learning approaches, and the proposed ET models achieved the best and extremely robust performances. Epigenetic profiles are highly predictive of R-loops genome-widely and they are strongly associated with R-loop formation. In addition, the presence of R-loops is significantly correlated with RNAPII transcription activity, H3K4me3 and open chromatin around the transcription start site, and H3K9me1 and H3K9me3 around the transcription termination site. RNAPII pausing defects were correlated with 5′R-loops accumulation, and transcriptional termination defects and read-throughs were correlated with 3′R-loops accumulation. Furthermore, we found driver genes with 5′R-loops and RNAPII pausing defects express significantly higher and genes with 3′R-loops and read-through transcription express significantly lower than genes without R-loops. These driver genes are enriched with chromosomal instability, Hippo–Merlin signaling Dysregulation, DNA damage response and TGF-β pathways, indicating R-loops accumulating at the 5′ end of genes play oncogenic roles, whereas at the 3′ end of genes play tumor-suppressive roles in tumorigenesis.

Epigenome-wide association study of incident type 2 diabetes: a meta-analysis of five prospective European cohorts

AIMS/HYPOTHESIS

Type 2 diabetes is a complex metabolic disease with increasing prevalence worldwide. Improving the prediction of incident type 2 diabetes using epigenetic markers could help tailor prevention efforts to those at the highest risk. The aim of this study was to identify predictive methylation markers for incident type 2 diabetes by combining epigenome-wide association study (EWAS) results from five prospective European cohorts.

METHODS

We conducted a meta-analysis of EWASs in blood collected 7-10 years prior to type 2 diabetes diagnosis. DNA methylation was measured with Illumina Infinium Methylation arrays. A total of 1250 cases and 1950 controls from five longitudinal cohorts were included: Doetinchem, ESTHER, KORA1, KORA2 and EPIC-Norfolk. Associations between DNA methylation and incident type 2 diabetes were examined using robust linear regression with adjustment for potential confounders. Inverse-variance fixed-effects meta-analysis of cohort-level individual CpG EWAS estimates was performed using METAL. The methylGSA R package was used for gene set enrichment analysis. Confirmation of genome-wide significant CpG sites was performed in a cohort of Indian Asians (LOLIPOP, UK).

RESULTS

The meta-analysis identified 76 CpG sites that were differentially methylated in individuals with incident type 2 diabetes compared with control individuals (p values <1.1 × 10-7). Sixty-four out of 76 (84.2%) CpG sites were confirmed by directionally consistent effects and p values <0.05 in an independent cohort of Indian Asians. However, on adjustment for baseline BMI only four CpG sites remained genome-wide significant, and addition of the 76 CpG methylation risk score to a prediction model including established predictors of type 2 diabetes (age, sex, BMI and HbA1c) showed no improvement (AUC 0.757 vs 0.753). Gene set enrichment analysis of the full epigenome-wide results clearly showed enrichment of processes linked to insulin signalling, lipid homeostasis and inflammation.

CONCLUSIONS/INTERPRETATION

By combining results from five European cohorts, and thus significantly increasing study sample size, we identified 76 CpG sites associated with incident type 2 diabetes. Replication of 64 CpGs in an independent cohort of Indian Asians suggests that the association between DNA methylation levels and incident type 2 diabetes is robust and independent of ethnicity. Our data also indicate that BMI partly explains the association between DNA methylation and incident type 2 diabetes. Further studies are required to elucidate the underlying biological mechanisms and to determine potential causal roles of the differentially methylated CpG sites in type 2 diabetes development.

The PAF1 complex promotes 3' processing of pervasive transcripts

The PAF1 complex (PAF1C) functions in multiple transcriptional processes involving RNA polymerase II (RNA Pol II). Enhancer RNAs (eRNAs) and promoter upstream transcripts (PROMPTs) are pervasive transcripts transcribed by RNA Pol II and degraded rapidly by the nuclear exosome complex after 3' endonucleolytic cleavage by the Integrator complex (Integrator). Here we show that PAF1C has a role in termination of eRNAs and PROMPTs that are cleaved 1-3 kb downstream of the transcription start site. Mechanistically, PAF1C facilitates recruitment of Integrator to sites of pervasive transcript cleavage, promoting timely cleavage and transcription termination. We also show that PAF1C recruits Integrator to coding genes, where PAF1C then dissociates from Integrator upon entry into processive elongation. Our results demonstrate a function of PAF1C in limiting the length and accumulation of pervasive transcripts that result from non-productive transcription.

Targeting of noncoding RNAs encoded by a novel MYC enhancers inhibits the proliferation of human hepatic carcinoma cells in vitro

The proto-oncogene MYC is important for development and cell growth, however, its abnormal regulation causes cancer. Recent studies identified distinct enhancers of MYC in various cancers, but any MYC enhancer(s) in hepatocellular carcinoma (HCC) remain(s) elusive. By analyzing H3K27ac enrichment and enhancer RNA (eRNA) expression in cultured HCC cells, we identified six putative MYC enhancer regions. Amongst these, two highly active enhancers, located ~ 800 kb downstream of the MYC gene, were identified by qRT-PCR and reporter assays. We functionally confirmed these enhancers by demonstrating a significantly reduced MYC expression and cell proliferation upon CRISPR/Cas9-based deletion and/or antisense oligonucleotide (ASO)-mediated inhibition. In conclusion, we identified potential MYC enhancers of HCC and propose that the associated eRNAs may be suitable targets for HCC treatment.

Identification of 90 NAFLD GWAS loci and establishment of NAFLD PRS and causal role of NAFLD in coronary artery disease

The prevalence of non-alcoholic fatty liver disease (NAFLD), now also known as metabolic dysfunction-associated fatty liver disease (MAFLD), is rapidly increasing worldwide due to the ongoing obesity epidemic. However, currently the NALFD diagnosis requires non-readily available imaging technologies or liver biopsy, which has drastically limited the sample sizes of NAFLD studies and hampered the discovery of its genetic component. Here we utilized the large UK Biobank (UKB) to accurately estimate the NAFLD status in UKB based on common serum traits and anthropometric measures. Scoring all individuals in UKB for NAFLD risk resulted in 28,396 NAFLD cases and 108,652 healthy individuals at a >90% confidence level. Using this imputed NAFLD status to perform the largest NAFLD genome-wide association study (GWAS) to date, we identified 94 independent (R² < 0.2) NAFLD GWAS loci, of which 90 have not been identified before; built a polygenic risk score (PRS) model to predict the genetic risk of NAFLD; and used the GWAS variants of imputed NAFLD for a tissue-aware Mendelian randomization analysis that discovered a significant causal effect of NAFLD on coronary artery disease (CAD). In summary, we accurately estimated the NAFLD status in UKB using common serum traits and anthropometric measures, which empowered us to identify 90 GWAS NAFLD loci, build NAFLD PRS, and discover a significant causal effect of NAFLD on CAD.

Hepatocellular carcinoma in nonalcoholic fatty liver disease: A growing challenge

Nonalcoholic fatty liver disease (NAFLD) is the most common cause of liver disease worldwide, and its prevalence increases continuously. As it predisposes to hepatocellular carcinoma both in the presence and in the absence of cirrhosis, it is not surprising that the incidence of NAFLD-related hepatocellular carcinoma would also rise. Some of the mechanisms involved in hepatocarcinogenesis are particular to individuals with fatty liver, and they help explain why liver cancer develops even in patients without cirrhosis. Genetic and immune-mediated mechanisms seem to play an important role in the development of hepatocellular carcinoma in this population. Currently, it is consensual that patients with NAFLD-related cirrhosis should be surveilled with ultrasonography every 6 mo (with or without alpha-fetoprotein), but it is known that they are less likely to follow this recommendation than individuals with other kinds of liver disease. Moreover, the performance of the methods of surveillance are lower in NAFLD than they are in other liver diseases. Furthermore, it is not clear which subgroups of patients without cirrhosis should undergo surveillance. Understanding the mechanisms of hepatocarcinogenesis in NAFLD could hopefully lead to the identification of biomarkers to be used in the surveillance for liver cancer in these individuals. By improving surveillance, tumors could be detected in earlier stages, amenable to curative treatments.

Animal-eRNAdb: a comprehensive animal enhancer RNA database

Enhancer RNAs (eRNAs) are a class of non-coding RNAs transcribed from enhancers. As the markers of active enhancers, eRNAs play important roles in gene regulation and are associated with various complex traits and characteristics. With increasing attention to eRNAs, numerous eRNAs have been identified in different human tissues. However, the expression landscape, regulatory network and potential functions of eRNAs in animals have not been fully elucidated. Here, we systematically characterized 185 177 eRNAs from 5085 samples across 10 species by mapping the RNA sequencing data to the regions of known enhancers. To explore their potential functions based on evolutionary conservation, we investigated the sequence similarity of eRNAs among multiple species. In addition, we identified the possible associations between eRNAs and transcription factors (TFs) or nearby genes to decipher their possible regulators and target genes, as well as characterized trait-related eRNAs to explore their potential functions in biological processes. Based on these findings, we further developed Animal-eRNAdb (http://gong_lab.hzau.edu.cn/Animal-eRNAdb/), a user-friendly database for data searching, browsing and downloading. With the comprehensive characterization of eRNAs in various tissues of different species, Animal-eRNAdb may greatly facilitate the exploration of functions and mechanisms of eRNAs.

Identification of TBX15 as an adipose master trans regulator of abdominal obesity genes

Background

Obesity predisposes individuals to multiple cardiometabolic disorders, including type 2 diabetes (T2D). As body mass index (BMI) cannot reliably differentiate fat from lean mass, the metabolically detrimental abdominal obesity has been estimated using waist-hip ratio (WHR). Waist-hip ratio adjusted for body mass index (WHRadjBMI) in turn is a well-established sex-specific marker for abdominal fat and adiposity, and a predictor of adverse metabolic outcomes, such as T2D. However, the underlying genes and regulatory mechanisms orchestrating the sex differences in obesity and body fat distribution in humans are not well understood.

Methods

We searched for genetic master regulators of WHRadjBMI by employing integrative genomics approaches on human subcutaneous adipose RNA sequencing (RNA-seq) data (n ~ 1400) and WHRadjBMI GWAS data (n ~ 700,000) from the WHRadjBMI GWAS cohorts and the UK Biobank (UKB), using co-expression network, transcriptome-wide association study (TWAS), and polygenic risk score (PRS) approaches. Finally, we functionally verified our genomic results using gene knockdown experiments in a human primary cell type that is critical for adipose tissue function.

Results

Here, we identified an adipose gene co-expression network that contains 35 obesity GWAS genes and explains a significant amount of polygenic risk for abdominal obesity and T2D in the UKB (n = 392,551) in a sex-dependent way. We showed that this network is preserved in the adipose tissue data from the Finnish Kuopio Obesity Study and Mexican Obesity Study. The network is controlled by a novel adipose master transcription factor (TF), TBX15, a WHRadjBMI GWAS gene that regulates the network in trans. Knockdown of TBX15 in human primary preadipocytes resulted in changes in expression of 130 network genes, including the key adipose TFs, PPARG and KLF15, which were significantly impacted (FDR < 0.05), thus functionally verifying the trans regulatory effect of TBX15 on the WHRadjBMI co-expression network.

Conclusions

Our study discovers a novel key function for the TBX15 TF in trans regulating an adipose co-expression network of 347 adipose, mitochondrial, and metabolically important genes, including PPARG, KLF15, PPARA, ADIPOQ, and 35 obesity GWAS genes. Thus, based on our converging genomic, transcriptional, and functional evidence, we interpret the role of TBX15 to be a main transcriptional regulator in the adipose tissue and discover its importance in human abdominal obesity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00939-2.

LXR directly regulates glycosphingolipid synthesis and affects human CD4+ T cell function

The liver X receptor (LXR) is a key transcriptional regulator of cholesterol, fatty acid, and phospholipid metabolism. Dynamic remodeling of immunometabolic pathways, including lipid metabolism, is a crucial step in T cell activation. Here, we explored the role of LXR-regulated metabolic processes in primary human CD4⁺ T cells and their role in controlling plasma membrane lipids (glycosphingolipids and cholesterol), which strongly influence T cell immune signaling and function. Crucially, we identified the glycosphingolipid biosynthesis enzyme glucosylceramide synthase as a direct transcriptional LXR target. LXR activation by agonist GW3965 or endogenous oxysterol ligands significantly altered the glycosphingolipid:cholesterol balance in the plasma membrane by increasing glycosphingolipid levels and reducing cholesterol. Consequently, LXR activation lowered plasma membrane lipid order (stability), and an LXR antagonist could block this effect. LXR stimulation also reduced lipid order at the immune synapse and accelerated activation of proximal T cell signaling molecules. Ultimately, LXR activation dampened proinflammatory T cell function. Finally, compared with responder T cells, regulatory T cells had a distinct pattern of LXR target gene expression corresponding to reduced lipid order. This suggests LXR-driven lipid metabolism could contribute to functional specialization of these T cell subsets. Overall, we report a mode of action for LXR in T cells involving the regulation of glycosphingolipid and cholesterol metabolism and demonstrate its relevance in modulating T cell function.

Deep Insight Into Long Non-coding RNA and mRNA Transcriptome Profiling in HepG2 Cells Expressing Genotype IV Swine Hepatitis E Virus ORF3

Swine hepatitis E (swine HE) is a new type of zoonotic infectious disease caused by the swine hepatitis E virus (swine HEV). Open reading frame 3 (ORF3) is an important virulent protein of swine HEV, but its function still is mainly unclear. In this study, we generated adenoviruses ADV4-ORF3 and ADV4 negative control (ADV4-NC), which successfully mediated overexpression of enhanced green fluorescent protein (EGFP)-ORF3 and EGFP, respectively, in HepG2 cells. High-throughput sequencing was used to screen for differentially expressed long non-coding RNAs (lncRNAs) and messenger RNAs (mRNAs). The cis-target genes of lncRNAs were predicted, functional enrichment (Gene Ontology [GO] and Kyoto Encyclopedia of Genes and Genomes [KEGG]) was performed, and 12 lncRNAs with statistically significant different expressions (p ≤ 0.05 and q ≤ 1) were selected for further quantitative real-time reverse transcription (qRT-PCR) validation. In HepG2 cells, we identified 62 significantly differentially expressed genes (DEGs) (6,564 transcripts) and 319 lncRNAs (124 known lncRNAs and 195 novel lncRNAs) that were affected by ORF3, which were involved in systemic lupus erythematosus, Staphylococcus aureus infection, signaling pathways pluripotency regulation of stem cells, the peroxisome proliferator-activated receptor (PPAR) signaling pathway, and platinum drug resistance pathways. Cis-target gene prediction identified 45 lncRNAs corresponding to candidate mRNAs, among which eight were validated by qRT-PCR: LINC02476 (two transcripts), RAP2C-AS1, AC016526, AL139099, and ZNF337-AS1 (3 transcripts). Our results revealed that the lncRNA profile in host cells affected by ORF3, swine HEV ORF3, might affect the pentose and glucuronate interconversions and mediate the formation of obstructive jaundice by influencing bile secretion, which will help to determine the function of ORF3 and the infection mechanism and treatment of swine HE.

Serum aromatic and branched-chain amino acids associated with NASH demonstrate divergent associations with serum lipids

BACKGROUND & AIMS

Non-alcoholic fatty liver disease (NAFLD) has been associated with multiple metabolic abnormalities. By applying a non-targeted metabolomics approach, we aimed at investigating whether serum metabolite profile that associates with NAFLD would differ in its association with NAFLD-related metabolic risk factors.

METHODS & RESULTS

A total of 233 subjects (mean ± SD: 48.3 ± 9.3 years old; BMI: 43.1 ± 5.4 kg/m² ; 64 male) undergoing bariatric surgery were studied. Of these participants, 164 with liver histology could be classified as normal liver (n = 79), simple steatosis (SS, n = 40) or non-alcoholic steatohepatitis (NASH, n = 45). Among the identified fasting serum metabolites with higher levels in those with NASH when compared to those with normal phenotype were the aromatic amino acids (AAAs: tryptophan, tyrosine and phenylalanine), the branched-chain amino acids (BCAAs: leucine and isoleucine), a phosphatidylcholine (PC(16:0/16:1)) and uridine (all FDRp < 0.05). Only tryptophan was significantly higher in those with NASH compared to those with SS (FDRp < 0.05). Only the AAAs tryptophan and tyrosine correlated positively with serum total and LDL cholesterol (FDRp < 0.1), and accordingly, with liver LDLR at mRNA expression level. In addition, tryptophan was the single AA associated with liver DNA methylation of CpG sites known to be differentially methylated in those with NASH.

CONCLUSIONS

We found that serum levels of the NASH-related AAAs and BCAAs demonstrate divergent associations with serum lipids. The specific correlation of tryptophan with LDL-c may result from the molecular events affecting LDLR mRNA expression and NASH-associated methylation of genes in the liver.

Integrative analysis of liver-specific non-coding regulatory SNPs associated with the risk of coronary artery disease

Genetic factors underlying coronary artery disease (CAD) have been widely studied using genome-wide association studies (GWASs). However, the functional understanding of the CAD loci has been limited by the fact that a majority of GWAS variants are located within non-coding regions with no functional role. High cholesterol and dysregulation of the liver metabolism such as non-alcoholic fatty liver disease confer an increased risk of CAD. Here, we studied the function of non-coding single-nucleotide polymorphisms in CAD GWAS loci located within liver-specific enhancer elements by identifying their potential target genes using liver cis-eQTL analysis and promoter Capture Hi-C in HepG2 cells. Altogether, 734 target genes were identified of which 121 exhibited correlations to liver-related traits. To identify potentially causal regulatory SNPs, the allele-specific enhancer activity was analyzed by (1) sequence-based computational predictions, (2) quantification of allele-specific transcription factor binding, and (3) STARR-seq massively parallel reporter assay. Altogether, our analysis identified 1,277 unique SNPs that display allele-specific regulatory activity. Among these, susceptibility enhancers near important cholesterol homeostasis genes (APOB, APOC1, APOE, and LIPA) were identified, suggesting that altered gene regulatory activity could represent another way by which genetic variation regulates serum lipoprotein levels. Using CRISPR-based perturbation, we demonstrate how the deletion/activation of a single enhancer leads to changes in the expression of many target genes located in a shared chromatin interaction domain. Our integrative genomics approach represents a comprehensive effort in identifying putative causal regulatory regions and target genes that could predispose to clinical manifestation of CAD by affecting liver function.

HeRA: an atlas of enhancer RNAs across human tissues

Enhancer RNA (eRNA) is a type of long non-coding RNA transcribed from DNA enhancer regions. Despite critical roles of eRNA in gene regulation, the expression landscape of eRNAs in normal human tissue remains unexplored. Using numerous samples from the Genotype-Tissue Expression project, we characterized 45 411 detectable eRNAs and identified tens of thousands of associations between eRNAs and traits, including gender, race, and age. We constructed a co-expression network to identify millions of putative eRNA regulators and target genes across different tissues. We further constructed a user-friendly data portal, Human enhancer RNA Atlas (HeRA, https://hanlab.uth.edu/HeRA/). In HeRA, users can search, browse, and download the eRNA expression profile, trait-related eRNAs, and eRNA co-expression network by searching the eRNA ID, gene symbol, and genomic region in one or multiple tissues. HeRA is the first data portal to characterize eRNAs from 9577 samples across 54 human tissues and facilitates functional and mechanistic investigations of eRNAs.

Emerging roles and potential clinical applications of noncoding RNAs in hepatocellular carcinoma

Hepatocellular carcinoma(HCC) is one of the most common forms of cancer, and accounts for a high proportion of cancer-associated deaths. Growing evidences have demonstrated that non- protein-coding regions of the genome could give rise to transcripts, termed noncoding RNA (ncRNA), that form novel functional layers of the cellular activity. ncRNAs are implicated in different molecular mechanisms and functions at transcriptional, translational and post-translational levels. An increasing number of studies have demonstrated a complex array of molecular and cellular functions of ncRNAs in different stages of the HCC tumorigenesis, either in an oncogenic or tumor-suppressive manner. As a result, several pre-clinical studies have highlighted the great potentials of ncRNAs as novel biomarkers for cancer diagnosis or therapeutics in targeting HCC progression. In this review, we briefly described the characteristics of several representative ncRNAs and summarized the latest findings of their roles and mechanisms in the development of HCC, in order to better understand the cancer biology and their potential clinical applications in this malignancy.

Emerging risk factors for nonalcoholic fatty liver disease associated hepatocellular carcinoma

Worldwide, nonalcoholic fatty liver disease (NAFLD) has reached epidemic proportions and in parallel, hepatocellular carcinoma (HCC) has become one of the fastest growing cancers. Epidemiological studies have not only shed light on the prevalence and incidence of the disease but have also unmasked important environmental risk factors, including the role of diabetes and dyslipidemia in disease pathogenesis. Genetic association studies have identified single nucleotide polymorphisms implicated in NAFLD-HCC, many of which are part of lipid metabolism pathways. Through these clinical studies and subsequently, translational and basic research, the role of statins as a chemoprotective agent has also emerged with ongoing clinical trials assessing their utility in HCC prevention and treatment. In this review, we summarize the recent epidemiological studies describing the burden of NAFLD-HCC in different patient populations and countries. We discuss the genetic and environmental risk factors for NAFLD-HCC and highlight the chemoprotective role of statins and aspirin. We also summarize what is known about NAFLD-HCC in the cirrhosis and non-cirrhosis populations and briefly address the role of surveillance in NAFLD-HCC patients.

The Epigenetic Drug Discovery Landscape for Metabolic-associated Fatty Liver Disease

Despite decades of research, effective therapies for metabolic (dysfunction)-associated fatty liver disease (MAFLD) are lacking. An increasing body of evidence suggests that epigenetic dysregulation is frequent in MAFLD, and orchestrates many aspects of its development and progression. Furthermore, the high plasticity of epigenetic modifications in response to environmental cues renders epigenetics a novel area for therapeutic drug discovery. Over recent years, several epigenetics-based drugs and diagnostic biomarkers have entered clinical development and/or obtained regulatory approval. Here, we review recent advances in our understanding of epigenetic regulation and programming during MAFLD, including DNA methylation, histone modifications, chromatin remodelling, transcriptional control, and noncoding (nc)RNAs. We also discuss the potential translational implications and challenges of epigenetics in the context of MAFLD.