A time-resolved multi-omic atlas of the developing mouse liver

ESRP2-microRNA-122 axis promotes the postnatal onset of liver polyploidization and maturation

Hepatocyte polyploidy and maturity are critical to acquiring specialized liver functions. Multiple intracellular and extracellular factors influence ploidy, but how they cooperate temporally to steer liver polyploidization and maturation or how post-transcriptional mechanisms integrate into these paradigms is unknown. Here, we identified an important regulatory hierarchy in which postnatal activation of epithelial splicing regulatory protein 2 (ESRP2) stimulates processing of liver-specific microRNA (miR-122) to facilitate polyploidization, maturation, and functional competence of hepatocytes. By determining transcriptome-wide protein-RNA interactions in vivo and integrating them with single-cell and bulk hepatocyte RNA-seq data sets, we delineated an ESRP2-driven RNA processing program that drives sequential replacement of fetal-to-adult transcript isoforms. Specifically, ESRP2 binds the primary miR-122 host gene transcript to promote its processing/biogenesis. Combining constitutive and inducible ESRP2 gain- and loss-of-function mouse models with miR-122 rescue experiments, we demonstrated that timed activation of ESRP2 augments the miR-122-driven program of cytokinesis failure, ensuring the proper onset and extent of hepatocyte polyploidization.

The temporal protein signature analyses of developing human deciduous molar tooth germ

The tooth serves as an exemplary model for developmental studies, encompassing epithelial-mesenchymal transition and cell differentiation. The essential factors and pathways identified in tooth development will help understand the natural development process and the malformations of mineralized tissues such as skeleton. The time-dependent proteomic changes were investigated through the proteomics of healthy human molars during embryonic stages, ranging from the cap-to-early bell stage. A comprehensive analysis revealed 713 differentially expressed proteins (DEPs) exhibiting five distinct temporal expression patterns. Through the application of weighted gene co-expression network analysis (WGCNA), 24 potential driver proteins of tooth development were screened, including CHID1, RAP1GDS1, HAPLN3, AKAP12, WLS, GSS, DDAH1, CLSTN1, AFM, RBP1, AGO1, SET, HMGB2, HMGB1, ANP32A, SPON1, FREM1, C8B, PRPS2, FCHO2, PPP1R12A, GPALPP1, U2AF2, and RCC2. Then, the proteomics and transcriptomics expression patterns of these proteins were further compared, complemented by single-cell RNA-sequencing (scRNA-seq). In summary, this study not only offers a wealth of information regarding the molecular intricacies of human embryonic epithelial and mesenchymal cell differentiation but also serves as an invaluable resource for future mechanistic inquiries into tooth development.

Systemic delivery of AAV-GCDH ameliorates HLD-induced phenotype in a glutaric aciduria type I mouse model

Glutaric aciduria type 1 (GA1) is a rare inherited metabolic disorder caused by a deficiency of glutaryl-coenzyme A dehydrogenase (GCDH), with accumulation of neurotoxic metabolites, resulting in a complex movement disorder, irreversible brain damage, and premature death in untreated individuals. While early diagnosis and a lysine restricted diet can extend survival, they do not prevent neurological damage in approximately one-third of treated patients, and more effective therapies are required. Here we report the efficacy of adeno-associated virus 9 (AAV9)-mediated systemic delivery of human GCDH at preventing a high lysine diet (HLD)-induced phenotype in Gcdh -/- mice. Neonatal treatment with AAV-GCDH restores GCDH expression and enzyme activity in liver and striatum. This treatment protects the mice from HLD-aggressive phenotype with all mice surviving this exposure; in stark contrast, a lack of treatment on an HLD triggers very high accumulation of glutaric acid, 3-hydroxyglutaric acid, and glutarylcarnitine in tissues, with about 60% death due to brain accumulation of toxic lysine metabolites. AAV-GCDH significantly ameliorates the striatal neuropathology, minimizing neuronal dysfunction, gliosis, and alterations in myelination. Magnetic resonance imaging findings show protection against striatal injury. Altogether, these results provide preclinical evidence to support AAV-GCDH gene therapy for GA1.

