Next Article in Journal
SV40 miR-S1 and Cellular miR-1266 Sequester Each Other from Their Targets, Enhancing Telomerase Activity and Viral Replication
Next Article in Special Issue
Computational Methods to Study DNA:DNA:RNA Triplex Formation by lncRNAs
Previous Article in Journal
miR-218: A Stress-Responsive Epigenetic Modifier
Previous Article in Special Issue
DisiMiR: Predicting Pathogenic miRNAs Using Network Influence and miRNA Conservation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
ncRNA
Volume 8
Issue 4
10.3390/ncrna8040056
ncrna-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease

by
Mirolyuba Ilieva
1,*,
James Dao
2,
Henry E. Miller
2,3,4,
Jens Hedelund Madsen
1,
Alexander J. R. Bishop
3,4,5,
Sakari Kauppinen
1 and
Shizuka Uchida
1,2,*
1
Center for RNA Medicine, Department of Clinical Medicine, Aalborg University, DK-2450 Copenhagen SV, Denmark
2
Bioinformatics Research Network, Atlanta, GA 30317, USA
3
Department of Cell Systems and Anatomy, UT Health San Antonio, San Antonio, TX 78229, USA
4
Greehey Children’s Cancer Research Institute, UT Health San Antonio, San Antonio, TX 78229, USA
5
May’s Cancer Center, UT Health San Antonio, San Antonio, TX 78229, USA
*
Authors to whom correspondence should be addressed.
Non-Coding RNA 2022, 8(4), 56; https://doi.org/10.3390/ncrna8040056
Submission received: 11 July 2022 / Revised: 18 July 2022 / Accepted: 19 July 2022 / Published: 22 July 2022
(This article belongs to the Special Issue Methods and Tools in RNA Biology)
Browse Figures
Review Reports Versions Notes

Abstract

:
The largest solid organ in humans, the liver, performs a variety of functions to sustain life. When damaged, cells in the liver can regenerate themselves to maintain normal liver physiology. However, some damage is beyond repair, which necessitates liver transplantation. Increasing rates of obesity, Western diets (i.e., rich in processed carbohydrates and saturated fats), and cardiometabolic diseases are interlinked to liver diseases, including non-alcoholic fatty liver disease (NAFLD), which is a collective term to describe the excess accumulation of fat in the liver of people who drink little to no alcohol. Alarmingly, the prevalence of NAFLD extends to 25% of the world population, which calls for the urgent need to understand the disease mechanism of NAFLD. Here, we performed secondary analyses of published RNA sequencing (RNA-seq) data of NAFLD patients compared to healthy and obese individuals to identify long non-coding RNAs (lncRNAs) that may underly the disease mechanism of NAFLD. Similar to protein-coding genes, many lncRNAs are dysregulated in NAFLD patients compared to healthy and obese individuals, suggesting that understanding the functions of dysregulated lncRNAs may shed light on the pathology of NAFLD. To demonstrate the functional importance of lncRNAs in the liver, loss-of-function experiments were performed for one NAFLD-related lncRNA, LINC01639, which showed that it is involved in the regulation of genes related to apoptosis, TNF/TGF, cytokine signaling, and growth factors as well as genes upregulated in NAFLD. Since there is no lncRNA database focused on the liver, especially NAFLD, we built a web database, LiverDB, to further facilitate functional and mechanistic studies of hepatic lncRNAs.

1. Introduction

The liver is the largest solid organ in the human body, which performs over 500 essential functions related to digestion, immunity, and metabolism [1]. Liver dysfunction leads to a multitude of health risks and diseases, including liver cancer, cirrhotic liver disease, non-alcoholic fatty liver disease, and acute liver failure [2,3]. Unlike most other organs (e.g., brain, heart), the liver possesses a regenerative capacity [4]. However, upon prolonged insults (e.g., excessive consumption of alcohol, drugs), the damage to the liver becomes irreparable, necessitating liver transplantation [5]. Globally, approximately two million deaths result from liver diseases annually [6]. Thus, intensive research is conducted in identifying the causes of liver diseases.
It is now well recognized that most of the mammalian genome is transcribed as RNA, yet only a few percentages of these transcripts encode for exons of protein-coding genes [7]. All of those that do not code for proteins are collectively called non-coding RNAs (ncRNAs), which include ribosomal RNAs (rRNAs), transfer RNAs (tRNAs), microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) [8]. Among these ncRNAs, lncRNAs are of particular interest to researchers in different fields because lncRNAs are regarded as missing links for understanding many cellular activities and signaling pathways. Because of their broad definition–any ncRNAs longer than 200 nucleotides–it is speculated that the number of lncRNAs well surpasses that of protein-coding genes [7,9]. Indeed, more and more lncRNAs are being discovered from next-generation sequencing (NGS), especially RNA sequencing (RNA-seq), although only limited numbers of lncRNAs are functionally studied. The amount of RNA-seq data deposited in public databases [e.g., ArrayExpress, Gene Expression Omnibus (GEO), Sequence Read Archive (SRA)] is increasing rapidly [10]. However, most of these RNA-seq data are analyzed only for protein-coding genes and not for lncRNAs. Although this is especially true for the hepatology field, recent studies indicate that the dysregulation of lncRNAs in the liver contribute to liver diseases [11,12].
Among different types of liver diseases, the incidence of non-alcoholic fatty liver disease (NAFLD) is increasing rapidly due to lack of exercise and obesity (which are two major causes of type 2 diabetes, another risk factor for NAFLD) [13,14]. NAFLD (also called “metabolic associated fatty liver disease” (MAFLD) [15]) is a collective term to describe the excess accumulation of fat in the liver of people who drink little to no alcohol [16]. It is divided into two subtypes: simple steatosis, non-alcoholic fatty liver (NAFL); and a more aggressive form, nonalcoholic steatohepatitis (NASH) [17,18]. Sadly, NAFLD occurs in about 25% of the world population [19], which highlights the urgent need to understand the mechanism of this life-threatening disease. To uncover the underlying disease mechanism of NAFLD, high-throughput methods (e.g., microarrays and RNA-seq) have been employed [20,21,22,23,24,25,26,27,28,29,30,31,32]. However, these prior studies are mainly focused on protein-coding genes. Given that lncRNAs are increasingly implicated in a variety of diseases, including liver disease and NAFLD (e.g., GAS5 [33], HOTAIR [34], LINC01260 [35], PVT1 [36]), a secondary analysis of published high-throughput data may reveal novel insights into NAFLD. For this purpose, RNA-seq data of previous studies investigating dysregulated protein-coding genes in NAFLD patients compared to healthy donors are of great interest as unlike microarrays, RNA-seq data contain expression information about lncRNAs.
To address the above point, we conducted secondary analyses of published RNA-seq data of NAFLD patients to provide a comprehensive overview of lncRNA dysregulation in this disease context. Since protein-coding genes were already analyzed and further experimentally confirmed in the original studies, here, we omitted the detailed analysis of protein-coding genes (e.g., gene ontology and pathway analyses). Instead, our analyses focus on the expression profiling of lncRNAs, which were not provided in the original studies. By performing a systematic analysis of RNA-seq data for lncRNAs using the latest gene annotation, we uncover the differential regulation of lncRNAs in various patient samples. To confirm the findings, we conducted loss-of-function experiments of dysregulated lncRNA to supplement the lack of knowledge regarding lncRNAs in the pathogenesis of NAFLD.

