Single-cell RNA and T-cell receptor sequencing unveil mycosis fungoides heterogeneity and a possible gene signature.

Background: Mycosis fungoides (MF) is the most common subtype of cutaneous T-cell lymphoma (CTCL). Comprehensive analysis of MF cells in situ and ex vivo is complicated by the fact that is challenging to distinguish malignant from reactive T cells with certainty. Methods: To overcome this limitation, we performed combined single-cell RNA (scRNAseq) and T-cell receptor TCR sequencing (scTCRseq) of skin lesions of cutaneous MF lesions from 12 patients. A sufficient quantity of living T cells was obtained from 9 patients, but 2 had to be excluded due to unclear diagnoses (coexisting CLL or revision to a fixed toxic drug eruption). Results: From the remaining patients we established single-cell mRNA expression profiles and the corresponding TCR repertoire of 18,630 T cells. TCR clonality unequivocally identified 13,592 malignant T cells. Reactive T cells of all patients clustered together, while malignant cells of each patient formed a unique cluster expressing genes typical of naive/memory, such as CD27, CCR7 and IL7R, or cytotoxic T cells, e.g., GZMA, NKG7 and GNLY. Genes encoding classic CTCL markers were not detected in all clusters, consistent with the fact that mRNA expression does not correlate linearly with protein expression. Nevertheless, we successfully pinpointed distinctive gene signatures differentiating reactive malignant from malignant T cells: keratins (KRT81, KRT86), galectins (LGALS1, LGALS3) and S100 genes (S100A4, S100A6) being overexpressed in malignant cells. Conclusions: Combined scRNAseq and scTCRseq not only allows unambiguous identification of MF cells, but also revealed marked heterogeneity between and within patients with unexpected functional phenotypes. While the correlation between mRNA and protein abundance was limited with respect to established MF markers, we were able to identify a single-cell gene expression signature that distinguishes malignant from reactive T cells.

Wnt/beta-catenin activated non pilomatrical carcinoma of the skin: a case series.

Among skin epithelial tumors, recurrent mutations in the APC/CTNNB1 genes resulting in activation of the Wnt/ß-catenin pathway have been reported predominantly in neoplasms with matrical differentiation. In the present study, we describe the morphologic, immunohistochemical, and genetic features of 16 primary cutaneous carcinomas harboring mutations activating the Wnt/ß-catenin pathway without evidence of matrical differentiation, as well as four combined tumors in which a similar Wnt/ß-catenin activated carcinoma component was associated with Merkel cell carcinoma or pilomatrical carcinoma. Among the pure tumor cases, 6/16 patients were female with a median age of 80 years (range: 58-98). Tumors were located on the head and neck (n=7, 44%), upper limb (n=4, 25%), trunk (n=3, 18%), and leg (n=2, 13%). Metastatic spread was observed in 4 cases resulting in death from disease in one patient. Microscopically, all cases were poorly differentiated neoplasms infiltrating the dermis and/or subcutaneous tissue. In 13 cases, solid "squamoid" areas were associated with a basophilic component characterized by rosette/pseudoglandular formation resulting in a biphasic appearance. Three specimens consisted only of poorly differentiated carcinoma lacking rosette formation. Immunohistochemical studies showed frequent expression of EMA (100%), BerEP4 (100%), cytokeratin 7 (94%), chromogranin A (44%), synaptophysin (82%) and cytokeratin 20 (69%). Complete loss of Rb expression was observed in all but one case. Nuclear ß-catenin and CDX2 expressions were detected in all cases. Recurrent pathogenic somatic mutations were observed in APC (60%), CTNNB1 (40%) and RB1 (n=47%). Global methylation analysis confirmed that cases with rosette formation constituted a homogenous tumor group distinct from established skin tumor entities (pilomatrical carcinoma, Merkel cell carcinoma and squamous cell carcinoma) while the 3 other cases lacking such morphologic features did not. In addition, we identified four combined neoplasms in which there was a component showing a similar poorly differentiated rosette forming carcinoma demonstrating Rb loss and beta-catenin activation associated with either Merkel cell carcinoma (n=3) or pilomatrical carcinoma (n=1). In conclusion, we describe a distinctive neoplasm, for which we propose the term "Wnt/ß-catenin activated rosette-forming carcinoma", morphologically characterized by the association of rosette formation, squamous and/or neuroendocrine differentiation, diffuse CDX2 expression, Rb loss, and mutations in CTNNB1/APC genes.