The mouse multi-organ proteome from infancy to adulthood

The early-life organ development and maturation shape the fundamental blueprint for later-life phenotype. However, a multi-organ proteome atlas from infancy to adulthood is currently not available. Herein, we present a comprehensive proteomic analysis of ten mouse organs (brain, heart, lung, liver, kidney, spleen, stomach, intestine, muscle and skin) at three crucial developmental stages (1-, 4- and 8-weeks after birth) acquired using data-independent acquisition mass spectrometry. We detect and quantify 11,533 protein groups across the ten organs and obtain 115 age-related differentially expressed protein groups that are co-expressed in all organs from infancy to adulthood. We find that spliceosome proteins prevalently play crucial regulatory roles in the early-life development of multiple organs, and detect organ-specific expression patterns and sexual dimorphism. This multi-organ proteome atlas provides a fundamental resource for understanding the molecular mechanisms underlying early-life organ development and maturation.

ESRP2-microRNA-122 axis directs the postnatal onset of liver polyploidization and maturation

Hepatocyte polyploidy and maturity are critical to acquiring specialized liver functions. Multiple intra- and extracellular factors influence ploidy, but how they cooperate temporally to steer liver polyploidization and maturation or how post-transcriptional mechanisms integrate into these paradigms is unknown. Here, we identified an important regulatory hierarchy in which postnatal activation of Epithelial-Splicing-Regulatory-Protein-2 (ESRP2) stimulates biogenesis of liver-specific microRNA (miR-122), thereby facilitating polyploidization, maturation, and functional competence of hepatocytes. By determining transcriptome-wide protein-RNA interactions in vivo and integrating them with single-cell and bulk hepatocyte RNA-seq datasets, we delineate an ESRP2-driven RNA processing program that drives sequential replacement of fetal-to-adult transcript isoforms. Specifically, ESRP2 binds the primary miR-122 host gene transcript to promote its processing/biogenesis. Combining constitutive and inducible ESRP2 gain- and loss-of-function mice models with miR-122 rescue experiments, we demonstrate that timed activation of ESRP2 augments miR-122-driven program of cytokinesis failure, ensuring proper onset and extent of hepatocyte polyploidization.

Generation of in vivo-like multicellular liver organoids by mimicking developmental processes: A review

Liver is involved in metabolic reactions, ammonia detoxification, and immunity. Multicellular liver tissue cultures are more desirable for drug screening, disease modeling, and researching transplantation therapy, than hepatocytes monocultures. Hepatocytes monocultures are not stable for long. Further, hepatocyte-like cells induced from pluripotent stem cells and in vivo hepatocytes are functionally dissimilar. Organoid technology circumvents these issues by generating functional ex vivo liver tissue from intrinsic liver progenitor cells and extrinsic stem cells, including pluripotent stem cells. To function as in vivo liver tissue, the liver organoid cells must be arranged precisely in the 3-dimensional space, closely mimicking in vivo liver tissue. Moreover, for long term functioning, liver organoids must be appropriately vascularized and in contact with neighboring epithelial tissues (e.g., bile canaliculi and intrahepatic bile duct, or intrahepatic and extrahepatic bile ducts). Recent discoveries in liver developmental biology allows one to successfully induce liver component cells and generate organoids. Thus, here, in this review, we summarize the current state of knowledge on liver development with a focus on its application in generating different liver organoids. We also cover the future prospects in creating (functionally and structurally) in vivo-like liver organoids using the current knowledge on liver development.

Endothelial Slc35a1 Deficiency Causes Loss of LSEC Identity and Exacerbates Neonatal Lipid Deposition in the Liver in Mice

BACKGROUND & AIMS

