Identification of Colorectal Cancer Cell Stemness from Single-Cell RNA Sequencing.

Cancer stem cells (CSC) play a critical role in metastasis, relapse, and therapy resistance in colorectal cancer. While characterization of the normal lineage of cell development in the intestine has led to the identification of many genes involved in the induction and maintenance of pluripotency, recent studies suggest significant heterogeneity in CSC populations. Moreover, while many canonical colorectal cancer CSC marker genes have been identified, the ability to use these classical markers to annotate stemness at the single-cell level is limited. In this study, we performed single-cell RNA sequencing on a cohort of 6 primary colon, 9 liver metastatic tumors, and 11 normal (nontumor) controls to identify colorectal CSCs at the single-cell level. Finding poor alignment of the 11 genes most used to identify colorectal CSC, we instead extracted a single-cell stemness signature (SCS_sig) that robustly identified "gold-standard" colorectal CSCs that expressed all marker genes. Using this SCS_sig to quantify stemness, we found that while normal epithelial cells show a bimodal distribution, indicating distinct stem and differentiated states, in tumor epithelial cells stemness is a continuum, suggesting greater plasticity in these cells. The SCS_sig score was quite variable between different tumors, reflective of the known transcriptomic heterogeneity of CRC. Notably, patients with higher SCS_sig scores had significantly shorter disease-free survival time after curative intent surgical resection, suggesting stemness is associated with relapse. IMPLICATIONS: This study reveals significant heterogeneity of expression of genes commonly used to identify colorectal CSCs, and identifies a novel stemness signature to identify these cells from scRNA-seq data.

Resolving the Spatial and Cellular Architecture of Lung Adenocarcinoma by Multiregion Single-Cell Sequencing.

Little is known of the geospatial architecture of individual cell populations in lung adenocarcinoma (LUAD) evolution. Here, we perform single-cell RNA sequencing of 186,916 cells from five early-stage LUADs and 14 multiregion normal lung tissues of defined spatial proximities from the tumors. We show that cellular lineages, states, and transcriptomic features geospatially evolve across normal regions to LUADs. LUADs also exhibit pronounced intratumor cell heterogeneity within single sites and transcriptional lineage-plasticity programs. T regulatory cell phenotypes are increased in normal tissues with proximity to LUAD, in contrast to diminished signatures and fractions of cytotoxic CD8+ T cells, antigen-presenting macrophages, and inflammatory dendritic cells. We further find that the LUAD ligand-receptor interactome harbors increased expression of epithelial CD24, which mediates protumor phenotypes. These data provide a spatial atlas of LUAD evolution, and a resource for identification of targets for its treatment. SIGNIFICANCE: The geospatial ecosystem of the peripheral lung and early-stage LUAD is not known. Our multiregion single-cell sequencing analyses unravel cell populations, states, and phenotypes in the spatial and ecologic evolution of LUAD from the lung that comprise high-potential targets for early interception.This article is highlighted in the In This Issue feature, p. 2355.

Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients.

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 pandemic. Severely symptomatic COVID-19 is associated with lung inflammation, pneumonia, and respiratory failure, thereby raising concerns of elevated risk of COVID-19-associated mortality among lung cancer patients. Angiotensin-converting enzyme 2 (ACE2) is the major receptor for SARS-CoV-2 entry into lung cells. The single-cell expression landscape of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung cancer patients remains unknown. We sought to delineate single-cell expression profiles of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung adenocarcinoma (LUAD) patients. We examined the expression levels and cellular distribution of ACE2 and SARS-CoV-2-priming proteases TMPRSS2 and TMPRSS4 in 5 LUADs and 14 matched normal tissues by single-cell RNA-sequencing (scRNA-seq) analysis. scRNA-seq of 186,916 cells revealed epithelial-specific expression of ACE2, TMPRSS2, and TMPRSS4. Analysis of 70,030 LUAD- and normal-derived epithelial cells showed that ACE2 levels were highest in normal alveolar type 2 (AT2) cells and that TMPRSS2 was expressed in 65% of normal AT2 cells. Conversely, the expression of TMPRSS4 was highest and most frequently detected (75%) in lung cells with malignant features. ACE2-positive cells co-expressed genes implicated in lung pathobiology, including COPD-associated HHIP, and the scavengers CD36 and DMBT1. Notably, the viral scavenger DMBT1 was significantly positively correlated with ACE2 expression in AT2 cells. We describe normal and tumor lung epithelial populations that express SARS-CoV-2 receptor and proteases, as well as major host defense genes, thus comprising potential treatment targets for COVID-19 particularly among lung cancer patients.