Merkel cell polyomavirus pan-T antigen knockdown reduces cancer cell stemness and promotes neural differentiation independent of RB1.

Merkel cell carcinoma (MCC) is a highly aggressive skin cancer associated with integration of Merkel cell polyomavirus (MCPyV). MCPyV-encoded T-antigens (TAs) are pivotal for sustaining MCC's oncogenic phenotype, i.e., repression of TAs results in reactivation of the RB pathway and subsequent cell cycle arrest. However, the MCC cell line LoKe, characterized by a homozygous loss of the RB1 gene, exhibits uninterrupted cell cycle progression after shRNA-mediated TA repression. This unique feature allows an in-depth analysis of the effects of TAs beyond inhibition of the RB pathway, revealing the decrease in expression of stem cell-related genes upon panTA-knockdown. Analysis of gene regulatory networks identified members of the E2F family (E2F1, E2F8, TFDP1) as key transcriptional regulators that maintain stem cell properties in TA-expressing MCC cells. Furthermore, minichromosome maintenance (MCM) genes, which encodes DNA-binding licensing proteins essential for stem cell maintenance, were suppressed upon panTA-knockdown. The decline in stemness occurred simultaneously with neural differentiation, marked by the increased expression of neurogenesis-related genes such as neurexins, BTG2, and MYT1L. This upregulation can be attributed to heightened activity of PBX1 and BPTF, crucial regulators of neurogenesis pathways. The observations in LoKe were confirmed in an additional MCPyV-positive MCC cell line in which RB1 was silenced before panTA-knockdown. Moreover, spatially resolved transcriptomics demonstrated reduced TA expression in situ in a part of a MCC tumor characterized by neural differentiation. In summary, TAs are critical for maintaining stemness of MCC cells and suppressing neural differentiation, irrespective of their impact on the RB-signaling pathway.

Analyses of combined Merkel cell carcinomas with neuroblastic components suggests that loss of T antigen expression in Merkel cell carcinoma may result in cell cycle arrest and neuroblastic transdifferentiation.

Merkel cell carcinoma (MCC) is an aggressive skin cancer frequently caused by genomic integration of the Merkel cell polyomavirus (MCPyV). MCPyV-negative cases often present as combined MCCs, which represent a distinctive subset of tumors characterized by association of an MCC with a second tumor component, mostly squamous cell carcinoma. Up to now, only exceptional cases of combined MCC with neuroblastic differentiation have been reported. Herein we describe two additional combined MCCs with neuroblastic differentiation and provide comprehensive morphologic, immunohistochemical, transcriptomic, genetic and epigenetic characterization of these tumors, which both arose in elderly men and appeared as an isolated inguinal adenopathy. Microscopic examination revealed biphasic tumors combining a poorly differentiated high-grade carcinoma with a poorly differentiated neuroblastic component lacking signs of proliferation. Immunohistochemical investigation revealed keratin 20 and MCPyV T antigen (TA) in the MCC parts, while neuroblastic differentiation was confirmed in the other component in both cases. A clonal relation of the two components can be deduced from 20 and 14 shared acquired point mutations detected by whole exome analysis in both combined tumors, respectively. Spatial transcriptomics demonstrated a lower expression of stem cell marker genes such as SOX2 and MCM2 in the neuroblastic component. Interestingly, although the neuroblastic part lacked TA expression, the same genomic MCPyV integration and the same large T-truncating mutations were observed in both tumor parts. Given that neuronal transdifferentiation upon TA repression has been reported for MCC cell lines, the most likely scenario for the two combined MCC/neuroblastic tumors is that neuroblastic transdifferentiation resulted from loss of TA expression in a subset of MCC cells. Indeed, DNA methylation profiling suggests an MCC-typical cellular origin for the combined MCC/neuroblastomas. © 2024 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