The functional maturation of the liver largely occurs after birth. In the early stages of life, the liver of a newborn encounters enormous high-fat metabolic stress caused by the consumption of breast milk. It is unclear how the maturing liver adapts to high lipid metabolism. Liver sinusoidal endothelial cells (LSECs) play a fundamental role in establishing liver vasculature and are decorated with many glycoproteins on their surface. The Slc35a1 gene encodes a cytidine-5'-monophosphate (CMP)-sialic acid transporter responsible for transporting CMP-sialic acids between the cytoplasm and the Golgi apparatus for protein sialylation. This study aimed to determine whether endothelial sialylation plays a role in hepatic vasculogenesis and functional maturation.

METHODS

Endothelial-specific Slc35a1 knockout mice were generated. Liver tissues were collected for histologic analysis, lipidomic profiling, RNA sequencing, confocal immunofluorescence, and immunoblot analyses.

RESULTS

Endothelial Slc35a1-deficient mice exhibited excessive neonatal hepatic lipid deposition, severe liver damage, and high mortality. Endothelial deletion of Slc35a1 led to sinusoidal capillarization and disrupted hepatic zonation. Mechanistically, vascular endothelial growth factor receptor 2 (VEGFR2) in LSECs was desialylated and VEGFR2 signaling was enhanced in Slc35a1-deficient mice. Inhibition of VEGFR2 signaling by SU5416 alleviated lipid deposition and restored hepatic vasculature in Slc35a1-deficient mice.

CONCLUSIONS

Our findings suggest that sialylation of LSECs is critical for maintaining hepatic vascular development and lipid homeostasis. Targeting VEGFR2 signaling may be a new strategy to prevent liver disorders associated with abnormal vasculature and lipid deposition.

Age-, sex- and proximal-distal-resolved multi-omics identifies regulators of intestinal aging in non-human primates

The incidence of intestinal diseases increases with age, yet the mechanisms governing gut aging and its link to diseases, such as colorectal cancer (CRC), remain elusive. In this study, while considering age, sex and proximal-distal variations, we used a multi-omics approach in non-human primates (Macaca fascicularis) to shed light on the heterogeneity of intestinal aging and identify potential regulators of gut aging. We explored the roles of several regulators, including those from tryptophan metabolism, in intestinal function and lifespan in Caenorhabditis elegans. Suggesting conservation of region specificity, tryptophan metabolism via the kynurenine and serotonin (5-HT) pathways varied between the proximal and distal colon, and, using a mouse colitis model, we observed that distal colitis was more sensitive to 5-HT treatment. Additionally, using proteomics analysis of human CRC samples, we identified links between gut aging and CRC, with high HPX levels predicting poor prognosis in older patients with CRC. Together, this work provides potential targets for preventing gut aging and associated diseases.

Multisite and multitimepoint proteomics reveal that patent foramen ovale closure improves migraine and epilepsy by reducing right-to-left shunt-induced hypoxia

Patent foramen ovale (PFO) is a congenital defect in the partition between two atria, which may cause right-to-left shunt (RLS), leading to neurological chronic diseases with episodic manifestations (NCDEMs), such as migraine and epilepsy. However, whether PFO closure was effective in improving NCDEMs and the mechanism were unclear. Twenty-eight patients with migraine or epilepsy who underwent PFO closure were recruited. Notably, approximately half of patients received 50% or more reduction in seizure or headache attacks. Meanwhile, the postoperative blood oxygen partial pressure and oxygen saturation were elevated after PFO closure. Multisite (peripheral, right, and left atrial) and multitimepoint (before and after surgery) plasma proteomics from patients showed that the levels of free hemoglobin and cell adhesion molecules (CAMs) were significantly increased after PFO closure, which may be related to the relief of the hypoxic state. Furtherly, the omics data from multiple brain regions of mice revealed that a large number of proteins were differentially expressed in the occipital region in response to PFO, including redox molecules and CAMs, suggesting PFO-caused hypoxia may have great impacts on occipital region. Collectively, PFO may cause NCDEMs due to RLS-induced hypoxia, and PFO closure could prevent RLS to improve migraine and epilepsy.

Two Hippo signaling modules orchestrate liver size and tumorigenesis