Alpha1 soluble guanylyl cyclase (sGC) splice forms as potential regulators of human sGC activity.

Soluble guanylyl cyclase (sGC), a key protein in the NO/cGMP signaling pathway, is an obligatory heterodimeric protein composed of one alpha- and one beta-subunit. The alpha(1)/beta(1) sGC heterodimer is the predominant form expressed in various tissues and is regarded as the major isoform mediating NO-dependent effects such as vasodilation. We have identified three new alpha(1) sGC protein variants generated by alternative splicing. The 363 residue N1-alpha(1) sGC splice variant contains the regulatory domain, but lacks the catalytic domain. The shorter N2-alpha(1) sGC maintains 126 N-terminal residues and gains an additional 17 unique residues. The C-alpha(1) sGC variant lacks 240 N-terminal amino acids, but maintains a part of the regulatory domain and the entire catalytic domain. Q-PCR of N1-alpha(1), N2-alpha(1) sGC mRNA levels together with RT-PCR analysis for C-alpha(1) sGC demonstrated that the expression of the alpha(1) sGC splice forms vary in different human tissues indicative of tissue-specific regulation. Functional analysis of the N1-alpha(1) sGC demonstrated that this protein has a dominant-negative effect on the activity of sGC when coexpressed with the alpha(1)/beta(1) heterodimer. The C-alpha(1) sGC variant heterodimerizes with the beta(1) subunit and produces a fully functional NO- and BAY41-2272-sensitive enzyme. We also found that despite identical susceptibility to inhibition by ODQ, intracellular levels of the 54-kDa C-alpha(1) band did not change in response to ODQ treatments, while the level of 83 kDa alpha(1) band was significantly affected by ODQ. These studies suggest that modulation of the level and diversity of splice forms may represent novel mechanisms modulating the function of sGC in different human tissues.

Ligand selectivity of soluble guanylyl cyclase: effect of the hydrogen-bonding tyrosine in the distal heme pocket on binding of oxygen, nitric oxide, and carbon monoxide.

Although soluble guanylyl cyclase (sGC) functions in an environment in which O(2), NO, and CO are potential ligands for its heme moiety, the enzyme displays a high affinity for only its physiological ligand, NO, but has a limited affinity for CO and no affinity for O(2). Recent studies of a truncated version of the sGC beta(1)-subunit containing the heme-binding domain (Boon, E. M., Huang, S H., and Marletta, M. A. (2005) Nat. Chem. Biol., 1, 53-59) showed that introduction of the hydrogen-bonding tyrosine into the distal heme pocket changes the ligand specificity of the heme moiety and results in an oxygen-binding sGC. The hypothesis that the absence of hydrogen-bonding residues in the distal heme pocket is sufficient to provide oxygen discrimination by sGC was put forward. We tested this hypothesis in a context of a complete sGC heterodimer containing both the intact alpha(1)- and beta(1)-subunits. We found that the I145Y substitution in the full-length beta-subunit of the sGC heterodimer did not produce an oxygen-binding enzyme. However, this substitution impeded the association of NO and destabilized the NO.heme complex. The tyrosine in the distal heme pocket also impeded both the binding and dissociation of the CO ligand. We propose that the mechanism of oxygen exclusion by sGC not only involves the lack of hydrogen bonding in the distal heme pocket, but also depends on structural elements from other domains of sGC.