Single-Cell RNA Sequencing Reveals Repair Features of Human Umbilical Cord Mesenchymal Stromal Cells

RATIONALE

The chronic lung disease bronchopulmonary dysplasia (BPD) is the most severe complication of extreme prematurity. BPD results in impaired lung alveolar and vascular development and long-term respiratory morbidity, for which only supportive therapies exist. Umbilical cord-derived mesenchymal stromal cells (UC-MSCs) improve lung structure and function in experimental BPD. Results of clinical trials with MSCs for many disorders do not yet match the promising preclinical studies. A lack of specific criteria to define functionally distinct MSCs persists.

OBJECTIVES

To determine and correlate single-cell UC-MSC transcriptomic profile with therapeutic potential.

METHODS

UC-MSCs from five term donors and human neonatal dermal fibroblasts (HNDFs, control cells of mesenchymal origin) transcriptomes were investigated by single-cell RNA sequencing analysis (scRNA-seq). The lung-protective effect of UC-MSCs with a distinct transcriptome and control HNDFs was tested in vivo in hyperoxia-induced neonatal lung injury in rats.

MEASUREMENTS AND MAIN RESULTS

UC-MSCs showed limited transcriptomic heterogeneity, but were different from HNDFs. Gene ontology enrichment analysis revealed distinct - progenitor-like and fibroblast-like - UC-MSC subpopulations. Only the treatment with progenitor-like UC-MSCs improved lung function and structure and attenuated pulmonary hypertension in hyperoxia-exposed rat pups. Moreover, scRNA-seq identified major histocompatibility complex class I as a molecular marker of non-therapeutic cells and associated with decreased lung retention.

CONCLUSIONS

UC-MSCs with a progenitor-like transcriptome, but not with a fibroblast-like transcriptome, provide lung protection in experimental BPD. High expression of major histocompatibility complex class I is associated with reduced therapeutic benefit. scRNA-seq may be useful to identify subsets of MSCs with superior repair capacity for clinical application.

FGL2 promotes tumour growth and attenuates infiltration of activated immune cells in melanoma and ovarian cancer models

The tumour microenvironment is infiltrated by immunosuppressive cells, such as regulatory T cells (Tregs), which contribute to tumour escape and impede immunotherapy outcomes. Soluble fibrinogen-like protein 2 (sFGL2), a Treg effector protein, inhibits immune cell populations, via receptors FcγRIIB and FcγRIII, leading to downregulation of CD86 in antigen presenting cells and limiting T cell activation. Increased FGL2 expression is associated with tumour progression and poor survival in several different cancers, such as glioblastoma multiforme, lung, renal, liver, colorectal, and prostate cancer. Querying scRNA-seq human cancer data shows FGL2 is produced by cells in the tumour microenvironment (TME), particularly monocytes and macrophages as well as T cells and dendritic cells (DCs), while cancer cells have minimal expression of FGL2. We studied the role of FGL2 exclusively produced by cells in the TME, by leveraging Fgl2 knockout mice. We tested two murine models of cancer in which the role of FGL2 has not been previously studied: epithelial ovarian cancer and melanoma. We show that absence of FGL2 leads to a more activated TME, including activated DCs (CD86+, CD40+) and T cells (CD25+, TIGIT+), as well as demonstrating for the first time that the absence of FGL2 leads to more activated natural killer cells (DNAM-1+, NKG2D+) in the TME. Furthermore, the absence of FGL2 leads to prolonged survival in the B16F10 melanoma model, while the absence of FGL2 synergizes with oncolytic virus to prolong survival in the ID8-p53-/-Brca2-/- ovarian cancer model. In conclusion, targeting FGL2 is a promising cancer treatment strategy alone and in combination immunotherapies.

Mesenchymal ovarian cancer cells promote CD8+ T cell exhaustion through the LGALS3-LAG3 axis

Cancer cells often metastasize by undergoing an epithelial-mesenchymal transition (EMT). Although abundance of CD8+ T-cells in the tumor microenvironment correlates with improved survival, mesenchymal cancer cells acquire greater resistance to antitumor immunity in some cancers. We hypothesized the EMT modulates the immune response to ovarian cancer. Here we show that cancer cells from infiltrated/inflamed tumors possess more mesenchymal cells, than excluded and desert tumors. We also noted high expression of LGALS3 is associated with EMT in vivo, a finding validated with in vitro EMT models. Dissecting the cellular communications among populations in the tumor revealed that mesenchymal cancer cells in infiltrated tumors communicate through LGALS3 to LAG3 receptor expressed by CD8+ T cells. We found CD8+ T cells express high levels of LAG3, a marker of T cell exhaustion. The results indicate that EMT in ovarian cancer cells promotes interactions between cancer cells and T cells through the LGALS3 - LAG3 axis, which could increase T cell exhaustion in infiltrated tumors, dampening antitumor immunity.

Comparative analysis of syngeneic mouse models of high-grade serous ovarian cancer

Ovarian cancers exhibit high rates of recurrence and poor treatment response. Preclinical models that recapitulate human disease are critical to develop new therapeutic approaches. Syngeneic mouse models allow for the generation of tumours comprising the full repertoire of non-malignant cell types but have expanded in number, varying in the cell type of origin, method for transformation, and ultimately, the properties of the tumours they produce. Here we have performed a comparative analysis of high-grade serous ovarian cancer models based on transcriptomic profiling of 22 cell line models, and intrabursal and intraperitoneal tumours from 12. Among cell lines, we identify distinct signalling activity, such as elevated inflammatory signalling in STOSE and OVE16 models, and MAPK/ERK signalling in ID8 and OVE4 models; metabolic differences, such as reduced glycolysis-associated expression in several engineered ID8 subclones; and relevant functional properties, including differences in EMT activation, PD-L1 and MHC class I expression, and predicted chemosensitivity. Among tumour samples, we observe increased variability and stromal content among intrabursal tumours. Finally, we predict differences in the microenvironment of ID8 models engineered with clinically relevant mutations. We anticipate that this work will serve as a valuable resource, providing new insight to help select models for specific experimental objectives.