The Hippo pathway is a central regulator of organ size and tumorigenesis and is commonly depicted as a kinase cascade, with an increasing number of regulatory and adaptor proteins linked to its regulation over recent years. Here, we propose that two Hippo signaling modules, MST1/2-SAV1-WWC1-3 (HPO1) and MAP4K1-7-NF2 (HPO2), together regulate the activity of LATS1/2 kinases and YAP/TAZ transcriptional co-activators. In mouse livers, the genetic inactivation of either HPO1 or HPO2 module results in partial activation of YAP/TAZ, bile duct hyperplasia, and hepatocellular carcinoma (HCC). On the contrary, inactivation of both HPO1 and HPO2 modules results in full activation of YAP/TAZ, rapid development of intrahepatic cholangiocarcinoma (iCCA), and early lethality. Interestingly, HPO1 has a predominant role in regulating organ size. HPO1 inactivation causes a homogenous YAP/TAZ activation and cell proliferation across the whole liver, resulting in a proportional and rapid increase in liver size. Thus, this study has reconstructed the order of the Hippo signaling network and suggests that LATS1/2 and YAP/TAZ activities are finetuned by HPO1 and HPO2 modules to cause different cell fates, organ size changes, and tumorigenesis trajectories.

The origins and functional effects of postzygotic mutations throughout the human life span

Postzygotic mutations (PZMs) begin to accrue in the human genome immediately after fertilization, but how and when PZMs affect development and lifetime health remain unclear. To study the origins and functional consequences of PZMs, we generated a multitissue atlas of PZMs spanning 54 tissue and cell types from 948 donors. Nearly half the variation in mutation burden among tissue samples can be explained by measured technical and biological effects, and 9% can be attributed to donor-specific effects. Through phylogenetic reconstruction of PZMs, we found that their type and predicted functional impact vary during prenatal development, across tissues, and through the germ cell life cycle. Thus, methods for interpreting effects across the body and the life span are needed to fully understand the consequences of genetic variants.

Changes of development from childhood to late adulthood in rats tracked by urinary proteome

To date, studies of development have mainly focused on the embryonic stage and a short time thereafter. There has been little research on the whole life of an individual from childhood to aging and death. For the first time, we used non-invasive urinary proteome technology to track changes in several important developmental timepoints in a group of rats, covering 10 timepoints from childhood, adolescence, young adulthood, middle adulthood, and near death in old age. Similar to previous studies on puberty, proteins were detected involved in sexual or reproductive maturation, mature spermatozoa in seminiferous tubules (first seen), gonadal hormones, decline of oestradiol, brain growth, and central nervous system myelination, and our differential protein enrichment pathways also included reproductive system development, tube development, response to hormone, response to oestradiol, brain development, and neuron development. Similar to previous studies in young adults, proteins were detected involved in musculoskeletal maturity, peak bone mass, development of the immune system, and growth and physical development, and our differential protein enrichment pathways also included skeletal system development, bone regeneration, system development, immune system processes, myeloid leukocyte differentiation, growth, and developmental growth. Studies on aging-related changes in neurons and neurogenesis have been reported, and we also found relevant pathways in aged rats, such as regulation of neuronal synaptic plasticity and positive regulation of long-term neuronal synaptic plasticity. However, at all timepoints throughout life, there were many biological pathways revealed by differential urinary protein enrichment involving multiple organs, tissues, systems, etc., that have not been mentioned in existing studies. This study shows comprehensive and detailed changes in rat lifetime development through the urinary proteome, helping to fill the gap in development research. Moreover, it provides a new approach to monitoring changes in human health and diseases of aging using the urinary proteome.

Time-resolved transcriptomics of mouse gastric pit cells during postnatal development reveals features distinct from whole stomach development