The HDAC inhibitor domatinostat induces type I interferon α in Merkel cell carcinoma by HES1 repression.

BACKGROUND: Class I selective histone deacetylase inhibitors (HDACi) have been previously demonstrated to not only increase major histocompatibility complex class I surface expression in Merkel cell carcinoma (MCC) cells by restoring the antigen processing and presentation machinery, but also exert anti-tumoral effect by inducing apoptosis. Both phenomena could be due to induction of type I interferons (IFN), as has been described for HDACi. However, the mechanism of IFN induction under HDACi is not fully understood because the expression of IFNs is regulated by both activating and inhibitory signaling pathways. Our own preliminary observations suggest that this may be caused by suppression of HES1. METHODS: The effect of the class I selective HDACi domatinostat and IFNα on cell viability and the apoptosis of MCPyV-positive (WaGa, MKL-1) and -negative (UM-MCC 34) MCC cell lines, as well as, primary fibroblasts were assessed by colorimetric methods or measuring mitochondrial membrane potential and intracellular caspase-3/7, respectively. Next, the impact of domatinostat on IFNA and HES1 mRNA expression was measured by RT-qPCR; intracellular IFNα production was detected by flow cytometry. To confirm that the expression of IFNα induced by HDACi was due to the suppression of HES1, it was silenced by RNA interference and then mRNA expression of IFNA and IFN-stimulated genes was assessed. RESULTS: Our studies show that the previously reported reduction in viability of MCC cell lines after inhibition of HDAC by domatinostat is accompanied by an increase in IFNα expression, both of mRNA and at the protein level. We confirmed that treatment of MCC cells with external IFNα inhibited their proliferation and induced apoptosis. Re-analysis of existing single-cell RNA sequencing data indicated that induction of IFNα by domatinostat occurs through repression of HES1, a transcriptional inhibitor of IFNA; this was confirmed by RT-qPCR. Finally, siRNA-mediated silencing of HES1 in the MCC cell line WaGa not only increased mRNA expression of IFNA and IFN-stimulated genes but also decreased cell viability. CONCLUSION: Our results demonstrate that the direct anti-tumor effect of HDACi domatinostat on MCC cells is at least in part mediated via decreased HES1 expression allowing the induction of IFNα, which in turn causes apoptosis.

Integrative omics analysis identifies biomarkers of idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is an interstitial lung disease characterized by chronic progressive pulmonary fibrosis and a poor prognosis. Genetic studies, including transcriptomic and proteomics, have provided new insight into revealing mechanisms of IPF. Herein we provided a novel strategy to identify biomarkers by integrative analysis of transcriptomic and proteomic profiles of IPF patients. We examined the landscape of IPF patients' gene expression in the transcription and translation phases and investigated the expression and functions of two new potential biomarkers. Differentially expressed (DE) mRNAs were mainly enriched in pathways associated with immune system activities and inflammatory responses, while DE proteins are related to extracellular matrix production and wound repair. The upregulated genes in both phases are associated with wound repair and cell differentiation, while the downregulated genes in both phases are associated with reduced immune activities and the damage of the alveolar tissues. On this basis, we identified thirteen potential marker genes. Among them, we validated the expression changes of butyrophilin-like 9 (BTNL9) and plasmolipin (PLLP) and investigated their functional pathways in the IPF mechanism. Both genes are downregulated in the tissues of IPF patients and Bleomycin-induced mice, and co-expression analysis indicates that they have a protective effect by inhibiting extracellular matrix production and promoting wound repair in alveolar epithelial cells.

Dissecting the heterogeneity and tumorigenesis of BRCA1 deficient mammary tumors via single cell RNA sequencing.