2. Results

2.1. Expression Profiling of LncRNAs in NAFLD Patients Compared to Healthy Donors

The primary causes of NAFLD are dyslipidemia (abnormally high amount of lipids in the blood), insulin resistance, obesity, and type 2 diabetes [37]. All these primary causes are intertwined and are linked to other metabolic diseases (e.g., cardiovascular disease, stroke). Not surprisingly, both microarrays and RNA-seq techniques have been employed to profile transcriptomic changes in NAFLD patients compared to healthy patients, including obese individuals without NAFLD [20,21,22,23,24,25,26,27,28,29,30,31,32]. Although these studies uncovered important dysregulated signaling pathways in NAFLD patients, there is an overall lack of detailed profiling of lncRNAs in these studies. To address this lack of information, we first reanalyzed the RNA-seq data from a study which profiled protein-coding genes in patients with either subtype of NAFLD [NAFL (n = 15) or NASH (n = 16)] compared with healthy normal weight (n = 14) and obese individuals (n = 12) using liver biopsy samples [32]. The major findings in this original study are: (1) many differentially expressed protein-coding genes are identified when NAFLD patients are compared to both healthy normal weight and obese individuals; (2) these differentially expressed protein-coding genes are enriched in signaling pathways involved in lipid metabolism, immunity, extracellular matrix, and cell cycle control; and (3) there is a large overlap of differentially expressed protein-coding genes in both subtypes of NAFLD compared to healthy normal weight and obese individuals [32].
When a secondary analysis of this dataset (GEO accession number, GSE126848) was performed with a threshold of 2-fold change and false discovery rate (FDR)-adjusted p-value < 0.05, it was observed that over 1000 genes are differentially expressed in both subtypes of NAFLD compared to the control (i.e., both healthy normal weight and obese individuals) (Figure 1A), whereas there are few genes that are differentially expressed between NAFL and NASH as well as between healthy normal weight and obese individuals, as reported previously [32] (Figure 1B; Supplementary Tables S1–S6). Compared to protein-coding genes, the numbers of differentially expressed lncRNAs are generally lower in all conditions. As stated earlier, the original study [32] provided a detailed analysis of protein-coding genes but not for lncRNAs. Thus, we focused on lncRNAs here. When the differentially expressed lncRNAs are compared between NAFLD patients and the control individuals, 163 up- and 12 down-regulated lncRNAs are shared among all comparisons (Figure 1C; Supplementary Tables S7 and S8), suggesting that some differentially expressed lncRNAs are common between two subtypes of NAFLD, as is the case for protein-coding genes [32].

2.2. Dysregulated LncRNAs during the Progression of NASH