Whole-organ transcriptomic analyses have emerged as a common method for characterizing developmental transitions in mammalian organs. However, it is unclear if all cell types in an organ follow the whole-organ defined developmental trajectory. Recently, a postnatal two-stage developmental process was described for the mouse stomach. Here, using laser capture microdissection to obtain in situ transcriptomic data, we show that mouse gastric pit cells exhibit four postnatal developmental stages. Interestingly, early stages are characterized by the up-regulation of genes associated with metabolism, a functionality not typically associated with pit cells. Hence, beyond revealing that not all constituent cells develop according to the whole-organ determined pathway, these results broaden our understanding of the pit cell phenotypic landscape during stomach development.

The imprinted gene Zac1 regulates steatosis in developmental cadmium-induced nonalcoholic fatty liver disease

Cadmium (Cd) exposure in adulthood is associated with nonalcoholic fatty liver disease (NAFLD), characterized by steatosis, inflammation, and fibrosis. The prevalence of NAFLD in children is increasing, suggesting a role for the developmental environment in programming susceptibility. However, the role of developmental Cd exposure in programming NAFLD and the underlying mechanisms remain unclear. We have proposed that imprinted genes are strong candidates for connecting the early life environment and later life disease. In support of this, we previously identified roles for the Imprinted Gene Network (IGN) and its regulator Zac1 in programming NAFLD in response to maternal metabolic dysfunction. Here, we test the hypothesis that developmental Cd exposure is sufficient to program NAFLD, and further, that this process is mediated by Zac1 and the IGN. Using mice, we show that developmental cadmium chloride (CdCl2) exposure leads to histological, biochemical, and molecular signatures of steatosis and fibrosis in juveniles. Transcriptomic analyses comparing livers of CdCl2-exposed and control mice show upregulation of Zac1 and the IGN coincident with disease presentation. Increased hepatic Zac1 expression is independent of promoter methylation and imprinting statuses. Finally, we show that over-expression of Zac1 in cultured hepatocytes is sufficient to induce lipid accumulation in a Pparγ-dependent manner and demonstrate direct binding of Zac1 to the Pparγ promoter. Our findings demonstrate that developmental Cd exposure is sufficient to program NAFLD in later life, and with our previous work, establish Zac1 and the IGN as key regulators of prosteatotic and profibrotic pathways, two of the major pathological hallmarks of NAFLD.

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

Untargeted Metabolome Analysis Reveals Reductions in Maternal Hepatic Glucose and Amino Acid Content That Correlate with Fetal Organ Weights in a Mouse Model of Fetal Alcohol Spectrum Disorders

Prenatal alcohol exposure (PAE) causes fetal growth restrictions. A major driver of fetal growth deficits is maternal metabolic disruption; this is under-investigated following PAE. Untargeted metabolomics on the dam and fetus exposed to alcohol (ALC) revealed that the hepatic metabolome of ALC and control (CON) dams were distinct, whereas that of ALC and CON fetuses were similar. Alcohol reduced maternal hepatic glucose content and enriched essential amino acid (AA) catabolites, N-acetylated AA products, urea content, and free fatty acids. These alterations suggest an attempt to minimize the glucose gap by increasing gluconeogenesis using AA and glycerol. In contrast, ALC fetuses had unchanged glucose and AA levels, suggesting an adequate draw of maternal nutrients, despite intensified stress on ALC dams. Maternal metabolites including glycolytic intermediates, AA catabolites, urea, and one-carbon-related metabolites correlated with fetal liver and brain weights, whereas lipid metabolites correlated with fetal body weight, indicating they may be drivers of fetal weight outcomes. Together, these data suggest that ALC alters maternal hepatic metabolic activity to limit glucose availability, thereby switching to alternate energy sources to meet the high-energy demands of pregnancy. Their correlation with fetal phenotypic outcomes indicates the influence of maternal metabolism on fetal growth and development.

Proteome-centric cross-omics characterization and integrated network analyses of triple-negative breast cancer