Background: BRCA1 plays critical roles in mammary gland development and mammary tumorigenesis. And loss of BRCA1 induces mammary tumors in a stochastic manner. These tumors present great heterogeneity at both intertumor and intratumor levels. Methods: To comprehensively elucidate the heterogeneity of BRCA1 deficient mammary tumors and the underlying mechanisms for tumor initiation and progression, we conducted bulk and single cell RNA sequencing (scRNA-seq) on both mammary gland cells and mammary tumor cells isolated from Brca1 knockout mice. Results: We found the BRCA1 deficient tumors could be classified into four subtypes with distinct molecular features and different sensitivities to anti-cancer drugs at the intertumor level. Whereas within the tumors, heterogeneous subgroups were classified mainly due to the different activities of cell proliferation, DNA damage response/repair and epithelial-to-mesenchymal transition (EMT). Besides, we reconstructed the BRCA1 related mammary tumorigenesis to uncover the transcriptomes alterations during this process via pseudo-temporal analysis of the scRNA-seq data. Furthermore, from candidate markers for BRCA1 mutant tumors, we discovered and validated one oncogene Mrc2, whose loss could reduce mammary tumor growth in vitro and in vivo. Conclusion: Our study provides a useful resource for better understanding of mammary tumorigenesis induced by BRCA1 deficiency.

Gene fusion of IL7 involved in the regulation of idiopathic pulmonary fibrosis.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a rare form of immune-mediated interstitial lung disease characterized by progressive pulmonary fibrosis and scarring. The pathogenesis of IPF is still unclear. Gene fusion events exist universally during transcription and show alternated patterns in a variety of lung diseases. Therefore, the comprehension of the function of gene fusion in IPF might shed light on IPF pathogenesis research and facilitate treatment development. METHODS: In this study, we included 91 transcriptome datasets from the National Center for Biotechnology Information (NCBI), including 52 IPF patients and 39 healthy controls. We detected fusion events in these datasets and probed gene fusion-associated differential gene expression and functional pathways. To obtain robust results, we corrected the batch bias across different projects. RESULTS: We identified 1550 gene fusion events in all transcriptomes and studied the possible impacts of IL7 = AC083837.1 gene fusion. The two genes locate adjacently in chromosome 8 and share the same promoters. Their fusion is associated with differential expression of 282 genes enriched in six Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 35 functional gene sets. Gene ontology (GO) enrichment analysis shows that IL7 = AC083837.1 gene fusion is associated with the enrichment of 187 gene sets. The co-expression network of interleukin-7 (IL7) indicates that decreased IL7 expression is associated with many pathways that regulate IPF progress. CONCLUSION: Based on the results, we conclude that IL7 = AC083837.1 gene fusion might exacerbate fibrosis in IPF via enhancing activities of natural killer cell-mediated cytotoxicity, skin cell apoptosis, and vessel angiogenesis, the interaction of which contributes to the development of fibrosis and the deterioration of respiratory function of IPF patients. Our work unveils the possible roles of gene fusion in regulating IPF and demonstrates that gene fusion investigation is a valid approach in probing immunologic mechanisms and searching potential therapeutic targets for treating IPF.The reviews of this paper are available via the supplemental material section.

Conservation analysis of SARS-CoV-2 spike suggests complicated viral adaptation history from bat to human.