Single-cell transcriptomic atlas of lung microvascular regeneration after targeted endothelial cell ablation

We sought to define the mechanism underlying lung microvascular regeneration in a model of severe acute lung injury (ALI) induced by selective lung endothelial cell ablation. Intratracheal instillation of DT in transgenic mice expressing human diphtheria toxin (DT) receptor targeted to ECs resulted in ablation of >70% of lung ECs, producing severe ALI with near complete resolution by 7 days. Using single-cell RNA sequencing, eight distinct endothelial clusters were resolved, including alveolar aerocytes (aCap) ECs expressing apelin at baseline and general capillary (gCap) ECs expressing the apelin receptor. At 3 days post-injury, a novel gCap EC population emerged characterized by de novo expression of apelin, together with the stem cell marker, protein C receptor. These stem-like cells transitioned at 5 days to proliferative endothelial progenitor-like cells, expressing apelin receptor together with the pro-proliferative transcription factor, Foxm1, and were responsible for the rapid replenishment of all depleted EC populations by 7 days post-injury. Treatment with an apelin receptor antagonist prevented ALI resolution and resulted in excessive mortality, consistent with a central role for apelin signaling in EC regeneration and microvascular repair. The lung has a remarkable capacity for microvasculature EC regeneration which is orchestrated by newly emergent apelin-expressing gCap endothelial stem-like cells that give rise to highly proliferative, apelin receptor-positive endothelial progenitors responsible for the regeneration of the lung microvasculature.

Prenatal low-dose methylmercury exposure causes premature neuronal differentiation and autism-like behaviors in a rodent model

Aberrant neurodevelopment is a core deficit of autism spectrum disorder (ASD). Here we ask whether a non-genetic factor, prenatal exposure to the environmental pollutant methylmercury (MeHg), is a contributing factor in ASD onset. We showed that adult mice prenatally exposed to non-apoptotic MeHg exhibited key ASD characteristics, including impaired communication, reduced sociability, and increased restrictive repetitive behaviors, whereas in the embryonic cortex, prenatal MeHg exposure caused premature neuronal differentiation. Further single-cell RNA sequencing (scRNA-seq) analysis disclosed that prenatal exposure to MeHg resulted in cortical radial glial precursors (RGPs) favoring asymmetric differentiation to directly generate cortical neurons, omitting the intermediate progenitor stage. In addition, MeHg exposure in cultured RGPs increased CREB phosphorylation and enhanced the interaction between CREB and CREB binding protein (CBP). Intriguingly, metformin, an FDA-approved drug, can reverse MeHg-induced premature neuronal differentiation via CREB/CBP repulsion. These findings provide insights into ASD etiology, its underlying mechanism, and a potential therapeutic strategy.

Absolute scaling of single-cell transcriptomes identifies pervasive hypertranscription in adult stem and progenitor cells

Hypertranscription supports biosynthetically demanding cellular states through global transcriptome upregulation. Despite its potential widespread relevance, documented examples of hypertranscription remain few and limited to early development. Here, we demonstrate that absolute scaling of single-cell RNA-sequencing data enables the estimation of total transcript abundances per cell. We validate absolute scaling in known cases of developmental hypertranscription and apply it to adult cell types, revealing a remarkable dynamic range in transcriptional output. In adult organs, hypertranscription marks activated stem/progenitor cells with multilineage potential and is redeployed in conditions of tissue injury, where it precedes bursts of proliferation during regeneration. Our analyses identify a common set of molecular pathways associated with both adult and embryonic hypertranscription, including chromatin remodeling, DNA repair, ribosome biogenesis, and translation. These shared features across diverse cell contexts support hypertranscription as a general and dynamic cellular program that is pervasively employed during development, organ maintenance, and regeneration.

Metformin prevents age-associated ovarian fibrosis by modulating the immune landscape in female mice

Ovarian fibrosis is a pathological condition associated with aging and is responsible for a variety of ovarian dysfunctions. Given the known contributions of tissue fibrosis to tumorigenesis, it is anticipated that ovarian fibrosis may contribute to ovarian cancer risk. We recently reported that diabetic postmenopausal women using metformin had ovarian collagen abundance and organization that were similar to premenopausal ovaries from nondiabetic women. In this study, we investigated the effects of aging and metformin on mouse ovarian fibrosis at a single-cell level. We discovered that metformin treatment prevented age-associated ovarian fibrosis by modulating the proportion of fibroblasts, myofibroblasts, and immune cells. Senescence-associated secretory phenotype (SASP)-producing fibroblasts increased in aged ovaries, and a unique metformin-responsive subpopulation of macrophages emerged in aged mice treated with metformin. The results demonstrate that metformin can modulate specific populations of immune cells and fibroblasts to prevent age-associated ovarian fibrosis and offers a new strategy to prevent ovarian fibrosis.

The Tumor Immune Profile of Murine Ovarian Cancer Models: An Essential Tool For Ovarian Cancer Immunotherapy Research

Epithelial ovarian cancer (EOC) is the most lethal gynecologic cancer with an imperative need for new treatments. Immunotherapy has had marked success in some cancer types; however, clinical trials studying the efficacy of immune checkpoint inhibitors for the treatment of EOC benefited less than 15% of patients. Given that EOC develops from multiple tissues in the reproductive system and metastasizes widely throughout the peritoneal cavity, responses to immunotherapy are likely hindered by heterogeneous tumor microenvironments (TME) containing a variety of immune profiles. To fully characterize and compare syngeneic model systems that may reflect this diversity, we determined the immunogenicity of six ovarian tumor models in vivo, the T and myeloid profile of orthotopic tumors and the immune composition and cytokine profile of ascites, by single-cell RNA sequencing, flow cytometry and IHC. The selected models reflect the different cellular origins of EOC (ovarian and fallopian tube epithelium) and harbor mutations relevant to human disease, including Tp53 mutation, PTEN suppression, and constitutive KRAS activation. ID8-p53-/- and ID8-C3 tumors were most highly infiltrated by T cells, whereas STOSE and MOE-PTEN/KRAS tumors were primarily infiltrated by tumor associated macrophages and were unique in MHC class I and II expression. MOE-PTEN/KRAS tumors were capable of forming T cell clusters. This panel of well-defined murine EOC models reflects some of the heterogeneity found in human disease and can serve as a valuable resource for studies that aim to test immunotherapies, explore the mechanisms of immune response to therapy, and guide selection of treatments for patient populations.

When killers become thieves: Trogocytosed PD-1 inhibits NK cells in cancer

Trogocytosis modulates immune responses, with still unclear underlying molecular mechanisms. Using leukemia mouse models, we found that lymphocytes perform trogocytosis at high rates with tumor cells. While performing trogocytosis, both Natural Killer (NK) and CD8+ T cells acquire the checkpoint receptor PD-1 from leukemia cells. In vitro and in vivo investigation revealed that PD-1 on the surface of NK cells, rather than being endogenously expressed, was derived entirely from leukemia cells in a SLAM receptor-dependent fashion. PD-1 acquired via trogocytosis actively suppressed NK cell antitumor immunity. PD-1 trogocytosis was corroborated in patients with clonal plasma cell disorders, where NK cells that stained for PD-1 also stained for tumor cell markers. Our results, in addition to shedding light on a previously unappreciated mechanism underlying the presence of PD-1 on NK and cytotoxic T cells, reveal the immunoregulatory effect of membrane transfer occurring when immune cells contact tumor cells.

VIVA1: a more invasive subclone of MDA-MB-134VI invasive lobular carcinoma cells with increased metastatic potential in xenograft models

BACKGROUND

Invasive lobular carcinoma (ILC) is the second most common type of breast cancer. As few tools exist to study ILC metastasis, we isolated ILC cells with increased invasive properties to establish a spontaneously metastasising xenograft model.

METHODS

MDA-MB-134VI ILC cells were placed in transwells for 7 days. Migrated cells were isolated and expanded to create the VIVA1 cell line. VIVA1 cells were compared to parental MDA-MB-134VI cells in vitro for ILC marker expression and relative proliferative and invasive ability. An intraductally injected orthotopic xenograft model was used to assess primary and metastatic tumour growth in vivo.

RESULTS

Similar to MDA-MB-134VI, VIVA1 cells retained expression of oestrogen receptor (ER) and lacked expression of E-cadherin, however showed increased invasion in vitro. Following intraductal injection, VIVA1 and MDA-MB-134VI cells had similar primary tumour growth and survival kinetics. However, macrometastases were apparent in 7/10 VIVA1-injected animals. Cells from a primary orthotopic tumour (VIVA-LIG43) were isolated and showed similar proliferative rates but were also more invasive than parental cells. Upon re-injection intraductally, VIVA-LIG43 cells had more rapid tumour growth with similar metastatic incidence and location.

CONCLUSIONS

We generated a new orthotopic spontaneously metastasising xenograft model for ER+ ILC amenable for the study of ILC metastasis.

Single-Cell RNA Sequencing-Based Characterization of Resident Lung Mesenchymal Stromal Cells in Bronchopulmonary Dysplasia

Late lung development is a period of alveolar and microvascular formation, which is pivotal in ensuring sufficient and effective gas exchange. Defects in late lung development manifest in premature infants as a chronic lung disease named bronchopulmonary dysplasia (BPD). Numerous studies demonstrated the therapeutic properties of exogenous bone marrow and umbilical cord-derived mesenchymal stromal cells (MSCs) in experimental BPD. However, very little is known regarding the regenerative capacity of resident lung MSCs (L-MSCs) during normal development and in BPD. In this study we aimed to characterize the L-MSC population in homeostasis and upon injury. We used single-cell RNA sequencing (scRNA-seq) to profile in situ Ly6a+ L-MSCs in the lungs of normal and O2-exposed neonatal mice (a well-established model to mimic BPD) at 3 developmental timepoints (postnatal days 3, 7, and 14). Hyperoxia exposure increased the number and altered the expression profile of L-MSCs, particularly by increasing the expression of multiple pro-inflammatory, pro-fibrotic, and anti-angiogenic genes. In order to identify potential changes induced in the L-MSCs transcriptome by storage and culture, we profiled 15 000 Ly6a+ L-MSCs after in vitro culture. We observed great differences in expression profiles of in situ and cultured L-MSCs, particularly those derived from healthy lungs. Additionally, we have identified the location of Ly6a+/Col14a1+ L-MSCs in the developing lung and propose Serpinf1 as a novel, culture-stable marker of L-MSCs. Finally, cell communication analysis suggests inflammatory signals from immune and endothelial cells as main drivers of hyperoxia-induced changes in L-MSCs transcriptome.

Transcriptional census of epithelial-mesenchymal plasticity in cancer

Epithelial-mesenchymal plasticity (EMP) contributes to tumor progression, promoting therapy resistance and immune cell evasion. Definitive molecular features of this plasticity have largely remained elusive due to the limited scale of most studies. Leveraging single-cell RNA sequencing data from 266 tumors spanning eight different cancer types, we identify expression patterns associated with intratumoral EMP. Integrative analysis of these programs confirmed a high degree of diversity among tumors. These diverse programs are associated with combinations of various common regulatory mechanisms initiated from cues within the tumor microenvironment. We show that inferring regulatory features can inform effective therapeutics to restrict EMP.

Periostin gene expression in neu-positive breast cancer cells is regulated by a FGFR signaling cross talk with TGFβ/PI3K/AKT pathways

Background

Breast cancer is a highly heterogeneous disease with multiple drivers and complex regulatory networks. Periostin (Postn) is a matricellular protein involved in a plethora of cancer types and other diseases. Postn has been shown to be involved in various processes of tumor development, such as angiogenesis, invasion, cell survival and metastasis. The expression of Postn in breast cancer cells has been correlated with a more aggressive phenotype. Despite extensive research, it remains unclear how epithelial cancer cells regulate Postn expression.

Methods

Using murine tumor models and human TMAs, we have assessed the proportion of tumor samples that have acquired Postn expression in tumor cells. Using biochemical approaches and tumor cell lines derived from Neu+ murine primary tumors, we have identified major regulators of Postn gene expression in breast cancer cell lines.

Results

Here, we show that, while the stromal compartment typically always expresses Postn, about 50% of breast tumors acquire Postn expression in the epithelial tumor cells. Furthermore, using an in vitro model, we show a cross-regulation between FGFR, TGFβ and PI3K/AKT pathways to regulate Postn expression. In HER2-positive murine breast cancer cells, we found that basic FGF can repress Postn expression through a PKC-dependent pathway, while TGFβ can induce Postn expression in a SMAD-independent manner. Postn induction following the removal of the FGF-suppressive signal is dependent on PI3K/AKT signaling.

Conclusion

Overall, these results reveal a novel regulatory mechanism and shed light on how breast tumor cells acquire Postn expression. This complex regulation is likely to be cell type and cancer specific as well as have important therapeutic implications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13058-021-01487-8.

Identification of Rac guanine nucleotide exchange factors promoting Lgl1 phosphorylation in glioblastoma

The protein Lgl1 is a key regulator of cell polarity. We previously showed that Lgl1 is inactivated by hyperphosphorylation in glioblastoma as a consequence of PTEN tumour suppressor loss and aberrant activation of the PI 3-kinase pathway; this contributes to glioblastoma pathogenesis both by promoting invasion and repressing glioblastoma cell differentiation. Lgl1 is phosphorylated by atypical protein kinase C that has been activated by binding to a complex of the scaffolding protein Par6 and active, GTP-bound Rac. The specific Rac guanine nucleotide exchange factors that generate active Rac to promote Lgl1 hyperphosphorylation in glioblastoma are unknown. We used CRISPR/Cas9 to knockout PREX1, a PI 3-kinase pathway-responsive Rac guanine nucleotide exchange factor, in patient-derived glioblastoma cells. Knockout cells had reduced Lgl1 phosphorylation, which was reversed by re-expressing PREX1. They also had reduced motility and an altered phenotype suggestive of partial neuronal differentiation; consistent with this, RNA-seq analyses identified sets of PREX1-regulated genes associated with cell motility and neuronal differentiation. PREX1 knockout in glioblastoma cells from a second patient did not affect Lgl1 phosphorylation. This was due to overexpression of a short isoform of the Rac guanine nucleotide exchange factor TIAM1; knockdown of TIAM1 in these PREX1 knockout cells reduced Lgl1 phosphorylation. These data show that PREX1 links aberrant PI 3-kinase signaling to Lgl1 phosphorylation in glioblastoma, but that TIAM1 is also to fill this role in a subset of patients. This redundancy between PREX1 and TIAM1 is only partial, as motility was impaired in PREX1 knockout cells from both patients.

Transcriptional heterogeneity of stemness phenotypes in the ovarian epithelium

The ovarian surface epithelium (OSE) is a monolayer of epithelial cells surrounding the ovary that ruptures during each ovulation to allow release of the oocyte. This wound is quickly repaired, but mechanisms promoting repair are poorly understood. The contribution of tissue-resident stem cells in the homeostasis of several epithelial tissues is widely accepted, but their involvement in OSE is unclear. We show that traits associated with stem cells can be increased following exposure to the cytokine TGFB1, overexpression of the transcription factor Snai1, or deletion of Brca1. We find that stemness is often linked to mesenchymal-associated gene expression and higher activation of ERK signalling, but is not consistently dependent on their activation. Expression profiles of these populations are extremely context specific, suggesting that stemness may not be associated with a single, distinct population, but rather is a heterogeneous cell state that may emerge from diverse environmental cues. These findings support that the OSE may not require distinct stem cells for long-term maintenance, and may instead achieve this through transient dedifferentiation into a stem-like state.

Using spheroids as a model, the authors report the molecular signatures of ovarian surface epithelial (OSE) cells exhibiting stemness phenotype under various conditions in vitro, and found that there is a large degree of heterogeneity in the transcriptional profiles of stem cells induced under different conditions. They suggest that maintenance of the OSE may not require a single stem cell population, but heterogeneous stem cells that can be induced transiently under diverse environmental cues.

Identification of novel lipid droplet factors that regulate lipophagy and cholesterol efflux in macrophage foam cells

Macrophage autophagy is a highly anti-atherogenic process that promotes the catabolism of cytosolic lipid droplets (LDs) to maintain cellular lipid homeostasis. Selective autophagy relies on tags such as ubiquitin and a set of selectivity factors including selective autophagy receptors (SARs) to label specific cargo for degradation. Originally described in yeast cells, "lipophagy" refers to the degradation of LDs by autophagy. Yet, how LDs are targeted for autophagy is poorly defined. Here, we employed mass spectrometry to identify lipophagy factors within the macrophage foam cell LD proteome. In addition to structural proteins (e.g., PLIN2), metabolic enzymes (e.g., ACSL) and neutral lipases (e.g., PNPLA2), we found the association of proteins related to the ubiquitination machinery (e.g., AUP1) and autophagy (e.g., HMGB, YWHA/14-3-3 proteins). The functional role of candidate lipophagy factors (a total of 91) was tested using a custom siRNA array combined with high-content cholesterol efflux assays. We observed that knocking down several of these genes, including Hmgb1, Hmgb2, Hspa5, and Scarb2, significantly reduced cholesterol efflux, and SARs SQSTM1/p62, NBR1 and OPTN localized to LDs, suggesting a role for these in lipophagy. Using yeast lipophagy assays, we established a genetic requirement for several candidate lipophagy factors in lipophagy, including HSPA5, UBE2G2 and AUP1. Our study is the first to systematically identify several LD-associated proteins of the lipophagy machinery, a finding with important biological and therapeutic implications. Targeting these to selectively enhance lipophagy to promote cholesterol efflux in foam cells may represent a novel strategy to treat atherosclerosis.Abbreviations:ADGRL3: adhesion G protein-coupled receptor L3; agLDL: aggregated low density lipoprotein; AMPK: AMP-activated protein kinase; APOA1: apolipoprotein A1; ATG: autophagy related; AUP1: AUP1 lipid droplet regulating VLDL assembly factor; BMDM: bone-marrow derived macrophages; BNIP3L: BCL2/adenovirus E1B interacting protein 3-like; BSA: bovine serum albumin; CALCOCO2: calcium binding and coiled-coil domain 2; CIRBP: cold inducible RNA binding protein; COLGALT1: collagen beta(1-O)galactosyltransferase 1; CORO1A: coronin 1A; DMA: deletion mutant array; Faa4: long chain fatty acyl-CoA synthetase; FBS: fetal bovine serum; FUS: fused in sarcoma; HMGB1: high mobility group box 1; HMGB2: high mobility group box 2: HSP90AA1: heat shock protein 90: alpha (cytosolic): class A member 1; HSPA5: heat shock protein family A (Hsp70) member 5; HSPA8: heat shock protein 8; HSPB1: heat shock protein 1; HSPH1: heat shock 105kDa/110kDa protein 1; LDAH: lipid droplet associated hydrolase; LIPA: lysosomal acid lipase A; LIR: LC3-interacting region; MACROH2A1: macroH2A.1 histone; MAP1LC3: microtubule-associated protein 1 light chain 3; MCOLN1: mucolipin 1; NBR1: NBR1, autophagy cargo receptor; NPC2: NPC intracellular cholesterol transporter 2; OPTN: optineurin; P/S: penicillin-streptomycin; PLIN2: perilipin 2; PLIN3: perilipin 3; PNPLA2: patatin like phospholipase domain containing 2; RAB: RAB, member RAS oncogene family; RBBP7, retinoblastoma binding protein 7, chromatin remodeling factor; SAR: selective autophagy receptor; SCARB2: scavenger receptor class B, member 2; SGA: synthetic genetic array; SQSTM1: sequestosome 1; TAX1BP1: Tax1 (human T cell leukemia virus type I) binding protein 1; TFEB: transcription factor EB; TOLLIP: toll interacting protein; UBE2G2: ubiquitin conjugating enzyme E2 G2; UVRAG: UV radiation resistance associated gene; VDAC2: voltage dependent anion channel 2; VIM: vimentin.

COX2 is induced in the ovarian epithelium during ovulatory wound repair and promotes cell survival†

The ovarian surface epithelium (OSE) is a monolayer of cells surrounding the ovary that is ruptured during ovulation. After ovulation, the wound is repaired, however, this process is poorly understood. In epithelial tissues, wound repair is mediated by an epithelial-to-mesenchymal transition (EMT). Transforming Growth Factor Beta-1 (TGFβ1) is a cytokine commonly known to induce an EMT and is present throughout the ovarian microenvironment. We, therefore, hypothesized that TGFβ1 induces an EMT in OSE cells and activates signaling pathways important for wound repair. Treating primary cultures of mouse OSE cells with TGFβ1 induced an EMT mediated by TGFβRI signaling. The transcription factor Snail was the only EMT-associated transcription factor increased by TGFβ1 and, when overexpressed, was shown to increase OSE cell migration. A polymerase chain reaction array of TGFβ signaling targets determined Cyclooxygenase-2 (Cox2) to be most highly induced by TGFβ1. Constitutive Cox2 expression modestly increased migration and robustly enhanced cell survival, under stress conditions similar to those observed during wound repair. The increase in Snail and Cox2 expression with TGFβ1 was reproduced in human OSE cultures, suggesting these responses are conserved between mouse and human. Finally, the induction of Cox2 expression in OSE cells during ovulatory wound repair was shown in vivo, suggesting TGFβ1 increases Cox2 to promote wound repair by enhancing cell survival. These data support that TGFβ1 promotes ovulatory wound repair by induction of an EMT and activation of a COX2-mediated pro-survival pathway. Understanding ovulatory wound repair may give insight into why ovulation is the primary non-hereditary risk factor for ovarian cancer.

RIPK1 gene variants associate with obesity in humans and can be therapeutically silenced to reduce obesity in mice

Obesity is a major public health burden worldwide and is characterized by chronic low-grade inflammation driven by the cooperation of the innate immune system and dysregulated metabolism in adipose tissue and other metabolic organs. Receptor-interacting serine/threonine-protein kinase 1 (RIPK1) is a central regulator of inflammatory cell function that coordinates inflammation, apoptosis and necroptosis in response to inflammatory stimuli. Here we show that genetic polymorphisms near the human RIPK1 locus associate with increased RIPK1 gene expression and obesity. We show that one of these single nucleotide polymorphisms is within a binding site for E4BP4 and increases RIPK1 promoter activity and RIPK1 gene expression in adipose tissue. Therapeutic silencing of RIPK1 in vivo in a mouse model of diet-induced obesity dramatically reduces fat mass, total body weight and improves insulin sensitivity, while simultaneously reducing macrophage and promoting invariant natural killer T cell accumulation in adipose tissue. These findings demonstrate that RIPK1 is genetically associated with obesity, and reducing RIPK1 expression is a potential therapeutic approach to target obesity and related diseases.

Publisher Correction: RIPK1 gene variants associate with obesity in humans and can be therapeutically silenced to reduce obesity in mice

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Ovarian cancer and the evolution of subtype classifications using transcriptional profiling†

Ovarian cancer is a complex disease with multiple subtypes, each having distinct histopathologies and variable responses to treatment. This review highlights the technological milestones and the studies that have applied them to change our definitions of ovarian cancer. Over the past 50 years, technologies such as microarrays and next-generation sequencing have led to the discovery of molecular alterations that define each of the ovarian cancer subtypes and has enabled further subclassification of the most common subtype, high-grade serous ovarian cancer (HGSOC). Improvements in mutational profiling have provided valuable insight, such as the ubiquity of TP53 mutations in HGSOC tumors. However, the information derived from these technological advances has also revealed the immense heterogeneity of this disease, from variation between patients to compositional differences within single masses. In looking forward, the emerging technologies for single-cell and spatially resolved transcriptomics will allow us to better understand the cellular composition and structure of tumors and how these contribute to the molecular subtypes. Attempts to incorporate the complexities ovarian cancer has resulted in increasing sophistication of model systems, and the increased precision in molecular profiling of ovarian cancers has already led to the introduction of inhibitors of poly (ADP-ribose) polymerases as a new class of treatments for ovarian cancer with DNA repair deficiencies. Future endeavors to define increasingly accurate classification strategies for ovarian cancer subtypes will allow for confident prediction of disease progression and provide important insight into potentially targetable molecular mechanisms specific to each subtype.

The Canadian Cancer Research Conference 2019

The 5th Canadian Cancer Research Conference (ccrc) took place 3-5 November 2019 in Ottawa, Ontario. Nearly 1000 participants-scientists, oncologists, community members, and patients-gathered to share knowledge, foster collaboration, and fuel the future of cancer research in Canada. The scientific program included 3 plenary sessions, 26 concurrent sessions, and 2 poster sessions presenting research described in more than 600 submitted abstracts, giving participants the opportunity to share health research that collectively encompassed the 4 pillars recognized by the Canadian Institutes of Health Research. In addition to the breadth of topics addressed by Canadian and international experts, the highlights of the meeting included the integration of patients and patient advocates, new rapid-fire sessions for abstract presentation, and events that enhanced learning opportunities for trainees.

Abstract P2-04-05: Loss of Ste20-like kinase induces a basal/stem-like phenotype in HER2-positive breast cancers

Over the past two decades, the mammalian Ste20-like kinase (SLK) has been characterized for its role in regulating cellular migration, proliferation and apoptosis in fibroblasts and myoblasts. In mammary epithelial cells, SLK has been shown to be required for efficient epithelial-to-mesenchymal transition and to be activated downstream of the HER2/Neu-oncogene to control chemotactic cellular migration. Here, we assessed the role of SLK in HER2/Neu-induced mammary tumorigenesis in vivo. As SLK is activated downstream of HER2/Neu, we hypothesized that the loss of SLK would significantly delay tumor progression in a mouse model of HER2-positive breast cancer. As we have shown that global attenuation of SLK kinase activity results in embryonic lethality, a conditional SLK knockout mouse model was generated. To study the role of SLK in HER2-positive breast cancer, we crossed these conditional SLK knockout mice with mice expressing HER2/Neu linked to Cre recombinase in the mammary luminal epithelium. Unexpectedly, we have demonstrated that conditional deletion of SLK significantly accelerates Neu-induced mammary tumor onset and decreases overall survival. SLK deletion results in the induction of Sox10 which drives mammary stem/progenitor activity and cooperates with HER2/Neu to drive tumor growth. Using the Cancer Genome Atlas, we have supported previous findings and validated Sox10 as a potential biomarker of the Triple-negative Breast Cancer subtype. Furthermore, we have uncovered that SLK deletion results in enhanced activation of AKT. We provide evidence that Sox10 induction requires signaling through a novel AKT/Sox9-dependent pathway following SLK deletion. Taken together, our data suggests that SLK may be required to maintain cells in a fully differentiated state and that loss of SLK in HER2/Neu-induced breast cancer drives a more basal/stem-like phenotype through the induction of Sox10.