We report a comprehensive proteomic study of a 90-case cohort of paired samples of triple-negative breast cancer (TNBC) in quantification, phosphorylation, and DNA-binding capacity. Four integrative subtypes (iP-1-4) are stratified on the basis of global proteome and phosphoproteome, each of which exhibits distinct molecular and pathway features. Scaffold and co-expression network analyses of three proteomic datasets, integrated with those from genome and transcriptome of the same cohort, reveal key pathways and master regulators that, characteristic of TNBC subtypes, play important regulatory roles within and between scaffold sub-structures and co-expression communities. We find that NAE1 is a potential drug target for subtype iP-1, and a series of key molecules in fatty acid metabolism, such as AKT1/FASN, are plausible targets for subtype iP-2. Libraries of proteins, pathways and networks of TNBC provide a valuable molecular infrastructure for further clinical exploration and in-depth studies of the molecular mechanisms of the disease.

Multi-omics analyses reveal the mechanisms of Arsenic-induced male reproductive toxicity in mice

Arsenic (As), a widespread environmental contaminant, can induce serious male reproductive injury; however, the underlying mechanisms remain unclear. Multi-omics analyses, including transcriptome, proteome, and phosphoproteome could promote our understanding of As-induced male reproductive toxicity. Here, we established the reproductive injured mice model by intraperitoneal injection of NaAsO2 (8 mg/kg body weight), which was validated by reduced reproductive cells, sperm motility, and litter size. The followed multi-omics analyses of mice revealed that As exposure inhibited ATP production by decreasing the expression of proteins HK1, and GAPDHS, and the enzymatic activities of PDH and SDH. The inhibition of mitochondrial activity and increase in HDAC2 and MTA3 dysregulated the lysine acetylation levels of histone and global proteins. Specifically, the downregulated histones H4K5ac and H4K12ac and upregulated histone H3K9ac disordered the distribution of TP1 to interfere with spermatogenesis. Moreover, As could reduce the expression of COL1A1, RAB13, and LSR to disrupt the junctions between seminiferous tubules, and thereinto, the inhibition of RAB13 increased PKA-dependent phosphorylation. Our study reveals that As causes male reproductive toxicity through decreasing energy production, altering histone acetylation, and impairing cell junctions. Our findings provide basic data for further studies on As reproductive toxicity.

Using Cell Type-Specific Genes to Identify Cell-Type Transitions Between Different in vitro Culture Conditions

In vitro differentiation or expansion of stem and progenitor cells under chemical stimulation or genetic manipulation is used for understanding the molecular mechanisms of cell differentiation and self-renewal. However, concerns around the cell identity of in vitro-cultured cells exist. Bioinformatics methods, which rely heavily on signatures of cell types, have been developed to estimate cell types in bulk samples. The Tabula Muris Senis project provides an important basis for the comprehensive identification of signatures for different cell types. Here, we identified 46 cell type-specific (CTS) gene clusters for 83 mouse cell types. We conducted Gene Ontology term enrichment analysis on the gene clusters and revealed the specific functions of the relevant cell types. Next, we proposed a simple method, named CTSFinder, to identify different cell types between bulk RNA-Seq samples using the 46 CTS gene clusters. We applied CTSFinder on bulk RNA-Seq data from 17 organs and from developing mouse liver over different stages. We successfully identified the specific cell types between organs and captured the dynamics of different cell types during liver development. We applied CTSFinder with bulk RNA-Seq data from a growth factor-induced neural progenitor cell culture system and identified the dynamics of brain immune cells and nonimmune cells during the long-time cell culture. We also applied CTSFinder with bulk RNA-Seq data from reprogramming induced pluripotent stem cells and identified the stage when those cells were massively induced. Finally, we applied CTSFinder with bulk RNA-Seq data from in vivo and in vitro developing mouse retina and captured the dynamics of different cell types in the two development systems. The CTS gene clusters and CTSFinder method could thus serve as promising toolkits for assessing the cell identity of in vitro culture systems.

A MicroRNA-Based Network Provides Potential Predictive Signatures and Reveals the Crucial Role of PI3K/AKT Signaling for Hepatic Lineage Maturation

Background

