Comprehensive single-cell transcriptome analysis reveals heterogeneity in endometrioid adenocarcinoma tissues

Microcirculatory disturbance in acute liver injury is triggered by IFNγ-CD40 axis

BACKGROUND

Acute liver failure (ALF) is a life-threatening disorder that progresses from self-limiting acute liver injury (ALI). Microcirculatory disturbance characterized by sinusoidal hypercoagulation and subsequent massive hypoxic hepatocyte damage have been proposed to be the mechanism by which ALI deteriorates to ALF; however, the precise molecular pathway of the sinusoidal hypercoagulation remains unknown. Here, we analyzed ALI patients and mice models to uncover the pathogenesis of ALI with microcirculatory disturbance.

METHODS

We conducted a single-center retrospective study for ALI and blood samples and liver tissues were analyzed to evaluate the microcirculatory disturbance in ALI patients (n = 120). Single-cell RNA sequencing analysis (scRNA-seq) was applied to the liver from the concanavalin A (Con A)‑induced mouse model of ALI. Interferon-gamma (IFNγ) and tumor necrosis factor-alpha knockout mice, and primary human liver sinusoidal endothelial cells (LSECs) were used to assess the mechanism of microcirculatory disturbance.

RESULTS

The serum IFNγ concentrations were significantly higher in ALI patients with microcirculatory disturbance than in patients without microcirculatory disturbance, and the IFNγ was upregulated in the Con A mouse model which presented microcirculatory disturbance. Hepatic IFNγ expression was increased as early as 1 hour after Con A treatment prior to sinusoidal hypercoagulation and hypoxic liver damage. scRNA-seq revealed that IFNγ was upregulated in innate lymphoid cells and stimulated hepatic vascular endothelial cells at the early stage of liver injury. In IFNγ knockout mice treated with Con A, the sinusoidal hypercoagulation and liver damage were remarkably attenuated, concomitant with the complete inhibition of CD40 and tissue factor (TF) upregulation in vascular endothelial cells. By ligand-receptor analysis, CD40-CD40 ligand interaction was identified in vascular endothelial cells. In human LSECs, IFNγ upregulated CD40 expression and TF was further induced by increased CD40-CD40 ligand interaction. Consistent with these findings, hepatic CD40 expression was significantly elevated in human ALI patients with microcirculatory disturbance.

CONCLUSION

We identified the critical role of the IFNγ-CD40 axis as the molecular mechanism of microcirculatory disturbance in ALI. This finding may provide novel insights into the pathogenesis of ALI and potentially contribute to the emergence of new therapeutic strategies for ALI patients.

Clonal Origin and Lineage Ambiguity in Mixed Neuroendocrine Carcinoma of the Uterine Cervix

Small-cell neuroendocrine carcinoma (SCNEC) of the cervix is a rare disease characterized by a high incidence of mixed tumors with other types of cancer. The mechanism underlying this mixed phenotype is not well understood. This study established a panel of organoid lines from patients with SCNEC of the cervix and ultimately focused on one line, which retained a mixed tumor phenotype, both in vitro and in vivo. Histologically, both organoids and xenograft tumors showed distinct differentiation into either SCNEC or adenocarcinoma in some regions and ambiguous differentiation in others. Tracking single cells indicated the existence of cells with bipotential differentiation toward SCNEC and adenocarcinomas. Single-cell transcriptional analysis identified three distinct clusters: SCNEC-like, adenocarcinoma-like, and a cluster lacking specific differentiation markers. The expression of neuroendocrine markers was enriched in the SCNEC-like cluster but not exclusively. Human papillomavirus 18 E6 was enriched in the SCNEC-like cluster, which showed higher proliferation and lower levels of the p53 pathway. After treatment with anticancer drugs, the expression of adenocarcinoma markers increased, whereas that of SCNEC decreased. Using a reporter system for keratin 19 expression, changes in the differentiation of each cell were shown to be associated with the shift in differentiation induced by drug treatment. These data suggest that mixed SCNEC/cervical tumors have a clonal origin and are characterized by an ambiguous and flexible differentiation state.

Phospholipase A2 Group IIA Is Associated with Inflammatory Hepatocellular Adenoma

Although benign hepatocellular adenomas (HCA) are very rare, recent observations have shown their occurrence in patients with diabetes mellitus. Consequently, most of these cases are treated by resection due to concerns regarding their potential progression to hepatocarcinoma (HCC). This decision is largely driven by the limited number of studies on HCC subtyping and the lack of molecular and biological insights into the carcinogenic potential of benign tumors. This study aimed to comprehensively investigate the subtype classification of HCA and to compare and analyze gene expression profiling between HCA and HCC tissues. One fresh inflammatory HCA (I-HCA), three non-B non-C HCCs, two hepatitis B virus-HCCs, and one normal liver tissue sample were subjected to single-cell RNA sequencing (scRNA-seq). Comparative analysis of scRNA-seq among different tissues showed that phospholipase A2 group IIA (PLA2G2A) mRNA was specifically expressed in I-HCA, following RNA-seq analysis in formalin-fixed paraffin-embedded tissues from other HCAs. Immunohistochemistry using the PLA2G2A antibody in these tissues indicated that the positive reaction was mainly observed in hepatocytes of I-HCAs and stromal cells surrounding the tumor tissue in HCC were also stained. According to a clinical database, PLA2G2A expression in HCC does not correlate with poor prognosis. This finding may potentially help develop a new definition for I-HCA, resulting in a significant clinical contribution, but it requires validation with other fresh HCA samples.

A New Era of Integration between Multiomics and Spatio-Temporal Analysis for the Translation of EMT towards Clinical Applications in Cancer

Epithelial-mesenchymal transition (EMT) is crucial to metastasis by increasing cancer cell migration and invasion. At the cellular level, EMT-related morphological and functional changes are well established. At the molecular level, critical signaling pathways able to drive EMT have been described. Yet, the translation of EMT into efficient diagnostic methods and anti-metastatic therapies is still missing. This highlights a gap in our understanding of the precise mechanisms governing EMT. Here, we discuss evidence suggesting that overcoming this limitation requires the integration of multiple omics, a hitherto neglected strategy in the EMT field. More specifically, this work summarizes results that were independently obtained through epigenomics/transcriptomics while comprehensively reviewing the achievements of proteomics in cancer research. Additionally, we prospect gains to be obtained by applying spatio-temporal multiomics in the investigation of EMT-driven metastasis. Along with the development of more sensitive technologies, the integration of currently available omics, and a look at dynamic alterations that regulate EMT at the subcellular level will lead to a deeper understanding of this process. Further, considering the significance of EMT to cancer progression, this integrative strategy may enable the development of new and improved biomarkers and therapeutics capable of increasing the survival and quality of life of cancer patients.

Molecular insights of a CBP/β-catenin-signaling inhibitor on nonalcoholic steatohepatitis-induced liver fibrosis and disorder

Nonalcoholic steatohepatitis (NASH) is a progressive fibrotic disease associated with an increased risk of developing hepatocellular carcinoma; at present, no efficient therapeutic strategy has been established. Herein, we examined the efficacy of PRI-724, a potent inhibitor of CBP/β-catenin signaling, for treating NASH-related liver fibrosis and disorder and characterized its mechanism. Choline-deficient, L-amino acid-defined, high-fat diet (CDAHFD)-fed mice exhibited NASH-induced liver fibrosis that is characterized by steatosis, lobular inflammation, hepatocellular injury and collagen fibrils. To examine the therapeutic effect, CDAHFD-fed mice were administered PRI-724. Serum levels of ALT and pro-fibrotic molecule, i.e. Mac-2 bp, alpha smooth muscle actin, type I and type III collagens, decreased significantly. mRNA levels of the matrix metalloproteinases Mmp8 and Mmp9 in the liver were significantly increased, and increases in the abundance of MMP9-producing neutrophils and macrophages were observed. Marco+Mmp9+Cd68+ Kupffer cells were only observed in the livers of mice treated with PRI-724, and Mmp9 expression in Marco+Cd68+ Kupffer cells increased 4.3-fold. Moreover, hepatic expression of the lipid metabolism regulator, pyruvate dehydrogenase kinase 4 and liver lipid droplets also decreased significantly. PRI-724-treated NASH mice not only recovered from NASH-related liver fibrosis through the effect of PRI-724 down-regulating the expression of pro-fibrotic genes and up-regulating the expression of anti-fibrotic genes, but they also recovered from NASH-induced liver disorder. PRI-724, a selective CBP/β-catenin inhibitor, thus shows a potent therapeutic effect for NASH-related liver fibrosis and for decreasing adipose tissue in the liver.

Single-cell RNA sequencing to detect age-associated genes that identify senescent cells in the liver of aged mice

Senescent cells are predicted to occur and increase in animal tissues with aging. However, senescent cells in the tissues of aged animals remain to be identified. We refer to the marker genes to identify senescent cells in tissues as "age-associated genes". In this study, we searched for age-associated genes to identify senescent cells in the livers of aged animals. We performed single-cell RNA sequencing (scRNA-seq) to screen candidates for age-associated genes using young and aged rat primary hepatocytes. To remove animal species specificity, gene expression analyses in mouse livers were performed, confirming age-associated increases in the mRNA expression levels of Glipr1, Clec12a, and Phlda3. Moreover, the mRNA expression levels of Glipr1 and Phlda3 were increased by stress-induced premature senescence using doxorubicin in primary hepatocytes and livers of young mice. Transcriptome data of aged rat hepatocytes suggested that Glipr1, Clec12a, and Phlda3 were expressed in almost identical cells. Fluorescence in situ hybridization (FISH) confirmed the presence of cells with abundant Glipr1, Clec12a, and Phlda3 mRNA in 27-month-old mouse primary hepatocytes, which are considered to be senescent cells. This study is the first to identify Glipr1, Clec12a, and Phlda3 as age-associated genes in the mouse liver.

TGF-β generates a population of cancer cells residing in G1 phase with high motility and metastatic potential via KRTAP2-3

Transforming growth factor β (TGF-β) increases epithelial cancer cell migration and metastasis by inducing epithelial-mesenchymal transition (EMT). TGF-β also inhibits cell proliferation by inducing G1 phase cell-cycle arrest. However, the correlation between these tumor-promoting and -suppressing effects remains unclear. Here, we show that TGF-β confers higher motility and metastatic ability to oral cancer cells in G1 phase. Mechanistically, keratin-associated protein 2-3 (KRTAP2-3) is a regulator of these dual effects of TGF-β, and its expression is correlated with tumor progression in patients with head and neck cancer and migratory and metastatic potentials of oral cancer cells. Furthermore, single-cell RNA sequencing reveals that TGF-β generates two populations of mesenchymal cancer cells with differential cell-cycle status through two distinctive EMT pathways mediated by Slug/HMGA2 and KRTAP2-3. Thus, TGF-β-induced KRTAP2-3 orchestrates cancer cell proliferation and migration by inducing EMT, suggesting motile cancer cells arrested in G1 phase as a target to suppress metastasis.

Single-cell analysis of human primary prostate cancer reveals the heterogeneity of tumor-associated epithelial cell states

Prostate cancer is the second most common malignancy in men worldwide and consists of a mixture of tumor and non-tumor cell types. To characterize the prostate cancer tumor microenvironment, we perform single-cell RNA-sequencing on prostate biopsies, prostatectomy specimens, and patient-derived organoids from localized prostate cancer patients. We uncover heterogeneous cellular states in prostate epithelial cells marked by high androgen signaling states that are enriched in prostate cancer and identify a population of tumor-associated club cells that may be associated with prostate carcinogenesis. ERG-negative tumor cells, compared to ERG-positive cells, demonstrate shared heterogeneity with surrounding luminal epithelial cells and appear to give rise to common tumor microenvironment responses. Finally, we show that prostate epithelial organoids harbor tumor-associated epithelial cell states and are enriched with distinct cell types and states from their parent tissues. Our results provide diagnostically relevant insights and advance our understanding of the cellular states associated with prostate carcinogenesis.

Effect of Aberrant Long Noncoding RNA on the Prognosis of Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is a kind of lethal cancer. Although there are mature treatment methods, there is still a lack of rigorous and scientific means for cancer diagnosis. Long noncoding RNAs (lncRNAs) are a kind of noncoding RNA (ncRNA). Recent studies find that alteration of lncRNA expression is related to the occurrence of many cancers. In order to find lncRNAs which can effectively predict the prognosis of ccRCC, RNA-seq count data and clinical information were downloaded from TCGA-KIRC, and gene expression profiles from 530 patients were included. Then, K-means was used for clustering, and the number of clusters was determined to be 5. The R-package "edgeR" was used to perform differential expression analysis. Subsequently, a risk model composed of 10 lncRNA biomarkers significantly related to prognosis was identified via Cox and LASSO regression analyses. Then, patients were divided into two groups according to the model-based risk score, and then, GSEA pathway enrichment was performed. The results showed that metabolism- and mTOR-related pathways were activated while immune-related pathways were inhibited in the high-risk patients. Combined with previous studies, it is believed that these 10 lncRNAs are potential targets for the treatment of ccRCC. In addition, Cox regression analysis was used to verify the independence of the risk model, and as results revealed, the risk model can be used to independently predict the prognosis of patients. In conclusion, our study found 10 lncRNAs related to the prognosis of ccRCC and provided new ideas for clinical diagnosis and drug development.

Advances in understanding the molecular pathology of gynecological malignancies: the role and potential of RNA sequencing

For many years technological limitations restricted the progress of identifying the underlying genetic causes of gynecologicalcancers. However, during the past decade, high-throughput next-generation sequencing technologies have revolutionized cancer research. RNA sequencing has arisen as a very useful technique in expanding our understanding of genome changes in cancer. Cancer is characterized by the accumulation of genetic alterations affecting genes, including substitutions, insertions, deletions, translocations, gene fusions, and alternative splicing. If these aberrant genes become transcribed, aberrations can be detected by RNA sequencing, which will also provide information on the transcript abundance revealing the expression levels of the aberrant genes. RNA sequencing is considered the technique of choice when studying gene expression and identifying new RNA species. This is due to the quantitative and qualitative improvement that it has brought to transcriptome analysis, offering a resolution that allows research into different layers of transcriptome complexity. It has also been successful in identifying biomarkers, fusion genes, tumor suppressors, and uncovering new targets responsible for drug resistance in gynecological cancers. To illustrate that we here review the role of RNA sequencing in studies that enhanced our understanding of the molecular pathology of gynecological cancers.

Organoids of the Female Reproductive Tract: Innovative Tools to Study Desired to Unwelcome Processes

The pelviperineal organs of the female reproductive tract form an essential cornerstone of human procreation. The system comprises the ectodermal external genitalia, the Müllerian upper-vaginal, cervical, endometrial and oviductal derivatives, and the endodermal ovaries. Each of these organs presents with a unique course of biological development as well as of malignant degeneration. For many decades, various preclinical in vitro models have been employed to study female reproductive organ (patho-)biology, however, facing important shortcomings of limited expandability, loss of representativeness and inadequate translatability to the clinic. The recent emergence of 3D organoid models has propelled the field forward by generating powerful research tools that in vitro replicate healthy as well as diseased human tissues and are amenable to state-of-the-art experimental interventions. Here, we in detail review organoid modeling of the different female reproductive organs from healthy and tumorigenic backgrounds, and project perspectives for both scientists and clinicians.

DOCK11 and DENND2A play pivotal roles in the maintenance of hepatitis B virus in host cells

Human hepatitis B virus (HBV) infection remains a serious health problem worldwide. However, the mechanism for the maintenance of HBV in a latent state within host cells remains unclear. Here, using single-cell RNA sequencing analysis, we identified four genes linked to the maintenance of HBV in a liver cell line expressing HBV RNA at a low frequency. These genes included DOCK11 and DENND2A, which encode small GTPase regulators. In primary human hepatocytes infected with HBV, knockdown of these two genes decreased the amount of both HBV DNA and covalently closed circular DNA to below the limit of detection. Our findings reveal a role for DOCK11 and DENND2A in the maintenance of HBV.

Follicle-Stimulating Hormone Induces Lipid Droplets via Gαi/o and β-Arrestin in an Endometrial Cancer Cell Line

Follicle-stimulating hormone (FSH) and its G protein-coupled receptor, FSHR, represents a paradigm for receptor signaling systems that activate multiple and complex pathways. Classically, FSHR activates Gαs to increase intracellular levels of cAMP, but its ability to activate other G proteins, and β-arrestin-mediated signaling is well documented in many different cell systems. The pleiotropic signal capacity of FSHR offers a mechanism for how FSH drives multiple and dynamic downstream functions in both gonadal and non-gonadal cell types, including distinct diseases, and how signal bias may be achieved at a pharmacological and cell system-specific manner. In this study, we identify an additional mechanism of FSH-mediated signaling and downstream function in the endometrial adenocarcinoma Ishikawa cell line. While FSH did not induce increases in cAMP levels, this hormone potently activated pertussis toxin sensitive Gαi/o signaling. A selective allosteric FSHR ligand, B3, also activated Gαi/o signaling in these cells, supporting a role for receptor-mediated activation despite the low levels of FSHR mRNA. The low expression levels may attribute to the lack of Gαs/cAMP signaling as increasing FSHR expression resulted in FSH-mediated activation of the Gαs pathway. Unlike prior reports for FSH-mediated Gαs/cAMP signaling, FSH-mediated Gαi/o signaling was not affected by inhibition of dynamin-dependent receptor internalization. While chronic FSH did not alter cell viability, FSH was able to increase lipid droplet size. The β-arrestins are key adaptor proteins known to regulate FSHR signaling. Indeed, a rapid, FSH-dependent increase in interactions between β-arrestin1 and Gαi1 was observed via NanoBiT complementation in Ishikawa cells. Furthermore, both inhibition of Gαi/o signaling and siRNA knockdown of β-arrestin 1/2 significantly reduced FSH-induced lipid droplet accumulation, implying a role for a Gαi/o/β-arrestin complex in FSH functions in this cell type. As FSH/FSHR has been implicated in distinct hormone-dependent cancers, including endometrial cancer, analysis of the cancer genome database from 575 human endometrial adenocarcinoma tumors revealed that a subpopulation of samples expressed FSHR. Overall, this study highlights a novel mechanism for FSHR signal pleiotropy that may be exploited for future personalized therapeutic approaches.

SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data

BACKGROUND

Droplet-based single-cell RNA sequence analyses assume that all acquired RNAs are endogenous to cells. However, any cell-free RNAs contained within the input solution are also captured by these assays. This sequencing of cell-free RNA constitutes a background contamination that confounds the biological interpretation of single-cell transcriptomic data.

RESULTS

We demonstrate that contamination from this "soup" of cell-free RNAs is ubiquitous, with experiment-specific variations in composition and magnitude. We present a method, SoupX, for quantifying the extent of the contamination and estimating "background-corrected" cell expression profiles that seamlessly integrate with existing downstream analysis tools. Applying this method to several datasets using multiple droplet sequencing technologies, we demonstrate that its application improves biological interpretation of otherwise misleading data, as well as improving quality control metrics.

CONCLUSIONS

We present SoupX, a tool for removing ambient RNA contamination from droplet-based single-cell RNA sequencing experiments. This tool has broad applicability, and its application can improve the biological utility of existing and future datasets.

Age-associated changes in the transcriptomes of non-cultured adipose-derived stem cells from young and old mice assessed via single-cell transcriptome analysis

Adipose-derived stem cells (ASCs) exhibit self-renewal and pluripotency. The differentiation potency of ASCs has been reported to deteriorate with aging; however, relevant studies used ASCs that were isolated and subcultured several times. It is still unclear whether subcultured ASCs accurately reflect the in vivo state. To address this question, we used freshly isolated stromal vascular fractions (SVFs) and performed comprehensive single-cell transcriptome analysis. In this study, we identified three cell populations as putative ASC candidates in SVFs and three novel ASC-related genes: Adamts7, Snai2, and Tgfbr1, that are reported to be negative regulators of cell differentiation. Moreover, we identified age-associated high gene expression levels of Adamts7, Egfr, and Igfbp4 in the earliest differentiation stage of ASCs. These results suggest that aging may make it impossible to maintain the stringency of the regulation of the expression of some genes related to ASC differentiation.

Human-Derived Model Systems in Gynecological Cancer Research

The human female reproductive tract (FRT) is a complex system that combines series of organs, including ovaries, fallopian tubes, uterus, cervix, vagina, and vulva; each of which possesses unique cellular characteristics and functions. This versatility, in turn, allows for the development of a wide range of epithelial gynecological cancers with distinct features. Thus, reliable model systems are required to better understand the diverse mechanisms involved in the regional pathogenesis of the reproductive tract and improve treatment strategies. Here, we review the current human-derived model systems available to study the multitude of gynecological cancers, including ovarian, endometrial, cervical, vaginal, and vulvar cancer, and the recent advances in the push towards personalized therapy.

Endometrial Tumour Microenvironment

Endometrial cancer (EC) is the most common gynaecological tumour in developed countries, and its incidence is increasing in part due to the prevalence of obesity and its related hormone dysregulation. As described in this chapter, the tumour microenvironment plays a principal role in unopposed oestrogen stimulation promoting tumour cell proliferation. Factors and cytokines secreted by the different cell types defining the reactive tumour stroma also determine the invasive abilities of the tumour cells. Cancer-associated fibroblasts and tumour-associated macrophages actively participate through SDF-1, TGF-b or HGF to promote epithelial-to-mesenchymal transition or to generate an appropriate tumour niche. Likewise, endothelial cells facilitate lymph node and vascular infiltration through VEGF. Finally, the possibility to balance the immunosuppressive phenotypes in advanced endometrial cancer through the tumour microenvironment will probably represent a main therapeutic strategy in the near future.

New tools for pathology: a user's review of a highly multiplexed method for in situ analysis of protein and RNA expression in tissue

Tumor cell heterogeneity and tumor cell-stromal interactions are being explored as determinants of disease progression and treatment resistance in solid tumor and hematological malignancies. As such, tools simultaneously capable of highly multiplexed profiling of tissues' protein and RNA content, as well as interrogation of rare or single cells, are required to precisely characterize constituent tumor cell populations, infiltrating lymphocytes and stromal elements. Access to spatial relationships will enable more precise characterization of tumors, support patient stratification and may help to identify novel drug targets. Multiple platforms are being developed to address these critical unmet needs. The NanoString digital spatial profiling (DSP) platform enables highly multiplexed, spatial assessment of protein and/or RNA targets in tissues by detecting oligonucleotide barcodes conjugated via a photocleavable linker to primary antibodies or nucleic acid probes. Although this platform enables high-dimensional spatial interrogation of tissue protein and RNA expression, a detailed understanding of its composition, function and chemistry is advisable to guide experimental design and data interpretation. The purpose of this review is to provide an independent, comprehensive description of the DSP technology, including an overview of NanoString's capture and antibody barcode conjugation chemistries, experimental workflow, data output and analysis methods. The DSP technology will be discussed in the context of other highly multiplexed immunohistochemistry methods, including imaging mass cytometry and multiplexed ion beam imaging, to inform potential users of the advantages and limitations of each. Additional issues such as preanalytical variability, sampling and specimen adequacy will be considered with respect to the platforms to inform potential experimental design. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Nx1-Seq (Well Based Single-Cell Analysis System)

Research on the hierarchical nature of cell differentiation and heterogeneity in tissues has been performed by isolating and identifying cells by the use of monoclonal antibodies, cell sorting, microdissection, and functional assays. However, it is difficult to analyze continuous changes in cell differentiation and the identification of cells for which cell markers are unclear. Furthermore, cell populations considered identical were shown to be diverse. Recently, single cell gene expression analysis was performed to help understand the complexity of cell populations. Single-cell analysis can analyze the diversity of individual cell populations as well as the tissue microenvironment, and is extremely useful for research on intercellular interactions in diseases and identifying specific marker genes. Recent advances in technology have made it possible to analyze hundreds of single cells. In this paper, we introduce our newly developed well-based single-cell transcriptome method, which includes other methods.