Citation Format: Khalid N Al-Zahrani, John Abou-Hamad, Ben Pryce, David P Cook, Johnathan Hodgins, Cedrik Labreche, Pascale Robineau-Charette, Christiano de Souza, John Bell, Rebecca Auer, Michele Ardolino, Barbara Vanderhyden, Daniel Schramek, Luc Sabourin. Loss of Ste20-like kinase induces a basal/stem-like phenotype in HER2-positive breast cancers [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P2-04-05.

Metformin Abrogates Age-Associated Ovarian Fibrosis

PURPOSE

The ovarian cancer risk factors of age and ovulation are curious because ovarian cancer incidence increases in postmenopausal women, long after ovulations have ceased. To determine how age and ovulation underlie ovarian cancer risk, we assessed the effects of these risk factors on the ovarian microenvironment.

EXPERIMENTAL DESIGN

Aged C57/lcrfa mice (0-33 months old) were generated to assess the aged ovarian microenvironment. To expand our findings into human aging, we assembled a cohort of normal human ovaries (n = 18, 21-71 years old). To validate our findings, an independent cohort of normal human ovaries was assembled (n = 9, 41-82 years old).

RESULTS

We first validated the presence of age-associated murine ovarian fibrosis. Using interdisciplinary methodologies, we provide novel evidence that ovarian fibrosis also develops in human postmenopausal ovaries across two independent cohorts (n = 27). Fibrotic ovaries have an increased CD206⁺:CD68⁺ cell ratio, CD8⁺ T-cell infiltration, and profibrotic DPP4⁺αSMA⁺ fibroblasts. Metformin use was associated with attenuated CD8⁺ T-cell infiltration and reduced CD206⁺:CD68⁺ cell ratio.

CONCLUSIONS

These data support a novel hypothesis that unifies the primary nonhereditary ovarian cancer risk factors through the development of ovarian fibrosis and the formation of a premetastatic niche, and suggests a potential use for metformin in ovarian cancer prophylaxis.See related commentary by Madariaga et al., p. 523.

Single-cell RNA-sequencing reveals transcriptional dynamics of estrogen-induced dysplasia in the ovarian surface epithelium

Estrogen therapy increases the risk of ovarian cancer and exogenous estradiol accelerates the onset of ovarian cancer in mouse models. Both in vivo and in vitro, ovarian surface epithelial (OSE) cells exposed to estradiol develop a subpopulation that loses cell polarity, contact inhibition, and forms multi-layered foci of dysplastic cells with increased susceptibility to transformation. Here, we use single-cell RNA-sequencing to characterize this dysplastic subpopulation and identify the transcriptional dynamics involved in its emergence. Estradiol-treated cells were characterized by up-regulation of genes associated with proliferation, metabolism, and survival pathways. Pseudotemporal ordering revealed that OSE cells occupy a largely linear phenotypic spectrum that, in estradiol-treated cells, diverges towards cell state consistent with the dysplastic population. This divergence is characterized by the activation of various cancer-associated pathways including an increase in Greb1 which was validated in fallopian tube epithelium and human ovarian cancers. Taken together, this work reveals possible mechanisms by which estradiol increases epithelial cell susceptibility to tumour initiation.

Assessing the efficacy of androgen receptor and Sox10 as independent markers of the triple-negative breast cancer subtype by transcriptome profiling

The Androgen Receptor (AR) has recently garnered a lot of attention as a potential biomarker and therapeutic target in hormone-dependent cancers, including breast cancer. However, several inconsistencies exist within the literature as to which subtypes of breast cancer express AR or whether it can be used to define its own unique subtype. Here, we analyze 1246 invasive breast cancer samples from the Cancer Genome Atlas and show that human breast cancers that have been subtyped based on their HER2, ESR1, or PGR expression contain four clusters of genes that are differentially expressed across all subtypes. We demonstrate that Sox10 is highly expressed in approximately one-third of all HER2/ESR1/PGR-low tumors and is a candidate biomarker of the triple-negative subtype. Although AR expression is acquired in many breast cancer cases, its expression could not define a unique subtype. Despite several reports stating that AR expression is acquired in HER2/ESR1/PGR triple-negative cancers, here we show that a low percentage of these cancers express AR (~20%). In contrast, AR is highly expressed in HER2-positive or ESR1/PGR-positive cancers (> 95%). Although AR expression cannot be used as an independent subtype biomarker, our analysis shows that routine evaluation of AR expression in tumors which express HER2, ESR1 and/or PGR may identify a unique subset of tumors which would benefit from anti-androgen based therapies.

The biomechanical characteristics of wearing FitFlop™ sandals highlight significant alterations in gait pattern: a comparative study

BACKGROUND

The net contribution of all muscles that act about a joint can be represented as an internal joint moment profile. This approach may be advantageous when studying footwear-induced perturbations during walking since the contribution of the smaller deeper muscles that cross the ankle joint cannot be evaluated with surface electromyography. Therefore, the present study aimed to advance the understanding of FitFlop™ footwear interaction by investigating lower extremity joint moment, and kinematic and centre of pressure profiles during gait.

METHODS

28 healthy participants performed 5 walking trials in 3 conditions: a FitFlop™ sandal, a conventional sandal and an athletic trainer. Three-dimensional ankle joint, and sagittal plane knee and hip joint moments, as well as corresponding kinematics and centre of pressure trajectories were evaluated.

FINDINGS

FitFlop™ differed significantly to both the conventional sandal and athletic trainer in: average anterior position of centre of pressure trajectory (P<0.0001) and peak hip extensor moment (P=0.001) during early stance; average medial position of centre of pressure trajectory during late stance; peak ankle dorsiflexion and corresponding range of motion; peak plantarflexor moment and total negative work performed at the ankle (all P<0.0001).

INTERPRETATION

The present findings demonstrate that FitFlop™ footwear significantly alters the gait pattern of wearers. An anterior displacement of the centre of pressure trajectory during early stance is the primary response to the destabilising effect of the mid-sole technology, and this leads to reductions in sagittal plane ankle joint range of motion and corresponding kinetics. Future investigations should consider the clinical implications of these findings.

Low-frequency accelerations over-estimate impact-related shock during walking

During gait, a failure to acknowledge the low-frequency component of a segmental acceleration signal will result in an overestimation of impact-related shock and may lead to inappropriately drawn conclusions. The present study was undertaken to investigate the significance of this low-frequency component in two distinctly different modalities of gait: barefoot (BF) and shod (SHOD) walking. Twenty-seven participants performed five walking trials at self-selected speed in each condition. Peak positive accelerations (PPA) at the shank and spine were first derived from the time-domain signal. The raw acceleration signals were then resolved in the frequency-domain and the active (low-frequency) and impact-related components of the power spectrum density (PSD) were quantified. PPA was significantly higher at the shank (P<0.0001) and spine (P=0.0007) in the BF condition. In contrast, no significant differences were apparent between conditions for shank (P=0.979) or spine (P=0.178) impact-related PSD when the low-frequency component was considered. This disparity between approaches was due to a significantly higher active PSD in both signals in the BF condition (P<0.0001; P=0.008, respectively), due to kinematic differences between conditions (P<0.05). These results indicate that the amplitude of the low-frequency component of an acceleration signal during gait is dependent on knee and ankle joint coordination behaviour, and highlight that impact-related shock is more accurately quantified in the frequency-domain following subtraction of this component.

Wide-pulse electrical stimulation to an intrinsic foot muscle induces acute functional changes in forefoot-rearfoot coupling behaviour during walking

Interventions for strengthening intrinsic foot muscles may be beneficial for rehabilitation from overuse injuries. In this study the acute effects of high-frequency, low-intensity wide-pulse electrical stimulation (WPS) over an intrinsic muscle on subsequent foot function during walking was assessed in healthy participants. WPS was delivered to the m. abductor hallucis (m.AH) of the non-dominant foot during relaxed standing. 3-dimensional forefoot (FF)--rearfoot (RF) coordination was quantified with a vector coding technique within separate periods of the stance phase to study WPS functional effects on foot motion. 4 types of coordinative strategies between the FF and RF were interpreted and compared PRE-to-POST-WPS for both the experimental and control feet. Bilateral electromyography (EMG) from m.AH was analysed during the intervention period for evidence of acute neuromuscular adaptation. The results showed that WPS significantly modulated FF-RF coordination during mid-stance, indicative of a more stable foot. Specifically, a statistically significant increase in FF eversion with concomitant RF inversion in the frontal plane and RF-dominated adduction in the transverse plane was observed. Subject-specific increases in post-stimulus m.AH EMG activation were observed but this was not reflected in an overall group effect. It is concluded that the structural integrity of the foot during walking is enhanced following an acute session of WPS and that this mechanical effect is most likely due to stimulation induced post-tetanic potentiation of synaptic transmission.

Effects of Olympic-style taekwondo kicks on an instrumented head-form and resultant injury measures

OBJECTIVE

The objective of this study was to assess the effect of taekwondo kicks and peak foot velocity (FVEL) on resultant head linear acceleration (RLA), head injury criterion (HIC15) and head velocity (HVEL).

METHODS

Each subject (n=12) randomly performed five repetitions of the turning kick (TK), clench axe kick (CA), front leg axe kick, jump back kick (JB) and jump spinning hook kick (JH) at the average standing head height for competitors in their weight division. A Hybrid II Crash Test Dummy head was fitted with a protective taekwondo helmet and instrumented with a triaxial accelerometer and fixed to a height-adjustable frame. Resultant head linear acceleration, HVEL, FVEL data were captured and processed using Qualysis Track Manager.

RESULTS

The TK (130.11 ± 51.67 g) produced a higher RLA than the CA (54.95 ± 20.08 g, p<0.001, d=1.84) and a higher HIC15 than the JH (672.74 ± 540.89 vs 300.19 ± 144.35, p<0.001, ES=0.97). There was no difference in HVEL of the TK (4.73 ± 1.67 m/s) and that of the JB (4.43 ± 0.78 m/s; p=0.977, ES<0.01).

CONCLUSIONS

The TK is of concern because it is the most common technique and cause of concussion in taekwondo. Future studies should aim to understand rotational accelerations of the head.

Throwing in cricket

This paper considers the kinematic characteristics of overarm throwing with particular emphasis on the techniques of throwing and pitching in baseball. The technique is subdivided into: (1) sequential pattern of throwing, (2) lead foot contact, (3) preparatory phase, (4) arm acceleration and (5) instant of ball release. Specific biomechanical principles that underpin throwing and their application within baseball are identified. The paper also presents a case study of the three-dimensional characteristics of throwing technique in cricket. The aim was to compare the skill in cricket to that previously researched in baseball. The findings for throwing in cricket are similar to those reported for baseball, indicating that there is a definite crossover in the rationale of how an individual should throw specific to the demands of cricket and baseball. The differences noted--greater elbow flexion at lead foot contact and less external rotation during the preparation phase--can be attributed to the demands placed on the fielder and pitcher specific to their respective sports.