Functions of miRNAs involved in tumorigenesis are well reported, yet, their roles in normal cell lineage commitment remain ambiguous. Here, we investigated a specific "transcription factor (TF)-miRNA-Target" regulatory network during the lineage maturation of biliary tree stem cells (BTSCs) into adult hepatocytes (hAHeps).

Method

Bioinformatic analysis was conducted based on our RNA-seq and microRNA-seq datasets with four human hepatic-lineage cell lines, including hBTSCs, hepatic stem cells (hHpSCs), hepatoblasts (hHBs), and hAHeps. Short time-series expression miner (STEM) analysis was performed to reveal the time-dependent dynamically changed miRNAs and mRNAs. GO and KEGG analyses were applied to reveal the potential function of key miRNAs and mRNAs. Then, the miRDB, miRTarBase, TargetScan, miRWalk, and DIANA-microT-CDS databases were adopted to predict the potential targets of miRNAs while the TransmiR v2.0 database was used to obtain the experimentally supported TFs that regulate miRNAs. The TCGA, Kaplan-Meier Plotter, and human protein atlas (HPA) databases and more than 10 sequencing data, including bulk RNA-seq, microRNA-seq, and scRNA-seq data related to hepatic development or lineage reprogramming, were obtained from both our or other published studies for validation.

Results

STEM analysis showed that during the maturation from hBTSCs to hAHeps, 52 miRNAs were downwardly expressed and 928 mRNA were upwardly expressed. Enrichment analyses revealed that those 52 miRNAs acted as pluripotency regulators for stem cells and participated in various novel signaling pathways, including PI3K/AKT, MAPK, and etc., while 928 mRNAs played important roles in liver-functional metabolism. With an extensive sorting of those key miRNAs and mRNAs based on the target prediction results, 23 genes were obtained which not only functioned as the targets of 17 miRNAs but were considered critical for the hepatic lineage commitment. A "TF-miRNA-Target" regulatory network for hepatic lineage commitment was therefore established and had been well validated by various datasets. The network revealed that the PI3K/AKT pathway was gradually suppressed during the hepatic commitment.

Conclusion

A total of 17 miRNAs act as suppressors during hepatic maturation mainly by regulating 23 targets and modulating the PI3K/AKT signaling pathway. The regulatory network uncovers possible signatures and guidelines enabling us to identify or obtain the functional hepatocytes for future study.

In silico analysis suggests disruption of interactions between HAMP from hepatocytes and SLC40A1 from macrophages in hepatocellular carcinoma

BACKGROUND

Identification of factors associated with proliferation in the hepatocellular carcinoma (HCC) microenvironment aids in understanding the mechanisms of disease progression and provides druggable targets. Gene expression profiles of individual cells in HCC and para-carcinoma tissues can be effectively obtained using the single-cell RNA sequencing (scRNA-Seq) technique. Here, we aimed to identify proliferative hepatocytes from HCC and para-carcinoma tissues, detect differentially expressed genes between the two types of proliferative hepatocytes, and investigate their potential roles in aberrant proliferation.

RESULTS

Two respective gene signatures for proliferative cells and hepatocytes were established and used to identify proliferative hepatocytes from HCC and para-carcinoma tissues based on scRNA-Seq data. Gene expression profiles between the two types of proliferative hepatocytes were compared. Overall, 40 genes were upregulated in proliferative hepatocytes from para-carcinoma tissue, whereas no upregulated genes were detected in those from HCC tissue. Twelve of the genes, including HAMP, were specifically expressed in the liver tissue. Based on previous reports, we found that HAMP modulates cell proliferation through interaction with its receptor SLC40A1. Comprehensive analysis of cells in HCC and para-carcinoma tissues revealed that: (1) HAMP is specifically expressed in hepatocytes and significantly downregulated in malignant hepatocytes; (2) a subset of macrophages expressing SLC40A1 and genes reacting to various infections is present in para-carcinoma but not in HCC tissue. We independently validated the findings with scRNA-Seq and large-scale tissue bulk RNA-Seq/microarray analyses.

CONCLUSION

HAMP was significantly downregulated in malignant hepatocytes. In addition, a subset of macrophages expressing SLC40A1 and genes reacting to various infections was absent in HCC tissue. These findings support the involvement of HAMP-SLC40A1 signaling in aberrant hepatocyte proliferation in the HCC microenvironment. The collective data from our in silico analysis provide novel insights into the mechanisms underlying HCC progression and require further validation with wet laboratory experiments.

Jonckheere-Terpstra-Kendall-based non-parametric analysis of temporal differential gene expression

Time-course experiments using parallel sequencers have the potential to uncover gradual changes in cells over time that cannot be observed in a two-point comparison. An essential step in time-series data analysis is the identification of temporal differentially expressed genes (TEGs) under two conditions (e.g. control versus case). Model-based approaches, which are typical TEG detection methods, often set one parameter (e.g. degree or degree of freedom) for one dataset. This approach risks modeling of linearly increasing genes with higher-order functions, or fitting of cyclic gene expression with linear functions, thereby leading to false positives/negatives. Here, we present a Jonckheere-Terpstra-Kendall (JTK)-based non-parametric algorithm for TEG detection. Benchmarks, using simulation data, show that the JTK-based approach outperforms existing methods, especially in long time-series experiments. Additionally, application of JTK in the analysis of time-series RNA-seq data from seven tissue types, across developmental stages in mouse and rat, suggested that the wave pattern contributes to the TEG identification of JTK, not the difference in expression levels. This result suggests that JTK is a suitable algorithm when focusing on expression patterns over time rather than expression levels, such as comparisons between different species. These results show that JTK is an excellent candidate for TEG detection.

Time-resolved multi-omics analysis reveals the role of nutrient stress-induced resource reallocation for TAG accumulation in oleaginous fungus Mortierella alpina

BACKGROUND

Global resource reallocation is an established critical strategy through which organisms deal with environmental stress. The regulation of intracellular lipid storage or utilization is one of the most important strategies for maintaining energy homeostasis and optimizing growth. Oleaginous microorganisms respond to nitrogen deprivation by inducing lipid hyper accumulation; however, the associations between resource allocation and lipid accumulation are poorly understood.

RESULTS

Here, the time-resolved metabolomics, lipidomics, and proteomics data were generated in response to nutrient availability to examine how metabolic alternations induced by nitrogen deprivation drive the triacylglycerols (TAG) accumulation in M. alpina. The subsequent accumulation of TAG under nitrogen deprivation was a consequence of the reallocation of carbon, nitrogen sources, and lipids, rather than an up-regulation of TAG biosynthesis genes. On one hand, nitrogen deprivation induced the down-regulation of isocitrate dehydrogenase level in TCA cycle and redirected glycolytic flux of carbon from amino acid biosynthesis into fatty acids' synthesis; on the other hand, nitrogen deprivation induced the up-regulation of cell autophagy and ubiquitin-mediated protein proteolysis which resulted in a recycling of preformed protein nitrogen and carbon. Combining with the up-regulation of glutamate decarboxylase and succinic semialdehyde dehydrogenase in GABA shunt, and the phosphoenolpyruvate carboxykinase in the central hub involving pyruvate/phosphoenolpyruvate/oxaloacetate, the products from nitrogen-containing compounds degradation were recycled to be intermediates of TCA cycle and be shunted toward de novo biosynthesis of fatty acids. We found that nitrogen deprivation increased the protein level of phospholipase C/D that contributes to degradation of phosphatidylcholine and phosphatidylethanolamine, and supplied acyl chains for TAG biosynthesis pathway. In addition, ATP from substrate phosphorylation was presumed to be a critical factor regulation of the global resource allocation and fatty acids' synthesis rate.

CONCLUSIONS

The present findings offer a panoramic view of resource allocation by M. alpina in response to nutrient stress and revealed a set of intriguing associations between resource reallocation and TAG accumulation. This system-level insight provides a rich resource with which to explore in-depth functional characterization and gain information about the strategic combination of strain development and process integration to achieve optimal lipid productivity under nutrient stress.