BACKGROUND: The current coronavirus disease 2019 (COVID-19) pandemic, caused by severe acute respiratory syndrome (SARS)-CoV-2, has become the most devastating public health emergency in the 21st century and one of the most influential plagues in history. Studies on the origin of SARS-CoV-2 have generally agreed that the virus probably comes from bat, closely related to a bat CoV named BCoV-RaTG13 taken from horseshoe bat (Rhinolophus affinis), with Malayan pangolin (Manis javanica) being a plausible intermediate host. However, due to the relatively low number of SARS-CoV-2-related strains available in public domain, the evolutionary history remains unclear. METHODOLOGY: Nine hundred ninety-five coronavirus sequences from NCBI Genbank and GISAID were obtained and multiple sequence alignment was carried out to categorize SARS-CoV-2 related groups. Spike sequences were analyzed using similarity analysis and conservation analyses. Mutation analysis was used to identify variations within receptor-binding domain (RBD) in spike for SARS-CoV-2-related strains. RESULTS: We identified a family of SARS-CoV-2-related strains, including the closest relatives, bat CoV RaTG13 and pangolin CoV strains. Sequence similarity analysis and conservation analysis on spike sequence identified that N-terminal domain, RBD and S2 subunit display different degrees of conservation with several coronavirus strains. Mutation analysis on contact sites in SARS-CoV-2 RBD reveals that human-susceptibility probably emerges in pangolin. CONCLUSION AND IMPLICATION: We conclude that the spike sequence of SARS-CoV-2 is the result of multiple recombination events during its transmission from bat to human, and we propose a framework of evolutionary history that resolve the relationship of BCoV-RaTG13 and pangolin coronaviruses with SARS-CoV-2. LAY SUMMARY: This study analyses whole-genome and spike sequences of coronavirus from NCBI using phylogenetic and conservation analyses to reconstruct the evolutionary history of severe acute respiratory syndrome (SARS)-CoV-2 and proposes an evolutionary history of spike in the progenitors of SARS-CoV-2 from bat to human through mammal hosts before they recombine into the current form.

Co-expression network analysis of lncRNAs and mRNAs in rat liver tissue reveals the complex interactions in response to pathogenic cytotoxicity.

Liver is one of the most vital organs to maintain homeostasis because of its peculiar detoxification functionalities to detoxify chemicals and metabolize drugs and toxins. Due to its crucial functions, the liver is also prone to various diseases, i.e., hepatitis, cirrhosis and hepatoma, etc. Additionally, long non-coding RNAs (lncRNAs) has emerged as key regulators which are found to play important roles in transcription, splicing, translation, replication, chromatin shaping and post translational modification of proteins in living cells. However, the underlying mechanisms of biological processes mediated by lncRNA remain unclear. Here, with the aim of disclosing potential lncRNAs implicated in the biological processes in liver in response to cytotoxicity, we performed a co-expression network analysis based on the transcriptome data of the damaged liver tissue of Rattus norvegicus induced by three cytotoxic compounds (carbon tetrachloride, chloroform and thioacetamide). Our analysis unveils that many biological processes and pathways were collectively affected by the three cytotoxic compounds, including drug metabolism, oxidation-reduction process, oxidative stress, glucuronidation, liver development and flavonoid biosynthetic process, etc. Also, our network analysis has identified several highly conserved lncRNA-mRNA interactions participating in those correlated processes and pathways, implying their potential roles in response to the induced cytotoxicity in liver. Our study provides new insights into lncRNA-mRNA regulatory mechanisms in response to pathogenic cytotoxic damaging in liver and facilitates the development of lncRNA-oriented therapies for hepatic diseases in the future.

Transcriptome analysis reveals lncRNA-mediated complex regulatory network response to DNA damage in the liver tissue of Rattus norvegicus.

DNA is prone to damages, which would result in genetic disorders and enhance risk of tumorigenesis. Hence, understanding the molecular mechanisms of DNA damage and repair will provide deep insights into tumorigenesis, carcinogenesis as well as the corresponding treatments. Aiming at investigating potential long noncoding RNAs (lncRNAs) response against DNA damage, we performed a comprehensive transcriptomic analysis based on RNA sequencing data of the liver tissue from Rattus norvegicus, in which DNA damage was induced using aflatoxin B1, ifosfamide and N-nitrosodimethylamine. Through our analyses, numerous novel lncRNAs are identified for the first time, and differential network analysis discloses lncRNA-mediated regulatory networks related to DNA-damage response. The result shows that these DNA-damage-inducing chemicals might disrupt many lncRNA-mediated interactions involved in diverse biological processes and pathways, for example, immune function and cell adhesion. In contrast, the host might also activate a few RNA interactions in response to DNA damage, involving response to drug and regulation of cell cycle.