Single-cell RNA sequencing of cystic fibrosis liver disease explants reveals endothelial complement activation.

BACKGROUND & AIMS: Cystic fibrosis (CF) is considered a multisystemic disorder in which CF-associated liver disease (CFLD) is the third most common cause of mortality. Currently, no effective treatment is available for CFLD because its pathophysiology is still unclear. Interestingly, CFLD exhibits identical vascular characteristics as non-cirrhotic portal hypertension, recently classified as porto-sinusoidal vascular disorders (PSVD). METHODS: Since endothelial cells (ECs) are an important component in PSVD, we performed single-cell RNA sequencing (scRNA-seq) on four explant livers from CFLD patients to identify differential endothelial characteristics which could contribute to the disease. We comprehensively characterized the endothelial compartment and compared it with publicly available scRNA-seq datasets from cirrhotic and healthy livers. Key gene signatures were validated ex vivo on patient tissues. RESULTS: We found that ECs from CF liver explants are more closely related to healthy than cirrhotic patients. In CF patients we also discovered a distinct population of liver sinusoidal ECs-coined CF LSECs-upregulating genes involved in the complement cascade and coagulation. Finally, our immunostainings further validated the predominant periportal location of CF LSECs. CONCLUSIONS: Our work showed novel aspects of human liver ECs at the single-cell level thereby supporting endothelial involvement in CFLD, and reinforcing the hypothesis that ECs could be a driver of PSVD. Therefore, considering the vascular compartment in CF and CFLD may help developing new therapeutic approaches for these diseases.

The pulmonary vasculature in lethal COVID-19 and idiopathic pulmonary fibrosis at single-cell resolution.

AIMS: Severe acute respiratory syndrome coronavirus-2 infection causes COVID-19, which in severe cases evokes life-threatening acute respiratory distress syndrome (ARDS). Transcriptome signatures and the functional relevance of non-vascular cell types (e.g. immune and epithelial cells) in COVID-19 are becoming increasingly evident. However, despite its known contribution to vascular inflammation, recruitment/invasion of immune cells, vascular leakage, and perturbed haemostasis in the lungs of severe COVID-19 patients, an in-depth interrogation of the endothelial cell (EC) compartment in lethal COVID-19 is lacking. Moreover, progressive fibrotic lung disease represents one of the complications of COVID-19 pneumonia and ARDS. Analogous features between idiopathic pulmonary fibrosis (IPF) and COVID-19 suggest partial similarities in their pathophysiology, yet, a head-to-head comparison of pulmonary cell transcriptomes between both conditions has not been implemented to date. METHODS AND RESULTS: We performed single-nucleus RNA-sequencing on frozen lungs from 7 deceased COVID-19 patients, 6 IPF explant lungs, and 12 controls. The vascular fraction, comprising 38 794 nuclei, could be subclustered into 14 distinct EC subtypes. Non-vascular cell types, comprising 137 746 nuclei, were subclustered and used for EC-interactome analyses. Pulmonary ECs of deceased COVID-19 patients showed an enrichment of genes involved in cellular stress, as well as signatures suggestive of dampened immunomodulation and impaired vessel wall integrity. In addition, increased abundance of a population of systemic capillary and venous ECs was identified in COVID-19 and IPF. COVID-19 systemic ECs closely resembled their IPF counterparts, and a set of 30 genes was found congruently enriched in systemic ECs across studies. Receptor-ligand interaction analysis of ECs with non-vascular cell types in the pulmonary micro-environment revealed numerous previously unknown interactions specifically enriched/depleted in COVID-19 and/or IPF. CONCLUSIONS: This study uncovered novel insights into the abundance, expression patterns, and interactomes of EC subtypes in COVID-19 and IPF, relevant for future investigations into the progression and treatment of both lethal conditions.

Single cell atlas identifies lipid-processing and immunomodulatory endothelial cells in healthy and malignant breast.

Since a detailed inventory of endothelial cell (EC) heterogeneity in breast cancer (BC) is lacking, here we perform single cell RNA-sequencing of 26,515 cells (including 8433 ECs) from 9 BC patients and compare them to published EC taxonomies from lung tumors. Angiogenic ECs are phenotypically similar, while other EC subtypes are different. Predictive interactome analysis reveals known but also previously unreported receptor-ligand interactions between ECs and immune cells, suggesting an involvement of breast EC subtypes in immune responses. We also identify a capillary EC subtype (LIPEC (Lipid Processing EC)), which expresses genes involved in lipid processing that are regulated by PPAR-γ and is more abundant in peri-tumoral breast tissue. Retrospective analysis of 4648 BC patients reveals that treatment with metformin (an indirect PPAR-γ signaling activator) provides long-lasting clinical benefit and is positively associated with LIPEC abundance. Our findings warrant further exploration of this LIPEC/PPAR-γ link for BC treatment.

Generation of vessel co-option lung metastases mouse models for single-cell isolation of metastases-derived cells and endothelial cells.

Tumor vessel co-option, a process in which cancer cells "hijack" pre-existing blood vessels to grow and invade healthy tissue, is poorly understood but is a proposed resistance mechanism against anti-angiogenic therapy (AAT). Here, we describe protocols for establishing murine renal (RENCA) and breast (4T1) cancer lung vessel co-option metastases models. Moreover, we outline a reproducible protocol for single-cell isolation from murine lung metastases using magnetic-activated cell sorting as well as immunohistochemical stainings to distinguish vessel co-option from angiogenesis. For complete details on the use and execution of this protocol, please refer to Teuwen et al. (2021).

Protocols for endothelial cell isolation from mouse tissues: kidney, spleen, and testis.

Endothelial cells (ECs) exhibit phenotypic and functional tissue specificities, critical for studies in the vascular field and beyond. Thus, tissue-specific methods for isolation of highly purified ECs are necessary. Kidney, spleen, and testis ECs are relevant players in health and diseases such as chronic kidney disease, acute kidney injury, myelofibrosis, and cancer. Here, we provide tailored protocols for rapid and reproducible EC purification established for scRNA sequencing from these adult murine tissues using the combination of magnetic- and fluorescence-activated cell sorting. For complete details on the use and execution of these protocols, please refer to Kalucka et al. (2020) and Dumas et al. (2020).

Arhgef2 regulates mitotic spindle orientation in hematopoietic stem cells and is essential for productive hematopoiesis.

How hematopoietic stem cells (HSCs) coordinate their divisional axis and whether this orientation is important for stem cell-driven hematopoiesis is poorly understood. Single-cell RNA sequencing data from patients with Shwachman-Diamond syndrome (SDS), an inherited bone marrow failure syndrome, show that ARHGEF2, a RhoA-specific guanine nucleotide exchange factor and determinant of mitotic spindle orientation, is specifically downregulated in SDS hematopoietic stem and progenitor cells (HSPCs). We demonstrate that transplanted Arhgef2-/- fetal liver and bone marrow cells yield impaired hematopoietic recovery and a production deficit from long-term HSCs, phenotypes that are not the result of differences in numbers of transplanted HSCs, their cell cycle status, level of apoptosis, progenitor output, or homing ability. Notably, these defects are functionally restored in vivo by overexpression of ARHGEF2 or its downstream activated RHOA GTPase. By using live imaging of dividing HSPCs, we show an increased frequency of misoriented divisions in the absence of Arhgef2. ARHGEF2 knockdown in human HSCs also impairs their ability to regenerate hematopoiesis, culminating in significantly smaller xenografts. Together, these data demonstrate a conserved role for Arhgef2 in orienting HSPC division and suggest that HSCs may divide in certain orientations to establish hematopoiesis, the loss of which could contribute to HSC dysfunction in bone marrow failure.

Tumor vessel co-option probed by single-cell analysis.

Tumor vessel co-option is poorly understood, yet it is a resistance mechanism against anti-angiogenic therapy (AAT). The heterogeneity of co-opted endothelial cells (ECs) and pericytes, co-opting cancer and myeloid cells in tumors growing via vessel co-option, has not been investigated at the single-cell level. Here, we use a murine AAT-resistant lung tumor model, in which VEGF-targeting induces vessel co-option for continued growth. Single-cell RNA sequencing (scRNA-seq) of 31,964 cells reveals, unexpectedly, a largely similar transcriptome of co-opted tumor ECs (TECs) and pericytes as their healthy counterparts. Notably, we identify cell types that might contribute to vessel co-option, i.e., an invasive cancer-cell subtype, possibly assisted by a matrix-remodeling macrophage population, and another M1-like macrophage subtype, possibly involved in keeping or rendering vascular cells quiescent.

Single-Cell RNA Sequencing Maps Endothelial Metabolic Plasticity in Pathological Angiogenesis.

Endothelial cell (EC) metabolism is an emerging target for anti-angiogenic therapy in tumor angiogenesis and choroidal neovascularization (CNV), but little is known about individual EC metabolic transcriptomes. By single-cell RNA sequencing 28,337 murine choroidal ECs (CECs) and sprouting CNV-ECs, we constructed a taxonomy to characterize their heterogeneity. Comparison with murine lung tumor ECs (TECs) revealed congruent marker gene expression by distinct EC phenotypes across tissues and diseases, suggesting similar angiogenic mechanisms. Trajectory inference predicted that differentiation of venous to angiogenic ECs was accompanied by metabolic transcriptome plasticity. ECs displayed metabolic transcriptome heterogeneity during cell-cycle progression and in quiescence. Hypothesizing that conserved genes are important, we used an integrated analysis, based on congruent transcriptome analysis, CEC-tailored genome-scale metabolic modeling, and gene expression meta-analysis in cross-species datasets, followed by in vitro and in vivo validation, to identify SQLE and ALDH18A1 as previously unknown metabolic angiogenic targets.

An Integrated Gene Expression Landscape Profiling Approach to Identify Lung Tumor Endothelial Cell Heterogeneity and Angiogenic Candidates.

Heterogeneity of lung tumor endothelial cell (TEC) phenotypes across patients, species (human/mouse), and models (in vivo/in vitro) remains poorly inventoried at the single-cell level. We single-cell RNA (scRNA)-sequenced 56,771 endothelial cells from human/mouse (peri)-tumoral lung and cultured human lung TECs, and detected 17 known and 16 previously unrecognized phenotypes, including TECs putatively regulating immune surveillance. We resolved the canonical tip TECs into a known migratory tip and a putative basement-membrane remodeling breach phenotype. Tip TEC signatures correlated with patient survival, and tip/breach TECs were most sensitive to vascular endothelial growth factor blockade. Only tip TECs were congruent across species/models and shared conserved markers. Integrated analysis of the scRNA-sequenced data with orthogonal multi-omics and meta-analysis data across different human tumors, validated by functional analysis, identified collagen modification as a candidate angiogenic pathway.

Post-transcriptional regulation in hematopoiesis: RNA binding proteins take control 1.

Normal hematopoiesis is sustained through a carefully orchestrated balance between hematopoietic stem cell (HSC) self-renewal and differentiation. The functional importance of this axis is underscored by the severity of disease phenotypes initiated by abnormal HSC function, including myelodysplastic syndromes and hematopoietic malignancies. Major advances in the understanding of transcriptional regulation of primitive hematopoietic cells have been achieved; however, the post-transcriptional regulatory layer that may impinge on their behavior remains underexplored by comparison. Key players at this level include RNA-binding proteins (RBPs), which execute precise and highly coordinated control of gene expression through modulation of RNA properties that include its splicing, polyadenylation, localization, degradation, or translation. With the recent identification of RBPs having essential roles in regulating proliferation and cell fate decisions in other systems, there has been an increasing appreciation of the importance of post-transcriptional control at the stem cell level. Here we discuss our current understanding of RBP-driven post-transcriptional regulation in HSCs, its implications for normal, perturbed, and malignant hematopoiesis, and the most recent technological innovations aimed at RBP-RNA network characterization at the systems level. Emerging evidence highlights RBP-driven control as an underappreciated feature of primitive hematopoiesis, the greater understanding of which has important clinical implications.

Poor-prognosis colon cancer is defined by a molecularly distinct subtype and develops from serrated precursor lesions.

Colon cancer is a clinically diverse disease. This heterogeneity makes it difficult to determine which patients will benefit most from adjuvant therapy and impedes the development of new targeted agents. More insight into the biological diversity of colon cancers, especially in relation to clinical features, is therefore needed. We demonstrate, using an unsupervised classification strategy involving over 1,100 individuals with colon cancer, that three main molecularly distinct subtypes can be recognized. Two subtypes have been previously identified and are well characterized (chromosomal-instable and microsatellite-instable cancers). The third subtype is largely microsatellite stable and contains relatively more CpG island methylator phenotype-positive carcinomas but cannot be identified on the basis of characteristic mutations. We provide evidence that this subtype relates to sessile-serrated adenomas, which show highly similar gene expression profiles, including upregulation of genes involved in matrix remodeling and epithelial-mesenchymal transition. The identification of this subtype is crucial, as it has a very unfavorable prognosis and, moreover, is refractory to epidermal growth factor receptor-targeted therapy.