More recently, an RNA-seq experiment using a large cohort of European NAFLD patients was conducted [21]. In this original study, the liver biopsy samples from 206 NAFLD patients with different fibrosis stages were compared to those from 10 healthy obese controls (GEO accession number, GSE135251). The NAFLD patients were divided into two groups based on the progression of steatohepatitis (i.e., NASH), which are denoted as early (n = 138) and moderate (n = 68). Using this new data set, we attempted to further narrow down the list of differentially expressed lncRNAs identified in the previous data set of NAFLD patients compared to healthy donors. When the same threshold (2-fold and FDR-adjusted p-value < 0.05) was applied to this new dataset, we observed that several hundred protein-coding and lncRNA genes are up- and down-regulated in NAFLD patients compared to the controls (healthy obese individuals), even in the case of early and moderate NASH patients, where more lncRNAs are down-regulated than up-regulated compared to protein-coding genes (Figure 2A; Supplementary Tables S9–S11). When the previously identified differentially expressed lncRNAs (i.e., 163 up- and 12 down-regulated lncRNAs) were compared within this new dataset, surprisingly, only very few lncRNAs were noted as having the same trend of up- or down-regulation (Figure 2B). This results likely reflects the difficulty of comparing RNA-seq data from different studies and possibly reflects the heterogeneity of biopsy samples [38,39]. To check whether such difficulty also exists in protein-coding genes, the same set of analyses as for lncRNAs was performed, which resulted in an analogous discrepancy of differentially expressed genes between the two studies (Supplementary Figure S1).
Although the overlaps among different studies and conditions are not big, those differentially expressed genes that displayed the same trend (i.e., up-regulation in all conditions and data) should represent the most likely candidates to be validated in other laboratories and studies, regardless of experimental biases. Indeed, among 24 commonly up- and 14 down-regulated protein-coding genes, well-known markers of advanced NAFLD [e.g., aldo-keto reductase family 1 member B10 (AKR1B10) [40,41], cell death inducing DFFA like effector c (CIDEC) [42,43,44], lipoprotein lipase (LPL) [45,46], platelet derived growth factor subunit A (PDGFA) [47], triggering receptor expressed on myeloid cells 2 (TREM2) [48,49,50,51], and thymidylate synthetase (TYMS) [52] are consistently up-regulated. In the case of lncRNAs, there are four commonly up- and two down-regulated lncRNAs in all conditions and data above (Figure 3). The up-regulated lncRNAs are AJ009632.2 (Ensembl Gene ID, ENSG00000229425), long intergenic non-protein coding RNA 1639 (LINC01639; ENSG00000236117), MIRLET7B host gene (MIRLET7BHG; ENSG00000197182), and SNAP25 antisense RNA 1 (SNAP25-AS1; ENSG00000227906). The down-regulated lncRNAs are AC110995.1 (ENSG00000236120) and mitotically associated long non-coding RNA (MANCR; ENSG00000231298).
Among the six differentially expressed lncRNAs, SNAP25-AS1 is shown to be up-regulated in lung cancer [53], whereas the expression of MANCR is up-regulated in breast [54,55], gastric [56], and thyroid cancers [57]. Functionally, MANCR downregulates miR-122a in hepatocellular carcinoma [58], and a recent study shows that bromodomain protein 4 (BRD4), a member of the bromodomain and extra-terminal domain (BET) family of nuclear proteins, binds the super-enhancer region of MANCR locus, which in turn, influences the migration and invasion of prostate cancer cell line PC3 [59]. Mechanistically, MANCR acts as a miRNA sponge to sequester miR-218, which targets RUNX2, in mantle cell lymphoma [60]. In addition, MIRLET7BHG is shown to act as a miRNA sponge to sequester miR-330-5p, which targets the key signal transducer of the hedgehog signaling pathway, smoothened (Smo), in the exosomes of activated hepatic stellate cells [61]. However, none of these differentially expressed lncRNAs are studied in detail within the NAFLD context, which indicates the novelty in this topic.

2.3. Loss-of-Function Experiments to Uncover the Roles of Differentially Expressed LncRNAs

Since only a handful of lncRNAs are functionally characterized in the liver, especially related to NAFLD, we conducted loss-of-function experiments for one commonly up-regulated lncRNA, LINC01639. This lncRNA was chosen for further analysis because of its high expression levels in the human hepatocyte-derived carcinoma cell line, Huh-7, and the presence of only one transcript (thus, no isoforms for this lncRNA gene), according to the latest annotation provided by the Ensembl database.
When silenced in Huh-7 cells by siRNAs (Figure 4A), a differential expression of several genes was observed in LINC01639 silenced cells compared to the control (cells treated with siRNA against scrambled sequence), including significant downregulation in the expression of NAFLD-linked (TYMS, DUSP2, SRF) TNF/TGF, and cytokine signaling (CXCL3, TNFRSF10D) and apoptosis-related genes (BAD), as well as growth factors FGF21 and IGF1 (Figure 4B).
To further understand the impact of silencing LINC01639, we developed a cellular model for NAFLD by treating Huh-7 cells with a mixture of fatty acids (FA) consisting of oleic and palmitic acids (Figure 4C). Interestingly, compared to the control, up-regulation of the cytokine signaling genes CCL2 and IL-6, as well as the Rho GTPase RAC2, was recorded in LINC01639 silenced cells (Figure 4D), suggesting that LINC01639 might function in cytokine signaling.

2.4. A Web Database, LiverDB, for Exploration of NAFLD-Related Genes

To facilitate further understanding of lncRNAs in the liver and NAFLD, we built a web database, LiverDB, to catalog NAFLD-related genes (Figure 5A). LiverDB is an easy-to-use web database that displays expression changes of protein-coding and lncRNA genes in counts per million (CPM), reads per kilobase of transcript per million mapped reads (RPKM), and transcripts per kilobase million (TPM) (Figure 5B). Each study is linked to the data information provided by the GEO [62], and the hyperlink to GeneCards [63] is provided for each gene to obtain the known information related to the target gene. To enable exploration of differentially expressed genes, in the Explore tab of LiverDB, a volcano plot was generated (Figure 5C). In this volcano plot, the gene of interest selected from the Results Table on the left of the screen can be dynamically highlighted to visually inspect the ratio and FDR values of the target gene. These differentially expressed genes can be visually examined further by heatmaps (Figure 5D) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment (Figure 5E). The displayed information should help to understand the genes and signaling pathways being affected by the two conditions being compared. Moreover, UpSet plots can be generated to display the numbers of differentially expressed genes (DEGs) that are found across study/contrast pairs in the Comparison tab (Figure 5F). Compared to a Venn diagram, an Upset plot provides a more efficient way of visualizing the intersections of all comparisons of the dataset, including in LiverDB. To allow further analysis, all data included in LiverDB can be downloaded as text files from the Download Table. In addition, the commands and nextflow [64] pipelines used for the content of LiverDB are available via the GitHub repository to allow for a further analysis of similar RNA-seq data.

3. Discussion

The salient findings of this study are: (1) many genes (both protein-coding and lncRNA genes) are differentially expressed in both subtypes of NAFLD patients compared to the healthy normal weight and obese individuals; (2) the overlap of differentially expressed genes between two independent studies are low, possibly due to the heterogeneity of affected cells in biopsy samples used for RNA-seq experiments; (3) the NAFLD-related lncRNA, LINC01639, is involved in the regulation of genes related to apoptosis, TNF/TGF and cytokine signaling, growth factors, as well as genes upregulated in NAFLD.
In this study, a large number of RNA-seq data were analyzed to interrogate NAFLD-related lncRNAs. Yet, there are several limitations to our study. First, all the RNA-seq data are from the strand-specific mRNA sequencing using the Illumina NextSeq 500 machine. Thus, lncRNAs without poly A tails were not detected in this study, which underestimates the impact of lncRNAs without poly A tails on the pathogenesis of NAFLD. Second, the RNA-seq data analysis pipeline used the annotation provided by the Ensembl database as a GTF file. Thus, only known and registered lncRNAs were analyzed in this study, which omits the novel lncRNAs. Third, all the RNA-seq data used in this study were generated from biopsy samples. Thus, the data variability might be high due to the heterogeneity of NALFD-affected livers [38,39]. To overcome this problem, single-cell RNA-seq (scRNA-seq) might be a solution; yet, only highly expressed lncRNAs without poly A tails can be recorded by the current method of scRNA-seq [65,66]. Fourth, only a molecular analysis of NAFLD-related lncRNA, LINC01639, was provided in this study. Thus, functional and mechanistic studies are needed to understand the contribution of LINC01639 to the pathogenesis of NAFLD.

4. Materials and Methods

4.1. RNA-Seq Data Analysis

The following two datasets were used in this study: (1) GSE126848 [32] and (2) GSE135251 [21]. RNA-Seq data preprocessing pipeline was managed via Nextflow [64] (Version 21.10.6.5660). The raw data from the Sequence Read Archive (SRA) database was downloaded using the prefetch command within the SRA Toolkit (https://trace.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?view=software (accessed on 17 February 2022) (Version 2.10.0). The raw data was then preprocessed by the fastq-dump command within the SRA Toolkit to generate fastq files. The fastq files were then preprocessed by fastp [67] (Version 0.23.2) to generate trimmed fastq files. The trimmed fastq files were aligned by STAR [68] (Version 2.7.10a) to align reads to the reference genome GRCh38.103 in order to generate raw read counts. From the raw read counts, the R packages, edgeR [69] (Version 3.38.1) and GenomicFeatures (Version 1.48.1) [70] were used to generate read count metrics (CPM, RPKM, TPM) and differentially expressed genes with fold change in the logarithmic of base two scale and FDR-adjusted p-values. Pathway analysis on select differentially expressed genes was conducted via Enrichr [71]. All codes are available at the LiverDB GitHub repository (https://github.com/Bishop-Laboratory/LiverDB (accessed on 23 June 2022) in the “preprocess” directory.

4.2. Data analysis and Visualization

Volcano plots were generated using the R-package ggplot2 [72]. Venn diagrams were plotted via http://bioinformatics.psb.ugent.be/webtools/Venn/ (accessed on 2 April 2022). To generate heat maps, MultiExperiment Viewer (MeV) [73] was used. The gene ontology terms were analyzed via the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 [74,75].

4.3. Interactive Web Database, LiverDB

The interactive web database LiverDB was built with R Shiny [76] (Version 1.7.1) and is hosted via shinyapps.io. The R packages used for plotting are ggplot2 [72] (Version 3.3.6) and ComplexHeatmap [77] (Version 2.12.0). For the displayed gene names, Ensembl IDs were mapped to the HUGO Gene Nomenclature Committee (HGNC) gene symbols via the R package biomaRt [78] (Version 2.52.0).

4.4. Cell Culture

Huh-7 cells were cultured in Minimum Essential Medium Eagle (Sigma-Aldrich, St. Louis, MO, USA, #M2279) supplemented with 10% fetal bovine serum (Sigma-Aldrich, #F4135), 1% MEM Non-essential Amino Acid Solution (Sigma-Aldrich, #M7145), 1% L-Glutamine solution (Sigma-Aldrich, #G7513), and 1% Penicillin-Streptomycin (Sigma-Aldrich, #P4333) at 37 °C with 5% CO2.

4.5. Transfection with siRNA and Treatment with Fatty Acids (FA)

The cells at 80% confluence were transfected with 50 nM of siRNA against LINC01639 (sense: CAGACAUAGCAGGAUUUAA[dT][dT]/antisense: UUAAAUCCUGCUAUGUCUG[dT][dT]) using Lipofectamine RNAiMAX transfection reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s protocol. Mission siRNA Universal Negative control #2 (Merck, Kenilworth, NJ, USA, SIC002) served as a control. The samples for RNA purification were collected 48 h after transfection. For development of NAFLD model, 24 h after the transfection, the growth media was exchanged to a starvation media (growth media without FBS) for another 24 h. Then, cells were treated with a mixture of fatty acids (FA) consisting of 2 mM oleic acid (OA; oleic acid-albumin from bovine serum, Merck O3008) and palmitic acid (PA; Sigma) conjugated to Bovine Serum Albumin (BSA) Fraction V, fatty acid free (Millipore, Burlington, MA, USA) in starvation media containing 10% BSA. Cells grown in starvation media containing 10% BSA served as a control. 48 h after treatment with FA cells were harvest for RNA isolation.

4.6. Lipid Accumulation Assay

Cells were seeded at density 1 × 104 in a 96-well plate and grown overnight. After 24 h starvation, cells were treated with an FA mixture (2 mM OA, 2 mM PA, in 10% BSA media) or grown in 10% BSA media alone for 48 h. Wells without cells served as blank controls. Each treatment was performed in triplicates. The lipid accumulation was measured using Hepatic lipid accumulation/Steatosis assay kit (Abcam, Cambridge, UK, ab133131) according to the manufacturer’s protocol. Briefly, cells were fixed for 15 min, washed, and stained with Oil Red O solution for 20 min. The quantification of lipid accumulation was performed after dye extraction for 30 min and absorbance read at 490–520 nm.

4.7. Isolation of Total RNA and RT-PCR

RNeasy Mini Kit (Qiagen, Hilden, Germany, #74104) was used to isolate total RNA from cells and purified following the manufacturer’s protocol. SuperScript IV VILO Master Mix with ezDNase™ Enzyme (Thermo Fisher Scientific, Waltham, MA, USA, #11766500) was used to digest the genomic DNA and reverse transcribe one μg of total RNA for each sample to synthesize the first-strand complementary DNA (cDNA). After reverse transcription, the first-strand cDNA was diluted with DNase/RNase-free water to the concentration of 1 ng/μL. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) reaction was performed with 1 ng of cDNA template per reaction using PowerUp SYBR Green Master Mix (Thermo Fisher Scientific, #A25777) via QuantStudio 6 Flex Real-Time PCR System (Thermo Fisher Scientific) with the annealing temperature at 60 °C. Relative fold expression was calculated by 2-DDCt using ribosomal protein lateral stalk subunit P0 (PRLP0) as an internal control. The primer pairs were designed using Primer3 (http://bioinfo.ut.ee/primer3-0.4.0/ (accessed on 21 March 2022) [79] and in silico validated with the UCSC In-Silico PCR tool (https://genome.ucsc.edu/cgi-bin/hgPcr (accessed on 21 March 2022) before extensive testing by conventional RT-PCR reaction followed by running the PCR product on an agarose gel to examine for a single band of the expected size for each primer pair. The primer sequences are provided in Supplementary Table S12.

4.8. Statistics

Data are presented as mean ± S.E.M. unless otherwise indicated. Two-sample, two-tail, heteroscedastic Student’s t-test was performed to calculate a p-value via Microsoft Excel.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ncrna8040056/s1, Supplementary Table S1: List of differentially expressed genes in obese individuals compared to healthy normal weight individuals; Supplementary Table S2: List of differentially expressed genes in NAFL patients compared to healthy normal weight individuals; Supplementary Table S3: List of differentially expressed genes in NASH patients compared to healthy normal weight individuals; Supplementary Table S4: List of differentially expressed genes in NAFL patients compared to obese individuals; Supplementary Table S5: List of differentially expressed genes in NASH patients compared to obese individuals; Supplementary Table S6: List of differentially expressed genes in NASH patients compared to NAFL patients; Supplementary Table S7: List of shared up-regulated lncRNAs; Supplementary Table S8: List of shared down-regulated lncRNAs; Supplementary Table S9: List of differentially expressed genes in NAFLD patients compared to healthy obese individuals (control); Supplementary Table S10: List of differentially expressed genes in early NASH patients compared to healthy obese individuals (control); Supplementary Table S11: List of differentially expressed genes in moderate NASH patients compared to healthy obese individuals (control); Supplementary Table S12: List of primers used in this study; Supplementary Figure S1: Differentially expressed protein-coding genes.

Author Contributions

Conceptualization, S.U. and M.I.; methodology, S.U., M.I. and H.E.M.; database, H.E.M. and J.D.; validation, M.I. and J.H.M., resources, A.J.R.B., S.K. and S.U.; writing—original draft preparation, S.U. and M.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the Novo Nordisk Foundation (NNF18OC0033438 to S.K., M.I. and S.U.); the National Institutes of Health (1R01CA241554 to A.J.R.B., F31AG072902 to H.E.M.); and Stand Up to Cancer—Cancer Research UK (RT6187 to A.J.R.B.).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Aalborg University (230 337 on 17 March 2017).

Informed Consent Statement

Not applicable as established human cell line was used.

Data Availability Statement

The Supplementary Materials can be found on the GitHub repository: https://github.com/heartlncrna/Analysis_of_NAFLD_Studies (accessed on 23 June 2022). All codes used to generate LiverDB are available on the GitHub repository: https://github.com/Bishop-Laboratory/LiverDB (accessed on 23 June 2022).

Acknowledgments

We would like to thank the members of Bioinformatics Research Network for their continued support to further advance the development of bioinformatics and bioinformatics tools.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zabaleta, N.; Hommel, M.; Salas, D.; Gonzalez-Aseguinolaza, G. Genetic-Based Approaches to Inherited Metabolic Liver Diseases. Hum. Gene Ther. 2019, 30, 1190–1203. [Google Scholar] [CrossRef] [PubMed]
  2. Llovet, J.M.; Montal, R.; Sia, D.; Finn, R.S. Molecular therapies and precision medicine for hepatocellular carcinoma. Nat. Rev. Clin. Oncol. 2018, 15, 599–616. [Google Scholar] [CrossRef] [PubMed]
  3. El-Serag, H.B. Hepatocellular carcinoma. N. Engl. J. Med. 2011, 365, 1118–1127. [Google Scholar] [CrossRef] [PubMed]
  4. Michalopoulos, G.K. Liver regeneration. J. Cell. Physiol. 2007, 213, 286–300. [Google Scholar] [CrossRef] [PubMed]
  5. Trotter, J.F. Current Issues in Liver Transplantation. Gastroenterol. Hepatol. 2016, 12, 214–219. [Google Scholar]
  6. Asrani, S.K.; Devarbhavi, H.; Eaton, J.; Kamath, P.S. Burden of liver diseases in the world. J. Hepatol. 2019, 70, 151–171. [Google Scholar] [CrossRef]
  7. Weirick, T.; Militello, G.; Muller, R.; John, D.; Dimmeler, S.; Uchida, S. The identification and characterization of novel transcripts from RNA-seq data. Brief. Bioinform. 2016, 17, 678–685. [Google Scholar] [CrossRef]
  8. Uchida, S.; Adams, J.C. Physiological roles of non-coding RNAs. Am. J. Physiol. Cell. Physiol. 2019, 317, C1–C2. [Google Scholar] [CrossRef]
  9. Rinn, J.L.; Chang, H.Y. Genome regulation by long noncoding RNAs. Annu. Rev. Biochem. 2012, 81, 145–166. [Google Scholar] [CrossRef] [Green Version]
  10. Katz, K.; Shutov, O.; Lapoint, R.; Kimelman, M.; Brister, J.R.; O’Sullivan, C. The Sequence Read Archive: A decade more of explosive growth. Nucleic Acids Res. 2022, 50, D387–D390. [Google Scholar] [CrossRef]
  11. De Vincentis, A.; Rahmani, Z.; Muley, M.; Vespasiani-Gentilucci, U.; Ruggiero, S.; Zamani, P.; Jamialahmadi, T.; Sahebkar, A. Long noncoding RNAs in nonalcoholic fatty liver disease and liver fibrosis: State-of-the-art and perspectives in diagnosis and treatment. Drug Discov. Today 2020, 25, 1277–1286. [Google Scholar] [CrossRef] [PubMed]
  12. Kim, Y.A.; Park, K.K.; Lee, S.J. LncRNAs Act as a Link between Chronic Liver Disease and Hepatocellular Carcinoma. Int. J. Mol. Sci. 2020, 21, 2883. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pfister, D.; Nunez, N.G.; Pinyol, R.; Govaere, O.; Pinter, M.; Szydlowska, M.; Gupta, R.; Qiu, M.; Deczkowska, A.; Weiner, A.; et al. NASH limits anti-tumour surveillance in immunotherapy-treated HCC. Nature 2021, 592, 450–456. [Google Scholar] [CrossRef] [PubMed]
  14. Parry, S.A.; Hodson, L. Managing NAFLD in Type 2 Diabetes: The Effect of Lifestyle Interventions, a Narrative Review. Adv. Ther. 2020, 37, 1381–1406. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Eslam, M.; Sanyal, A.J.; George, J. MAFLD: A Consensus-Driven Proposed Nomenclature for Metabolic Associated Fatty Liver Disease. Gastroenterology 2020, 158, 1999–2014.e1991. [Google Scholar] [CrossRef]
  16. Byrne, C.D.; Targher, G. NAFLD: A multisystem disease. J. Hepatol. 2015, 62, S47–S64. [Google Scholar] [CrossRef] [Green Version]
  17. Diehl, A.M.; Day, C. Cause, Pathogenesis, and Treatment of Nonalcoholic Steatohepatitis. N. Engl. J Med. 2017, 377, 2063–2072. [Google Scholar] [CrossRef]
  18. Rinella, M.E. Nonalcoholic fatty liver disease: A systematic review. JAMA 2015, 313, 2263–2273. [Google Scholar] [CrossRef]
  19. Younossi, Z.M. Non-alcoholic fatty liver disease—A global public health perspective. J. Hepatol. 2019, 70, 531–544. [Google Scholar] [CrossRef] [Green Version]
  20. Huang, P.; Tu, B.; Liao, H.J.; Huang, F.Z.; Li, Z.Z.; Zhu, K.Y.; Dai, F.; Liu, H.Z.; Zhang, T.Y.; Sun, C.Z. Elevation of plasma tRNA fragments as a promising biomarker for liver fibrosis in nonalcoholic fatty liver disease. Sci. Rep. 2021, 11, 5886. [Google Scholar] [CrossRef]
  21. Govaere, O.; Cockell, S.; Tiniakos, D.; Queen, R.; Younes, R.; Vacca, M.; Alexander, L.; Ravaioli, F.; Palmer, J.; Petta, S.; et al. Transcriptomic profiling across the nonalcoholic fatty liver disease spectrum reveals gene signatures for steatohepatitis and fibrosis. Sci. Transl. Med. 2020, 12, eaba4448. [Google Scholar] [CrossRef] [PubMed]
  22. Luo, Z.; Hu, H.; Liu, S.; Zhang, Z.; Li, Y.; Zhou, L. Comprehensive analysis of the translatome reveals the relationship between the translational and transcriptional control in high fat diet-induced liver steatosis. RNA Biol. 2020, 18, 863–874. [Google Scholar] [CrossRef] [PubMed]
  23. Veyel, D.; Wenger, K.; Broermann, A.; Bretschneider, T.; Luippold, A.H.; Krawczyk, B.; Rist, W.; Simon, E. Biomarker discovery for chronic liver diseases by multi-omics—A preclinical case study. Sci. Rep. 2020, 10, 1314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Haas, J.T.; Vonghia, L.; Mogilenko, D.A.; Verrijken, A.; Molendi-Coste, O.; Fleury, S.; Deprince, A.; Nikitin, A.; Woitrain, E.; Ducrocq-Geoffroy, L.; et al. Transcriptional Network Analysis Implicates Altered Hepatic Immune Function in NASH development and resolution. Nat. Metab. 2019, 1, 604–614. [Google Scholar] [CrossRef] [PubMed]
  25. Jia, X.; Zhai, T. Integrated Analysis of Multiple Microarray Studies to Identify Novel Gene Signatures in Non-alcoholic Fatty Liver Disease. Front. Endocrinol. 2019, 10, 599. [Google Scholar] [CrossRef] [Green Version]
  26. Lefebvre, P.; Lalloyer, F.; Bauge, E.; Pawlak, M.; Gheeraert, C.; Dehondt, H.; Vanhoutte, J.; Woitrain, E.; Hennuyer, N.; Mazuy, C.; et al. Interspecies NASH disease activity whole-genome profiling identifies a fibrogenic role of PPARalpha-regulated dermatopontin. JCI Insight 2017, 2, e92264. [Google Scholar] [CrossRef]
  27. Teufel, A.; Itzel, T.; Erhart, W.; Brosch, M.; Wang, X.Y.; Kim, Y.O.; von Schonfels, W.; Herrmann, A.; Bruckner, S.; Stickel, F.; et al. Comparison of Gene Expression Patterns Between Mouse Models of Nonalcoholic Fatty Liver Disease and Liver Tissues From Patients. Gastroenterology 2016, 151, 513–525.E0. [Google Scholar] [CrossRef]
  28. Arendt, B.M.; Comelli, E.M.; Ma, D.W.; Lou, W.; Teterina, A.; Kim, T.; Fung, S.K.; Wong, D.K.; McGilvray, I.; Fischer, S.E.; et al. Altered hepatic gene expression in nonalcoholic fatty liver disease is associated with lower hepatic n-3 and n-6 polyunsaturated fatty acids. Hepatology 2015, 61, 1565–1578. [Google Scholar] [CrossRef] [Green Version]
  29. Zhao, S.; Fung-Leung, W.P.; Bittner, A.; Ngo, K.; Liu, X. Comparison of RNA-Seq and microarray in transcriptome profiling of activated T cells. PLoS ONE 2014, 9, e78644. [Google Scholar] [CrossRef]
  30. Moylan, C.A.; Pang, H.; Dellinger, A.; Suzuki, A.; Garrett, M.E.; Guy, C.D.; Murphy, S.K.; Ashley-Koch, A.E.; Choi, S.S.; Michelotti, G.A.; et al. Hepatic gene expression profiles differentiate presymptomatic patients with mild versus severe nonalcoholic fatty liver disease. Hepatology 2014, 59, 471–482. [Google Scholar] [CrossRef]
  31. Starmann, J.; Falth, M.; Spindelbock, W.; Lanz, K.L.; Lackner, C.; Zatloukal, K.; Trauner, M.; Sultmann, H. Gene expression profiling unravels cancer-related hepatic molecular signatures in steatohepatitis but not in steatosis. PLoS ONE 2012, 7, e46584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Suppli, M.P.; Rigbolt, K.T.G.; Veidal, S.S.; Heeboll, S.; Eriksen, P.L.; Demant, M.; Bagger, J.I.; Nielsen, J.C.; Oro, D.; Thrane, S.W.; et al. Hepatic transcriptome signatures in patients with varying degrees of nonalcoholic fatty liver disease compared with healthy normal-weight individuals. Am. J. Physiol. Gastrointest. Liver Physiol. 2019, 316, G462–G472. [Google Scholar] [CrossRef] [PubMed]
  33. Cui, J.; Wang, Y.; Xue, H. Long non-coding RNA GAS5 contributes to the progression of nonalcoholic fatty liver disease by targeting the microRNA-29a-3p/NOTCH2 axis. Bioengineered 2022, 13, 8370–8381. [Google Scholar] [CrossRef] [PubMed]
  34. Guo, B.; Cheng, Y.; Yao, L.; Zhang, J.; Lu, J.; Qi, H.; Chen, H. LncRNA HOTAIR regulates the lipid accumulation in non-alcoholic fatty liver disease via miR-130b-3p/ROCK1 axis. Cell. Signal. 2022, 90, 110190. [Google Scholar] [CrossRef]
  35. Ge, X.; Sun, T.; Zhang, Y.; Li, Y.; Gao, P.; Zhang, D.; Zhang, B.; Wang, P.; Ma, W.; Lu, S. The role and possible mechanism of the long noncoding RNA LINC01260 in nonalcoholic fatty liver disease. Nutr. Metab. 2022, 19, 3. [Google Scholar] [CrossRef]
  36. Zhang, H.; Niu, Q.; Liang, K.; Li, X.; Jiang, J.; Bian, C. Effect of LncPVT1/miR-20a-5p on Lipid Metabolism and Insulin Resistance in NAFLD. Diabetes Metab. Syndr. Obes. 2021, 14, 4599–4608. [Google Scholar] [CrossRef]
  37. Kneeman, J.M.; Misdraji, J.; Corey, K.E. Secondary causes of nonalcoholic fatty liver disease. Ther. Adv. Gastroenterol. 2012, 5, 199–207. [Google Scholar] [CrossRef] [Green Version]
  38. Pal, P.; Palui, R.; Ray, S. Heterogeneity of non-alcoholic fatty liver disease: Implications for clinical practice and research activity. World J. Hepatol. 2021, 13, 1584–1610. [Google Scholar] [CrossRef]
  39. Arun, J.; Jhala, N.; Lazenby, A.J.; Clements, R.; Abrams, G.A. Influence of liver biopsy heterogeneity and diagnosis of nonalcoholic steatohepatitis in subjects undergoing gastric bypass. Obes. Surg. 2007, 17, 155–161. [Google Scholar] [CrossRef]
  40. Shi, L.; Guo, S.; Zhang, S.; Gao, X.; Liu, A.; Wang, Q.; Zhang, T.; Zhang, Y.; Wen, A. Glycyrrhetinic acid attenuates disturbed vitamin a metabolism in non-alcoholic fatty liver disease through AKR1B10. Eur. J. Pharmacol. 2020, 883, 173167. [Google Scholar] [CrossRef]
  41. Kanno, M.; Kawaguchi, K.; Honda, M.; Horii, R.; Takatori, H.; Shimakami, T.; Kitamura, K.; Arai, K.; Yamashita, T.; Sakai, Y.; et al. Serum aldo-keto reductase family 1 member B10 predicts advanced liver fibrosis and fatal complications of nonalcoholic steatohepatitis. J. Gastroenterol. 2019, 54, 549–557. [Google Scholar] [CrossRef] [PubMed]
  42. Sans, A.; Bonnafous, S.; Rousseau, D.; Patouraux, S.; Canivet, C.M.; Leclere, P.S.; Tran-Van-Nhieu, J.; Luci, C.; Bailly-Maitre, B.; Xu, X.; et al. The Differential Expression of Cide Family Members is Associated with Nafld Progression from Steatosis to Steatohepatitis. Sci. Rep. 2019, 9, 7501. [Google Scholar] [CrossRef] [PubMed]
  43. Sahini, N.; Borlak, J. Genomics of human fatty liver disease reveal mechanistically linked lipid droplet-associated gene regulations in bland steatosis and nonalcoholic steatohepatitis. Transl. Res. 2016, 177, 41–69. [Google Scholar] [CrossRef] [PubMed]
  44. Tanaka, N.; Takahashi, S.; Matsubara, T.; Jiang, C.; Sakamoto, W.; Chanturiya, T.; Teng, R.; Gavrilova, O.; Gonzalez, F.J. Adipocyte-specific disruption of fat-specific protein 27 causes hepatosteatosis and insulin resistance in high-fat diet-fed mice. J. Biol. Chem. 2015, 290, 3092–3105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bruce, K.D.; Dobrinskikh, E.; Wang, H.; Rudenko, I.; Gao, H.; Libby, A.E.; Gorkhali, S.; Yu, T.; Zsombok, A.; Eckel, R.H. Neuronal Lipoprotein Lipase Deficiency Alters Neuronal Function and Hepatic Metabolism. Metabolites 2020, 10, 385. [Google Scholar] [CrossRef]
  46. Pardina, E.; Baena-Fustegueras, J.A.; Llamas, R.; Catalan, R.; Galard, R.; Lecube, A.; Fort, J.M.; Llobera, M.; Allende, H.; Vargas, V.; et al. Lipoprotein lipase expression in livers of morbidly obese patients could be responsible for liver steatosis. Obes. Surg. 2009, 19, 608–616. [Google Scholar] [CrossRef]
  47. Abderrahmani, A.; Yengo, L.; Caiazzo, R.; Canouil, M.; Cauchi, S.; Raverdy, V.; Plaisance, V.; Pawlowski, V.; Lobbens, S.; Maillet, J.; et al. Increased Hepatic PDGF-AA Signaling Mediates Liver Insulin Resistance in Obesity-Associated Type 2 Diabetes. Diabetes 2018, 67, 1310–1321. [Google Scholar] [CrossRef] [Green Version]
  48. Hou, J.; Zhang, J.; Cui, P.; Zhou, Y.; Liu, C.; Wu, X.; Ji, Y.; Wang, S.; Cheng, B.; Ye, H.; et al. TREM2 sustains macrophage-hepatocyte metabolic coordination in nonalcoholic fatty liver disease and sepsis. J. Clin. Investig. 2021, 131, e135197. [Google Scholar] [CrossRef]
  49. Seidman, J.S.; Troutman, T.D.; Sakai, M.; Gola, A.; Spann, N.J.; Bennett, H.; Bruni, C.M.; Ouyang, Z.; Li, R.Z.; Sun, X.; et al. Niche-Specific Reprogramming of Epigenetic Landscapes Drives Myeloid Cell Diversity in Nonalcoholic Steatohepatitis. Immunity 2020, 52, 1057–1074.e7. [Google Scholar] [CrossRef]
  50. Xiong, X.; Kuang, H.; Ansari, S.; Liu, T.; Gong, J.; Wang, S.; Zhao, X.Y.; Ji, Y.; Li, C.; Guo, L.; et al. Landscape of Intercellular Crosstalk in Healthy and NASH Liver Revealed by Single-Cell Secretome Gene Analysis. Mol. Cell 2019, 75, 644–660.e5. [Google Scholar] [CrossRef]
  51. Maguire, M.; Bushkofsky, J.R.; Larsen, M.C.; Foong, Y.H.; Tanumihardjo, S.A.; Jefcoate, C.R. Diet-dependent retinoid effects on liver gene expression include stellate and inflammation markers and parallel effects of the nuclear repressor Shp. J. Nutr. Biochem. 2017, 47, 63–74. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Zhang, D.; Liu, J.; Xie, T.; Jiang, Q.; Ding, L.; Zhu, J.; Ye, Q. Oleate acid-stimulated HMMR expression by CEBPalpha is associated with nonalcoholic steatohepatitis and hepatocellular carcinoma. Int. J. Biol. Sci. 2020, 16, 2812–2827. [Google Scholar] [CrossRef] [PubMed]
  53. Cao, Q.; Dong, Z.; Liu, S.; An, G.; Yan, B.; Lei, L. Construction of a metastasis-associated ceRNA network reveals a prognostic signature in lung cancer. Cancer Cell Int. 2020, 20, 208. [Google Scholar] [CrossRef] [PubMed]
  54. Tahmouresi, F.; Razmara, E.; Pakravan, K.; Mossahebi-Mohammadi, M.; Rouhollah, F.; Montazeri, M.; Sarrafzadeh, A.; Fahimi, H.; Babashah, S. Upregulation of the long noncoding RNAs DSCAM-AS1 and MANCR is a potential diagnostic marker for breast carcinoma. Biotechnol. Appl. Biochem. 2020, 68, 1250–1256. [Google Scholar] [CrossRef]
  55. Tracy, K.M.; Tye, C.E.; Ghule, P.N.; Malaby, H.L.H.; Stumpff, J.; Stein, J.L.; Stein, G.S.; Lian, J.B. Mitotically-Associated lncRNA (MANCR) Affects Genomic Stability and Cell Division in Aggressive Breast Cancer. Mol. Cancer Res. 2018, 16, 587–598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Yao, L.; Yan, J.; Gan, L.; Huang, S.; Cheng, F.; Fang, N. Upregulation of MANCR predicts poor survival in patients with gastric cancer. Oncol. Lett. 2019, 18, 6801–6806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Huang, N.S.; Lei, B.W.; Tan, L.C.; Yu, P.C.; Shi, X.; Wang, Y.; Ji, Q.H.; Wei, W.J.; Lu, Z.W.; Wang, Y.L. Mitotically associated long non-coding RNA is a tumor promoter in anaplastic thyroid cancer. Ann. Transl. Med. 2020, 8, 1226. [Google Scholar] [CrossRef]
  58. Zhang, X.; Li, H.; Mao, M.; Wang, X.; Zheng, J.; Yang, S. Mitotically-associated long non-coding RNA promotes cancer cell proliferation in hepatocellular carcinoma by downregulating miR-122a. Oncol. Lett. 2019, 18, 6237–6242. [Google Scholar] [CrossRef] [Green Version]
  59. Nagasawa, M.; Tomimatsu, K.; Terada, K.; Kondo, K.; Miyazaki, K.; Miyazaki, M.; Motooka, D.; Okuzaki, D.; Yoshida, T.; Kageyama, S.; et al. Long non-coding RNA MANCR is a target of BET bromodomain protein BRD4 and plays a critical role in cellular migration and invasion abilities of prostate cancer. Biochem. Biophys. Res. Commun. 2020, 526, 128–134. [Google Scholar] [CrossRef]
  60. Wen, S.; Zeng, M.; Li, Y.; Hu, X.; Li, S.; Liang, X.; Zhu, L.; Yang, S. Downregulation of MANCR inhibits cancer cell proliferation in mantle cell lymphoma possibly by interacting with RUNX2. Acta Biochim. Biophys. Sin. 2019, 51, 1142–1147. [Google Scholar] [CrossRef]
  61. Xia, Y.; Zhen, L.; Li, H.; Wang, S.; Chen, S.; Wang, C.; Yang, X. MIRLET7BHG promotes hepatocellular carcinoma progression by activating hepatic stellate cells through exosomal SMO to trigger Hedgehog pathway. Cell Death Dis. 2021, 12, 326. [Google Scholar] [CrossRef] [PubMed]
  62. Barrett, T.; Wilhite, S.E.; Ledoux, P.; Evangelista, C.; Kim, I.F.; Tomashevsky, M.; Marshall, K.A.; Phillippy, K.H.; Sherman, P.M.; Holko, M.; et al. NCBI GEO: Archive for functional genomics data sets—Update. Nucleic Acids Res. 2013, 41, D991–D995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinform. 2016, 54, 1–30. [Google Scholar] [CrossRef]
  64. Di Tommaso, P.; Chatzou, M.; Floden, E.W.; Barja, P.P.; Palumbo, E.; Notredame, C. Nextflow enables reproducible computational workflows. Nat. Biotechnol. 2017, 35, 316–319. [Google Scholar] [CrossRef] [PubMed]
  65. Chen, G.; Ning, B.; Shi, T. Single-Cell RNA-Seq Technologies and Related Computational Data Analysis. Front. Genet. 2019, 10, 317. [Google Scholar] [CrossRef] [PubMed]
  66. Lahnemann, D.; Koster, J.; Szczurek, E.; McCarthy, D.J.; Hicks, S.C.; Robinson, M.D.; Vallejos, C.A.; Campbell, K.R.; Beerenwinkel, N.; Mahfouz, A.; et al. Eleven grand challenges in single-cell data science. Genome Biol. 2020, 21, 31. [Google Scholar] [CrossRef] [PubMed]
  67. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  68. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef]
  69. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  70. Lawrence, M.; Huber, W.; Pages, H.; Aboyoun, P.; Carlson, M.; Gentleman, R.; Morgan, M.T.; Carey, V.J. Software for computing and annotating genomic ranges. PLoS Comput. Biol. 2013, 9, e1003118. [Google Scholar] [CrossRef]
  71. Xie, Z.; Bailey, A.; Kuleshov, M.V.; Clarke, D.J.B.; Evangelista, J.E.; Jenkins, S.L.; Lachmann, A.; Wojciechowicz, M.L.; Kropiwnicki, E.; Jagodnik, K.M.; et al. Gene Set Knowledge Discovery with Enrichr. Curr. Protoc. 2021, 1, e90. [Google Scholar] [CrossRef] [PubMed]
  72. Wickham, H. ggplot2: Elegant Graphics for Data Analysis. In Use R! 2nd ed.; Springer International Publishing: Cham, Switzerland, 2016. [Google Scholar] [CrossRef]
  73. Howe, E.A.; Sinha, R.; Schlauch, D.; Quackenbush, J. RNA-Seq analysis in MeV. Bioinformatics 2011, 27, 3209–3210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  75. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [Green Version]
  76. Chang, W.; Cheng, J.; Allaire, J.; Xie, Y.; McPherson, J. Shiny: Web application framework for R. In R Package Version; R Foundation for Statistical Computing: Vienna, Austria, 2017; Volume 1. [Google Scholar]
  77. Gu, Z.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [Green Version]
  78. Durinck, S.; Spellman, P.T.; Birney, E.; Huber, W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat. Protoc. 2009, 4, 1184–1191. [Google Scholar] [CrossRef] [Green Version]
  79. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3—New capabilities and interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef] [Green Version]
Ncrna 08 00056 g001 550
Figure 1. Gene expression profiling of NAFL and NASH patients compared to healthy normal weight and obese individuals. (A) Volcano plots comparing different conditions. The numbers of samples are NAFL (n = 15), NASH (n = 16), healthy normal weight (n = 14), and obese individuals (n = 12). Up-regulated genes are colored in red, and down-regulated genes are in green. (B) Tables of differentially expressed genes. The Others category includes any genes other than protein-coding or lncRNAs based on the biotypes provided by the Ensembl database. (C) Venn diagrams of up- and down-regulated lncRNAs.
Figure 1. Gene expression profiling of NAFL and NASH patients compared to healthy normal weight and obese individuals. (A) Volcano plots comparing different conditions. The numbers of samples are NAFL (n = 15), NASH (n = 16), healthy normal weight (n = 14), and obese individuals (n = 12). Up-regulated genes are colored in red, and down-regulated genes are in green. (B) Tables of differentially expressed genes. The Others category includes any genes other than protein-coding or lncRNAs based on the biotypes provided by the Ensembl database. (C) Venn diagrams of up- and down-regulated lncRNAs.
Ncrna 08 00056 g001
Ncrna 08 00056 g002 550
Figure 2. Gene expression changes in different stages of NASH compared to obese individuals. (A) Tables of differentially expressed genes. The Others category includes any genes other than protein-coding or lncRNAs based on the biotypes provided by the Ensembl database. (B) Volcano plot of early (n = 138) and moderate NASH patients (n = 68) compared to the control (obese individuals; n = 10). Only 163 up- and 12 down-regulated lncRNAs identified in the previous dataset (GEO accession number, GSE126848) are shown in each volcano plot, as indicated in the figure. Up-regulated genes are colored in red, and down-regulated genes are in green.
Figure 2. Gene expression changes in different stages of NASH compared to obese individuals. (A) Tables of differentially expressed genes. The Others category includes any genes other than protein-coding or lncRNAs based on the biotypes provided by the Ensembl database. (B) Volcano plot of early (n = 138) and moderate NASH patients (n = 68) compared to the control (obese individuals; n = 10). Only 163 up- and 12 down-regulated lncRNAs identified in the previous dataset (GEO accession number, GSE126848) are shown in each volcano plot, as indicated in the figure. Up-regulated genes are colored in red, and down-regulated genes are in green.
Ncrna 08 00056 g002
Ncrna 08 00056 g003 550
Figure 3. Expression profiles of commonly up- or down-regulated lncRNAs. Violin plots show counts per million (CPM) values of each sample group for each dataset, as indicated in the figure. The number of samples for each category is indicated in the figure.
Figure 3. Expression profiles of commonly up- or down-regulated lncRNAs. Violin plots show counts per million (CPM) values of each sample group for each dataset, as indicated in the figure. The number of samples for each category is indicated in the figure.
Ncrna 08 00056 g003
Ncrna 08 00056 g004 550
Figure 4. Role of LINC01639 in regulation of gene expression in hepatocytes. (A) Silencing of LINC01639 in Huh-7 cells, Error bars represent mean ± S.E.M. *** p < 0.001. (B) Gene expression profile after silencing of LINC01639 in Huh-7 cells. n = three biological replicates. * p < 0.05; ** p < 0.01. (C) NAFLD cellular model. Huh-7 cells treated with a mixture of fatty acids (FA) consisting of 2 mM oleic acid and 2 mM palmitic acid at the ratio of 1:1. The schematic experimental timeline for NAFLD cellular model generation is presented. (D) Gene expression profile of NAFLD cellular model upon silencing of LINC01639 n = three biological replicates. * p < 0.05; ** p < 0.01. nd = not detected.
Figure 4. Role of LINC01639 in regulation of gene expression in hepatocytes. (A) Silencing of LINC01639 in Huh-7 cells, Error bars represent mean ± S.E.M. *** p < 0.001. (B) Gene expression profile after silencing of LINC01639 in Huh-7 cells. n = three biological replicates. * p < 0.05; ** p < 0.01. (C) NAFLD cellular model. Huh-7 cells treated with a mixture of fatty acids (FA) consisting of 2 mM oleic acid and 2 mM palmitic acid at the ratio of 1:1. The schematic experimental timeline for NAFLD cellular model generation is presented. (D) Gene expression profile of NAFLD cellular model upon silencing of LINC01639 n = three biological replicates. * p < 0.05; ** p < 0.01. nd = not detected.
Ncrna 08 00056 g004
Ncrna 08 00056 g005 550
Figure 5. LiverDB. (A) The top page of LiverDB. (B) An interactive box plot showing the distribution of gene expression values across samples within each biological condition in the selected study. (C) A volcano plot of selected comparative conditions. (D) A heatmap of differentially expressed genes. (E) A pathway enrichment plot displaying the top results from KEGG pathway analysis as a heatmap. (F) An UpSet plot showing the number of differentially expressed genes that are found across study/contrast pairs.
Figure 5. LiverDB. (A) The top page of LiverDB. (B) An interactive box plot showing the distribution of gene expression values across samples within each biological condition in the selected study. (C) A volcano plot of selected comparative conditions. (D) A heatmap of differentially expressed genes. (E) A pathway enrichment plot displaying the top results from KEGG pathway analysis as a heatmap. (F) An UpSet plot showing the number of differentially expressed genes that are found across study/contrast pairs.
Ncrna 08 00056 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ilieva, M.; Dao, J.; Miller, H.E.; Madsen, J.H.; Bishop, A.J.R.; Kauppinen, S.; Uchida, S. Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease. Non-Coding RNA 2022, 8, 56. https://doi.org/10.3390/ncrna8040056

AMA Style

Ilieva M, Dao J, Miller HE, Madsen JH, Bishop AJR, Kauppinen S, Uchida S. Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease. Non-Coding RNA. 2022; 8(4):56. https://doi.org/10.3390/ncrna8040056

Chicago/Turabian Style

Ilieva, Mirolyuba, James Dao, Henry E. Miller, Jens Hedelund Madsen, Alexander J. R. Bishop, Sakari Kauppinen, and Shizuka Uchida. 2022. "Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease" Non-Coding RNA 8, no. 4: 56. https://doi.org/10.3390/ncrna8040056

APA Style

Ilieva, M., Dao, J., Miller, H. E., Madsen, J. H., Bishop, A. J. R., Kauppinen, S., & Uchida, S. (2022). Systematic Analysis of Long Non-Coding RNA Genes in Nonalcoholic Fatty Liver Disease. Non-Coding RNA, 8(4), 56. https://doi.org/10.3390/ncrna8040056

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop