The mechanisms that drive nonalcoholic fatty liver disease (NAFLD) remain incompletely understood. This large multicenter study characterized the transcriptional changes that occur in liver tissue across the NAFLD spectrum as disease progresses to cirrhosis to identify potential circulating markers. We performed high-throughput RNA sequencing on a discovery cohort comprising histologically characterized NAFLD samples from 206 patients. Unsupervised clustering stratified NAFLD on the basis of disease activity and fibrosis stage with differences in age, aspartate aminotransferase (AST), type 2 diabetes mellitus, and carriage of PNPLA3 rs738409, a genetic variant associated with NAFLD. Relative to early disease, we consistently identified 25 differentially expressed genes as fibrosing steatohepatitis progressed through stages F2 to F4. This 25-gene signature was independently validated by logistic modeling in a separate replication cohort (n = 175), and an integrative analysis with publicly available single-cell RNA sequencing data elucidated the likely relative contribution of specific intrahepatic cell populations. Translating these findings to the protein level, SomaScan analysis in more than 300 NAFLD serum samples confirmed that circulating concentrations of proteins AKR1B10 and GDF15 were strongly associated with disease activity and fibrosis stage. Supporting the biological plausibility of these data, in vitro functional studies determined that endoplasmic reticulum stress up-regulated expression of AKR1B10, GDF15, and PDGFA, whereas GDF15 supplementation tempered the inflammatory response in macrophages upon lipid loading and lipopolysaccharide stimulation. This study provides insights into the pathophysiology of progressive fibrosing steatohepatitis, and proof of principle that transcriptomic changes represent potentially tractable and clinically relevant markers of disease progression.

Genetic factors associated with non-alcoholic fatty liver disease (NAFLD) remain incompletely understood. To date, most genome-wide association studies (GWASs) have adopted radiologically assessed hepatic triglyceride content as the reference phenotype and so cannot address steatohepatitis or fibrosis. We describe a GWAS encompassing the full spectrum of histologically characterised NAFLD. The GWAS involved 1,483 European NAFLD cases and 17,781 genetically matched controls. A replication cohort of 559 NAFLD cases and 945 controls was genotyped to confirm signals showing genome-wide or close to genome-wide significance. Case-control analysis identified signals showing p values ≤5 × 10-8 at 4 locations (chromosome [chr] 2 GCKR/C2ORF16; chr4 HSD17B13; chr19 TM6SF2; chr22 PNPLA3) together with 2 other signals with p <1 × 10-7 (chr1 near LEPR and chr8 near IDO2/TC1). Case-only analysis of quantitative traits showed that the PNPLA3 signal (rs738409) had genome-wide significance for steatosis, fibrosis and NAFLD activity score and a new signal (PYGO1 rs62021874) had close to genome-wide significance for steatosis (p = 8.2 × 10-8). Subgroup case-control analysis for NASH confirmed the PNPLA3 signal. The chr1 LEPR single nucleotide polymorphism also showed genome-wide significance for this phenotype. Considering the subgroup with advanced fibrosis (≥F3), the signals on chr2, chr19 and chr22 maintained their genome-wide significance. Except for GCKR/C2ORF16, the genome-wide significance signals were replicated. This study confirms PNPLA3 as a risk factor for the full histological spectrum of NAFLD at genome-wide significance levels, with important contributions from TM6SF2 and HSD17B13. PYGO1 is a novel steatosis modifier, suggesting that Wnt signalling pathways may be relevant in NAFLD pathogenesis. Non-alcoholic fatty liver disease is a common disease where excessive fat accumulates in the liver and may result in cirrhosis. To understand who is at risk of developing this disease and suffering liver damage, we undertook a genetic study to compare the genetic profiles of people suffering from fatty liver disease with genetic profiles seen in the general population. We found that particular sequences in 4 different areas of the human genome were seen at different frequencies in the fatty liver disease cases. These sequences may help predict an individual's risk of developing advanced disease. Some genes where these sequences are located may also be good targets for future drug treatments.

Non-alcoholic fatty liver disease (NAFLD) is a common, progressive liver disease strongly associated with the metabolic syndrome. It is unclear how progression of NAFLD towards cirrhosis translates into systematic changes in circulating proteins. Here, we provide a detailed proteo-transcriptomic map of steatohepatitis and fibrosis during progressive NAFLD. In this multicentre proteomic study, we characterize 4,730 circulating proteins in 306 patients with histologically characterized NAFLD and integrate this with transcriptomic analysis in paired liver tissue. We identify circulating proteomic signatures for active steatohepatitis and advanced fibrosis, and correlate these with hepatic transcriptomics to develop a proteo-transcriptomic signature of 31 markers. Deconvolution of this signature by single-cell RNA sequencing reveals the hepatic cell types likely to contribute to proteomic changes with disease progression. As an exemplar of use as a non-invasive diagnostic, logistic regression establishes a composite model comprising four proteins (ADAMTSL2, AKR1B10, CFHR4 and TREM2), body mass index and type 2 diabetes mellitus status, to identify at-risk steatohepatitis.

The rapid improvements in single cell sequencing technologies and analyses afford greater scope for dissecting organoid cultures composed of multiple cell types and create an opportunity to interrogate these models to understand tissue biology, cellular behavior and interactions. To this end, retinal organoids generated from human embryonic stem cells (hESCs) were analyzed by single cell RNA-sequencing (scRNA-Seq) at three time points of differentiation. Combinatorial data from all time points revealed the presence of nine clusters, five of which corresponded to key retinal cell types: retinal pigment epithelium (RPE), retinal ganglion cells (RGCs), cone and rod photoreceptors, and Müller glia. The remaining four clusters expressed genes typical of mitotic cells, extracellular matrix components and those involved in homeostasis. The cell clustering analysis revealed the decreasing presence of mitotic cells and RGCs, formation of a distinct RPE cluster, the emergence of cone and rod photoreceptors from photoreceptor precursors, and an increasing number of Müller glia cells over time. Pseudo-time analysis resembled the order of cell birth during retinal development, with the mitotic cluster commencing the trajectory and the large majority of Müller glia completing the time line. Together, these data demonstrate the feasibility and potential of scRNA-Seq to dissect the inherent complexity of retinal organoids and the orderly birth of key retinal cell types. Stem Cells 2019;37:593–598

RNA sequencing in situ allows for whole-transcriptome characterization at high resolution, while retaining spatial information. These data present an analytical challenge for bioinformatics-how to leverage spatial information effectively? Properties of data with a spatial dimension require special handling, which necessitate a different set of statistical and inferential considerations when compared to non-spatial data. The geographical sciences primarily use spatial data and have developed methods to analye them. Here we discuss the challenges associated with spatial analysis and examine how we can take advantage of practice from the geographical sciences to realize the full potential of spatial information in transcriptomic datasets.

The progression of non-alcoholic steatohepatitis (NASH) to fibrosis and hepatocellular carcinoma (HCC) is aggravated by auto-aggressive T cells. The gut-liver axis contributes to NASH, but the mechanisms involved and the consequences for NASH-induced fibrosis and liver cancer remain unknown. We investigated the role of gastrointestinal B cells in the development of NASH, fibrosis and NASH-induced HCC. C57BL/6J wild-type (WT), B cell-deficient and different immunoglobulin-deficient or transgenic mice were fed distinct NASH-inducing diets or standard chow for 6 or 12 months, whereafter NASH, fibrosis, and NASH-induced HCC were assessed and analysed. Specific pathogen-free/germ-free WT and μMT mice (containing B cells only in the gastrointestinal tract) were fed a choline-deficient high-fat diet, and treated with an anti-CD20 antibody, whereafter NASH and fibrosis were assessed. Tissue biopsy samples from patients with simple steatosis, NASH and cirrhosis were analysed to correlate the secretion of immunoglobulins to clinicopathological features. Flow cytometry, immunohistochemistry and single-cell RNA-sequencing analysis were performed in liver and gastrointestinal tissue to characterise immune cells in mice and humans. Activated intestinal B cells were increased in mouse and human NASH samples and licensed metabolic T-cell activation to induce NASH independently of antigen specificity and gut microbiota. Genetic or therapeutic depletion of systemic or gastrointestinal B cells prevented or reverted NASH and liver fibrosis. IgA secretion was necessary for fibrosis induction by activating CD11b+CCR2+F4/80+CD11c-FCGR1+ hepatic myeloid cells through an IgA-FcR signalling axis. Similarly, patients with NASH had increased numbers of activated intestinal B cells; additionally, we observed a positive correlation between IgA levels and activated FcRg+ hepatic myeloid cells, as well the extent of liver fibrosis. Intestinal B cells and the IgA-FcR signalling axis represent potential therapeutic targets for the treatment of NASH. There is currently no effective treatment for non-alcoholic steatohepatitis (NASH), which is associated with a substantial healthcare burden and is a growing risk factor for hepatocellular carcinoma (HCC). We have previously shown that NASH is an auto-aggressive condition aggravated, amongst others, by T cells. Therefore, we hypothesized that B cells might have a role in disease induction and progression. Our present work highlights that B cells have a dual role in NASH pathogenesis, being implicated in the activation of auto-aggressive T cells and the development of fibrosis via activation of monocyte-derived macrophages by secreted immunoglobulins (e.g., IgA). Furthermore, we show that the absence of B cells prevented HCC development. B cell-intrinsic signalling pathways, secreted immunoglobulins, and interactions of B cells with other immune cells are potential targets for combinatorial NASH therapies against inflammation and fibrosis.