Autophagy mediates glucose starvation-induced glioblastoma cell quiescence and chemoresistance through coordinating cell metabolism, cell cycle, and survival.
Wang Lian,Shang Zhouchun,Zhou Yang,Hu Xinyu,Chen Yihong,Fan Yantao,Wei Xiaoyu,Wu Liang,Liang Qiujuan,Zhang Jun,Gao Zhengliang
Cell death & disease
Metabolic reprogramming is pivotal to sustain cancer growth and progression. As such dietary restriction therapy represents a promising approach to starve and treat cancers. Nonetheless, tumors are dynamic and heterogeneous populations of cells with metabolic activities modulated by spatial and temporal contexts. Autophagy is a major pathway controlling cell metabolism. It can downregulate cell metabolism, leading to cancer cell quiescence, survival, and chemoresistance. To understand treatment dynamics and provide rationales for better future therapeutic strategies, we investigated whether and how autophagy is involved in the chemo-cytotoxicity and -resistance using two commonly used human glioblastoma (GBM) cell lines U87 and U251 together with primary cancer cells from the GBM patients. Our results suggest that autophagy mediates chemoresistance through reprogramming cancer cell metabolism and promoting quiescence and survival. Further unbiased transcriptome profiling identified a number of clinically relevant pathways and genes, strongly correlated with TCGA data. Our analyses have not only reported many well-known tumor players, but also uncovered a number of genes that were not previously implicated in cancers and/or GBM. The known functions of these genes are highly suggestive. It would be of high interest to investigate their potential involvement in GBM tumorigenesis, progression, and/or drug resistance. Taken together, our results suggest that autophagy inhibition could be a viable approach to aid GBM chemotherapy and combat drug resistance.
10.1038/s41419-017-0242-x
Low-Dose Mixture Hypothesis of Carcinogenesis Workshop: Scientific Underpinnings and Research Recommendations.
Miller Mark F,Goodson William H,Manjili Masoud H,Kleinstreuer Nicole,Bisson William H,Lowe Leroy
Environmental health perspectives
BACKGROUND:The current single-chemical-as-carcinogen risk assessment paradigm might underestimate or miss the cumulative effects of exposure to chemical mixtures, as highlighted in recent work from the Halifax Project. This is particularly important for chemical exposures in the low-dose range that may be affecting crucial cancer hallmark mechanisms that serve to enable carcinogenesis. OBJECTIVE:Could ongoing low-dose exposures to a mixture of commonly encountered environmental chemicals produce effects in concert that lead to carcinogenesis? A workshop held at the NIEHS in August 2015 evaluated the scientific support for the low-dose mixture hypothesis of carcinogenesis and developed a research agenda. Here we describe the science that supports this novel theory, identify knowledge gaps, recommend future methodologies, and explore preventative risk assessment and policy decision-making that incorporates cancer biology, environmental health science, translational toxicology, and clinical epidemiology. DISCUSSION AND CONCLUSIONS:The theoretical merits of the low-dose carcinogenesis hypothesis are well founded with clear biological relevance, and therefore, the premise warrants further investigation. Expert recommendations include the need for better insights into the ways in which noncarcinogenic constituents might combine to uniquely affect the process of cellular transformation (in vitro) and environmental carcinogenesis (in vivo), including investigations of the role of key defense mechanisms in maintaining transformed cells in a dormant state. The scientific community will need to acknowledge limitations of animal-based models in predicting human responses; evaluate biological events leading to carcinogenesis both spatially and temporally; examine the overlap between measurable cancer hallmarks and characteristics of carcinogens; incorporate epigenetic biomarkers, in silico modelling, high-performance computing and high-resolution imaging, microbiome, metabolomics, and transcriptomics into future research efforts; and build molecular annotations of network perturbations. The restructuring of many existing regulatory frameworks will require adequate testing of relevant environmental mixtures to build a critical mass of evidence on which to base policy decisions. Citation: Miller MF, Goodson WH III, Manjili MH, Kleinstreuer N, Bisson WH, Lowe L. 2017. Low-Dose Mixture Hypothesis of Carcinogenesis Workshop: scientific underpinnings and research recommendations. Environ Health Perspect 125:163-169; http://dx.doi.org/10.1289/EHP411.
10.1289/EHP411
Interfaces of Malignant and Immunologic Clonal Dynamics in Ovarian Cancer.
Zhang Allen W,McPherson Andrew,Milne Katy,Kroeger David R,Hamilton Phineas T,Miranda Alex,Funnell Tyler,Little Nicole,de Souza Camila P E,Laan Sonya,LeDoux Stacey,Cochrane Dawn R,Lim Jamie L P,Yang Winnie,Roth Andrew,Smith Maia A,Ho Julie,Tse Kane,Zeng Thomas,Shlafman Inna,Mayo Michael R,Moore Richard,Failmezger Henrik,Heindl Andreas,Wang Yi Kan,Bashashati Ali,Grewal Diljot S,Brown Scott D,Lai Daniel,Wan Adrian N C,Nielsen Cydney B,Huebner Curtis,Tessier-Cloutier Basile,Anglesio Michael S,Bouchard-Côté Alexandre,Yuan Yinyin,Wasserman Wyeth W,Gilks C Blake,Karnezis Anthony N,Aparicio Samuel,McAlpine Jessica N,Huntsman David G,Holt Robert A,Nelson Brad H,Shah Sohrab P
Cell
High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.
10.1016/j.cell.2018.03.073
Immune mechanisms linking metabolic injury to inflammation and fibrosis in fatty liver disease - novel insights into cellular communication circuits.
Journal of hepatology
Non-alcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver disease and is emerging as the leading cause of cirrhosis, liver transplantation and hepatocellular carcinoma (HCC). NAFLD is a metabolic disease that is considered the hepatic manifestation of the metabolic syndrome; however, during the evolution of NAFLD from steatosis to non-alcoholic steatohepatitis (NASH), to more advanced stages of NASH with liver fibrosis, the immune system plays an integral role. Triggers for inflammation are rooted in hepatic (lipid overload, lipotoxicity, oxidative stress) and extrahepatic (gut-liver axis, adipose tissue, skeletal muscle) systems, resulting in unique immune-mediated pathomechanisms in NAFLD. In recent years, the implementation of single-cell RNA-sequencing and high dimensional multi-omics (proteogenomics, lipidomics) and spatial transcriptomics have tremendously advanced our understanding of the complex heterogeneity of various liver immune cell subsets in health and disease. In NAFLD, several emerging inflammatory mechanisms have been uncovered, including profound macrophage heterogeneity, auto-aggressive T cells, the role of unconventional T cells and platelet-immune cell interactions, potentially yielding novel therapeutics. In this review, we will highlight the recent discoveries related to inflammation in NAFLD, discuss the role of immune cell subsets during the different stages of the disease (including disease regression) and integrate the multiple systems driving inflammation. We propose a refined concept by which the immune system contributes to all stages of NAFLD and discuss open scientific questions arising from this paradigm shift that need to be unravelled in the coming years. Finally, we discuss novel therapeutic approaches to target the multiple triggers of inflammation, including combination therapy via nuclear receptors (FXR agonists, PPAR agonists).
10.1016/j.jhep.2022.06.012
Helicase-like transcription factor (HLTF)-deleted CDX/TME model of colorectal cancer increased transcription of oxidative phosphorylation genes and diverted glycolysis to boost S-glutathionylation in lymphatic intravascular metastatic niches.
PloS one
Helicase-like transcription factor (HLTF) also known as SMARCA3, protects genome integrity. A tumor suppressor, HLTF is expressed in tumor cells but not in the tumor microenvironment (TME) in early-stage colorectal cancer (CRC). With disease progression, there is high concordance between epigenetic silencing of HLTF in CRC cells and negligible HLTF expression in the TME. We developed a cell line-derived xenograft (CDX) model and show for the first time that HLTF-deletion in cancer cells and the TME results in metabolic reprogramming that mitigates oxidative stress in lymphatic intravascular metastatic niches. The two metabolic pathways that derive energy from glucose-glycolysis and oxidative phosphorylation (OXPHOS)-are variously utilized by cancer cells depending upon the TME. HIF-1α, a master regulator of glycolysis, was eliminated from a role in reprogramming metabolism to satisfy CDX energetic requirements by RNAseq and spatial transcriptomics. Variability in the gut microbiome, with a putative role in altered metabolism, was also eliminated. HLTF-deleted cancer cells recovered from DNA damage at a transcriptomic level induction of DNA repair and OXPHOS genes linked to an amoeboid-associated phenotype at the tumor border (confocal microscopy). HLTF-deleted cancer and endothelial cells of lymphatic (PDPN) intravascular niches in the TME shared a site-specific protein S-glutathionylation signature (2D DIGE, MALDI-TOF/TOF mass spectrometry) for three glycolytic enzymes (PGK1 Cys379/380, PGAM1 Cys55, ENOA1 Cys119) that diverted glycolysis in support of continued glutathione biosynthesis. The collective absence of HLTF/Hltf from tumor and TME achieved redox homeostasis throughout the CDX and promoted metastasis.
10.1371/journal.pone.0291023
Advanced progress of spatial metabolomics in head and neck cancer research.
Neoplasia (New York, N.Y.)
Head and neck cancer ranks as the sixth most prevalent malignancy, constituting 5 % of all cancer cases. Its inconspicuous onset often leads to advanced stage diagnoses, prompting the need for early detection to enhance patient prognosis. Currently, research into early diagnostic markers relies predominantly on genomics, proteomics, transcriptomics, and other methods, which, unfortunately, necessitate tumor tissue homogenization, resulting in the loss of temporal and spatial information. Emerging as a recent addition to the omics toolkit, spatial metabolomics stands out. This method conducts in situ mass spectrometry analyses on fresh tissue specimens while effectively preserving their spatiotemporal information. The utilization of spatial metabolomics in life science research offers distinct advantages. This article comprehensively reviews the progress of spatial metabolomics in head and neck cancer research, encompassing insights into cancer cell metabolic reprogramming. Various mass spectrometry imaging techniques, such as secondary ion mass spectrometry, stroma-assisted laser desorption/ionization, and desorption electrospray ionization, enable in situ metabolite analysis for head and neck cancer. Finally, significant emphasis is placed on the application of presently available techniques for early diagnosis, margin assessment, and prognosis of head and neck cancer.
10.1016/j.neo.2023.100958
Multi-dimensional characterization of immunological profiles in small cell lung cancer uncovers clinically relevant immune subtypes with distinct prognoses and therapeutic vulnerabilities.
Pharmacological research
Small-cell lung cancer (SCLC) is generally considered a 'homogenous' disease, with little documented inter-tumor heterogeneity in treatment guidance or prognosis evaluation. The precise identification of clinically relevant molecular subtypes remains incomplete and their translation into clinical practice is limited. In this retrospective cohort study, we comprehensively characterized the immune microenvironment in SCLC by integrating transcriptional and protein profiling of formalin-fixation-and-paraffin-embedded (FFPE) samples from 29 patients. We identified two distinct disease subtypes: immune-enriched (IE-subtype) and immune-deprived (ID-subtype), displaying heterogeneity in immunological, biological, and clinical features. The IE-subtype was characterized by abundant immune infiltrate and elevated levels of interferon-alpha/gamma (IFNα/IFNγ) and inflammatory response, while the ID-subtype featured a complete lack of immune infiltration and a more proliferative phenotype. These two immune subtypes are associated with clinical benefits in SCLC patients treated with adjuvant therapy, with the IE-subtype exhibiting a more favorable response leading to improved survival and reduced disease recurrence risk. Additionally, we identified and validated a personalized prognosticator of immunophenotyping, the CCL5/CXCL9 chemokine index (CCI), using machine learning. The CCI demonstrated superior predictive abilities for prognosis and clinical benefits in SCLC patients, validated in our institute immunohistochemistry cohort and multicenter bulk transcriptomic data cohorts. In conclusion, our study provides a comprehensive and multi-dimensional characterization of the immune architecture of SCLC using clinical FFPE samples and proposes a new immune subtyping conceptual framework enabling risk stratification and the appropriate selection of individualized therapy.
10.1016/j.phrs.2023.106844
CD8+ and FoxP3+ T-Cell Cellular Density and Spatial Distribution After Programmed Death-Ligand 1 Check Point Inhibition.
The Laryngoscope
OBJECTIVES:To analyze CD8+ and FoxP3+ T-cell cellular density (CD) and intercellular distances (ID) in head and neck squamous cell carcinoma (HNSCC) samples from a neoadjuvant trial of durvalumab +/- metformin. METHODS:Paired pre- and post-treatment primary HNSCC tumor samples were stained for CD8+ and FoxP3+. Digital image analysis was used to determine estimated mean CD8+ and FoxP3+ CDs and CD8+-FoxP3+ IDs in the leading tumor edge (LTE) and tumor adjacent stroma (TAS) stratified by treatment arm, human papillomavirus (HPV) status, and pathologic treatment response. A subset of samples was characterized for T-cell related signatures using digital spatial genomic profiling. RESULTS:Post-treatment analysis revealed a significant decrease in FoxP3+ CD and an increase in CD8+ CDs in the TAS between patients receiving durvalumab and metformin versus durvlaumab alone. Both treatment arms demonstrated significant post-treatment increases in ID. Although HPV+ and HPV- had similar immune cell CDs in the tumor microenvironment, HPV+ pre-treatment samples had 1.60 times greater ID compared with HPV- samples, trending toward significance (p = 0.05). At baseline, pathologic responders demonstrated a 1.16-fold greater CD8+ CDs in the LTE (p = 0.045) and 2.28-fold greater ID (p = 0.001) than non-responders. Digital spatial profiling revealed upregulation of FoxP3+ and cytotoxic T-lymphocyte-associated antigen 4 (CTLA-4) in the TAS (p = 0.006, p = 0.026) in samples from pathologic responders. CONCLUSIONS:Analysis of CD8+ and FoxP3+ detected population differences according to HPV status, pathologic response, and treatment. Greater CD8+-FoxP3+ ID was associated with pathologic response. CD8+ and FoxP3+ T-cell distributions may be predictive of response to immune checkpoint inhibition. CLINICALTRIALS:gov (Identifier NCT03618654). LEVEL OF EVIDENCE:3 Laryngoscope, 133:1875-1884, 2023.
10.1002/lary.30389
Neoadjuvant intratumoral influenza vaccine treatment in patients with proficient mismatch repair colorectal cancer leads to increased tumor infiltration of CD8+ T cells and upregulation of PD-L1: a phase 1/2 clinical trial.
Journal for immunotherapy of cancer
BACKGROUND:In colorectal cancer, the effects of immune checkpoint inhibitors are mostly limited to patients with deficient mismatch repair tumors, characterized by a high grade infiltration of CD8+T cells. Interventions aimed at increasing intratumoral CD8+T-cell infiltration in proficient mismatch repair tumors are lacking. METHODS:We conducted a proof of concept phase 1/2 clinical trial, where patients with non-metastasizing sigmoid or rectal cancer, scheduled for curative intended surgery, were treated with an endoscopic intratumorally administered neoadjuvant influenza vaccine. Blood and tumor samples were collected before the injection and at the time of surgery. The primary outcome was safety of the intervention. Evaluation of pathological tumor regression grade, immunohistochemistry, flow cytometry of blood, tissue bulk transcriptional analyses, and spatial protein profiling of tumor regions were all secondary outcomes. RESULTS:A total of 10 patients were included in the trial. Median patient age was 70 years (range 54-78), with 30% women. All patients had proficient mismatch repair Union of International Cancer Control stage I-III tumors. No endoscopic safety events occurred, with all patients undergoing curative surgery as scheduled (median 9 days after intervention). Increased CD8+T-cell tumor infiltration was evident after vaccination (median 73 vs 315 cells/mm, p<0.05), along with significant downregulation of messenger RNA gene expression related to neutrophils and upregulation of transcripts encoding cytotoxic functions. Spatial protein analysis showed significant local upregulation of programmed death-ligand 1 (PD-L1) (adjusted p value<0.05) and downregulation of FOXP3 (adjusted p value<0.05). CONCLUSIONS:Neoadjuvant intratumoral influenza vaccine treatment in this cohort was demonstrated to be safe and feasible, and to induce CD8+T-cell infiltration and upregulation of PD-L1 proficient mismatch repair sigmoid and rectal tumors. Definitive conclusions regarding safety and efficacy can only be made in larger cohorts. TRIAL REGISTRATION NUMBER:NCT04591379.
10.1136/jitc-2023-006774
Integrative Bayesian tensor regression for imaging genetics applications.
Frontiers in neuroscience
Identifying biomarkers for Alzheimer's disease with a goal of early detection is a fundamental problem in clinical research. Both medical imaging and genetics have contributed informative biomarkers in literature. To further improve the performance, recently, there is an increasing interest in developing analytic approaches that combine data across modalities such as imaging and genetics. However, there are limited methods in literature that are able to systematically combine high-dimensional voxel-level imaging and genetic data for accurate prediction of clinical outcomes of interest. Existing prediction models that integrate imaging and genetic features often use region level imaging summaries, and they typically do not consider the spatial configurations of the voxels in the image or incorporate the dependence between genes that may compromise prediction ability. We propose a novel integrative Bayesian scalar-on-image regression model for predicting cognitive outcomes based on high-dimensional spatially distributed voxel-level imaging data, along with correlated transcriptomic features. We account for the spatial dependencies in the imaging voxels via a tensor approach that also enables massive dimension reduction to address the curse of dimensionality, and models the dependencies between the transcriptomic features via a Graph-Laplacian prior. We implement this approach via an efficient Markov chain Monte Carlo (MCMC) computation strategy. We apply the proposed method to the analysis of longitudinal ADNI data for predicting cognitive scores at different visits by integrating voxel-level cortical thickness measurements derived from T1w-MRI scans and transcriptomics data. We illustrate that the proposed imaging transcriptomics approach has significant improvements in prediction compared to prediction using a subset of features from only one modality (imaging or genetics), as well as when using imaging and transcriptomics features but ignoring the inherent dependencies between the features. Our analysis is one of the first to conclusively demonstrate the advantages of prediction based on combining voxel-level cortical thickness measurements along with transcriptomics features, while accounting for inherent structural information.
10.3389/fnins.2023.1212218
DNA methylation atlas of the mouse brain at single-cell resolution.
Nature
Mammalian brain cells show remarkable diversity in gene expression, anatomy and function, yet the regulatory DNA landscape underlying this extensive heterogeneity is poorly understood. Here we carry out a comprehensive assessment of the epigenomes of mouse brain cell types by applying single-nucleus DNA methylation sequencing to profile 103,982 nuclei (including 95,815 neurons and 8,167 non-neuronal cells) from 45 regions of the mouse cortex, hippocampus, striatum, pallidum and olfactory areas. We identified 161 cell clusters with distinct spatial locations and projection targets. We constructed taxonomies of these epigenetic types, annotated with signature genes, regulatory elements and transcription factors. These features indicate the potential regulatory landscape supporting the assignment of putative cell types and reveal repetitive usage of regulators in excitatory and inhibitory cells for determining subtypes. The DNA methylation landscape of excitatory neurons in the cortex and hippocampus varied continuously along spatial gradients. Using this deep dataset, we constructed an artificial neural network model that precisely predicts single neuron cell-type identity and brain area spatial location. Integration of high-resolution DNA methylomes with single-nucleus chromatin accessibility data enabled prediction of high-confidence enhancer-gene interactions for all identified cell types, which were subsequently validated by cell-type-specific chromatin conformation capture experiments. By combining multi-omic datasets (DNA methylation, chromatin contacts, and open chromatin) from single nuclei and annotating the regulatory genome of hundreds of cell types in the mouse brain, our DNA methylation atlas establishes the epigenetic basis for neuronal diversity and spatial organization throughout the mouse cerebrum.
10.1038/s41586-020-03182-8
Single-Cell and Spatial Transcriptome Analysis of Dermal Fibroblast Development in Perinatal Mouse Skin: Dynamic Lineage Differentiation and Key Driver Genes.
The Journal of investigative dermatology
Several single-cell RNA studies of developing mouse skin have elucidated the molecular and cellular processes involved in skin development. However, they have primarily focused on either the fetal or early postnatal period, leaving a gap in our understanding of skin development. In this study, we conducted a comprehensive time-series analysis by combining single-cell RNA-sequencing datasets collected at different stages of development (embryonic days 13.5, 14.5, and 16.5 and postnatal days 0, 2, and 4) and validated our findings through multipanel in situ spatial transcriptomics. Our analysis indicated that embryonic fibroblasts exhibit heterogeneity from a very early stage and that the rapid determination of each lineage occurs within days after birth. The expression of putative key driver genes, including Hey1, Ebf1, Runx3, and Sox11 for the dermal papilla trajectory; Lrrc15 for the dermal sheath trajectory; Zfp536 and Nrn1 for the papillary fibroblast trajectory; and Lrrn4cl and Mfap5 for the fascia fibroblast trajectory, was detected in the corresponding, spatially identified cell types. Finally, cell-to-cell interaction analysis indicated that the dermal papilla lineage is the primary source of the noncanonical Wnt pathway during skin development. Together, our study provides a transcriptomic reference for future research in the field of skin development and regeneration.
10.1016/j.jid.2023.11.008
First-in-human Intravesical Delivery of Pembrolizumab Identifies Immune Activation in Bladder Cancer Unresponsive to Bacillus Calmette-Guérin.
European urology
BACKGROUND:Intravenous immune checkpoint inhibition is an effective anticancer strategy for bacillus Calmette-Guérin (BCG)-unresponsive non-muscle-invasive bladder cancer (NMIBC) but may be associated with greater systemic toxicity compared with localized therapies. OBJECTIVE:We assessed the safety and antitumor activity of intravesical pembrolizumab combined with BCG. DESIGN, SETTING, AND PARTICIPANTS:A 3 + 3 phase 1 trial of pembrolizumab + BCG was conducted in patients with BCG-unresponsive NMIBC (NCT02808143). INTERVENTION:Pembrolizumab was given intravesically (1-5 mg/kg for 2 h) beginning 2 weeks prior to BCG induction until recurrence. Urine profiling during treatment and spatial transcriptomic profiling of pre- and post-treatment tumors were conducted to identify biomarkers that correlated with response. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS:Safety and tolerability of immune checkpoint inhibition were assessed, and Kaplan-Meier survival analysis was performed. RESULTS AND LIMITATIONS:Nine patients completed therapy. Median follow-up was 35 months for five patients still alive at the end of the trial. The trial was closed due to the COVID-19 pandemic. Grade 1-2 urinary symptoms were common. The maximum tolerated dose was not reached; however, one dose-limiting toxicity was reported (grade 2 diarrhea) in the only patient who reached 52 weeks without recurrence. One death occurred from myasthenia gravis that was deemed potentially related to treatment. The 6-mo and 1-yr recurrence-free rates were 67% (95% confidence interval [CI]: 42-100%) and 22% (95% CI: 6.5-75%), respectively. Pembrolizumab was detected in the urine and not in blood. CD4 T cells were significantly increased in the urine after treatment, and a transcriptomic analysis identified decreased expression of T-cell exhaustion markers in late recurrences. CONCLUSIONS:We demonstrate that intravesical pembrolizumab is safe, feasible, and capable of eliciting strong immune responses in a clinical setting and should be investigated further. PATIENT SUMMARY:Direct application of pembrolizumab to the bladder is a promising alternative for non-muscle-invasive bladder cancer unresponsive to Bacillus Calmette-Guérin and should be investigated further.
10.1016/j.eururo.2022.08.004
Enhanced antitumour immunity following neoadjuvant chemoradiotherapy mediates a favourable prognosis in women with resected pancreatic cancer.
Gut
BACKGROUND:This study investigates sex disparities in clinical outcomes and tumour immune profiles in patients with pancreatic ductal adenocarcinoma (PDAC) who underwent upfront resection or resection preceded by gemcitabine-based neoadjuvant chemoradiotherapy (nCRT). METHODS:Patients originated from the PREOPANC randomised controlled trial. Upfront surgery was performed in 82 patients, and 66 received nCRT before resection. The impact of sex on overall survival (OS) was investigated using Cox proportional hazards models. The immunological landscape within the tumour microenvironment (TME) was mapped using transcriptomic and spatial proteomic profiling. RESULTS:The 5-year OS rate differed between the sexes following resection preceded by nCRT, with 43% for women compared with 22% for men. In multivariate analysis, the female sex was a favourable independent prognostic factor for OS only in the nCRT group (HR 0.19; 95% CI 0.07 to 0.52). Multivariate heterogeneous treatment effects analysis revealed a significant interaction between sex and treatment, implying increased nCRT efficacy among women with resected PDAC. The TME of women contained fewer protumoural CD163+MRC1+M2 macrophages than that of men after nCRT, as indicated by transcriptomic and validated using spatial proteomic profiling. CONCLUSION:PDAC tumours of women are more sensitive to gemcitabine-based nCRT, resulting in longer OS after resection compared with men. This may be due to enhanced immunity impeding the infiltration of protumoral M2 macrophages into the TME. Our findings highlight the importance of considering sex disparities and mitigating immunosuppressive macrophage polarisation for personalised PDAC treatment.
10.1136/gutjnl-2023-330480
Stress Keratin 17 Is a Predictive Biomarker Inversely Associated with Response to Immune Check-Point Blockade in Head and Neck Squamous Cell Carcinomas and Beyond.
Cancers
Low response rates in immune check-point blockade (ICB)-treated head and neck squamous cell carcinoma (HNSCC) drive a critical need for robust, clinically validated predictive biomarkers. Our group previously showed that stress keratin 17 (CK17) suppresses macrophage-mediated CXCL9/CXCL10 chemokine signaling involved in attracting activated CD8+ T cells into tumors, correlating with decreased response rate to pembrolizumab-based therapy in a pilot cohort of ICB-treated HNSCC ( = 26). Here, we performed an expanded analysis of the predictive value of CK17 in ICB-treated HNSCC according to the REMARK criteria and investigated the gene expression profiles associated with high CK17 expression. Pretreatment samples from pembrolizumab-treated HNSCC patients were stained via immunohistochemistry using a CK17 monoclonal antibody ( = 48) and subjected to spatial transcriptomic profiling ( = 8). Our findings were validated in an independent retrospective cohort ( = 22). CK17 RNA expression in pembrolizumab-treated patients with various cancer types was investigated for predictive significance. Of the 48 patients (60% male, median age of 61.5 years), 21 (44%) were CK17 high, and 27 (56%) were CK17 low. A total of 17 patients (35%, 77% CK17 low) had disease control, while 31 patients (65%, 45% CK17 low) had progressive disease. High CK17 expression was associated with a lack of disease control ( = 0.037), shorter time to treatment failure ( = 0.025), and progression-free survival (PFS, = 0.004), but not overall survival (OS, = 0.06). A high CK17 expression was associated with lack of disease control in an independent validation cohort ( = 0.011). PD-L1 expression did not correlate with CK17 expression or clinical outcome. CK17 RNA expression was predictive of PFS and OS in 552 pembrolizumab-treated cancer patients. Our findings indicate that high CK17 expression may predict resistance to ICB in HNSCC patients and beyond.
10.3390/cancers15194905
Site-Specific Cytosol Sampling from a Single Cell in an Intact Tumor Spheroid Using an Electrochemical Syringe.
Nashimoto Yuji,Echigo Masakuni,Ino Kosuke,Shiku Hitoshi
Analytical chemistry
A multicellular tumor aggregate, known as a spheroid, is an indispensable tool to study cancer biology. Owing to its three-dimensional organization, a spheroid exhibits an inherent gradient of nutrients, oxygen, and metabolites within itself. The spheroid provides culture conditions that resemble the microenvironment of certain cancer cells and causes these cells to acquire characteristics relevant to tumors in our body. However, site-specific gene expression analysis in an intact spheroid with single-cell resolution has not been explored. Recently, some types of electrochemical syringes were developed to extract cellular materials from living single cells for transcriptomic analysis. Here, we investigated whether an electrochemical syringe could be used to evaluate site-specific gene expression in a spheroid. A small amount of cytosol (roughly 540-1480 fL, less than the volume of a single cell) was successfully collected from the first, second, and third layers of the spheroid using an electrochemical syringe without causing damage to the spheroid architecture. We found that the and expression levels were different between the surface and the average of the entire spheroid, indicating that there are heterogeneous cellular functions across different regions of the spheroid. This method provides opportunities to improve our understanding of spatial gene expression of single cells in a three-dimensional environment.
10.1021/acs.analchem.9b02062
Decidualization induces a secretome switch in perivascular niche cells of the human endometrium.
Murakami Keisuke,Lee Yie Hou,Lucas Emma S,Chan Yi-Wah,Durairaj Ruban Peter,Takeda Satoru,Moore Jonathan D,Tan Bee K,Quenby Siobhan,Chan Jerry K Y,Gargett Caroline E,Brosens Jan J
Endocrinology
The endometrial perivascular microenvironment is rich in mesenchymal stem-like cells that express type 1 integral membrane protein Sushi domain containing 2 (SUSD2) but the role of these cells in the decidual transformation of this tissue in pregnancy is unknown. We used an antibody directed against SUSD2 (W5C5) to isolate perivascular (W5C5(+)) and nonperivascular (W5C5(-)) fibroblasts from mid-luteal biopsies. We show that SUSD2 expression, and hence the ratio of W5C5(+):W5C5(-) cells, changes in culture depending on cell-cell contact and activation of the Notch signaling pathway. RNA sequencing revealed that cultures derived from W5C5(+) progenitor cells remain phenotypically distinct by the enrichment of novel and established endometrial perivascular signature genes. In an undifferentiated state, W5C5(+)-derived cells produced lower levels of various chemokines and inflammatory modulators when compared with their W5C5(-) counterparts. This divergence in secretomes was switched and became more pronounced upon decidualization, which transformed perivascular W5C5(+) cells into the dominant source of a range of chemokines and cytokines, including leukemia inhibitory factor and chemokine (C-C motif) ligand 7. Our findings suggest that the decidual response is spatially organized at the embryo-maternal interface with differentiating perivascular cells establishing distinct cytokine and chemokine profiles that could potentially direct trophoblast toward maternal vessels and govern local immune responses in pregnancy.
10.1210/en.2014-1370
Tissue niches and immunopathology through the lens of spatial tissue profiling techniques.
European journal of immunology
Spatial organization plays a fundamental role in biology, influencing the function of biological structures at various levels. The immune system, in particular, relies on the orchestrated interactions of immune cells with their microenvironment to mount protective or pathogenic immune responses. The COVID-19 pandemic has underscored the significance of studying immunity within target organs to understand disease progression and severity. To achieve this, multiplex histology and spatial transcriptomics have proven indispensable in providing a spatial context to protein and gene expression patterns. By combining these techniques, researchers gain a more comprehensive understanding of the complex interactions at the cellular and molecular level in distinct tissue niches, key functional units modulating health and disease. In this review, we discuss recent advances in spatial tissue profiling techniques, highlighting their advantages over traditional histopathology studies. The insights gained from these approaches have the potential to revolutionize the diagnosis and treatment of various diseases including cancer, autoimmune disorders, and infectious diseases. However, we also acknowledge their challenges and limitations. Despite these, spatial tissue profiling offers promising opportunities to improve our understanding of how tissue niches direct regional immunity, and their relevance in tissue immunopathology, as a basis for novel therapeutic strategies and personalized medicine.
10.1002/eji.202350484
Single-cell profiling of human dura and meningioma reveals cellular meningeal landscape and insights into meningioma immune response.
Genome medicine
BACKGROUND:Recent investigations of the meninges have highlighted the importance of the dura layer in central nervous system immune surveillance beyond a purely structural role. However, our understanding of the meninges largely stems from the use of pre-clinical models rather than human samples. METHODS:Single-cell RNA sequencing of seven non-tumor-associated human dura samples and six primary meningioma tumor samples (4 matched and 2 non-matched) was performed. Cell type identities, gene expression profiles, and T cell receptor expression were analyzed. Copy number variant (CNV) analysis was performed to identify putative tumor cells and analyze intratumoral CNV heterogeneity. Immunohistochemistry and imaging mass cytometry was performed on selected samples to validate protein expression and reveal spatial localization of select protein markers. RESULTS:In this study, we use single-cell RNA sequencing to perform the first characterization of both non-tumor-associated human dura and primary meningioma samples. First, we reveal a complex immune microenvironment in human dura that is transcriptionally distinct from that of meningioma. In addition, we characterize a functionally diverse and heterogenous landscape of non-immune cells including endothelial cells and fibroblasts. Through imaging mass cytometry, we highlight the spatial relationship among immune cell types and vasculature in non-tumor-associated dura. Utilizing T cell receptor sequencing, we show significant TCR overlap between matched dura and meningioma samples. Finally, we report copy number variant heterogeneity within our meningioma samples. CONCLUSIONS:Our comprehensive investigation of both the immune and non-immune cellular landscapes of human dura and meningioma at single-cell resolution builds upon previously published data in murine models and provides new insight into previously uncharacterized roles of human dura.
10.1186/s13073-022-01051-9
gain-of-function mutation modulates the immunosuppressive microenvironment in non-HPV-associated oral squamous cell carcinoma.
Journal for immunotherapy of cancer
BACKGROUND:, the most mutated gene in solid cancers, has a profound impact on most hallmarks of cancer. Somatic mutations occur in high frequencies in head and neck cancers, including oral squamous cell carcinoma (OSCC). Our study aims to understand the role of gain-of-function mutation in modulating the tumor immune microenvironment (TIME) in OSCC. METHODS:Short hairpin RNA knockdown of mutant in syngeneic oral tumors demonstrated changes in tumor growth between immunocompetent and immunodeficient mice. HTG EdgeSeq targeted messenger RNA sequencing was used to analyze cytokine and immune cell markers in tumors with inactivated mutant . Flow cytometry and multiplex immunofluorescence (mIF) confirmed the role of mutant in the TIME. The gene expression of patients with OSCC was analyzed by CIBERSORT and mIF was used to validate the immune landscape at the protein level. RESULTS:Mutant contributes to a cytokine transcriptome network that inhibits the infiltration of cytotoxic CD8 T cells and promotes intratumoral recruitment of regulatory T cells and M2 macrophages. Moreover, also regulates the spatial distribution of immunocyte populations, and their distribution between central and peripheral intratumoral locations. Interestingly, -mutated tumors are infiltrated with CD8 and CD4 T cells expressing programmed cell death protein 1, and these tumors responded to immune checkpoint inhibitor and stimulator of interferon gene 1 agonist therapy. CIBERSORT analysis of human OSCC samples revealed associations between immune cell populations and the mutation, which paralleled the findings from our syngeneic mouse tumor model. CONCLUSIONS:These findings demonstrate that syngeneic tumors bearing the gain-of-function mutation modulate the TIME to evade tumor immunity, leading to tumor progression and decreased survival.
10.1136/jitc-2023-006666
Nerve growth factor scales endocannabinoid signaling by regulating monoacylglycerol lipase turnover in developing cholinergic neurons.
Keimpema Erik,Tortoriello Giuseppe,Alpár Alán,Capsoni Simona,Arisi Ivan,Calvigioni Daniela,Hu Sherry Shu-Jung,Cattaneo Antonino,Doherty Patrick,Mackie Kenneth,Harkany Tibor
Proceedings of the National Academy of Sciences of the United States of America
Endocannabinoid, particularly 2-arachidonoyl glycerol (2-AG), signaling has recently emerged as a molecular determinant of neuronal migration and synapse formation during cortical development. However, the cell type specificity and molecular regulation of spatially and temporally confined morphogenic 2-AG signals remain unexplored. Here, we demonstrate that genetic and pharmacological manipulation of CB(1) cannabinoid receptors permanently alters cholinergic projection neuron identity and hippocampal innervation. We show that nerve growth factor (NGF), implicated in the morphogenesis and survival of cholinergic projection neurons, dose-dependently and coordinately regulates the molecular machinery for 2-AG signaling via tropomyosine kinase A receptors in vitro. In doing so, NGF limits the sorting of monoacylglycerol lipase (MGL), rate limiting 2-AG bioavailability, to proximal neurites, allowing cell-autonomous 2-AG signaling at CB(1) cannabinoid receptors to persist at atypical locations to induce superfluous neurite extension. We find that NGF controls MGL degradation in vitro and in vivo and identify the E3 ubiquitin ligase activity of breast cancer type 1 susceptibility protein (BRCA1) as a candidate facilitating MGL's elimination from motile neurite segments, including growth cones. BRCA1 inactivation by cisplatin or genetically can rescue and reposition MGL, arresting NGF-induced growth responses. These data indicate that NGF can orchestrate endocannabinoid signaling to promote cholinergic differentiation and implicate BRCA1 in determining neuronal morphology.
10.1073/pnas.1212563110
A map of the spatial distribution and tumour-associated macrophage states in glioblastoma and grade 4 IDH-mutant astrocytoma.
The Journal of pathology
Tumour-associated macrophages (TAMs) abundantly infiltrate high-grade gliomas and orchestrate immune response, but their diversity in isocitrate dehydrogenase (IDH)-differential grade 4 gliomas remains largely unknown. This study aimed to dissect the transcriptional states, spatial distribution, and clinicopathological significance of distinct monocyte-derived TAM (Mo-TAM) and microglia-derived TAM (Mg-TAM) clusters across glioblastoma-IDH-wild type and astrocytoma-IDH-mutant-grade 4 (Astro-IDH-mut-G4). Single-cell RNA sequencing was performed on four cases of human glioblastoma and three cases of Astro-IDH-mut-G4. Cell clustering, single-cell regulatory network inference, and gene set enrichment analysis were performed to characterize the functional states of myeloid clusters. The spatial distribution of TAM subsets was determined in human glioma tissues using multiplex immunostaining. The prognostic value of different TAM-cluster specific gene sets was evaluated in the TCGA glioma cohort. Profiling and unbiased clustering of 24,227 myeloid cells from glioblastoma and Astro-IDH-mut-G4 identified nine myeloid cell clusters including monocytes, six Mo/Mg-TAM subsets, dendritic cells, and proliferative myeloid clusters. Different Mo/Mg-TAM clusters manifest functional and transcriptional diversity controlled by specific regulons. Multiplex immunostaining of subset-specific markers identified spatial enrichment of distinct TAM clusters at peri-vascular/necrotic areas in tumour parenchyma or at the tumour-brain interface. Glioblastoma harboured a substantially higher number of monocytes and Mo-TAM-inflammatory clusters, whereas Astro-IDH-mut-G4 had a higher proportion of TAM subsets mediating antigen presentation. Glioblastomas with a higher proportion of monocytes exhibited a mesenchymal signature, increased angiogenesis, and worse patient outcome. Our findings provide insight into myeloid cell diversity and its clinical relevance in IDH-differential grade 4 gliomas, and may serve as a resource for immunotherapy development. © 2022 The Pathological Society of Great Britain and Ireland.
10.1002/path.5984
The co-evolution of the genome and epigenome in colorectal cancer.
Nature
Colorectal malignancies are a leading cause of cancer-related deathand have undergone extensive genomic study. However, DNA mutations alone do not fully explain malignant transformation. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.
10.1038/s41586-022-05202-1
Opposing immune and genetic mechanisms shape oncogenic programs in synovial sarcoma.
Nature medicine
Synovial sarcoma (SyS) is an aggressive neoplasm driven by the SS18-SSX fusion, and is characterized by low T cell infiltration. Here, we studied the cancer-immune interplay in SyS using an integrative approach that combines single-cell RNA sequencing (scRNA-seq), spatial profiling and genetic and pharmacological perturbations. scRNA-seq of 16,872 cells from 12 human SyS tumors uncovered a malignant subpopulation that marks immune-deprived niches in situ and is predictive of poor clinical outcomes in two independent cohorts. Functional analyses revealed that this malignant cell state is controlled by the SS18-SSX fusion, is repressed by cytokines secreted by macrophages and T cells, and can be synergistically targeted with a combination of HDAC and CDK4/CDK6 inhibitors. This drug combination enhanced malignant-cell immunogenicity in SyS models, leading to induced T cell reactivity and T cell-mediated killing. Our study provides a blueprint for investigating heterogeneity in fusion-driven malignancies and demonstrates an interplay between immune evasion and oncogenic processes that can be co-targeted in SyS and potentially in other malignancies.
10.1038/s41591-020-01212-6
Bexmarilimab-induced macrophage activation leads to treatment benefit in solid tumors: The phase I/II first-in-human MATINS trial.
Cell reports. Medicine
Macrophage Clever-1 contributes to impaired antigen presentation and suppression of anti-tumor immunity. This first-in-human trial investigates the safety and tolerability of Clever-1 blockade with bexmarilimab in patients with treatment-refractory solid tumors and assesses preliminary anti-tumor efficacy, pharmacodynamics, and immunologic correlates. Bexmarilimab shows no dose-limiting toxicities in part I (n = 30) and no additional safety signals in part II (n = 108). Disease control (DC) rates of 25%-40% are observed in cutaneous melanoma, gastric, hepatocellular, estrogen receptor-positive breast, and biliary tract cancers. DC associates with improved survival in a landmark analysis and correlates with high pre-treatment intratumoral Clever-1 positivity and increasing on-treatment serum interferon γ (IFNγ) levels. Spatial transcriptomics profiling of DC and non-DC tumors demonstrates bexmarilimab-induced macrophage activation and stimulation of IFNγ and T cell receptor signaling selectively in DC patients. These data suggest that bexmarilimab therapy is well tolerated and show that macrophage targeting can promote immune activation and tumor control in late-stage cancer.
10.1016/j.xcrm.2023.101307
Immune Landscape of Viral- and Carcinogen-Driven Head and Neck Cancer.
Cillo Anthony R,Kürten Cornelius H L,Tabib Tracy,Qi Zengbiao,Onkar Sayali,Wang Ting,Liu Angen,Duvvuri Umamaheswar,Kim Seungwon,Soose Ryan J,Oesterreich Steffi,Chen Wei,Lafyatis Robert,Bruno Tullia C,Ferris Robert L,Vignali Dario A A
Immunity
Head and neck squamous cell carcinoma (HNSCC) arises through exposure to environmental carcinogens or malignant transformation by human papillomavirus (HPV). Here, we assessed the transcriptional profiles of 131,224 single cells from peripheral and intra-tumoral immune populations from patients with HPV and HPV HNSCC and healthy donors. Immune cells within tumors of HPV and HPV HNSCC displayed a spectrum of transcriptional signatures, with helper CD4 T cells and B cells being relatively divergent and CD8+ T cells and CD4+ regulatory T cells being relatively similar. Transcriptional results were contextualized through multispectral immunofluorescence analyses and evaluating putative cell-cell communication based on spatial proximity. These analyses defined a gene expression signature associated with CD4 T follicular helper cells that is associated with longer progression-free survival in HNSCC patients. The datasets and analytical approaches herein provide a resource for the further study of the impact of immune cells on viral- and carcinogen-induced cancers.
10.1016/j.immuni.2019.11.014
The planarian TCF/LEF factor Smed-tcf1 is required for the regeneration of dorsal-lateral neuronal subtypes.
Brown David D R,Molinaro Alyssa M,Pearson Bret J
Developmental biology
The adult brain of the planarian Schmidtea mediterranea (a freshwater flatworm) is a dynamic structure with constant cell turnover as well as the ability to completely regenerate de novo. Despite this, function and pattern is achieved in a reproducible manner from individual to individual in terms of the correct spatial and temporal production of specific neuronal subtypes. Although several signaling molecules have been found to be key to scaling and cell turnover, the mechanisms by which specific neural subtypes are specified remain largely unknown. Here we performed a 6 day RNAseq time course on planarians that were regenerating either 0, 1, or 2 heads in order to identify novel regulators of brain regeneration. Focusing on transcription factors, we identified a TCF/LEF factor, Smed-tcf1, which was required to correctly pattern the dorsal-lateral cell types of the regenerating brain. The most severely affected neurons in Smed-tcf1(RNAi) animals were the dorsal GABAergic neurons, which failed to regenerate, leading to an inability of the animals to phototaxis away from light. Together, Smed-tcf1 is a critical regulator, required to pattern the dorsal-lateral region of the regenerating planarian brain.
10.1016/j.ydbio.2017.08.024
Distinct Transcriptional Programs in Ascitic and Solid Cancer Cells Induce Different Responses to Chemotherapy in High-Grade Serous Ovarian Cancer.
Molecular cancer research : MCR
High-grade serous ovarian cancer (HGSOC) is responsible for the largest number of ovarian cancer deaths. The frequent therapy-resistant relapses necessitate a better understanding of mechanisms driving therapy resistance. Therefore, we mapped more than a hundred thousand cells of HGSOC patients in different phases of the disease, using single-cell RNA sequencing. Within patients, we compared chemonaive with chemotreated samples. As such, we were able to create a single-cell atlas of different HGSOC lesions and their treatment. This revealed a high intrapatient concordance between spatially distinct metastases. In addition, we found remarkable baseline differences in transcriptomics of ascitic and solid cancer cells, resulting in a different response to chemotherapy. Moreover, we discovered different robust subtypes of cancer-associated fibroblasts (CAF) in all patients. Besides inflammatory CAFs, vascular CAFs, and matrix CAFs, we identified a new CAF subtype that was characterized by high expression of STAR, TSPAN8, and ALDH1A1 and clearly enriched after chemotherapy. Together, tumor heterogeneity in both cancer and stromal cells contributes to therapy resistance in HGSOC and could form the basis of novel therapeutic strategies that differentiate between ascitic and solid disease. IMPLICATIONS:The newly characterized differences between ascitic and solid cancer cells before and after chemotherapy could inform novel treatment strategies for metastatic HGSOC.
10.1158/1541-7786.MCR-21-0565
Molecular and immune signatures, and pathological trajectories of fatal COVID-19 lungs defined by in situ spatial single-cell transcriptome analysis.
Journal of medical virology
Despite intensive studies during the last 3 years, the pathology and underlying molecular mechanism of coronavirus disease 2019 (COVID-19) remain poorly defined. In this study, we investigated the spatial single-cell molecular and cellular features of postmortem COVID-19 lung tissues using in situ sequencing (ISS). We detected 10 414 863 transcripts of 221 genes in whole-slide tissues and segmented them into 1 719 459 cells that were mapped to 18 major parenchymal and immune cell types, all of which were infected by SARS-CoV-2. Compared with the non-COVID-19 control, COVID-19 lungs exhibited reduced alveolar cells (ACs) and increased innate and adaptive immune cells. We also identified 19 differentially expressed genes in both infected and uninfected cells across the tissues, which reflected the altered cellular compositions. Spatial analysis of local infection rates revealed regions with high infection rates that were correlated with high cell densities (HIHD). The HIHD regions expressed high levels of SARS-CoV-2 entry-related factors including ACE2, FURIN, TMPRSS2 and NRP1, and co-localized with organizing pneumonia (OP) and lymphocytic and immune infiltration, which exhibited increased ACs and fibroblasts but decreased vascular endothelial cells and epithelial cells, mirroring the tissue damage and wound healing processes. Sparse nonnegative matrix factorization (SNMF) analysis of niche features identified seven signatures that captured structure and immune niches in COVID-19 tissues. Trajectory inference based on immune niche signatures defined two pathological routes. Trajectory A primarily progressed with increased NK cells and granulocytes, likely reflecting the complication of microbial infections. Trajectory B was marked by increased HIHD and OP, possibly accounting for the increased immune infiltration. The OP regions were marked by high numbers of fibroblasts expressing extremely high levels of COL1A1 and COL1A2. Examination of single-cell RNA-seq data (scRNA-seq) from COVID-19 lung tissues and idiopathic pulmonary fibrosis (IPF) identified similar cell populations consisting mainly of myofibroblasts. Immunofluorescence staining revealed the activation of IL6-STAT3 and TGF-β-SMAD2/3 pathways in these cells, likely mediating the upregulation of COL1A1 and COL1A2 and excessive fibrosis in the lung tissues. Together, this study provides a spatial single-cell atlas of cellular and molecular signatures of fatal COVID-19 lungs, which reveals the complex spatial cellular heterogeneity, organization, and interactions that characterized the COVID-19 lung pathology.
10.1002/jmv.29009
Large-Scale Topological Changes Restrain Malignant Progression in Colorectal Cancer.
Cell
Widespread changes to DNA methylation and chromatin are well documented in cancer, but the fate of higher-order chromosomal structure remains obscure. Here we integrated topological maps for colon tumors and normal colons with epigenetic, transcriptional, and imaging data to characterize alterations to chromatin loops, topologically associated domains, and large-scale compartments. We found that spatial partitioning of the open and closed genome compartments is profoundly compromised in tumors. This reorganization is accompanied by compartment-specific hypomethylation and chromatin changes. Additionally, we identify a compartment at the interface between the canonical A and B compartments that is reorganized in tumors. Remarkably, similar shifts were evident in non-malignant cells that have accumulated excess divisions. Our analyses suggest that these topological changes repress stemness and invasion programs while inducing anti-tumor immunity genes and may therefore restrain malignant progression. Our findings call into question the conventional view that tumor-associated epigenomic alterations are primarily oncogenic.
10.1016/j.cell.2020.07.030
Identification of quiescent and spatially restricted mammary stem cells that are hormone responsive.
Fu Nai Yang,Rios Anne C,Pal Bhupinder,Law Charity W,Jamieson Paul,Liu Ruijie,Vaillant François,Jackling Felicity,Liu Kevin He,Smyth Gordon K,Lindeman Geoffrey J,Ritchie Matthew E,Visvader Jane E
Nature cell biology
Despite accumulating evidence for a mammary differentiation hierarchy, the basal compartment comprising stem cells remains poorly characterized. Through gene expression profiling of Lgr5 basal epithelial cells, we identify a new marker, Tetraspanin8 (Tspan8). Fractionation based on Tspan8 and Lgr5 expression uncovered three distinct mammary stem cell (MaSC) subsets in the adult mammary gland. These exist in a largely quiescent state but differ in their reconstituting ability, spatial localization, and their molecular and epigenetic signatures. Interestingly, the deeply quiescent MaSC subset (Lgr5Tspan8) resides within the proximal region throughout life, and has a transcriptome strikingly similar to that of claudin-low tumours. Lgr5Tspan8 cells appear to originate from the embryonic mammary primordia before switching to a quiescent state postnatally but can be activated by ovarian hormones. Our findings reveal an unexpected degree of complexity within the adult MaSC compartment and identify a dormant subset poised for activation in response to physiological stimuli.
10.1038/ncb3471
PLAC8 Localizes to the Inner Plasma Membrane of Pancreatic Cancer Cells and Regulates Cell Growth and Disease Progression through Critical Cell-Cycle Regulatory Pathways.
Kaistha Brajesh P,Lorenz Holger,Schmidt Harald,Sipos Bence,Pawlak Michael,Gierke Berthold,Kreider Ramona,Lankat-Buttgereit Brigitte,Sauer Melanie,Fiedler Lisa,Krattenmacher Anja,Geisel Bettina,Kraus Johann M,Frese Kristopher K,Kelkenberg Sabine,Giese Nathalia A,Kestler Hans A,Gress Thomas M,Buchholz Malte
Cancer research
Pancreatic ductal adenocarcinoma (PDAC) carries the most dismal prognosis of all solid tumors and is generally strongly resistant to currently available chemo- and/or radiotherapy regimens, including targeted molecular therapies. Therefore, unraveling the molecular mechanisms underlying the aggressive behavior of pancreatic cancer is a necessary prerequisite for the development of novel therapeutic approaches. We previously identified the protein placenta-specific 8 (PLAC8, onzin) in a genome-wide search for target genes associated with pancreatic tumor progression and demonstrated that PLAC8 is strongly ectopically expressed in advanced preneoplastic lesions and invasive human PDAC. However, the molecular function of PLAC8 remained unclear, and accumulating evidence suggested its role is highly dependent on cellular and physiologic context. Here, we demonstrate that in contrast to other cellular systems, PLAC8 protein localizes to the inner face of the plasma membrane in pancreatic cancer cells, where it interacts with specific membranous structures in a temporally and spatially stable manner. Inhibition of PLAC8 expression strongly inhibited pancreatic cancer cell growth by attenuating cell-cycle progression, which was associated with transcriptional and/or posttranslational modification of the central cell-cycle regulators CDKN1A, retinoblastoma protein, and cyclin D1 (CCND1), but did not impact autophagy. Moreover, Plac8 deficiency significantly inhibited tumor formation in genetically engineered mouse models of pancreatic cancer. Together, our findings establish PLAC8 as a central mediator of tumor progression in PDAC and as a promising candidate gene for diagnostic and therapeutic targeting.
10.1158/0008-5472.CAN-15-0216
Whole-genome sequencing reveals the molecular implications of the stepwise progression of lung adenocarcinoma.
Nature communications
The mechanism underlying the development of tumors, particularly at early stages, still remains mostly elusive. Here, we report whole-genome long and short read sequencing analysis of 76 lung cancers, focusing on very early-stage lung adenocarcinomas such as adenocarcinoma in situ (AIS) and minimally invasive adenocarcinoma. The obtained data is further integrated with bulk and spatial transcriptomic data and epigenomic data. These analyses reveal key events in lung carcinogenesis. Minimal somatic mutations in pivotal driver mutations and essential proliferative factors are the only detectable somatic mutations in the very early-stage of AIS. These initial events are followed by copy number changes and global DNA hypomethylation. Particularly, drastic changes are initiated at the later AIS stage, i.e., in Noguchi type B tumors, wherein cancer cells are exposed to the surrounding microenvironment. This study sheds light on the pathogenesis of lung adenocarcinoma from integrated pathological and molecular viewpoints.
10.1038/s41467-023-43732-y
Deep learning-based phenotyping reclassifies combined hepatocellular-cholangiocarcinoma.
Nature communications
Primary liver cancer arises either from hepatocytic or biliary lineage cells, giving rise to hepatocellular carcinoma (HCC) or intrahepatic cholangiocarcinoma (ICCA). Combined hepatocellular- cholangiocarcinomas (cHCC-CCA) exhibit equivocal or mixed features of both, causing diagnostic uncertainty and difficulty in determining proper management. Here, we perform a comprehensive deep learning-based phenotyping of multiple cohorts of patients. We show that deep learning can reproduce the diagnosis of HCC vs. CCA with a high performance. We analyze a series of 405 cHCC-CCA patients and demonstrate that the model can reclassify the tumors as HCC or ICCA, and that the predictions are consistent with clinical outcomes, genetic alterations and in situ spatial gene expression profiling. This type of approach could improve treatment decisions and ultimately clinical outcome for patients with rare and biphenotypic cancers such as cHCC-CCA.
10.1038/s41467-023-43749-3
Longitudinal dynamics of the tumor hypoxia response: From enzyme activity to biological phenotype.
Science advances
Poor oxygenation (hypoxia) is a common spatially heterogeneous feature of human tumors. Biological responses to tumor hypoxia are orchestrated by the decreased activity of oxygen-dependent enzymes. The affinity of these enzymes for oxygen positions them along a continuum of oxygen sensing that defines their roles in launching reactive and adaptive cellular responses. These responses encompass regulation of all steps in the central dogma, with rapid perturbation of the metabolome and proteome followed by more persistent reprogramming of the transcriptome and epigenome. Core hypoxia response genes and pathways are commonly regulated at multiple inflection points, fine-tuning the dependencies on oxygen concentration and hypoxia duration. Ultimately, shifts in the activity of oxygen-sensing enzymes directly or indirectly endow cells with intrinsic hypoxia tolerance and drive processes that are associated with aggressive phenotypes in cancer including angiogenesis, migration, invasion, immune evasion, epithelial mesenchymal transition, and stemness.
10.1126/sciadv.adj6409
Myc Expression Drives Aberrant Lipid Metabolism in Lung Cancer.
Cancer research
MYC-mediated pathogenesis in lung cancer continues to attract interest for new therapeutic strategies. In this study, we describe a transgenic mouse model of KRAS-driven lung adenocarcinoma that affords reversible activation of MYC, used here as a tool for lipidomic profiling of MYC-dependent lung tumors formed in this model. Advanced mass spectrometric imaging and surface analysis techniques were used to characterize the spatial and temporal changes in lipid composition in lung tissue. We found that normal lung tissue was characterized predominantly by saturated phosphatidylcholines and phosphatidylglycerols, which are major lipid components of pulmonary surfactant. In contrast, tumor tissues displayed an increase in phosphatidylinositols and arachidonate-containing phospholipids that can serve as signaling precursors. Deactivating MYC resulted in a rapid and dramatic decrease in arachidonic acid and its eicosanoid metabolites. In tumors with high levels of MYC, we found an increase in cytosolic phospholipase A2 (cPLA2) activity with a preferential release of membrane-bound arachidonic acid, stimulating the lipoxygenase (LOX) and COX pathways also amplified by MYC at the level of gene expression. Deactivating MYC lowered cPLA2 activity along with COX2 and 5-LOX mRNA levels. Notably, inhibiting the COX/5-LOX pathways in vivo reduced tumor burden in a manner associated with reduced cell proliferation. Taken together, our results show how MYC drives the production of specific eicosanoids critical for lung cancer cell survival and proliferation, with possible implications for the use of COX and LOX pathway inhibitors for lung cancer therapy. Cancer Res; 76(16); 4608-18. ©2016 AACR.
10.1158/0008-5472.CAN-15-3403
Spatial Regulation of Mitochondrial Heterogeneity by Stromal Confinement in Micropatterned Tumor Models.
Scientific reports
Heterogeneity of mitochondrial activities in cancer cells exists across different disease stages and even in the same patient, with increased mitochondrial activities associated with invasive cancer phenotypes and circulating tumor cells. Here, we use a micropatterned tumor-stromal assay (μTSA) comprised of MCF-7 breast cancer cells and bone marrow stromal cells (BMSCs) as a model to investigate the role of stromal constraints in altering the mitochondrial activities of cancer cells within the tumor microenvironment (TME). Using microdissection and RNA sequencing, we revealed a differentially regulated pattern of gene expression related to mitochondrial activities and metastatic potential at the tumor-stromal interface. Gene expression was confirmed by immunostaining of mitochondrial mass, and live microscopic imaging of mitochondrial membrane potential (ΔΨ) and optical redox ratio. We demonstrated that physical constraints by the stromal cells play a major role in ΔΨ heterogeneity, which was positively associated with nuclear translocation of the YAP/TAZ transcriptional co-activators. Importantly, inhibiting actin polymerization and Rho-associated protein kinase disrupted the differential ΔΨ pattern. In addition, we showed a positive correlation between ΔΨ level and metastatic burden in vivo in mice injected with MDA-MB-231 breast cancer cells. This study supports a new regulatory role for the TME in mitochondrial heterogeneity and metastatic potential.
10.1038/s41598-019-47593-8
High-Plex Spatial RNA Profiling Reveals Cell Type‒Specific Biomarker Expression during Melanoma Development.
The Journal of investigative dermatology
Early diagnosis of melanoma is critical for improved survival. However, the biomarkers of early melanoma evolution and their origin within the tumor and its microenvironment, including the keratinocytes, are poorly defined. To address this, we used spatial transcript profiling that maintains the morphological tumor context to measure the expression of >1,000 RNAs in situ in patient-derived formalin-fixed, paraffin-embedded tissue sections in primary melanoma and melanocytic nevi. We profiled 134 regions of interest (each 200 μm in diameter) enriched in melanocytes, neighboring keratinocytes, or immune cells. This approach captured distinct expression patterns across cell types and tumor types during melanoma development. Unexpectedly, we discovered that S100A8 is expressed by keratinocytes within the tumor microenvironment during melanoma growth. Immunohistochemistry of 252 tumors showed prominent keratinocyte-derived S100A8 expression in melanoma but not in benign tumors and confirmed the same pattern for S100A8's binding partner S100A9, suggesting that injury to the epidermis may be an early and readily detectable indicator of melanoma development. Together, our results establish a framework for high-plex, spatial, and cell type‒specific resolution of gene expression in archival tissue applicable to the development of biomarkers and characterization of tumor microenvironment interactions in tumor evolution.
10.1016/j.jid.2021.06.041
Human Ductal Carcinoma In Situ: Advances and Future Perspectives.
Cold Spring Harbor perspectives in medicine
Due to widespread adoption of screening mammography, there has been a significant increase in new diagnoses of ductal carcinoma in situ (DCIS). However, DCIS outcomes remain unclear. A large fraction of human DCIS (>50%) may not need the multimodality treatment options currently offered to all DCIS patients. More importantly, while we may be overtreating many, we cannot identify those most at risk of invasion or metastasis following a DCIS diagnosis. This review summarizes the studies that have furthered our understanding of DCIS pathology and mechanisms of invasive progression by using advanced technologies including spatial genomics, transcriptomics, and multiplex proteomics. This review also highlights a need for rethinking DCIS with a more focused view on epithelial states and programs and their cross talk with the microenvironment.
10.1101/cshperspect.a041319
Deciphering Intratumoral Molecular Heterogeneity in Clear Cell Renal Cell Carcinoma with a Radiogenomics Platform.
Udayakumar Durga,Zhang Ze,Xi Yin,Dwivedi Durgesh K,Fulkerson Michael,Haldeman Sydney,McKenzie Tiffani,Yousuf Qurratulain,Joyce Allison,Hajibeigi Asghar,Notgrass Hollis,de Leon Alberto Diaz,Yuan Qing,Lewis Matthew A,Madhuranthakam Ananth J,Sibley Robert C,Elias Roy,Guo Junyu,Christie Alana,McKay Renée M,Cadeddu Jeffrey A,Bagrodia Aditya,Margulis Vitaly,Brugarolas James,Wang Tao,Kapur Payal,Pedrosa Ivan
Clinical cancer research : an official journal of the American Association for Cancer Research
PURPOSE:Intratumoral heterogeneity (ITH) challenges the molecular characterization of clear cell renal cell carcinoma (ccRCC) and is a confounding factor for therapy selection. Most approaches to evaluate ITH are limited by two-dimensional tissue analyses. Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) can noninvasively assess the spatial landscape of entire tumors in their natural milieu. To assess the potential of DCE-MRI, we developed a vertically integrated radiogenomics colocalization approach for multi-region tissue acquisition and analyses. We investigated the potential of spatial imaging features to predict molecular subtypes using histopathologic and transcriptome correlatives. EXPERIMENTAL DESIGN:We report the results of a prospective study of 49 patients with ccRCC who underwent DCE-MRI prior to nephrectomy. Surgical specimens were sectioned to match the MRI acquisition plane. RNA sequencing data from multi-region tumor sampling (80 samples) were correlated with percent enhancement on DCE-MRI in spatially colocalized regions of the tumor. Independently, we evaluated clinical applicability of our findings in 19 patients with metastatic RCC (39 metastases) treated with first-line antiangiogenic drugs or checkpoint inhibitors. RESULTS:DCE-MRI identified tumor features associated with angiogenesis and inflammation, which differed within and across tumors, and likely contribute to the efficacy of antiangiogenic drugs and immunotherapies. Our vertically integrated analyses show that angiogenesis and inflammation frequently coexist and spatially anti-correlate in the same tumor. Furthermore, MRI contrast enhancement identifies phenotypes with better response to antiangiogenic therapy among patients with metastatic RCC. CONCLUSIONS:These findings have important implications for decision models based on biopsy samples and highlight the potential of more comprehensive imaging-based approaches.
10.1158/1078-0432.CCR-21-0706
Single cell transcriptomic landscape of diabetic foot ulcers.
Nature communications
Diabetic foot ulceration (DFU) is a devastating complication of diabetes whose pathogenesis remains incompletely understood. Here, we profile 174,962 single cells from the foot, forearm, and peripheral blood mononuclear cells using single-cell RNA sequencing. Our analysis shows enrichment of a unique population of fibroblasts overexpressing MMP1, MMP3, MMP11, HIF1A, CHI3L1, and TNFAIP6 and increased M1 macrophage polarization in the DFU patients with healing wounds. Further, analysis of spatially separated samples from the same patient and spatial transcriptomics reveal preferential localization of these healing associated fibroblasts toward the wound bed as compared to the wound edge or unwounded skin. Spatial transcriptomics also validates our findings of higher abundance of M1 macrophages in healers and M2 macrophages in non-healers. Our analysis provides deep insights into the wound healing microenvironment, identifying cell types that could be critical in promoting DFU healing, and may inform novel therapeutic approaches for DFU treatment.
10.1038/s41467-021-27801-8
Harnessing Tumor Evolution to Circumvent Resistance.
Pogrebniak Katherine L,Curtis Christina
Trends in genetics : TIG
High-throughput sequencing can be used to measure changes in tumor composition across space and time. Specifically, comparisons of pre- and post-treatment samples can reveal the underlying clonal dynamics and resistance mechanisms. Here, we discuss evidence for distinct modes of tumor evolution and their implications for therapeutic strategies. In addition, we consider the utility of spatial tissue sampling schemes, single-cell analysis, and circulating tumor DNA to track tumor evolution and the emergence of resistance, as well as approaches that seek to forestall resistance by targeting tumor evolution. Ultimately, characterization of the (epi)genomic, transcriptomic, and phenotypic changes that occur during tumor progression coupled with computational and mathematical modeling of tumor evolutionary dynamics may inform personalized treatment strategies.
10.1016/j.tig.2018.05.007
Tumor Niches: Perspectives for Targeted Therapies in Glioblastoma.
Antioxidants & redox signaling
Glioblastoma (GBM), the most common and lethal primary brain tumor with a median survival rate of only 15 months and a 5-year survival rate of only 6.8%, remains largely incurable despite the intensive multimodal treatment of surgical resection and radiochemotherapy. Developing effective new therapies is an unmet need for patients with GBM. Targeted therapies, such as antiangiogenesis therapy and immunotherapy, show great promise in treating GBM based upon increasing knowledge about brain tumor biology. Single-cell transcriptomics reveals the plasticity, heterogeneity, and dynamics of tumor cells during GBM development and progression. While antiangiogenesis therapy and immunotherapy have been highly effective in some types of cancer, the disappointing results from clinical trials represent continued challenges in applying these treatments to GBM. Molecular and cellular heterogeneity of GBM is developed temporally and spatially, which profoundly contributes to therapeutic resistance and tumor recurrence. Deciphering mechanisms of tumor heterogeneity and mapping tumor niche trajectories and functions will provide a foundation for the development of more effective therapies for GBM patients. In this review, we discuss five different tumor niches and the intercellular and intracellular communications among these niches, including the perivascular, hypoxic, invasive, immunosuppressive, and glioma-stem cell niches. We also highlight the cellular and molecular biology of these niches and discuss potential strategies to target these tumor niches for GBM therapy. 39, 904-922.
10.1089/ars.2022.0187
Circadian Expression of Migratory Factors Establishes Lineage-Specific Signatures that Guide the Homing of Leukocyte Subsets to Tissues.
He Wenyan,Holtkamp Stephan,Hergenhan Sophia Martina,Kraus Kerstin,de Juan Alba,Weber Jasmin,Bradfield Paul,Grenier Julien Martin Pierre,Pelletier Jeoffrey,Druzd David,Chen Chien-Sin,Ince Louise Madeleine,Bierschenk Susanne,Pick Robert,Sperandio Markus,Aurrand-Lions Michel,Scheiermann Christoph
Immunity
The number of leukocytes present in circulation varies throughout the day, reflecting bone marrow output and emigration from blood into tissues. Using an organism-wide circadian screening approach, we detected oscillations in pro-migratory factors that were distinct for specific vascular beds and individual leukocyte subsets. This rhythmic molecular signature governed time-of-day-dependent homing behavior of leukocyte subsets to specific organs. Ablation of BMAL1, a transcription factor central to circadian clock function, in endothelial cells or leukocyte subsets demonstrated that rhythmic recruitment is dependent on both microenvironmental and cell-autonomous oscillations. These oscillatory patterns defined leukocyte trafficking in both homeostasis and inflammation and determined detectable tumor burden in blood cancer models. Rhythms in the expression of pro-migratory factors and migration capacities were preserved in human primary leukocytes. The definition of spatial and temporal expression profiles of pro-migratory factors guiding leukocyte migration patterns to organs provides a resource for the further study of the impact of circadian rhythms in immunity.
10.1016/j.immuni.2018.10.007
Whole-organ analysis of TGF-β-mediated remodelling of the tumour microenvironment by tissue clearing.
Communications biology
Tissue clearing is one of the most powerful strategies for a comprehensive analysis of disease progression. Here, we established an integrated pipeline that combines tissue clearing, 3D imaging, and machine learning and applied to a mouse tumour model of experimental lung metastasis using human lung adenocarcinoma A549 cells. This pipeline provided the spatial information of the tumour microenvironment. We further explored the role of transforming growth factor-β (TGF-β) in cancer metastasis. TGF-β-stimulated cancer cells enhanced metastatic colonization of unstimulated-cancer cells in vivo when both cells were mixed. RNA-sequencing analysis showed that expression of the genes related to coagulation and inflammation were up-regulated in TGF-β-stimulated cancer cells. Further, whole-organ analysis revealed accumulation of platelets or macrophages with TGF-β-stimulated cancer cells, suggesting that TGF-β might promote remodelling of the tumour microenvironment, enhancing the colonization of cancer cells. Hence, our integrated pipeline for 3D profiling will help the understanding of the tumour microenvironment.
10.1038/s42003-021-01786-y
Single-cell transcriptomic mapping of intestinal epithelium that undergoes 3D morphogenesis and mechanodynamic stimulation in a gut-on-a-chip.
iScience
The human gut-on-a-chip has demonstrated -relevant cellular fidelity and physiological functions of Caco-2 intestinal epithelium compared to its static cultures. However, transcriptomic dynamics that controls the morphogenic and mechanodynamic perturbation of Caco-2 epithelium in a microphysiological culture remain elusive. Single-cell transcriptomic analysis revealed that a gut-on-a-chip culture drives three clusters that illustrate distinct gene expressions and their spatial representation in three-dimensional (3D) epithelial layers. A pseudotemporal trajectory analysis elucidated the evolutionary transition from a homogeneous ancestral genotype in Transwell cultures into heterogeneous transcriptomes in gut-on-a-chip cultures, verified in cell cycle perturbation, cytodifferentiation, and intestinal functions in digestion, transport, and drug metabolism. Furthermore, the inversed transcriptomic signature of oncogenes and tumor-suppressor genes of Caco-2 cells confirmed that gut-on-a-chip cultures induce post-mitotic reprogramming of cancer-associated genes. Our study suggests that a physiological gut-on-a-chip culture induces the transcriptomic perturbation of Caco-2 epithelium to elicit -relevant morphogenesis and restoration of normal physiological functions.
10.1016/j.isci.2022.105521
Investigation of Inter- and Intratumoral Heterogeneity of Glioblastoma Using TOF-SIMS.
Molecular & cellular proteomics : MCP
Glioblastoma (GBM) is one of the most aggressive human cancers with a median survival of less than two years. A distinguishing pathological feature of GBM is a high degree of inter- and intratumoral heterogeneity. Intertumoral heterogeneity of GBM has been extensively investigated on genomic, methylomic, transcriptomic, proteomic and metabolomics levels, however only a few studies describe intratumoral heterogeneity because of the lack of methods allowing to analyze GBM samples with high spatial resolution. Here, we applied TOF-SIMS (Time-of-flight secondary ion mass spectrometry) for the analysis of single cells and clinical samples such as paraffin and frozen tumor sections obtained from 57 patients. We developed a technique that allows us to simultaneously detect the distribution of proteins and metabolites in glioma tissue with 800 nm spatial resolution. Our results demonstrate that according to TOF-SIMS data glioma samples can be subdivided into clinically relevant groups and distinguished from the normal brain tissue. In addition, TOF-SIMS was able to elucidate differences between morphologically distinct regions of GBM within the same tumor. By staining GBM sections with gold-conjugated antibodies against Caveolin-1 we could visualize border between zones of necrotic and cellular tumor and subdivide glioma samples into groups characterized by different survival of the patients. Finally, we demonstrated that GBM contains cells that are characterized by high levels of Caveolin-1 protein and cholesterol. This population may partly represent a glioma stem cells. Collectively, our results show that the technique described here allows to analyze glioma tissues with a spatial resolution beyond reach of most of other omics approaches and the obtained data may be used to predict clinical behavior of the tumor.
10.1074/mcp.RA120.001986
Promiscuous gene expression in medullary thymic epithelial cells mirrors the peripheral self.
Derbinski J,Schulte A,Kyewski B,Klein L
Nature immunology
Expression of peripheral antigens in the thymus has been implicated in T cell tolerance and autoimmunity. Here we identified medullary thymic epithelial cells as being a unique cell type that expresses a diverse range of tissue-specific antigens. We found that this promiscuous gene expression was a cell-autonomous property of medullary epithelial cells and was maintained during the entire period of thymic T cell output. It may facilitate tolerance induction to self-antigens that would otherwise be temporally or spatially secluded from the immune system. However, the array of promiscuously expressed self-antigens appeared random rather than selected and was not confined to secluded self-antigens.
10.1038/ni723
Sequential Binding of MEIS1 and NKX2-5 on the Popdc2 Gene: A Mechanism for Spatiotemporal Regulation of Enhancers during Cardiogenesis.
Dupays Laurent,Shang Catherine,Wilson Robert,Kotecha Surendra,Wood Sophie,Towers Norma,Mohun Timothy
Cell reports
The homeobox transcription factors NKX2-5 and MEIS1 are essential for vertebrate heart development and normal physiology of the adult heart. We show that, during cardiac differentiation, the two transcription factors have partially overlapping expression patterns, with the result that as cardiac progenitors from the anterior heart field differentiate and migrate into the cardiac outflow tract, they sequentially experience high levels of MEIS1 and then increasing levels of NKX2-5. Using the Popdc2 gene as an example, we also show that a significant proportion of target genes for NKX2-5 contain a binding motif recognized by NKX2-5, which overlaps with a binding site for MEIS1. Binding of the two factors to such overlapping sites is mutually exclusive, and this provides a simple regulatory mechanism for spatial and temporal synchronization of a common pool of targets between NKX2-5 and MEIS1.
10.1016/j.celrep.2015.08.065
High expression of STAT3 within the tumour-associated stroma predicts poor outcome in breast cancer patients.
Cancer medicine
INTRODUCTION:Triple-negative breast cancer (TNBC) patients have the poorest clinical outcomes compared to other molecular subtypes of breast cancer. IL6/JAK/STAT3 signalling is upregulated in breast cancer; however, there is limited evidence for its role in TNBC. This study aimed to assess the expression of IL6/JAK/STAT3 in TNBC as a prognostic biomarker. METHODS:Tissue microarrays consisting of breast cancer specimens from a retrospective cohort (n = 850) were stained for IL6R, JAK1, JAK2 and STAT3 via immunohistochemistry. Staining intensity was assessed by weighted histoscore and analysed for association with survival/clinical characteristics. In a subset of patients (n = 14) bulk transcriptional profiling was performed using TempO-Seq. Nanostring GeoMx® digital spatial profiling was utilised to establish the differential spatial gene expression in high STAT3 tumours. RESULTS:In TNBC patients, high expression of stromal STAT3 was associated with reduced cancer-specific survival (HR = 2.202, 95% CI: 1.148-4.224, log rank p = 0.018). TNBC patients with high stromal STAT3 had reduced CD4 T-cell infiltrates within the tumour (p = 0.001) and higher tumour budding (p = 0.003). Gene set enrichment analysis (GSEA) of bulk RNA sequencing showed high stromal STAT3 tumours were characterised by enrichment of IFNγ, upregulation of KRAS signalling and inflammatory signalling Hallmark pathways. GeoMx™ spatial profiling showed high stromal STAT3 samples. Pan cytokeratin (panCK)-negative regions were enriched for CD27 (p < 0.001), CD3 (p < 0.05) and CD8 (p < 0.001). In panCK-positive regions, high stromal STAT3 regions had higher expression of VEGFA (p < 0.05). CONCLUSION:High expression of IL6/JAK/STAT3 proteins was associated with poor prognosis and characterised by distinct underlying biology in TNBC.
10.1002/cam4.6014
Small chromosomal regions position themselves autonomously according to their chromatin class.
van de Werken Harmen J G,Haan Josien C,Feodorova Yana,Bijos Dominika,Weuts An,Theunis Koen,Holwerda Sjoerd J B,Meuleman Wouter,Pagie Ludo,Thanisch Katharina,Kumar Parveen,Leonhardt Heinrich,Marynen Peter,van Steensel Bas,Voet Thierry,de Laat Wouter,Solovei Irina,Joffe Boris
Genome research
The spatial arrangement of chromatin is linked to the regulation of nuclear processes. One striking aspect of nuclear organization is the spatial segregation of heterochromatic and euchromatic domains. The mechanisms of this chromatin segregation are still poorly understood. In this work, we investigated the link between the primary genomic sequence and chromatin domains. We analyzed the spatial intranuclear arrangement of a human artificial chromosome (HAC) in a xenospecific mouse background in comparison to an orthologous region of native mouse chromosome. The two orthologous regions include segments that can be assigned to three major chromatin classes according to their gene abundance and repeat repertoire: (1) gene-rich and SINE-rich euchromatin; (2) gene-poor and LINE/LTR-rich heterochromatin; and (3) gene-depleted and satellite DNA-containing constitutive heterochromatin. We show, using fluorescence in situ hybridization (FISH) and 4C-seq technologies, that chromatin segments ranging from 0.6 to 3 Mb cluster with segments of the same chromatin class. As a consequence, the chromatin segments acquire corresponding positions in the nucleus irrespective of their chromosomal context, thereby strongly suggesting that this is their autonomous property. Interactions with the nuclear lamina, although largely retained in the HAC, reveal less autonomy. Taken together, our results suggest that building of a functional nucleus is largely a self-organizing process based on mutual recognition of chromosome segments belonging to the major chromatin classes.
10.1101/gr.213751.116
Better Efficacy of Synchrotron Spatially Microfractionated Radiation Therapy Than Uniform Radiation Therapy on Glioma.
Bouchet Audrey,Bräuer-Krisch Elke,Prezado Yolanda,El Atifi Michèle,Rogalev Léonid,Le Clec'h Céline,Laissue Jean Albert,Pelletier Laurent,Le Duc Géraldine
International journal of radiation oncology, biology, physics
PURPOSE:Synchrotron microbeam radiation therapy (MRT) is based on the spatial fractionation of the incident, highly focused synchrotron beam into arrays of parallel microbeams, typically a few tens of microns wide and depositing several hundred grays. This irradiation modality was shown to have a high therapeutic impact on tumors, especially in intracranial locations. However, mechanisms responsible for such a property are not fully understood. METHODS AND MATERIALS:Thanks to recent progress in dosimetry, we compared the effect of MRT and synchrotron broad beam (BB) radiation therapy delivered at comparable doses (equivalent to MRT valley dose) on tumor growth control and on classical radiobiological functions by histologic evaluation and/or transcriptomic analysis. RESULTS:MRT significantly improved survival of rats bearing 9L intracranial glioma compared with BB radiation therapy delivered at a comparable dose (P<.001); the efficacy of MRT and BB radiation therapy was similar when the MRT dose was half that of BB. The greater efficacy of MRT was not correlated with a difference in cell proliferation (Mki67 and proliferating cell nuclear antigen) or in transcriptomic stimulation of angiogenesis (vascular endothelial growth factor A or tyrosine kinase with immunoglobulin-like and epidermal growth factor-like domains 2) but was correlated with a higher cell death rate (factor for apoptosis signals) and higher recruitment of macrophages (tyrosine kinase with immunoglobulin-like and epidermal growth factor-like domains 1 and CD68 transcripts) a few days after MRT. CONCLUSIONS:These results show the superiority of MRT over BB radiation therapy when applied at comparable doses, suggesting that spatial fractionation is responsible for a specific and particularly efficient tissue response. The higher induction of cell death and immune cell activation in brain tumors treated by MRT may be involved in such responses.
10.1016/j.ijrobp.2016.03.040
The presence of large focal lesions is a strong independent prognostic factor in multiple myeloma.
Rasche Leo,Angtuaco Edgardo J,Alpe Terri L,Gershner Grant H,McDonald James E,Samant Rohan S,Kumar Manoj,Van Hemert Rudy,Epstein Joshua,Deshpande Shayu,Tytarenko Ruslana,Yaccoby Shmuel,Hillengass Jens,Thanendrarajan Sharmilan,Schinke Carolina,van Rhee Frits,Zangari Maurizio,Walker Brian A,Barlogie Bart,Morgan Gareth J,Davies Faith E,Weinhold Niels
Blood
Spatial intratumor heterogeneity is frequently seen in multiple myeloma (MM) and poses a significant challenge for risk classifiers, which rely on tumor samples from the iliac crest. Because biopsy-based assessment of multiple skeletal sites is difficult, alternative strategies for risk stratification are required. Recently, the size of focal lesions (FLs) was shown to be a surrogate marker for spatial heterogeneity, suggesting that data from medical imaging could be used to improve risk stratification approaches. Here, we investigated the prognostic value of FL size in 404 transplant-eligible, newly diagnosed MM patients. Using diffusion-weighted magnetic resonance imaging with background suppression, we identified the presence of multiple large FLs as a strong prognostic factor. Patients with at least 3 large FLs with a product of the perpendicular diameters >5 cm were associated with poor progression-free survival (PFS) and overall survival (OS; median, 2.3 and 3.6 years, respectively). This pattern, seen in 13.8% of patients, was independent of the Revised International Staging System (RISS), gene expression profiling (GEP)-based risk score, gain(1q), or extramedullary disease (hazard ratio, 2.7 and 2.2 for PFS and OS in multivariate analysis, respectively). The number of FLs lost its negative impact on outcome after adjusting for FL size. In conclusion, the presence of at least 3 large FL is a feature of high risk, which can be used to refine the diagnosis of this type of disease behavior and as an entry criterion for risk-stratified trials.
10.1182/blood-2018-04-842880
AP-2γ regulates oestrogen receptor-mediated long-range chromatin interaction and gene transcription.
Tan Si Kee,Lin Zhen Hua,Chang Cheng Wei,Varang Vipin,Chng Kern Rei,Pan You Fu,Yong Eu Leong,Sung Wing Kin,Sung Win King,Cheung Edwin
The EMBO journal
Oestrogen receptor α (ERα) is key player in the progression of breast cancer. Recently, the cistrome and interactome of ERα were mapped in breast cancer cells, revealing the importance of spatial organization in oestrogen-mediated transcription. However, the underlying mechanism of this process is unclear. Here, we show that ERα binding sites (ERBS) identified from the Chromatin Interaction Analysis-Paired End DiTag of ERα are enriched for AP-2 motifs. We demonstrate the transcription factor, AP-2γ, which has been implicated in breast cancer oncogenesis, binds to ERBS in a ligand-independent manner. Furthermore, perturbation of AP-2γ expression impaired ERα DNA binding, long-range chromatin interactions, and gene transcription. In genome-wide analyses, we show that a large number of AP-2γ and ERα binding events converge together across the genome. The majority of these shared regions are also occupied by the pioneer factor, FoxA1. Molecular studies indicate there is functional interplay between AP-2γ and FoxA1. Finally, we show that most ERBS associated with long-range chromatin interactions are colocalized with AP-2γ and FoxA1. Together, our results suggest AP-2γ is a novel collaborative factor in ERα-mediated transcription.
10.1038/emboj.2011.151
The Spatiotemporal Evolution of Lymph Node Spread in Early Breast Cancer.
Clinical cancer research : an official journal of the American Association for Cancer Research
The most significant prognostic factor in early breast cancer is lymph node involvement. This stage between localized and systemic disease is key to understanding breast cancer progression; however, our knowledge of the evolution of lymph node malignant invasion remains limited, as most currently available data are derived from primary tumors. In 11 patients with treatment-naïve node-positive early breast cancer without clinical evidence of distant metastasis, we investigated lymph node evolution using spatial multiregion sequencing ( = 78 samples) of primary and lymph node deposits and genomic profiling of matched longitudinal circulating tumor DNA (ctDNA). Linear evolution from primary to lymph node was rare (1/11), whereas the majority of cases displayed either early divergence between primary and nodes (4/11) or no detectable divergence (6/11), where both primary and nodal cells belonged to a single recent expansion of a metastatic clone. Divergence of metastatic subclones was driven in part by APOBEC. Longitudinal ctDNA samples from 2 of 7 subjects with evaluable plasma taken perioperatively reflected the two major evolutionary patterns and demonstrate that private mutations can be detected even from early metastatic nodal deposits. Moreover, node removal resulted in disappearance of private lymph node mutations in ctDNA. This study sheds new light on a crucial evolutionary step in the natural history of breast cancer, demonstrating early establishment of axillary lymph node metastasis in a substantial proportion of patients. .
10.1158/1078-0432.CCR-17-3374
Olig2 and Hes regulatory dynamics during motor neuron differentiation revealed by single cell transcriptomics.
PLoS biology
During tissue development, multipotent progenitors differentiate into specific cell types in characteristic spatial and temporal patterns. We addressed the mechanism linking progenitor identity and differentiation rate in the neural tube, where motor neuron (MN) progenitors differentiate more rapidly than other progenitors. Using single cell transcriptomics, we defined the transcriptional changes associated with the transition of neural progenitors into MNs. Reconstruction of gene expression dynamics from these data indicate a pivotal role for the MN determinant Olig2 just prior to MN differentiation. Olig2 represses expression of the Notch signaling pathway effectors Hes1 and Hes5. Olig2 repression of Hes5 appears to be direct, via a conserved regulatory element within the Hes5 locus that restricts expression from MN progenitors. These findings reveal a tight coupling between the regulatory networks that control patterning and neuronal differentiation and demonstrate how Olig2 acts as the developmental pacemaker coordinating the spatial and temporal pattern of MN generation.
10.1371/journal.pbio.2003127
Mammary Tumor-Derived Transplants as Breast Cancer Models to Evaluate Tumor-Immune Interactions and Therapeutic Responses.
Cancer research
In breast cancer, the type and distribution of infiltrating immune cells are associated with clinical outcome. Moreover, cancers with abundant tumor-infiltrating lymphocytes (TIL) are more likely to respond to immunotherapy, whereas those in which CD8 T cells are completely absent (deserts) or excluded are less likely to respond. Detailed understanding of this biology is limited by a lack of preclinical breast cancer models that recapitulate TIL distributions and their associated biology. Here we established mammary tumor-derived transplants (mTDT) from 12 Trp53-null mammary tumors in syngeneic BALB/cJ mice and examined the stability of their growth rate, TIL distribution, and transcriptomic profiles. All mTDTs were estrogen receptor negative. Half of the parental tumors were classified as infiltrated, and the rest were divided between excluded and desert phenotypes. After two orthotopic passages, most (70%) mTDT from infiltrated parents recapitulated this pattern, whereas the desert or excluded parental patterns were maintained in about half of daughter mTDT. Approximately 30% of mTDT gave rise to lung or liver metastases, although metastasis was not associated with a TIL phenotype. Unsupervised transcriptomic analysis clustered mTDT according to their TIL spatial patterns. Infiltrated mTDT transplanted subcutaneously or orthotopically were resistant to anti-PD-L1. Profiling implicated prolonged antigen stimulation and loss of effector function of lymphocytes rather than T-cell exhaustion in the lack of response of infiltrated mTDT to checkpoint blockade. In summary, the molecular diversity and immune complexity of mTDT should facilitate the dissection of mechanisms of breast cancer response to immunotherapies. SIGNIFICANCE: A set of diverse preclinical models of breast cancer is characterized to enable mechanistic dissection of tumor-immune interactions and to improve the efficacy of immunotherapies.
10.1158/0008-5472.CAN-21-0253
The Pediatric Cell Atlas: Defining the Growth Phase of Human Development at Single-Cell Resolution.
Developmental cell
Single-cell gene expression analyses of mammalian tissues have uncovered profound stage-specific molecular regulatory phenomena that have changed the understanding of unique cell types and signaling pathways critical for lineage determination, morphogenesis, and growth. We discuss here the case for a Pediatric Cell Atlas as part of the Human Cell Atlas consortium to provide single-cell profiles and spatial characterization of gene expression across human tissues and organs. Such data will complement adult and developmentally focused HCA projects to provide a rich cytogenomic framework for understanding not only pediatric health and disease but also environmental and genetic impacts across the human lifespan.
10.1016/j.devcel.2019.03.001
Cell-laden hydrogel/titanium microhybrids: Site-specific cell delivery to metallic implants for improved integration.
Koenig Geraldine,Ozcelik Hayriye,Haesler Lisa,Cihova Martina,Ciftci Sait,Dupret-Bories Agnes,Debry Christian,Stelzle Martin,Lavalle Philippe,Vrana Nihal Engin
Acta biomaterialia
Porous titanium implants are widely used in dental, orthopaedic and otorhinolaryngology fields to improve implant integration to host tissue. A possible step further to improve the integration with the host is the incorporation of autologous cells in porous titanium structures via cell-laden hydrogels. Fast gelling hydrogels have advantageous properties for in situ applications such as localisation of specific cells and growth factors at a target area without dispersion. The ability to control the cell types in different regions of an implant is important in applications where the target tissue (i) has structural heterogeneity (multiple cell types with a defined spatial configuration with respect to each other); (ii) has physical property gradients essential for its function (such as in the case of osteochondral tissue transition). Due to their near immediate gelation, such gels can also be used for site-specific modification of porous titanium structures, particularly for implants which would face different tissues at different locations. Herein, we describe a step by step design of a model system: the model cell-laden gel-containing porous titanium implants in the form of titanium microbead/hydrogel (maleimide-dextran or maleimide-PVA based) microhybrids. These systems enable the determination of the effect of titanium presence on gel properties and encapsulated cell behaviour as a miniaturized version of full-scale implants, providing a system compatible with conventional analysis methods. We used a fibroblast/vascular endothelial cell co-cultures as our model system and by utilising single microbeads we have quantified the effect of gel microenvironment (degradability, presence of RGD peptides within gel formulation) on cell behaviour and the effect of the titanium presence on cell behaviour and gel formation. Titanium presence slightly changed gel properties without hindering gel formation or affecting cell viability. Cells showed a preference to move towards the titanium beads and fibroblast proliferation was significantly higher in hybrids compared to gel only controls. The MMP (Matrix Metalloproteinase)-sensitive hydrogels induced sprouting by cells in co-culture configuration which was quantified by fluorescence microscopy, confocal microscopy and qRT-PCR (Quantitative Reverse transcription polymerase chain reaction). When the microhybrid up-scaled to 3D thick structures, cellular localisation in specific areas of the 3D titanium structures was achieved, without decreasing overall cell proliferation compared to titanium only scaffolds. Microhybrids of titanium and hydrogels are useful models for deciding the necessary modifications of metallic implants and they can be used as a modelling system for the study of tissue/titanium implant interactions. STATEMENT OF SIGNIFICANCE:This article demonstrates a method to apply cell-laden hydrogels to porous titanium implants and a model of titanium/hydrogel interaction at micro-level using titanium microbeads. The feasibility of site-specific modification of titanium implants with cell-laden microgels has been demonstrated. Use of titanium microbeads in combination with hydrogels with conventional analysis techniques as described in the article can facilitate the characterisation of surface modification of titanium in a relevant model system.
10.1016/j.actbio.2016.01.023
Altered MicroRNA expression confined to specific epithelial cell subpopulations in breast cancer.
Sempere Lorenzo F,Christensen Mette,Silahtaroglu Asli,Bak Mads,Heath Catherine V,Schwartz Gary,Wells Wendy,Kauppinen Sakari,Cole Charles N
Cancer research
MicroRNAs (miRNAs) are a new class of short noncoding regulatory RNAs (18-25 nucleotides) that are involved in diverse developmental and pathologic processes. Altered miRNA expression has been associated with several types of human cancer. However, most studies did not establish whether miRNA expression changes occurred within cells undergoing malignant transformation. To obtain insight into miRNA deregulation in breast cancer, we implemented an in situ hybridization (ISH) method to reveal the spatial distribution of miRNA expression in archived formalin-fixed, paraffin-embedded specimens representing normal and tumor tissue from >100 patient cases. Here, we report that expression of miR-145 and miR-205 was restricted to the myoepithelial/basal cell compartment of normal mammary ducts and lobules, whereas their accumulation was reduced or completely eliminated in matching tumor specimens. Conversely, expression of other miRNAs was detected at varying levels predominantly within luminal epithelial cells in normal tissue; expression of miR-21 was frequently increased, whereas that of let-7a was decreased in malignant cells. We also analyzed the association of miRNA expression with that of epithelial markers; prognostic indicators such as estrogen receptor, progesterone receptor, and HER2; as well as clinical outcome data. This ISH approach provides a more direct and informative assessment of how altered miRNA expression contributes to breast carcinogenesis compared with miRNA expression profiling in gross tissue biopsies. Most significantly, early manifestation of altered miR-145 expression in atypical hyperplasia and carcinoma in situ lesions suggests that this miRNA may have a potential clinical application as a novel biomarker for early detection.
10.1158/0008-5472.CAN-07-5019
Towards deciphering glioblastoma intra-tumoral heterogeneity: The importance of integrating multidimensional models.
Proteomics
Glioblastoma (GBM) is the most common and severe form of brain cancer among adults. Its aggressiveness is largely attributed to its complex and heterogeneous biology that despite maximal surgery and multimodal chemoradiation treatment, inevitably recurs. Traditional large-scale profiling approaches have contributed substantially to the understanding of patient-to-patient inter-tumoral differences in GBM. However, it is now clear that biological differences within an individual (intra-tumoral heterogeneity) are also a prominent factor in treatment resistance and recurrence of GBM and will likely require integration of data from multiple recently developed omics platforms to fully unravel. Here we dissect the growing geospatial model of GBM, which layers intra-tumoral heterogeneity on a GBM stem cell (GSC) precursor, single cell, and spatial level. We discuss potential unique and inter-dependant aspects of the model including potential discordances between observed genotypes and phenotypes in GBM.
10.1002/pmic.202200401
Commensal-driven immune zonation of the liver promotes host defence.
Nature
The liver connects the intestinal portal vasculature with the general circulation, using a diverse array of immune cells to protect from pathogens that translocate from the gut. In liver lobules, blood flows from portal triads that are situated in periportal lobular regions to the central vein via a polarized sinusoidal network. Despite this asymmetry, resident immune cells in the liver are considered to be broadly dispersed across the lobule. This differs from lymphoid organs, in which immune cells adopt spatially biased positions to promote effective host defence. Here we used quantitative multiplex imaging, genetic perturbations, transcriptomics, infection-based assays and mathematical modelling to reassess the relationship between the localization of immune cells in the liver and host protection. We found that myeloid and lymphoid resident immune cells concentrate around periportal regions. This asymmetric localization was not developmentally controlled, but resulted from sustained MYD88-dependent signalling induced by commensal bacteria in liver sinusoidal endothelial cells, which in turn regulated the composition of the pericellular matrix involved in the formation of chemokine gradients. In vivo experiments and modelling showed that this immune spatial polarization was more efficient than a uniform distribution in protecting against systemic bacterial dissemination. Together, these data reveal that liver sinusoidal endothelial cells sense the microbiome, actively orchestrating the localization of immune cells, to optimize host defence.
10.1038/s41586-020-2977-2
Spatial intratumoral heterogeneity and temporal clonal evolution in esophageal squamous cell carcinoma.
Hao Jia-Jie,Lin De-Chen,Dinh Huy Q,Mayakonda Anand,Jiang Yan-Yi,Chang Chen,Jiang Ye,Lu Chen-Chen,Shi Zhi-Zhou,Xu Xin,Zhang Yu,Cai Yan,Wang Jin-Wu,Zhan Qi-Min,Wei Wen-Qiang,Berman Benjamin P,Wang Ming-Rong,Koeffler H Phillip
Nature genetics
Esophageal squamous cell carcinoma (ESCC) is among the most common malignancies, but little is known about its spatial intratumoral heterogeneity (ITH) and temporal clonal evolutionary processes. To address this, we performed multiregion whole-exome sequencing on 51 tumor regions from 13 ESCC cases and multiregion global methylation profiling for 3 of these 13 cases. We found an average of 35.8% heterogeneous somatic mutations with strong evidence of ITH. Half of the driver mutations located on the branches of tumor phylogenetic trees targeted oncogenes, including PIK3CA, NFE2L2 and MTOR, among others. By contrast, the majority of truncal and clonal driver mutations occurred in tumor-suppressor genes, including TP53, KMT2D and ZNF750, among others. Interestingly, phyloepigenetic trees robustly recapitulated the topological structures of the phylogenetic trees, indicating a possible relationship between genetic and epigenetic alterations. Our integrated investigations of spatial ITH and clonal evolution provide an important molecular foundation for enhanced understanding of tumorigenesis and progression in ESCC.
10.1038/ng.3683
Potential Molecular Markers Related to Lymph Node Metastasis and Stalk Resection Margins in Pedunculated T1 Colorectal Cancers Using Digital Spatial Profiling: A Pilot Study with a Small Case Series.
International journal of molecular sciences
There is a debate regarding the prediction of lymph node metastasis (LNM) in pedunculated T1 colorectal cancer (CRC). In this study with four cases of pedunculated T1 CRCs, we aimed to investigate gene expression variations based on the distance from the Haggitt line (HL) and identify potential molecular risk factors for LNM. By leveraging the Cancer Transcriptome Atlas and digital spatial profiling technology, we meticulously analyzed discrete regions, including the head, HL, proximal stalk region (300-1000 μm from HL), and distal stalk region (1500-2000 μm from HL) to identify spatially sequential molecular changes. Our findings showed significant overall gene expression variations among the head, proximal stalk, and distal stalk regions of pedunculated T1 CRCs compared to the control adenoma. Compared to LNM-negative T1 CRCs, LNM-positive T1 CRC showed that the expression of genes involved in immune-related pathways such as , -, and - were significantly downregulated in the distal stalk region compared to the proximal stalk region. In summary, our results may tentatively suggest considering endoscopic resection of the stalk with a minimum 2000 μm margin from the HL, taking into account the gene expression alterations related to immune-related pathways. However, we acknowledge the limitations of this pilot study, notably the small case series, which may restrict the depth of interpretation. Further validation is imperative to substantiate these findings.
10.3390/ijms25021103
Intratumor DNA methylation heterogeneity in glioblastoma: implications for DNA methylation-based classification.
Wenger Anna,Ferreyra Vega Sandra,Kling Teresia,Bontell Thomas Olsson,Jakola Asgeir Store,Carén Helena
Neuro-oncology
BACKGROUND:A feature of glioblastoma (GBM) is cellular and molecular heterogeneity, both within and between tumors. This variability causes a risk for sampling bias and potential tumor escape from future targeted therapy. Heterogeneous intratumor gene expression in GBM is well documented, but little is known regarding the epigenetic heterogeneity. Variability in DNA methylation within tumors would have implications for diagnostics, as methylation can be used for tumor classification, subtyping, and determination of the clinically used biomarker O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation. We therefore aimed to profile the intratumor DNA methylation heterogeneity in GBM and its effect on diagnostic properties. METHODS:Three to 4 spatially separated biopsies per tumor were collected from 12 GBM patients. We performed genome-wide DNA methylation analysis and investigated intratumor variation. RESULTS:All samples were classified as GBM isocitrate dehydrogenase (IDH) wild type (wt)/mutated by methylation profiling, but the subclass differed within 5 tumors. Some GBM samples exhibited higher DNA methylation differences within tumors than between, and many cytosine-phosphate-guanine (CpG) sites (mean: 17 000) had different methylation levels within the tumors. MGMT methylation status differed in IDH mutated patients (1/1). CONCLUSIONS:We demonstrated that intratumor DNA methylation heterogeneity is a feature of GBM. Although all biopsies were classified as GBM IDH wt/mutated by methylation analysis, the assigned subclass differed in samples from the same patient. The observed heterogeneity within tumors is important to consider for methylation-based biomarkers and future improvements in stratification of GBM patients.
10.1093/neuonc/noz011
The multiple myeloma microenvironment is defined by an inflammatory stromal cell landscape.
de Jong Madelon M E,Kellermayer Zoltán,Papazian Natalie,Tahri Sabrin,Hofste Op Bruinink Davine,Hoogenboezem Remco,Sanders Mathijs A,van de Woestijne Pieter C,Bos P Koen,Khandanpour Cyrus,Vermeulen Jessica,Moreau Philippe,van Duin Mark,Broijl Annemiek,Sonneveld Pieter,Cupedo Tom
Nature immunology
Progression and persistence of malignancies are influenced by the local tumor microenvironment, and future eradication of currently incurable tumors will, in part, hinge on our understanding of malignant cell biology in the context of their nourishing surroundings. Here, we generated paired single-cell transcriptomic datasets of tumor cells and the bone marrow immune and stromal microenvironment in multiple myeloma. These analyses identified myeloma-specific inflammatory mesenchymal stromal cells, which spatially colocalized with tumor cells and immune cells and transcribed genes involved in tumor survival and immune modulation. Inflammatory stromal cell signatures were driven by stimulation with proinflammatory cytokines, and analyses of immune cell subsets suggested interferon-responsive effector T cell and CD8 stem cell memory T cell populations as potential sources of stromal cell-activating cytokines. Tracking stromal inflammation in individuals over time revealed that successful antitumor induction therapy is unable to revert bone marrow inflammation, predicting a role for mesenchymal stromal cells in disease persistence.
10.1038/s41590-021-00931-3
Intratumor heterogeneity and branched evolution revealed by multiregion sequencing.
The New England journal of medicine
BACKGROUND:Intratumor heterogeneity may foster tumor evolution and adaptation and hinder personalized-medicine strategies that depend on results from single tumor-biopsy samples. METHODS:To examine intratumor heterogeneity, we performed exome sequencing, chromosome aberration analysis, and ploidy profiling on multiple spatially separated samples obtained from primary renal carcinomas and associated metastatic sites. We characterized the consequences of intratumor heterogeneity using immunohistochemical analysis, mutation functional analysis, and profiling of messenger RNA expression. RESULTS:Phylogenetic reconstruction revealed branched evolutionary tumor growth, with 63 to 69% of all somatic mutations not detectable across every tumor region. Intratumor heterogeneity was observed for a mutation within an autoinhibitory domain of the mammalian target of rapamycin (mTOR) kinase, correlating with S6 and 4EBP phosphorylation in vivo and constitutive activation of mTOR kinase activity in vitro. Mutational intratumor heterogeneity was seen for multiple tumor-suppressor genes converging on loss of function; SETD2, PTEN, and KDM5C underwent multiple distinct and spatially separated inactivating mutations within a single tumor, suggesting convergent phenotypic evolution. Gene-expression signatures of good and poor prognosis were detected in different regions of the same tumor. Allelic composition and ploidy profiling analysis revealed extensive intratumor heterogeneity, with 26 of 30 tumor samples from four tumors harboring divergent allelic-imbalance profiles and with ploidy heterogeneity in two of four tumors. CONCLUSIONS:Intratumor heterogeneity can lead to underestimation of the tumor genomics landscape portrayed from single tumor-biopsy samples and may present major challenges to personalized-medicine and biomarker development. Intratumor heterogeneity, associated with heterogeneous protein function, may foster tumor adaptation and therapeutic failure through Darwinian selection. (Funded by the Medical Research Council and others.).
10.1056/NEJMoa1113205
Spatial multiomics map of trophoblast development in early pregnancy.
Nature
The relationship between the human placenta-the extraembryonic organ made by the fetus, and the decidua-the mucosal layer of the uterus, is essential to nurture and protect the fetus during pregnancy. Extravillous trophoblast cells (EVTs) derived from placental villi infiltrate the decidua, transforming the maternal arteries into high-conductance vessels. Defects in trophoblast invasion and arterial transformation established during early pregnancy underlie common pregnancy disorders such as pre-eclampsia. Here we have generated a spatially resolved multiomics single-cell atlas of the entire human maternal-fetal interface including the myometrium, which enables us to resolve the full trajectory of trophoblast differentiation. We have used this cellular map to infer the possible transcription factors mediating EVT invasion and show that they are preserved in in vitro models of EVT differentiation from primary trophoblast organoids and trophoblast stem cells. We define the transcriptomes of the final cell states of trophoblast invasion: placental bed giant cells (fused multinucleated EVTs) and endovascular EVTs (which form plugs inside the maternal arteries). We predict the cell-cell communication events contributing to trophoblast invasion and placental bed giant cell formation, and model the dual role of interstitial EVTs and endovascular EVTs in mediating arterial transformation during early pregnancy. Together, our data provide a comprehensive analysis of postimplantation trophoblast differentiation that can be used to inform the design of experimental models of the human placenta in early pregnancy.
10.1038/s41586-023-05869-0
Reconfigurable open microfluidics for studying the spatiotemporal dynamics of paracrine signalling.
Yu Jiaquan,Berthier Erwin,Craig Alexandria,de Groot Theodorus E,Sparks Sidney,Ingram Patrick N,Jarrard David F,Huang Wei,Beebe David J,Theberge Ashleigh B
Nature biomedical engineering
The study of intercellular signalling networks requires organotypic microscale systems that facilitate the culture, conditioning and manipulation of cells. Here, we describe a reconfigurable microfluidic cell-culture system that facilitates the assembly of three-dimensional tissue models by stacking layers that contain preconditioned microenvironments. By using principles of open and suspended microfluidics, the Stacks system is easily assembled or disassembled to provide spatial and temporal manoeuvrability in two-dimensional and three-dimensional assays of multiple cell types, enabling the modelling of sequential paracrine-signalling events, such as tumour-cell-mediated differentiation of macrophages and macrophage-facilitated angiogenesis. We used Stacks to recapitulate the in vivo observation that different prostate cancer tissues polarize macrophages with distinct gene-expression profiles of pro-inflammatory and anti-inflammatory cytokines. Stacks also enabled us to show that these two types of macrophages signal distinctly to endothelial cells, leading to blood vessels with different morphologies. Our proof-of-concept experiments exemplify how Stacks can efficiently model multicellular interactions and highlight the importance of spatiotemporal specificity in intercellular signalling.
10.1038/s41551-019-0421-4
Biomarker Technologies to Support Early Clinical Immuno-oncology Development: Advances and Interpretation.
Cannarile Michael A,Gomes Bruno,Canamero Marta,Reis Bernhard,Byrd Allyson,Charo Jehad,Yadav Mahesh,Karanikas Vaios
Clinical cancer research : an official journal of the American Association for Cancer Research
Today, there is a huge effort to develop cancer immunotherapeutics capable of combating cancer cells as well as the biological environment in which they can grow, adapt, and survive. For such treatments to benefit more patients, there is a great need to dissect the complex interplays between tumor cells and the host's immune system. Monitoring mechanisms of resistance to immunotherapeutics can delineate the evolution of key players capable of driving an efficacious antitumor immune response. In doing so, simultaneous and systematic interrogation of multiple biomarkers beyond single biomarker approaches needs to be undertaken. Zooming into cell-to-cell interactions using technological advancements with unprecedented cellular resolution such as single-cell spatial transcriptomics, advanced tissue histology approaches, and new molecular immune profiling tools promises to provide a unique level of molecular granularity of the tumor environment and may support better decision-making during drug development. This review will focus on how such technological tools are applied in clinical settings, to inform the underlying tumor-immune biology of patients and offer a deeper understanding of cancer immune responsiveness to immuno-oncology treatments.
10.1158/1078-0432.CCR-20-2345
Insights into Immune Escape During Tumor Evolution and Response to Immunotherapy Using a Rat Model of Breast Cancer.
Cancer immunology research
Animal models are critical for the preclinical validation of cancer immunotherapies. Unfortunately, mouse breast cancer models do not faithfully reproduce the molecular subtypes and immune environment of the human disease. In particular, there are no good murine models of estrogen receptor-positive (ER+) breast cancer, the predominant subtype in patients. Here, we show that Nitroso-N-methylurea-induced mammary tumors in outbred Sprague-Dawley rats recapitulate the heterogeneity for mutational profiles, ER expression, and immune evasive mechanisms observed in human breast cancer. We demonstrate the utility of this model for preclinical studies by dissecting mechanisms of response to immunotherapy using combination TGFBR inhibition and PD-L1 blockade. Short-term treatment of early-stage tumors induced durable responses. Gene expression profiling and spatial mapping classified tumors as inflammatory and noninflammatory, and identified IFNγ, T-cell receptor (TCR), and B-cell receptor (BCR) signaling, CD74/MHC II, and epithelium-interacting CD8+ T cells as markers of response, whereas the complement system, M2 macrophage phenotype, and translation in mitochondria were associated with resistance. We found that the expression of CD74 correlated with leukocyte fraction and TCR diversity in human breast cancer. We identified a subset of rat ER+ tumors marked by expression of antigen-processing genes that had an active immune environment and responded to treatment. A gene signature characteristic of these tumors predicted disease-free survival in patients with ER+ Luminal A breast cancer and overall survival in patients with metastatic breast cancer receiving anti-PD-L1 therapy. We demonstrate the usefulness of this preclinical model for immunotherapy and suggest examination to expand immunotherapy to a subset of patients with ER+ disease. See related Spotlight by Roussos Torres, p. 672.
10.1158/2326-6066.CIR-21-0804
Dark Side of Cancer Therapy: Cancer Treatment-Induced Cardiopulmonary Inflammation, Fibrosis, and Immune Modulation.
Boopathi Ettickan,Thangavel Chellappagounder
International journal of molecular sciences
Advancements in cancer therapy increased the cancer free survival rates and reduced the malignant related deaths. Therapeutic options for patients with thoracic cancers include surgical intervention and the application of chemotherapy with ionizing radiation. Despite these advances, cancer therapy-related cardiopulmonary dysfunction (CTRCPD) is one of the most undesirable side effects of cancer therapy and leads to limitations to cancer treatment. Chemoradiation therapy or immunotherapy promote acute and chronic cardiopulmonary damage by inducing reactive oxygen species, DNA damage, inflammation, fibrosis, deregulation of cellular immunity, cardiopulmonary failure, and non-malignant related deaths among cancer-free patients who received cancer therapy. CTRCPD is a complex entity with multiple factors involved in this pathogenesis. Although the mechanisms of cancer therapy-induced toxicities are multifactorial, damage to the cardiac and pulmonary tissue as well as subsequent fibrosis and organ failure seem to be the underlying events. The available biomarkers and treatment options are not sufficient and efficient to detect cancer therapy-induced early asymptomatic cell fate cardiopulmonary toxicity. Therefore, application of cutting-edge multi-omics technology, such us whole-exome sequencing, DNA methylation, whole-genome sequencing, metabolomics, protein mass spectrometry and single cell transcriptomics, and 10 X spatial genomics, are warranted to identify early and late toxicity, inflammation-induced carcinogenesis response biomarkers, and cancer relapse response biomarkers. In this review, we summarize the current state of knowledge on cancer therapy-induced cardiopulmonary complications and our current understanding of the pathological and molecular consequences of cancer therapy-induced cardiopulmonary fibrosis, inflammation, immune suppression, and tumor recurrence, and possible treatment options for cancer therapy-induced cardiopulmonary toxicity.
10.3390/ijms221810126
Wnt3a/beta-catenin signaling controls posterior body development by coordinating mesoderm formation and segmentation.
Dunty William C,Biris Kristin K,Chalamalasetty Ravindra B,Taketo Makoto M,Lewandoski Mark,Yamaguchi Terry P
Development (Cambridge, England)
Somitogenesis is thought to be controlled by a segmentation clock, which consists of molecular oscillators in the Wnt3a, Fgf8 and Notch pathways. Using conditional alleles of Ctnnb1 (beta-catenin), we show that the canonical Wnt3a/beta-catenin pathway is necessary for molecular oscillations in all three signaling pathways but does not function as an integral component of the oscillator. Small, irregular somites persist in abnormally posterior locations in the absence of beta-catenin and cycling clock gene expression. Conversely, Notch pathway genes continue to oscillate in the presence of stabilized beta-catenin but boundary formation is delayed and anteriorized. Together, these results suggest that the Wnt3a/beta-catenin pathway is permissive but not instructive for oscillating clock genes and that it controls the anterior-posterior positioning of boundary formation in the presomitic mesoderm (PSM). The Wnt3a/beta-catenin pathway does so by regulating the activation of the segment boundary determination genes Mesp2 and Ripply2 in the PSM through the activation of the Notch ligand Dll1 and the mesodermal transcription factors T and Tbx6. Spatial restriction of Ripply2 to the anterior PSM is ensured by the Wnt3a/beta-catenin-mediated repression of Ripply2 in posterior PSM. Thus, Wnt3a regulates somitogenesis by activating a network of interacting target genes that promote mesodermal fates, activate the segmentation clock, and position boundary determination genes in the anterior PSM.
10.1242/dev.009266
Thyroid cancer under the scope of emerging technologies.
Tarabichi Maxime,Demetter Pieter,Craciun Ligia,Maenhaut Carine,Detours Vincent
Molecular and cellular endocrinology
The vast majority of thyroid cancers originate from follicular cells. We outline outstanding issues at each step along the path of cancer patient care, from prevention to post-treatment follow-up and highlight how emerging technologies will help address them in the coming years. Three directions will dominate the coming technological landscape. Genomics will reveal tumoral evolutionary history and shed light on how these cancers arise from the normal epithelium and the genomics alteration driving their progression. Transcriptomics will gain cellular and spatial resolution providing a full account of intra-tumor heterogeneity and opening a window on the microenvironment supporting thyroid tumor growth. Artificial intelligence will set morphological analysis on an objective quantitative ground laying the foundations of a systematic thyroid tumor classification system. It will also integrate into unified representations the molecular and morphological perspectives on thyroid cancer.
10.1016/j.mce.2021.111491
Landscape of Noncoding RNA in Prostate Cancer.
Hua Junjie T,Chen Sujun,He Housheng H
Trends in genetics : TIG
The transcriptome of prostate cancer is highly heterogeneous, with noncoding transcripts being essential players. Long noncoding RNAs (lncRNAs) and circular RNAs (circRNAs) are two unique classes of noncoding RNA drawing increasing attention. Biologically, they have intriguing properties with important regulatory functions. Clinically, they present as promising biomarkers and therapeutic targets. Recent advancements in technologies have opened up new directions for noncoding RNA research, which include RNA-protein interaction, RNA secondary structure, and spatial transcriptomics. Furthermore, recent work has also evaluated the clinical applications of these noncoding RNAs in noninvasive liquid biopsy and RNA-based therapies. In this review, we summarize recent findings on lncRNAs and circRNAs in prostate cancer, discuss their clinical utilities, and highlight these exciting areas of research.
10.1016/j.tig.2019.08.004
Integrated spatial genomics reveals global architecture of single nuclei.
Nature
Identifying the relationships between chromosome structures, nuclear bodies, chromatin states and gene expression is an overarching goal of nuclear-organization studies. Because individual cells appear to be highly variable at all these levels, it is essential to map different modalities in the same cells. Here we report the imaging of 3,660 chromosomal loci in single mouse embryonic stem (ES) cells using DNA seqFISH+, along with 17 chromatin marks and subnuclear structures by sequential immunofluorescence and the expression profile of 70 RNAs. Many loci were invariably associated with immunofluorescence marks in single mouse ES cells. These loci form 'fixed points' in the nuclear organizations of single cells and often appear on the surfaces of nuclear bodies and zones defined by combinatorial chromatin marks. Furthermore, highly expressed genes appear to be pre-positioned to active nuclear zones, independent of bursting dynamics in single cells. Our analysis also uncovered several distinct mouse ES cell subpopulations with characteristic combinatorial chromatin states. Using clonal analysis, we show that the global levels of some chromatin marks, such as H3 trimethylation at lysine 27 (H3K27me3) and macroH2A1 (mH2A1), are heritable over at least 3-4 generations, whereas other marks fluctuate on a faster time scale. This seqFISH+-based spatial multimodal approach can be used to explore nuclear organization and cell states in diverse biological systems.
10.1038/s41586-020-03126-2
Genetic, cellular, and connectomic characterization of the brain regions commonly plagued by glioma.
Brain : a journal of neurology
For decades, it has been known that gliomas follow a non-random spatial distribution, appearing more often in some brain regions (e.g. the insula) compared to others (e.g. the occipital lobe). A better understanding of the localization patterns of gliomas could provide clues to the origins of these types of tumours, and consequently inform treatment targets. Following hypotheses derived from prior research into neuropsychiatric disease and cancer, gliomas may be expected to localize to brain regions characterized by functional hubness, stem-like cells, and transcription of genetic drivers of gliomagenesis. We combined neuroimaging data from 335 adult patients with high- and low-grade glioma to form a replicable tumour frequency map. Using this map, we demonstrated that glioma frequency is elevated in association cortex and correlated with multiple graph-theoretical metrics of high functional connectedness. Brain regions populated with putative cells of origin for glioma, neural stem cells and oligodendrocyte precursor cells, exhibited a high glioma frequency. Leveraging a human brain atlas of post-mortem gene expression, we found that gliomas were localized to brain regions enriched with expression of genes associated with chromatin organization and synaptic signalling. A set of glioma proto-oncogenes was enriched among the transcriptomic correlates of glioma distribution. Finally, a regression model incorporating connectomic, cellular, and genetic factors explained 58% of the variance in glioma frequency. These results add to previous literature reporting the vulnerability of hub regions to neurological disease, as well as provide support for cancer stem cell theories of glioma. Our findings illustrate how factors of diverse scale, from genetic to connectomic, can independently influence the anatomic localization of brain dysfunction.
10.1093/brain/awaa277
Spatial and Temporal Changes in PD-L1 Expression in Cancer: The Role of Genetic Drivers, Tumor Microenvironment and Resistance to Therapy.
Shklovskaya Elena,Rizos Helen
International journal of molecular sciences
Immunotherapies blocking immune inhibitory receptors programmed cell death-1 (PD-1) and cytotoxic T-lymphocyte-associated protein-4 (CTLA-4) on T-cells have dramatically improved patient outcomes in a range of advanced cancers. However, the lack of response, and the development of resistance remain major obstacles to long-term improvements in patient outcomes. There is significant interest in the clinical use of biomarkers to improve patient selection, and the expression of PD-1 ligand 1 (PD-L1) is often reported as a potential biomarker of response. However, accumulating evidence suggests that the predictive value of PD-L1 expression in tumor biopsies is relatively low due, in part, to its complex biology. In this review, we discuss the biological consequences of PD-L1 expression by various cell types within the tumor microenvironment, and the complex mechanisms that regulate PD-L1 expression at the genomic, transcriptomic and proteomic levels.
10.3390/ijms21197139
Applications of tissue-specific and cancer-selective gene promoters for cancer diagnosis and therapy.
Advances in cancer research
Current treatment of solid tumors with standard of care chemotherapies, radiation therapy and/or immunotherapies are often limited by severe adverse toxic effects, resulting in a narrow therapeutic index. Cancer gene therapy represents a targeted approach that in principle could significantly reduce undesirable side effects in normal tissues while significantly inhibiting tumor growth and progression. To be effective, this strategy requires a clear understanding of the molecular biology of cancer development and evolution and developing biological vectors that can serve as vehicles to target cancer cells. The advent and fine tuning of omics technologies that permit the collective and spatial recognition of genes (genomics), mRNAs (transcriptomics), proteins (proteomics), metabolites (metabolomics), epiomics (epigenomics, epitranscriptomics, and epiproteomics), and their interactomics in defined complex biological samples provide a roadmap for identifying crucial targets of relevance to the cancer paradigm. Combining these strategies with identified genetic elements that control target gene expression uncovers significant opportunities for developing guided gene-based therapeutics for cancer. The purpose of this review is to overview the current state and potential limitations in developing gene promoter-directed targeted expression of key genes and highlights their potential applications in cancer gene therapy.
10.1016/bs.acr.2023.03.005
Pseudobudding: ruptured glands do not represent true tumor buds.
The Journal of pathology
Tumor budding (TB) is a strong biomarker of poor prognosis in colorectal cancer and other solid cancers. TB is defined as isolated single cancer cells or clusters of up to four cancer cells at the invasive tumor front. In areas with a large inflammatory response at the invasive front, single cells and cell clusters surrounding fragmented glands are observed appearing like TB. Occurrence of these small groups is referred to as pseudobudding (PsB), which arises due to external influences such as inflammation and glandular disruption. Using a combination of orthogonal approaches, we show that there are clear biological differences between TB and PsB. TB is representative of active invasion by presenting features of epithelial-mesenchymal transition and exhibiting increased deposition of extracellular matrix within the surrounding tumor microenvironment (TME), whereas PsB represents a reactive response to heavy inflammation where increased levels of granulocytes within the surrounding TME are observed. Our study provides evidence that areas with a strong inflammatory reaction should be avoided in the routine diagnostic assessment of TB. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.
10.1002/path.6146
Machine learning-assisted analysis of epithelial mesenchymal transition pathway for prognostic stratification and immune infiltration assessment in ovarian cancer.
Frontiers in endocrinology
Background:Ovarian cancer is the most lethal gynaecological malignancy, and serous ovarian cancer (SOC) is one of the more important pathological subtypes. Previous studies have reported a significant association of epithelial tomesenchymal transition (EMT) with invasive metastasis and immune modulation of SOC, however, there is a lack of prognostic and immune infiltration biomarkers reported for SOC based on EMT. Methods:Gene expression data for ovarian cancer and corresponding patient clinical data were collected from the TCGA database and the GEO database, and cell type annotation and spatial expression analysis were performed on single cell sequencing data from the GEO database. To understand the cell type distribution of EMT-related genes in SOC single-cell data and the enrichment relationships of biological pathways and tumour functions. In addition, GO functional annotation analysis and KEGG pathway enrichment analysis were performed on mRNAs predominantly expressed with EMT to predict the biological function of EMT in ovarian cancer. The major differential genes of EMT were screened to construct a prognostic risk prediction model for SOC patients. Data from 173 SOC patient samples obtained from the GSE53963 database were used to validate the prognostic risk prediction model for ovarian cancer. Here we also analysed the direct association between SOC immune infiltration and immune cell modulation and EMT risk score. and calculate drug sensitivity scores in the GDSC database.In addition, we assessed the specific relationship between GAS1 gene and SOC cell lines. Results:Single cell transcriptome analysis in the GEO database annotated the major cell types of SOC samples, including: T cell, Myeloid, Epithelial cell, Fibroblast, Endothelial cell, and Bcell. cellchat revealed several cell type interactions that were shown to be associated with EMT-mediated SOC invasion and metastasis. A prognostic stratification model for SOC was constructed based on EMT-related differential genes, and the Kapan-Meier test showed that this biomarker had significant prognostic stratification value for several independent SOC databases. The EMT risk score has good stratification and identification properties for drug sensitivity in the GDSC database. Conclusions:This study constructed a prognostic stratification biomarker based on EMT-related risk genes for immune infiltration mechanisms and drug sensitivity analysis studies in SOC. This lays the foundation for in-depth clinical studies on the role of EMT in immune regulation and related pathway alterations in SOC. It is also hoped to provide effective potential solutions for early diagnosis and clinical treatment of ovarian cancer.
10.3389/fendo.2023.1196094
Comprehensive Molecular Portraits of Invasive Lobular Breast Cancer.
Ciriello Giovanni,Gatza Michael L,Beck Andrew H,Wilkerson Matthew D,Rhie Suhn K,Pastore Alessandro,Zhang Hailei,McLellan Michael,Yau Christina,Kandoth Cyriac,Bowlby Reanne,Shen Hui,Hayat Sikander,Fieldhouse Robert,Lester Susan C,Tse Gary M K,Factor Rachel E,Collins Laura C,Allison Kimberly H,Chen Yunn-Yi,Jensen Kristin,Johnson Nicole B,Oesterreich Steffi,Mills Gordon B,Cherniack Andrew D,Robertson Gordon,Benz Christopher,Sander Chris,Laird Peter W,Hoadley Katherine A,King Tari A, ,Perou Charles M
Cell
Invasive lobular carcinoma (ILC) is the second most prevalent histologic subtype of invasive breast cancer. Here, we comprehensively profiled 817 breast tumors, including 127 ILC, 490 ductal (IDC), and 88 mixed IDC/ILC. Besides E-cadherin loss, the best known ILC genetic hallmark, we identified mutations targeting PTEN, TBX3, and FOXA1 as ILC enriched features. PTEN loss associated with increased AKT phosphorylation, which was highest in ILC among all breast cancer subtypes. Spatially clustered FOXA1 mutations correlated with increased FOXA1 expression and activity. Conversely, GATA3 mutations and high expression characterized luminal A IDC, suggesting differential modulation of ER activity in ILC and IDC. Proliferation and immune-related signatures determined three ILC transcriptional subtypes associated with survival differences. Mixed IDC/ILC cases were molecularly classified as ILC-like and IDC-like revealing no true hybrid features. This multidimensional molecular atlas sheds new light on the genetic bases of ILC and provides potential clinical options.
10.1016/j.cell.2015.09.033
Differential expression of the actin-binding proteins, alpha-actinin-2 and -3, in different species: implications for the evolution of functional redundancy.
Mills M,Yang N,Weinberger R,Vander Woude D L,Beggs A H,Easteal S,North K
Human molecular genetics
The alpha-actinins are a multigene family of four actin-binding proteins related to dystrophin. The two skeletal muscle isoforms of alpha-actinin (ACTN2 and ACTN3) are major structural components of the Z-line involved in anchoring the actin-containing thin filaments. In humans, ACTN2 is expressed in all muscle fibres, while ACTN3 expression is restricted to a subset of type 2 fibres. We have recently demonstrated that alpha-actinin-3 is absent in approximately 18% of individuals in a range of human populations, and that homozygosity for a premature stop codon (577X) accounts for most cases of true alpha-actinin-3 deficiency. Absence of alpha-actinin-3 is not associated with an obvious disease phenotype, raising the possibility that ACTN3 is functionally redundant in humans, and that alpha-actinin-2 is able to compensate for alpha-actinin-3 deficiency. We now present data concerning the expression of ACTN3 in other species. Genotyping of non-human primates indicates that the 577X null mutation has likely arisen in humans. The mouse genome contains four orthologues which all map to evolutionarily conserved syntenic regions for the four human genes. Murine Actn2 and Actn3 are differentially expressed, spatially and temporally, during embryonic development and, in contrast to humans, alpha-actinin-2 expression does not completely overlap alpha-actinin-3 in postnatal skeletal muscle, suggesting independent function. Furthermore, sequence comparison of human, mouse and chicken alpha-actinin genes demonstrates that ACTN3 has been conserved over a long period of evolutionary time, implying a constraint on evolutionary rate imposed by continued function of the gene. These observations provide a real framework in which to test theoretical models of genetic redundancy as they apply to human populations. In addition we highlight the need for caution in making conclusions about gene function from the phenotypic consequences of loss-of-function mutations in animal knockout models.
10.1093/hmg/10.13.1335
Nodal reactive proliferation of monocytoid B-cells may represent atypical memory B-cells.
Journal of microbiology, immunology, and infection = Wei mian yu gan ran za zhi
BACKGROUND:Reactive lymphadenopathies such as toxoplasmosis and cytomegalovirus lymphadenitis are associated with monocytoid cell proliferation. Monocytoid cells are B-lymphocytes with an undetermined subset. METHODS:Using digital spatial profiling whole transcriptome analyses, this study compared monocytoid and control B-cells. The B-cell subset of monocytoid cells was assigned according to gene expression profiles. RESULTS:This study identified 466 differentially expressed genes between monocytoid and control B-cells. The cellular deconvolution algorithm identified monocytoid cells as memory B-cells instead of as naïve B-cells. A comparison of the upregulated genes revealed that atypical memory B-cells had the largest number of genes overlapping with monocytoid cells compared with other memory B-cell subsets. Atypical memory B-cell markers, namely TBX21 (T-bet), FCRL4 (IRTA1), and ITGAX (CD11c), were all upregulated in monocytoid cells. Similar to atypical memory B-cells, monocytoid cells exhibited (1) upregulated transcription factors (TBX21, TOX), (2) upregulated genes associated with B-cell inhibition (FCRL5, FCRL4) and downregulated genes associated with B-cell activation (PIK3CG, NFKB1A, CD40), (3) downregulated cell cycle-related genes (CDK6, MYC), and (4) downregulated cytokine receptors (IL4R). This study also analyzed the expression of monocytoid cell signature genes in various memory B-cell subsets. Atypical memory B-cells exhibited a gene expression pattern similar to that of monocytoid cells, but other memory B-cell subsets did not. Furthermore, monocytoid cells and marginal zone lymphomas differed in gene expression profiles. CONCLUSION:Spatial transcriptomic analyses indicated that monocytoid cells may be atypical memory B-cells.
10.1016/j.jmii.2023.03.010
Non-thermal plasma activates human keratinocytes by stimulation of antioxidant and phase II pathways.
Schmidt Anke,Dietrich Stephan,Steuer Anna,Weltmann Klaus-Dieter,von Woedtke Thomas,Masur Kai,Wende Kristian
The Journal of biological chemistry
Non-thermal atmospheric pressure plasma provides a novel therapeutic opportunity to control redox-based processes, e.g. wound healing, cancer, and inflammatory diseases. By spatial and time-resolved delivery of reactive oxygen and nitrogen species, it allows stimulation or inhibition of cellular processes in biological systems. Our data show that both gene and protein expression is highly affected by non-thermal plasma. Nuclear factor erythroid-related factor 2 (NRF2) and phase II enzyme pathway components were found to act as key controllers orchestrating the cellular response in keratinocytes. Additionally, glutathione metabolism, which is a marker for NRF2-related signaling events, was affected. Among the most robustly increased genes and proteins, heme oxygenase 1, NADPH-quinone oxidoreductase 1, and growth factors were found. The roles of NRF2 targets, investigated by siRNA silencing, revealed that NRF2 acts as an important switch for sensing oxidative stress events. Moreover, the influence of non-thermal plasma on the NRF2 pathway prepares cells against exogenic noxae and increases their resilience against oxidative species. Via paracrine mechanisms, distant cells benefit from cell-cell communication. The finding that non-thermal plasma triggers hormesis-like processes in keratinocytes facilitates the understanding of plasma-tissue interaction and its clinical application.
10.1074/jbc.M114.603555
Single-Cell Transcriptomics Uncovers Zonation of Function in the Mesenchyme during Liver Fibrosis.
Cell reports
Iterative liver injury results in progressive fibrosis disrupting hepatic architecture, regeneration potential, and liver function. Hepatic stellate cells (HSCs) are a major source of pathological matrix during fibrosis and are thought to be a functionally homogeneous population. Here, we use single-cell RNA sequencing to deconvolve the hepatic mesenchyme in healthy and fibrotic mouse liver, revealing spatial zonation of HSCs across the hepatic lobule. Furthermore, we show that HSCs partition into topographically diametric lobule regions, designated portal vein-associated HSCs (PaHSCs) and central vein-associated HSCs (CaHSCs). Importantly we uncover functional zonation, identifying CaHSCs as the dominant pathogenic collagen-producing cells in a mouse model of centrilobular fibrosis. Finally, we identify LPAR1 as a therapeutic target on collagen-producing CaHSCs, demonstrating that blockade of LPAR1 inhibits liver fibrosis in a rodent NASH model. Taken together, our work illustrates the power of single-cell transcriptomics to resolve the key collagen-producing cells driving liver fibrosis with high precision.
10.1016/j.celrep.2019.10.024
The Identity of Human Tissue-Emigrant CD8 T Cells.
Cell
Lymphocyte migration is essential for adaptive immune surveillance. However, our current understanding of this process is rudimentary, because most human studies have been restricted to immunological analyses of blood and various tissues. To address this knowledge gap, we used an integrated approach to characterize tissue-emigrant lineages in thoracic duct lymph (TDL). The most prevalent immune cells in human and non-human primate efferent lymph were T cells. Cytolytic CD8 T cell subsets with effector-like epigenetic and transcriptional signatures were clonotypically skewed and selectively confined to the intravascular circulation, whereas non-cytolytic CD8 T cell subsets with stem-like epigenetic and transcriptional signatures predominated in tissues and TDL. Moreover, these anatomically distinct gene expression profiles were recapitulated within individual clonotypes, suggesting parallel differentiation programs independent of the expressed antigen receptor. Our collective dataset provides an atlas of the migratory immune system and defines the nature of tissue-emigrant CD8 T cells that recirculate via TDL.
10.1016/j.cell.2020.11.019
Spatial organization-dependent EphA2 transcriptional responses revealed by ligand nanocalipers.
Verheyen Toon,Fang Trixy,Lindenhofer Dominik,Wang Yang,Akopyan Karen,Lindqvist Arne,Högberg Björn,Teixeira Ana I
Nucleic acids research
Ligand binding induces extensive spatial reorganization and clustering of the EphA2 receptor at the cell membrane. It has previously been shown that the nanoscale spatial distribution of ligands modulates EphA2 receptor reorganization, activation and the invasive properties of cancer cells. However, intracellular signaling downstream of EphA2 receptor activation by nanoscale spatially distributed ligands has not been elucidated. Here, we used DNA origami nanostructures to control the positions of ephrin-A5 ligands at the nanoscale and investigated EphA2 activation and transcriptional responses following ligand binding. Using RNA-seq, we determined the transcriptional profiles of human glioblastoma cells treated with DNA nanocalipers presenting a single ephrin-A5 dimer or two dimers spaced 14, 40 or 100 nm apart. These cells displayed divergent transcriptional responses to the differing ephrin-A5 nano-organization. Specifically, ephrin-A5 dimers spaced 40 or 100 nm apart showed the highest levels of differential expressed genes compared to treatment with nanocalipers that do not present ephrin-A5. These findings show that the nanoscale organization of ephrin-A5 modulates transcriptional responses to EphA2 activation.
10.1093/nar/gkaa274
Caenorhabditis elegans chromosome arms are anchored to the nuclear membrane via discontinuous association with LEM-2.
Ikegami Kohta,Egelhofer Thea A,Strome Susan,Lieb Jason D
Genome biology
BACKGROUND:Although Caenorhabditis elegans was the first multicellular organism with a completely sequenced genome, how this genome is arranged within the nucleus is not known. RESULTS:We determined the genomic regions associated with the nuclear transmembrane protein LEM-2 in mixed-stage C. elegans embryos via chromatin immunoprecipitation. Large regions of several megabases on the arms of each autosome were associated with LEM-2. The center of each autosome was mostly free of such interactions, suggesting that they are largely looped out from the nuclear membrane. Only the left end of the X chromosome was associated with the nuclear membrane. At a finer scale, the large membrane-associated domains consisted of smaller subdomains of LEM-2 associations. These subdomains were characterized by high repeat density, low gene density, high levels of H3K27 trimethylation, and silent genes. The subdomains were punctuated by gaps harboring highly active genes. A chromosome arm translocated to a chromosome center retained its association with LEM-2, although there was a slight decrease in association near the fusion point. CONCLUSIONS:Local DNA or chromatin properties are the main determinant of interaction with the nuclear membrane, with position along the chromosome making a minor contribution. Genes in small gaps between LEM-2 associated regions tend to be highly expressed, suggesting that these small gaps are especially amenable to highly efficient transcription. Although our data are derived from an amalgamation of cell types in mixed-stage embryos, the results suggest a model for the spatial arrangement of C. elegans chromosomes within the nucleus.
10.1186/gb-2010-11-12-r120
Multi-omic and spatial dissection of immunotherapy response groups in non-small cell lung cancer.
Immunology
The composition and activation status of the cellular milieu contained within the tumour microenvironment (TME) is becoming increasingly recognized as a driving factor for immunotherapy response. Here, we employed multiplex immunohistochemistry (mIHC), and digital spatial profiling (DSP) to capture the targeted immune proteome and transcriptome of tumour and TME compartments from an immune checkpoint inhibitor (ICI)-treated (n = 41) non-small cell lung cancer (NSCLC) patient cohort. We demonstrate by mIHC that the interaction of CD68 macrophages with PD1 , FoxP3 cells is enriched in ICI refractory tumours (p = 0.012). Patients responsive to ICI therapy expressed higher levels of IL2 receptor alpha (CD25, p = 0.028) within their tumour compartments, which corresponded with increased IL2 mRNA (p = 0.001) within their stroma. In addition, stromal IL2 mRNA levels positively correlated with the expression of pro-apoptotic markers cleaved caspase 9 (p = 2e ) and BAD (p = 5.5e ) and negatively with levels of memory marker, CD45RO (p = 7e ). Immuno-inhibitory markers CTLA-4 (p = 0.021) and IDO-1 (p = 0.023) were suppressed in ICI-responsive patients. Tumour expression of CD44 was depleted in the responsive patients (p = 0.02), while higher stromal expression of one of its ligands, SPP1 (p = 0.008), was observed. Cox survival analysis also indicated tumour CD44 expression was associated with poorer prognosis (hazard ratio [HR] = 1.61, p = 0.01), consistent with its depletion in ICI-responsive patients. Through multi-modal approaches, we have dissected the characteristics of NSCLC immunotherapy treatment groups and provide evidence for the role of several markers including IL2, CD25, CD44 and SPP1 in the efficacy of current generations of ICI therapy.
10.1111/imm.13646
Simple oligonucleotide-based multiplexing of single-cell chromatin accessibility.
Molecular cell
Microdroplet single-cell ATAC-seq is widely used to measure chromatin accessibility, however, highly scalable and simple sample multiplexing procedures are not available. Here, we present a transposome-assisted single nucleus barcoding approach for ATAC-seq (SNuBar-ATAC) that utilizes a single oligonucleotide adaptor for multiplexing samples during the existing tagmentation step and does not require a pre-labeling procedure. The accuracy and scalability of SNuBar-ATAC was evaluated using cell line mixture experiments. We applied SNuBar-ATAC to investigate treatment-induced chromatin accessibility dynamics by multiplexing 28 mice with lung tumors that received different combinations of chemo, radiation, and targeted immunotherapy. We also applied SNuBar-ATAC to study spatial epigenetic heterogeneity by multiplexing 32 regions from a human breast tissue. Additionally, we show that SNuBar can multiplex single cell ATAC and RNA multiomic assays in cell lines and human breast tissue samples. Our data show that SNuBar is a highly accurate, easy-to-use, and scalable system for multiplexing scATAC-seq and scATAC and RNA co-assay experiments.
10.1016/j.molcel.2021.09.026
Defining and identifying cell sub-crosstalk pairs for characterizing cell-cell communication patterns.
Scientific reports
Current cell-cell communication analysis focuses on quantifying intercellular interactions at cell type level. In the tissue microenvironment, one type of cells could be divided into multiple cell subgroups that function differently and communicate with other cell types or subgroups via different ligand-receptor-mediated signaling pathways. Given two cell types, we define a cell sub-crosstalk pair (CSCP) as a combination of two cell subgroups with strong and similar intercellular crosstalk signals and identify CSCPs based on coupled non-negative matrix factorization. Using single-cell spatial transcriptomics data of mouse olfactory bulb and visual cortex, we find that cells of different types within CSCPs are significantly spatially closer with each other than those in the whole single-cell spatial map. To demonstrate the utility of CSCPs, we apply 13 cell-cell communication analysis methods to sampled single-cell transcriptomics datasets at CSCP level and reveal ligand-receptor interactions masked at cell type level. Furthermore, by analyzing single-cell transcriptomics data from 29 breast cancer patients with different immunotherapy responses, we find that CSCPs are useful predictive features to discriminate patients responding to anti-PD-1 therapy from non-responders. Taken together, partitioning a cell type pair into CSCPs enables fine-grained characterization of cell-cell communication in tissue and tumor microenvironments.
10.1038/s41598-023-42883-8
Spatial multi-omic map of human myocardial infarction.
Nature
Myocardial infarction is a leading cause of death worldwide. Although advances have been made in acute treatment, an incomplete understanding of remodelling processes has limited the effectiveness of therapies to reduce late-stage mortality. Here we generate an integrative high-resolution map of human cardiac remodelling after myocardial infarction using single-cell gene expression, chromatin accessibility and spatial transcriptomic profiling of multiple physiological zones at distinct time points in myocardium from patients with myocardial infarction and controls. Multi-modal data integration enabled us to evaluate cardiac cell-type compositions at increased resolution, yielding insights into changes of the cardiac transcriptome and epigenome through the identification of distinct tissue structures of injury, repair and remodelling. We identified and validated disease-specific cardiac cell states of major cell types and analysed them in their spatial context, evaluating their dependency on other cell types. Our data elucidate the molecular principles of human myocardial tissue organization, recapitulating a gradual cardiomyocyte and myeloid continuum following ischaemic injury. In sum, our study provides an integrative molecular map of human myocardial infarction, represents an essential reference for the field and paves the way for advanced mechanistic and therapeutic studies of cardiac disease.
10.1038/s41586-022-05060-x
Differential gene expression patterns in vein regions susceptible versus resistant to neointimal hyperplasia.
Physiological genomics
Arteriovenous hemodialysis graft (AVG) stenosis results in thrombosis and AVG failure, but prevention of stenosis has been unsuccessful due in large part to our limited understanding of the molecular processes involved in neointimal hyperplasia (NH) formation. AVG stenosis develops chiefly as a consequence of highly localized NH formation in the vein-graft anastomosis region. Surprisingly, the vein region just downstream of the vein-graft anastomosis (herein termed proximal vein region) is relatively resistant to NH. We hypothesized that the gene expression profiles of the NH-prone and NH-resistant regions will be different from each other after graft placement, and analysis of their genomic profiles may yield potential therapeutic targets to prevent AVG stenosis. To test this, we evaluated the vein-graft anastomosis (NH-prone) and proximal vein (NH-resistant) regions in a porcine model of AVG stenosis with a porcine microarray. Gene expression changes in these two distinct vein regions, relative to the gene expression in unoperated control veins, were examined at early (5 days) and later (14 days) time points following graft placement. Global genomic changes were much greater in the NH-prone region than in the NH-resistant region at both time points. In the NH-prone region, genes related to regulation of cell proliferation and osteo-/chondrogenic vascular remodeling were most enriched among the significantly upregulated genes, and genes related to smooth muscle phenotype were significantly downregulated. These results provide insights into the spatial and temporal genomic modulation underlying NH formation in AVG and suggest potential therapeutic strategies to prevent and/or limit AVG stenosis.
10.1152/physiolgenomics.00082.2017
Molecular logic of cellular diversification in the mouse cerebral cortex.
Nature
The mammalian cerebral cortex has an unparalleled diversity of cell types, which are generated during development through a series of temporally orchestrated events that are under tight evolutionary constraint and are critical for proper cortical assembly and function. However, the molecular logic that governs the establishment and organization of cortical cell types remains unknown, largely due to the large number of cell classes that undergo dynamic cell-state transitions over extended developmental timelines. Here we generate a comprehensive atlas of the developing mouse neocortex, using single-cell RNA sequencing and single-cell assay for transposase-accessible chromatin using sequencing. We sampled the neocortex every day throughout embryonic corticogenesis and at early postnatal ages, and complemented the sequencing data with a spatial transcriptomics time course. We computationally reconstruct developmental trajectories across the diversity of cortical cell classes, and infer their spatial organization and the gene regulatory programs that accompany their lineage bifurcation decisions and differentiation trajectories. Finally, we demonstrate how this developmental map pinpoints the origin of lineage-specific developmental abnormalities that are linked to aberrant corticogenesis in mutant mice. The data provide a global picture of the regulatory mechanisms that govern cellular diversification in the neocortex.
10.1038/s41586-021-03670-5
Spatial mapping of thymic stromal microenvironments reveals unique features influencing T lymphoid differentiation.
Griffith Ann V,Fallahi Mohammad,Nakase Hiroshi,Gosink Mark,Young Brandon,Petrie Howard T
Immunity
Interaction of hematopoietic progenitors with the thymic microenvironment induces them to proliferate, adopt the T lineage fate, and asymmetrically diverge into multiple functional lineages. Progenitors at various developmental stages are stratified within the thymus, implying that the corresponding microenvironments provide distinct sets of signals to progenitors migrating between them. These differences remain largely undefined. Here we used physical and computational approaches to generate a comprehensive spatial map of stromal gene expression in the thymus. Although most stromal regions were characterized by a unique gene expression signature, the central cortex lacked distinctive features. Instead, a key function of this region appears to be the sequestration of unique microenvironments found at the cortical extremities, thus modulating the relative proximity of progenitors moving between them. Our findings compel reexamination of how cell migration, lineage specification, and proliferation are controlled by thymic architecture and provide an in-depth resource for global characterization of this control.
10.1016/j.immuni.2009.09.024
Toxicogenomics of endoplasmic reticulum stress inducer tunicamycin in the small intestine and liver of Nrf2 knockout and C57BL/6J mice.
Nair Sujit,Xu Changjiang,Shen Guoxiang,Hebbar Vidya,Gopalakrishnan Avantika,Hu Rong,Jain Mohit Raja,Liew Celine,Chan Jefferson Y,Kong Ah-Ng Tony
Toxicology letters
This objective of this study was to investigate the toxicogenomics and the spatial regulation of global gene expression profiles elicited by endoplasmic reticulum (ER) stress inducer tunicamycin (TM) in mouse small intestine and liver as well as to identify TM-modulated nuclear factor-E2-related factor 2 (Nrf2)-dependent genes. Gene expression profiles were analyzed using 45,000 Affymetrix mouse genome 430 2.0 array and GeneSpring 7.2 software. Microarray results were validated by quantitative real-time reverse transcription-PCR analyses. Clusters of genes that were either induced or suppressed more than two-fold by TM treatment compared with vehicle in C57BL/6J/Nrf2 (-/-; knockout) and C57BL/6J Nrf2 (+/+; wildtype) mice genotypes were identified. Amongst these, in small intestine and liver, 1291 and 750 genes, respectively, were identified as Nrf2-dependent and upregulated, and 1370 and 943 genes, respectively, as Nrf2-dependent and downregulated. Based on their biological functions, these genes can be categorized into molecular chaperones and heat shock proteins, ubiquitination/proteolysis, apoptosis/cell cycle, electron transport, detoxification, cell growth/differentiation, signaling molecules/interacting partners, kinases and phosphatases, transport, biosynthesis/metabolism, nuclear assembly and processing, and genes related to calcium and glucose homeostasis. Phase II detoxification/antioxidant genes as well as putative interacting partners of Nrf2 such as nuclear corepressors and coactivators, were also identified as Nrf2-dependent genes. The identification of TM-regulated and Nrf2-dependent genes in the unfolded protein response to ER stress not only provides potential novel insights into the gestalt biological effects of TM on the toxicogenomics and spatial regulation of global gene expression profiles in cancer pharmacology and toxicology, but also points to the pivotal role of Nrf2 in these biological processes.
10.1016/j.toxlet.2006.10.012
Recent Advances in Pathology: the 2021 Annual Review Issue of The Journal of Pathology.
Herrington C Simon,Poulsom Richard,Koeppen Hartmut,Coates Philip J
The Journal of pathology
The 2021 Annual Review Issue of The Journal of Pathology contains 14 invited reviews on current research areas of particular importance in pathology. The subjects included here reflect the broad range of interests covered by the journal, including both basic and applied research fields but always with the aim of improving our understanding of human disease. This year, our reviews encompass the huge impact of the COVID-19 pandemic, the development and application of biomarkers for immune checkpoint inhibitors, recent advances in multiplexing antigen/nucleic acid detection in situ, the use of genomics to aid drug discovery, organoid methodologies in research, the microbiome in cancer, the role of macrophage-stroma interactions in fibrosis, and TGF-β as a driver of fibrosis in multiple pathologies. Other reviews revisit the p53 field and its lack of clinical impact to date, dissect the genetics of mitochondrial diseases, summarise the cells of origin and genetics of sarcomagenesis, provide new data on the role of TRIM28 in tumour predisposition, review our current understanding of cancer stem cell niches, and the function and regulation of p63. The reviews are authored by experts in their field from academia and industry, and provide comprehensive updates of the chosen areas, in which there has been considerable recent progress. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.
10.1002/path.5687
The promising role of new molecular biomarkers in prostate cancer: from coding and non-coding genes to artificial intelligence approaches.
Prostate cancer and prostatic diseases
BACKGROUND:Risk stratification or progression in prostate cancer is performed with the support of clinical-pathological data such as the sum of the Gleason score and serum levels PSA. For several decades, methods aimed at the early detection of prostate cancer have included the determination of PSA serum levels. The aim of this systematic review is to provide an overview about recent advances in the discovery of new molecular biomarkers through transcriptomics, genomics and artificial intelligence that are expected to improve clinical management of the prostate cancer patient. METHODS:An exhaustive search was conducted by Pubmed, Google Scholar and Connected Papers using keywords relating to the genetics, genomics and artificial intelligence in prostate cancer, it includes "biomarkers", "non-coding RNAs", "lncRNAs", "microRNAs", "repetitive sequence", "prognosis", "prediction", "whole-genome sequencing", "RNA-Seq", "transcriptome", "machine learning", and "deep learning". RESULTS:New advances, including the search for changes in novel biomarkers such as mRNAs, microRNAs, lncRNAs, and repetitive sequences, are expected to contribute to an earlier and accurate diagnosis for each patient in the context of precision medicine, thus improving the prognosis and quality of life of patients. We analyze several aspects that are relevant for prostate cancer including its new molecular markers associated with diagnosis, prognosis, and prediction to therapy and how bioinformatic approaches such as machine learning and deep learning can contribute to clinic. Furthermore, we also include current techniques that will allow an earlier diagnosis, such as Spatial Transcriptomics, Exome Sequencing, and Whole-Genome Sequencing. CONCLUSION:Transcriptomic and genomic analysis have contributed to generate knowledge in the field of prostate carcinogenesis, new information about coding and non-coding genes as biomarkers has emerged. Synergies created by the implementation of artificial intelligence to analyze and understand sequencing data have allowed the development of clinical strategies that facilitate decision-making and improve personalized management in prostate cancer.
10.1038/s41391-022-00537-2
Molecularly targeted photothermal ablation improves tumor specificity and immune modulation in a rat model of hepatocellular carcinoma.
Muñoz Nina M,Dupuis Crystal,Williams Malea,Dixon Katherine,McWatters Amanda,Avritscher Rony,Bouchard Richard,Kaseb Ahmed,Schachtschneider Kyle M,Rao Arvind,Sheth Rahul A
Communications biology
Thermal ablation is a standard therapy for patients with hepatocellular carcinoma (HCC). Contemporary ablation devices are imperfect, as they lack tumor specificity. An ideal ablation modality would generate thermal energy only within tumoral tissue. Furthermore, as hyperthermia is known to influence tumor immunity, such a tumor-specific ablation modality may have the ability to favorably modulate the tumor immune landscape. Here we show a clinically relevant thermal ablation modality that generates tumor-specific hyperthermia, termed molecularly targeted photothermal ablation (MTPA), that is based upon the excellent localization of indocyanine green to HCC. In a syngeneic rat model, we demonstrate the tumor-specific hyperthermia generated by MTPA. We also show through spatial and transcriptomic profiling techniques that MTPA favorably modulates the intratumoral myeloid population towards tumor immunogenicity and diminishes the systemic release of oncogenic cytokines relative to conventional ablation modalities.
10.1038/s42003-020-01522-y
Human pericentromeric tandemly repeated DNA is transcribed at the end of oocyte maturation and is associated with membraneless mitochondria-associated structures.
Scientific reports
Most of the human genome is non-coding. However, some of the non-coding part is transcriptionally active. In humans, the tandemly repeated (TR) pericentromeric non-coding DNA-human satellites 2 and 3 (HS2, HS3)-are transcribed in somatic cells. These transcripts are also found in pre- and post-implantation embryos. The aim of this study was to analyze HS2/HS3 transcription and cellular localization of transcripts in human maturating oocytes. The maternal HS2/HS3 TR transcripts transcribed from both strands were accumulated in the ooplasm in GV-MI oocytes as shown by DNA-RNA FISH (fluorescence in-situ hybridization). The transcripts' content was higher in GV oocytes than in somatic cumulus cells according to real-time PCR. Using bioinformatics analysis, we demonstrated the presence of polyadenylated HS2 and HS3 RNAs in datasets of GV and MII oocyte transcriptomes. The transcripts shared a high degree of homology with HS2, HS3 transcripts previously observed in cancer cells. The HS2/HS3 transcripts were revealed by a combination of FISH and immunocytochemical staining within membraneless RNP structures that contained DEAD-box helicases DDX5 and DDX4. The RNP structures were closely associated with mitochondria, and are therefore similar to membraneless bodies described previously only in oogonia. These membraneless structures may be a site for spatial sequestration of RNAs and proteins in both maturating oocytes and cancer cells.
10.1038/s41598-020-76628-8
A fast small-sample kernel independence test for microbiome community-level association analysis.
Zhan Xiang,Plantinga Anna,Zhao Ni,Wu Michael C
Biometrics
To fully understand the role of microbiome in human health and diseases, researchers are increasingly interested in assessing the relationship between microbiome composition and host genomic data. The dimensionality of the data as well as complex relationships between microbiota and host genomics pose considerable challenges for analysis. In this article, we apply a kernel RV coefficient (KRV) test to evaluate the overall association between host gene expression and microbiome composition. The KRV statistic can capture nonlinear correlations and complex relationships among the individual data types and between gene expression and microbiome composition through measuring general dependency. Testing proceeds via a similar route as existing tests of the generalized RV coefficients and allows for rapid p-value calculation. Strategies to allow adjustment for confounding effects, which is crucial for avoiding misleading results, and to alleviate the problem of selecting the most favorable kernel are considered. Simulation studies show that KRV is useful in testing statistical independence with finite samples given the kernels are appropriately chosen, and can powerfully identify existing associations between microbiome composition and host genomic data while protecting type I error. We apply the KRV to a microbiome study examining the relationship between host transcriptome and microbiome composition within the context of inflammatory bowel disease and are able to derive new biological insights and provide formal inference on prior qualitative observations.
10.1111/biom.12684
Trespassing cancer cells: 'fingerprinting' invasive protrusions reveals metastatic culprits.
Klemke Richard L
Current opinion in cell biology
Metastatic cancer cells produce invasive membrane protrusions called invadopodia and pseudopodia, which play a central role in driving cancer cell dissemination in the body. Malignant cells use these structures to attach to and degrade extracellular matrix proteins, generate force for cell locomotion, and to penetrate the vasculature. Recent work using unique subcellular fractionation methodologies combined with spatial genomic, proteomic, and phosphoproteomic profiling has provided insight into the invadopodiome and pseudopodiome signaling networks that control the protrusion of invasive membranes. Here I highlight how these powerful spatial 'omics' approaches reveal important signatures of metastatic cancer cells and possible new therapeutic targets aimed at treating metastatic disease.
10.1016/j.ceb.2012.08.005
Host-Microbiota Interactions in the Esophagus During Homeostasis and Allergic Inflammation.
Gastroenterology
BACKGROUND & AIMS:Microbiota composition and mechanisms of host-microbiota interactions in the esophagus are unclear. We aimed to uncover fundamental information about the esophageal microbiome and its potential significance to eosinophilic esophagitis (EoE). METHODS:Microbiota composition, transplantation potential, and antibiotic responsiveness in the esophagus were established via 16S ribosomal RNA sequencing. Functional outcomes of microbiota colonization were assessed by RNA sequencing analysis of mouse esophageal epithelium and compared with the human EoE transcriptome. The impact of dysbiosis was assessed using a preclinical model of EoE. RESULTS:We found that the murine esophagus is colonized with diverse microbial communities within the first month of life. The esophageal microbiota is distinct, dominated by Lactobacillales, and demonstrates spatial heterogeneity as the proximal and distal esophagus are enriched in Bifidobacteriales and Lactobacillales, respectively. Fecal matter transplantation restores the esophageal microbiota, demonstrating that the local environment drives diversity. Microbiota colonization modifies esophageal tissue morphology and gene expression that is enriched in pathways associated with epithelial barrier function and overlapping with genes involved in EoE, including POSTN, KLK5, and HIF1A. Finally, neonatal antibiotic treatment reduces the abundance of Lactobacillales and exaggerates type 2 inflammation in the esophagus. Clinical data substantiated loss of esophageal Lactobacillales in EoE compared with controls. CONCLUSIONS:The esophagus has a unique microbiome with notable differences between its proximal and distal regions. Fecal matter transplantation restores the esophageal microbiome. Antibiotic-induced dysbiosis exacerbates disease in a murine model of EoE. Collectively, these data establish the composition, transplantation potential, antibiotic responsiveness, and host-microbiota interaction in the esophagus and have implications for gastrointestinal health and disease.
10.1053/j.gastro.2021.10.002
Interleukin-17 governs hypoxic adaptation of injured epithelium.
Science (New York, N.Y.)
Mammalian cells autonomously activate hypoxia-inducible transcription factors (HIFs) to ensure survival in low-oxygen environments. We report here that injury-induced hypoxia is insufficient to trigger HIF1α in damaged epithelium. Instead, multimodal single-cell and spatial transcriptomics analyses and functional studies reveal that retinoic acid-related orphan receptor γt (RORγt) γδ T cell-derived interleukin-17A (IL-17A) is necessary and sufficient to activate HIF1α. Protein kinase B (AKT) and extracellular signal-regulated kinase 1/2 (ERK1/2) signaling proximal of IL-17 receptor C (IL-17RC) activates mammalian target of rapamycin (mTOR) and consequently HIF1α. The IL-17A-HIF1α axis drives glycolysis in wound front epithelia. Epithelial-specific loss of IL-17RC, HIF1α, or blockade of glycolysis derails repair. Our findings underscore the coupling of inflammatory, metabolic, and migratory programs to expedite epithelial healing and illuminate the immune cell-derived inputs in cellular adaptation to hypoxic stress during repair.
10.1126/science.abg9302
Brain Gene Expression Pattern Correlated with the Differential Brain Activation by Pain and Touch in Humans.
Tang Jie,Su Qian,Zhang Xue,Qin Wen,Liu Huaigui,Liang Meng,Yu Chunshui
Cerebral cortex (New York, N.Y. : 1991)
Genes involved in pain and touch sensations have been studied extensively, but very few studies have tried to link them with neural activities in the brain. Here, we aimed to identify genes preferentially correlated to painful activation patterns by linking the spatial patterns of gene expression of Allen Human Brain Atlas with the pain-elicited neural responses in the human brain, with a parallel, control analysis for identification of genes preferentially correlated to tactile activation patterns. We identified 1828 genes whose expression patterns preferentially correlated to painful activation patterns and 411 genes whose expression patterns preferentially correlated to tactile activation pattern at the cortical level. In contrast to the enrichment for astrocyte and inhibitory synaptic transmission of genes preferentially correlated to tactile activation, the genes preferentially correlated to painful activation were mainly enriched for neuron and opioid- and addiction-related pathways and showed significant overlap with pain-related genes identified in previous studies. These findings not only provide important evidence for the differential genetic architectures of specific brain activation patterns elicited by painful and tactile stimuli but also validate a new approach to studying pain- and touch-related genes more directly from the perspective of neural responses in the human brain.
10.1093/cercor/bhab028
The -Mutant Consensus Molecular Subtype 3 Reveals an Immunosuppressive Tumor Microenvironment in Colorectal Cancer.
Cancers
Colorectal cancers (CRC) with mutations () are frequently included in consensus molecular subtype 3 (CMS3) with profound metabolic deregulation. We explored the transcriptomic impact of , focusing on the tumor microenvironment (TME) and pathways beyond metabolic deregulation. The status of in patients with CRC was investigated and overall survival (OS) was compared with wild-type (). Next, we identified CMS, and further investigated differentially expressed genes (DEG) of and distinctive pathways. Lastly, we used spatially resolved gene expression profiling to define the effect of in the TME regions of CMS3-classified CRC tissues. CRC patients with were mainly enriched in CMS3. Their specific enrichments of immune gene signatures in immunosuppressive TME were associated with worse OS. Activation of TGFβ signaling by was related to reduced pro-inflammatory and cytokine gene signatures, leading to suppression of immune infiltration. Digital spatial profiling in TME regions of CMS3-classified tissues suggested up-regulated genes, , , , , , and , that could be characteristic of immune suppression in TME. This study may help to depict the complex transcriptomic profile of in immunosuppressive TME. Future studies and clinical trials in CRC patients with should consider these transcriptional landscapes.
10.3390/cancers15041098
A 3D physio-mimetic interpenetrating network-based platform to decode the pro and anti-tumorigenic properties of cancer-associated fibroblasts.
Cao Huan,Cheng Hong Sheng,Wang Jun Kit,Tan Nguan Soon,Tay Chor Yong
Acta biomaterialia
Three-dimensional (3D) biomaterials with physiologically relevant and experimentally tractable biomechanical features are important platforms to advance our understanding of the influence of tissue mechanics in disease progression. Herein, an interpenetrating network (IPN) of collagen and alginate 3D culture system with tunable extracellular microstructure and mechanics is exploited as a tumor stroma proxy to study phenotypic plasticity of colorectal cancer-associated fibroblasts (CAF). In combination with Next Generation Sequencing (NGS) data analysis, we demonstrated that tuning the storage modulus of the IPN hydrogel between 49 and 419 Pa can trigger a reversible switch between an inflammatory (i-state, α-SMAIL-6) and myofibroblastic (m-state, α-SMAIL-6) state in CAF that is dependent on the polymer network confinement effect and ROS-HIF1-α mechanotransduction signaling axis. Secretome from m-state CAF upregulated several epithelial-mesenchymal-transition (EMT) transcripts and induced robust scattering in DLD-1, HCT116, and SW480 human colorectal adenocarcinoma, while the EMT-inducing capacity is muted in i-state CAF, suggestive of an anti-tumorigenic role. Our findings were further validated through Gene Expression Profiling Interactive Analysis (GEPIA), which showed that cytokines secreted at higher levels by i-state CAF are correlated (p < 0.05) with good overall colorectal cancer patient survival. Therefore, 3D network density and spatial cellular confinement are critical biophysical determinants that can profoundly influence CAF states, paracrine signaling, and EMT-inducing potential. STATEMENT OF SIGNIFICANCE: The communication between cancer cells and cancer-associated fibroblasts (CAF) contributes to tumor metastasis. CAF represent a diverse population of cellular subsets that can either promote or restrain tumor progression. However, the origin and cause of CAF heterogeneity remain elusive, limiting CAF-directed therapies for clinical use. We studied the dynamic phenotypes of CAF using a 3D physio-mimetic culture platform consisting of an interpenetrating collagen-alginate network. Combined with transcriptomic stratification and correlative analysis using cancer patient dataset, we showed phenotypic interconversion between inflammatory and myofibroblastic states, with anti- and pro-tumorigenic functions, in human colorectal CAF. This multidisciplinary study reveals the functional diversity of colorectal CAF caused by biophysical cues. The finding will influence the development of new CAF biomarkers and cancer therapies.
10.1016/j.actbio.2021.03.037
Deciphering sex-specific miRNAs as heat-recorders in zebrafish.
Scientific reports
In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.
10.1038/s41598-022-21864-3
Proteomics reveals the importance of the dynamic redistribution of the subcellular location of proteins in breast cancer cells.
Pinto Gabriella,Alhaiek Abdulrab Ahmed M,Godovac-Zimmermann Jasminka
Expert review of proteomics
At the molecular level, living cells are enormously complicated complex adaptive systems in which intertwined genomic, transcriptomic, proteomic and metabolic networks all play a crucial role. At the same time, cells are spatially heterogeneous systems in which subcellular compartmentalization of different functions is ubiquitous and requires efficient cross-compartmental communication. Dynamic redistribution of multitudinous proteins to different subcellular locations in response to cellular functional state is increasingly recognized as a crucial characteristic of cellular function that seems to be at least as important as overall changes in protein abundance. Characterization of the subcellular spatial dynamics of protein distribution is a major challenge for proteomics and recent results with MCF7 breast cancer cells suggest that this may be of particular importance for cancer cells.
10.1586/14789450.2015.1002474
Bi-directional and shared epigenomic signatures following proton and Fe irradiation.
Impey Soren,Jopson Timothy,Pelz Carl,Tafessu Amanuel,Fareh Fatema,Zuloaga Damian,Marzulla Tessa,Riparip Lara-Kirstie,Stewart Blair,Rosi Susanna,Turker Mitchell S,Raber Jacob
Scientific reports
The brain's response to radiation exposure is an important concern for patients undergoing cancer therapy and astronauts on long missions in deep space. We assessed whether this response is specific and prolonged and is linked to epigenetic mechanisms. We focused on the response of the hippocampus at early (2-weeks) and late (20-week) time points following whole body proton irradiation. We examined two forms of DNA methylation, cytosine methylation (5mC) and hydroxymethylation (5hmC). Impairments in object recognition, spatial memory retention, and network stability following proton irradiation were observed at the two-week time point and correlated with altered gene expression and 5hmC profiles that mapped to specific gene ontology pathways. Significant overlap was observed between DNA methylation changes at the 2 and 20-week time points demonstrating specificity and retention of changes in response to radiation. Moreover, a novel class of DNA methylation change was observed following an environmental challenge (i.e. space irradiation), characterized by both increased and decreased 5hmC levels along the entire gene body. These changes were mapped to genes encoding neuronal functions including postsynaptic gene ontology categories. Thus, the brain's response to proton irradiation is both specific and prolonged and involves novel remodeling of non-random regions of the epigenome.
10.1038/s41598-017-09191-4
Single-cell reconstruction of developmental trajectories during zebrafish embryogenesis.
Farrell Jeffrey A,Wang Yiqun,Riesenfeld Samantha J,Shekhar Karthik,Regev Aviv,Schier Alexander F
Science (New York, N.Y.)
During embryogenesis, cells acquire distinct fates by transitioning through transcriptional states. To uncover these transcriptional trajectories during zebrafish embryogenesis, we sequenced 38,731 cells and developed URD, a simulated diffusion-based computational reconstruction method. URD identified the trajectories of 25 cell types through early somitogenesis, gene expression along them, and their spatial origin in the blastula. Analysis of Nodal signaling mutants revealed that their transcriptomes were canalized into a subset of wild-type transcriptional trajectories. Some wild-type developmental branch points contained cells that express genes characteristic of multiple fates. These cells appeared to trans-specify from one fate to another. These findings reconstruct the transcriptional trajectories of a vertebrate embryo, highlight the concurrent canalization and plasticity of embryonic specification, and provide a framework with which to reconstruct complex developmental trees from single-cell transcriptomes.
10.1126/science.aar3131
Yeast Synthetic Minimal Biosensors for Evaluating Protein Production.
Peng Kai,Kroukamp Heinrich,Pretorius Isak S,Paulsen Ian T
ACS synthetic biology
The unfolded protein response (UPR) is a highly conserved cellular response in eukaryotic cells to counteract endoplasmic reticulum (ER) stress, typically triggered by unfolded protein accumulation. In addition to its relevance to human diseases like cancer, the induction of the UPR has a significant impact on the recombinant protein production in eukaryotic cell factories, including the industrial workhorse. Being able to accurately detect and measure this ER stress response in single cells, enables the rapid optimization of protein production conditions and high-throughput strain selection strategies. Current methodologies to monitor the UPR in are often temporally and spatially removed from the cultivation stage or lack updated systematic evaluation. To this end, we constructed and systematically evaluated a series of high-throughput UPR sensors by different designs, incorporating either yeast native UPR promoters or novel synthetic minimal UPR promoters. The native promoters of and were identified to have suitable UPR biosensor properties and served as an expression level guide for orthogonal sensor benchmarking. Our best synthetic minimal sensor is only 98 bp in length, has minimal homology to other native yeast sequences and displayed superior sensor characteristics. The synthetic minimal UPR sensor was able to accurately distinguish between cells expressing different heterologous proteins and between the different secretion levels of the same protein. This work demonstrated the potential of synthetic UPR biosensors as high-throughput tools to predict the protein production capacity of strains, interrogate protein properties hampering their secretion, and guide rational engineering strategies for optimal heterologous protein production.
10.1021/acssynbio.0c00633
MBD2 is required for correct spatial gene expression in the gut.
Molecular and cellular biology
Gene expression in the gut is segmentally regulated, but little is known of the molecular origin of patterning. Analysis of gene expression in colons from mice lacking the methyl-CpG binding repressor MBD2 revealed frequent activation of genes that are normally only expressed in the exocrine pancreas and duodenum. Reduced DNA methylation activated the same gene set in the colon. No significant differences in DNA methylation between the colon and duodenum were detected, but MBD2 was significantly more abundant in the colon. The relevance of MBD2 concentration was tested in a human colon cancer cell line. Depletion of MBD2 was again found to activate exocrine pancreatic genes. Gene activation in this cell culture model was accompanied by loss of promoter-bound MBD2 and increased histone acetylation. The results suggest that modulation of MBD2 during gut development establishes a region-specific gene expression pattern that is essential for establishing correct segmental character.
10.1128/MCB.02023-06
Cancer cell employs a microenvironmental neural signal trans-activating nucleus-mitochondria coordination to acquire stemness.
Signal transduction and targeted therapy
Cancer cell receives extracellular signal inputs to obtain a stem-like status, yet how tumor microenvironmental (TME) neural signals steer cancer stemness to establish the hierarchical tumor architectures remains elusive. Here, a pan-cancer transcriptomic screening for 10852 samples of 33 TCGA cancer types reveals that cAMP-responsive element (CRE) transcription factors are convergent activators for cancer stemness. Deconvolution of transcriptomic profiles, specification of neural markers and illustration of norepinephrine dynamics uncover a bond between TME neural signals and cancer-cell CRE activity. Specifically, neural signal norepinephrine potentiates the stemness of proximal cancer cells by activating cAMP-CRE axis, where ATF1 serves as a conserved hub. Upon activation by norepinephrine, ATF1 potentiates cancer stemness by coordinated trans-activation of both nuclear pluripotency factors MYC/NANOG and mitochondrial biogenesis regulators NRF1/TFAM, thereby orchestrating nuclear reprograming and mitochondrial rejuvenating. Accordingly, single-cell transcriptomes confirm the coordinated activation of nuclear pluripotency with mitochondrial biogenesis in cancer stem-like cells. These findings elucidate that cancer cell acquires stemness via a norepinephrine-ATF1 driven nucleus-mitochondria collaborated program, suggesting a spatialized stemness acquisition by hijacking microenvironmental neural signals.
10.1038/s41392-023-01487-4
Expression of zebrafish nos2b surrounds oral cavity.
Poon Kar-Lai,Richardson Michael,Korzh Vladimir
Developmental dynamics : an official publication of the American Association of Anatomists
Inducible nitric oxide synthase (NOS2) catalyzes the production of nitric oxide (NO), and is one of the factors establishing innate immunity. In zebrafish, Nos2 is represented by nos2a and nos2b. Here, we report the cloning and expression pattern of the zebrafish nos2b gene, which does not seem to participate in induced immune response. nos2b was mapped to zebrafish linkage group 15. The spatial and temporal expression pattern of nos2b in embryonic zebrafish was analyzed by whole-mount in situ hybridization. nos2b is expressed constitutively in two primordia located along the ventral midline. The first group of cells contributes to the neurohypophysis. Initially at the level of the ventral hindbrain, the second group of cells migrates closely with the thyroid primordium to its final position at the basihyal by 3 dpf. Thus, the analysis of expression pattern of nos2b reveals complex morphogenetic movements resulting in its expression surrounding the oral cavity.
10.1002/dvdy.21566
Stromal cell diversity associated with immune evasion in human triple-negative breast cancer.
Wu Sunny Z,Roden Daniel L,Wang Chenfei,Holliday Holly,Harvey Kate,Cazet Aurélie S,Murphy Kendelle J,Pereira Brooke,Al-Eryani Ghamdan,Bartonicek Nenad,Hou Rui,Torpy James R,Junankar Simon,Chan Chia-Ling,Lam Chuan En,Hui Mun N,Gluch Laurence,Beith Jane,Parker Andrew,Robbins Elizabeth,Segara Davendra,Mak Cindy,Cooper Caroline,Warrier Sanjay,Forrest Alistair,Powell Joseph,O'Toole Sandra,Cox Thomas R,Timpson Paul,Lim Elgene,Liu X Shirley,Swarbrick Alexander
The EMBO journal
The tumour stroma regulates nearly all stages of carcinogenesis. Stromal heterogeneity in human triple-negative breast cancers (TNBCs) remains poorly understood, limiting the development of stromal-targeted therapies. Single-cell RNA sequencing of five TNBCs revealed two cancer-associated fibroblast (CAF) and two perivascular-like (PVL) subpopulations. CAFs clustered into two states: the first with features of myofibroblasts and the second characterised by high expression of growth factors and immunomodulatory molecules. PVL cells clustered into two states consistent with a differentiated and immature phenotype. We showed that these stromal states have distinct morphologies, spatial relationships and functional properties in regulating the extracellular matrix. Using cell signalling predictions, we provide evidence that stromal-immune crosstalk acts via a diverse array of immunoregulatory molecules. Importantly, the investigation of gene signatures from inflammatory-CAFs and differentiated-PVL cells in independent TNBC patient cohorts revealed strong associations with cytotoxic T-cell dysfunction and exclusion, respectively. Such insights present promising candidates to further investigate for new therapeutic strategies in the treatment of TNBCs.
10.15252/embj.2019104063
A Knock-In Mouse Model of Thymoma With the GTF2I L424H Mutation.
Journal of thoracic oncology : official publication of the International Association for the Study of Lung Cancer
INTRODUCTION:The pathogenesis of thymic epithelial tumors remains largely unknown. We previously identified GTF2I L424H as the most frequently recurrent mutation in thymic epithelial tumors. Nevertheless, the precise role of this mutation in tumorigenesis of thymic epithelial cells is unclear. METHODS:To investigate the role of GTF2I L424H mutation in thymic epithelial cells in vivo, we generated and characterized a mouse model in which the Gtf2i L424H mutation was conditionally knocked-in in the Foxn1+ thymic epithelial cells. Digital spatial profiling was performed on thymomas and normal thymic tissues with GeoMx-mouse whole transcriptome atlas. Immunohistochemistry staining was performed using both mouse tissues and human thymic epithelial tumors. RESULTS:We observed that the Gtf2i mutation impairs development of the thymic medulla and maturation of medullary thymic epithelial cells in young mice and causes tumor formation in the thymus of aged mice. Cell cycle-related pathways, such as E2F targets and MYC targets, are enriched in the tumor epithelial cells. Results of gene set variation assay analysis revealed that gene signatures of cortical thymic epithelial cells and thymic epithelial progenitor cells are also enriched in the thymomas of the knock-in mice, which mirrors the human counterparts in The Cancer Genome Atlas database. Immunohistochemistry results revealed similar expression pattern of epithelial cell markers between mouse and human thymomas. CONCLUSIONS:We have developed and characterized a novel thymoma mouse model. This study improves knowledge of the molecular drivers in thymic epithelial cells and provides a tool for further study of the biology of thymic epithelial tumors and for development of novel therapies.
10.1016/j.jtho.2022.08.008
Dynamics of the 4D genome during in vivo lineage specification and differentiation.
Nature communications
Mammalian gene expression patterns are controlled by regulatory elements, which interact within topologically associating domains (TADs). The relationship between activation of regulatory elements, formation of structural chromatin interactions and gene expression during development is unclear. Here, we present Tiled-C, a low-input chromosome conformation capture (3C) technique. We use this approach to study chromatin architecture at high spatial and temporal resolution through in vivo mouse erythroid differentiation. Integrated analysis of chromatin accessibility and single-cell expression data shows that regulatory elements gradually become accessible within pre-existing TADs during early differentiation. This is followed by structural re-organization within the TAD and formation of specific contacts between enhancers and promoters. Our high-resolution data show that these enhancer-promoter interactions are not established prior to gene expression, but formed gradually during differentiation, concomitant with progressive upregulation of gene activity. Together, these results provide new insight into the close, interdependent relationship between chromatin architecture and gene regulation during development.
10.1038/s41467-020-16598-7
Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions.
Nature communications
In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.
10.1038/s41467-021-21361-7
Differential expression of the angiogenic factor genes vascular endothelial growth factor (VEGF) and endocrine gland-derived VEGF in normal and polycystic human ovaries.
Ferrara Napoleone,Frantz Gretchen,LeCouter Jennifer,Dillard-Telm Lisa,Pham Thinh,Draksharapu Aparna,Giordano Thomas,Peale Franklin
The American journal of pathology
Angiogenesis is a key aspect of the dynamic changes occurring during the normal ovarian cycle. Hyperplasia and hypervascularity of the ovarian theca interna and stroma are also prominent features of the polycystic ovary syndrome (PCOS), a leading cause of infertility. Compelling evidence indicated that vascular endothelial growth factor (VEGF) is a key mediator of the cyclical corpus luteum angiogenesis. However, the nature of the factor(s) that mediate angiogenesis in PCOS is less clearly understood. Endocrine gland-derived (EG)-VEGF has been recently identified as an endothelial cell mitogen with selectivity for the endothelium of steroidogenic glands and is expressed in normal human ovaries. In the present study, we compared the expression of EG-VEGF and VEGF mRNA in a series of 13 human PCOS and 13 normal ovary specimens by in situ hybridization. EG-VEGF expression in normal ovaries is dynamic and generally complementary to VEGF expression in both follicles and corpora lutea. A particularly high expression of EG-VEGF was detected in the Leydig-like hilus cells found in the highly vascularized ovarian hilus. In PCOS ovaries, we found strong expression of EG-VEGF mRNA in theca interna and stroma in most of the specimens examined, thus spatially related to the new blood vessels. In contrast, VEGF mRNA expression was most consistently associated with the granulosa cell layer and sometimes the theca, but rarely with the stroma. These findings indicate that both EG-VEGF and VEGF are expressed in PCOS ovaries, but in different cell types at different stages of differentiation, thus suggesting complementary functions for the two factors in angiogenesis and possibly cyst formation.
10.1016/S0002-9440(10)64322-2
Mycobacterium tuberculosis Exploits a Molecular Off Switch of the Immune System for Intracellular Survival.
von Both Ulrich,Berk Maurice,Agapow Paul-Michael,Wright Joseph D,Git Anna,Hamilton Melissa Shea,Goldgof Greg,Siddiqui Nazneen,Bellos Evangelos,Wright Victoria J,Coin Lachlan J,Newton Sandra M,Levin Michael
Scientific reports
Mycobacterium tuberculosis (M. tuberculosis) survives and multiplies inside human macrophages by subversion of immune mechanisms. Although these immune evasion strategies are well characterised functionally, the underlying molecular mechanisms are poorly understood. Here we show that during infection of human whole blood with M. tuberculosis, host gene transcriptional suppression, rather than activation, is the predominant response. Spatial, temporal and functional characterisation of repressed genes revealed their involvement in pathogen sensing and phagocytosis, degradation within the phagolysosome and antigen processing and presentation. To identify mechanisms underlying suppression of multiple immune genes we undertook epigenetic analyses. We identified significantly differentially expressed microRNAs with known targets in suppressed genes. In addition, after searching regions upstream of the start of transcription of suppressed genes for common sequence motifs, we discovered novel enriched composite sequence patterns, which corresponded to Alu repeat elements, transposable elements known to have wide ranging influences on gene expression. Our findings suggest that to survive within infected cells, mycobacteria exploit a complex immune "molecular off switch" controlled by both microRNAs and Alu regulatory elements.
10.1038/s41598-017-18528-y
Molecular architecture of lineage allocation and tissue organization in early mouse embryo.
Peng Guangdun,Suo Shengbao,Cui Guizhong,Yu Fang,Wang Ran,Chen Jun,Chen Shirui,Liu Zhiwen,Chen Guoyu,Qian Yun,Tam Patrick P L,Han Jing-Dong J,Jing Naihe
Nature
During post-implantation development of the mouse embryo, descendants of the inner cell mass in the early epiblast transit from the naive to primed pluripotent state. Concurrently, germ layers are formed and cell lineages are specified, leading to the establishment of the blueprint for embryogenesis. Fate-mapping and lineage-analysis studies have revealed that cells in different regions of the germ layers acquire location-specific cell fates during gastrulation. The regionalization of cell fates preceding the formation of the basic body plan-the mechanisms of which are instrumental for understanding embryonic programming and stem-cell-based translational study-is conserved in vertebrate embryos. However, a genome-wide molecular annotation of lineage segregation and tissue architecture of the post-implantation embryo has yet to be undertaken. Here we report a spatially resolved transcriptome of cell populations at defined positions in the germ layers during development from pre- to late-gastrulation stages. This spatiotemporal transcriptome provides high-resolution digitized in situ gene-expression profiles, reveals the molecular genealogy of tissue lineages and defines the continuum of pluripotency states in time and space. The transcriptome further identifies the networks of molecular determinants that drive lineage specification and tissue patterning, supports a role of Hippo-Yap signalling in germ-layer development and reveals the contribution of visceral endoderm to the endoderm in the early mouse embryo.
10.1038/s41586-019-1469-8
Alterations in intranuclear localization of Runx2 affect biological activity.
Zaidi Sayyed K,Javed Amjad,Pratap Jitesh,Schroeder Tania M,J Westendorf Jennifer,Lian Jane B,van Wijnen Andre J,Stein Gary S,Stein Janet L
Journal of cellular physiology
The transcription factor Runx2 controls osteoblast proliferation and differentiation. Runx2 organizes and assembles gene-regulatory complexes in nuclear microenvironments where target genes are activated or suppressed in a context-dependent manner. Intranuclear localization of Runx2 is mediated by the nuclear matrix-targeting signal (NMTS), an autonomous motif with a loop (L1)-turn-loop (L2) structure that forms predicted protein-protein interaction surfaces. Here we examined the functional consequences of introducing mutations in the L1 and L2 loops of the NMTS. These mutant proteins enter the nucleus, interact with the hetero-dimeric partner Cbfbeta, and bind to DNA in vitro and in vivo. In addition, these mutants retain interaction with the carboxy-terminus interacting co-regulatory proteins that include TLE, YAP, and Smads. However, two critical mutations in the L2 domain of the NMTS decrease association of Runx2 with the nuclear matrix. These subnuclear targeting defective (STD) mutants of Runx2 compromise target gene activation or repression. The biological significance of these findings is reflected by decreased osteogenic differentiation of mesenchymal progenitors, concomitant with major changes in gene expression profiles, upon expression of the STD Runx2 mutant. Our results demonstrate that fidelity of temporal and spatial localization of Runx2 within the nucleus is functionally linked with biological activity.
10.1002/jcp.20791
SFRP2 induces a mesenchymal subtype transition by suppression of SOX2 in glioblastoma.
Guo Min,Goudarzi Kaveh M,Abedi Shiva,Pieber Melanie,Sjöberg Elin,Behnan Jinan,Zhang Xing-Mei,Harris Robert A,Bartek Jiri,Lindström Mikael S,Nistér Monica,Hägerstrand Daniel
Oncogene
Intratumoral heterogeneity is a characteristic of glioblastomas that contain an intermixture of cell populations displaying different glioblastoma subtype gene expression signatures. Proportions of these populations change during tumor evolution, but the occurrence and regulation of glioblastoma subtype transition is not well described. To identify regulators of glioblastoma subtypes we utilized a combination of in vitro experiments and in silico analyses, using experimentally generated as well as publicly available data. Through this combined approach SOX2 was identified to confer a proneural glioblastoma subtype gene expression signature. SFRP2 was subsequently identified as a SOX2-antagonist, able to induce a mesenchymal glioblastoma subtype signature. A subset of patient glioblastoma samples with high SFRP2 and low SOX2 expression was particularly enriched with mesenchymal subtype samples. Phenotypically, SFRP2 decreased tumor sphere formation, stemness as assessed by limiting dilution assay, and overall cell proliferation but increased cell motility, whereas SOX2 induced the opposite effects. Furthermore, an SFRP2/non-canonical-WNT/KLF4/PDGFR/phospho-AKT/SOX2 signaling axis was found to be involved in the mesenchymal transition. Analysis of human tumor tissue spatial gene expression patterns showed distinct expression of SFRP2- and SOX2-correlated genes in vascular and cellular areas, respectively. Finally, conditioned media from SFRP2 overexpressing cells increased CD206 on macrophages. Together, these findings present SFRP2 as a SOX2-antagonist with the capacity to induce a mesenchymal subtype transition in glioma cells located in vascular tumor areas, highlighting its role in glioblastoma tumor evolution and intratumoral heterogeneity.
10.1038/s41388-021-01825-2
Coupled Proliferation and Apoptosis Maintain the Rapid Turnover of Microglia in the Adult Brain.
Cell reports
Microglia play key roles in brain development, homeostasis, and function, and it is widely assumed that the adult population is long lived and maintained by self-renewal. However, the precise temporal and spatial dynamics of the microglial population are unknown. We show in mice and humans that the turnover of microglia is remarkably fast, allowing the whole population to be renewed several times during a lifetime. The number of microglial cells remains steady from late postnatal stages until aging and is maintained by the spatial and temporal coupling of proliferation and apoptosis, as shown by pulse-chase studies, chronic in vivo imaging of microglia, and the use of mouse models of dysregulated apoptosis. Our results reveal that the microglial population is constantly and rapidly remodeled, expanding our understanding of its role in the maintenance of brain homeostasis.
10.1016/j.celrep.2016.12.041
MPL expression on AML blasts predicts peripheral blood neutropenia and thrombocytopenia.
Rauch Philipp J,Ellegast Jana M,Widmer Corinne C,Fritsch Kristin,Goede Jeroen S,Valk Peter J M,Löwenberg Bob,Takizawa Hitoshi,Manz Markus G
Blood
Although the molecular pathways that cause acute myeloid leukemia (AML) are increasingly well understood, the pathogenesis of peripheral blood cytopenia, a major cause of AML mortality, remains obscure. A prevailing assumption states that AML spatially displaces nonleukemic hematopoiesis from the bone marrow. However, examining an initial cohort of 223 AML patients, we found no correlation between bone marrow blast content and cytopenia, questioning the displacement theory. Measuring serum concentration of thrombopoietin (TPO), a key regulator of hematopoietic stem cells and megakaryocytes, revealed loss of physiologic negative correlation with platelet count in AML cases with blasts expressing MPL, the thrombopoietin (scavenging) receptor. Mechanistic studies demonstrated that MPL blasts could indeed clear TPO, likely therefore leading to insufficient cytokine levels for nonleukemic hematopoiesis. Microarray analysis in an independent multicenter study cohort of 437 AML cases validated MPL expression as a central predictor of thrombocytopenia and neutropenia in AML. Moreover, t(8;21) AML cases demonstrated the highest average MPL expression and lowest average platelet and absolute neutrophil counts among subgroups. Our work thus explains the pathophysiology of peripheral blood cytopenia in a relevant number of AML cases.
10.1182/blood-2016-04-711986
Distinct skeletal stem cell types orchestrate long bone skeletogenesis.
eLife
Skeletal stem and progenitor cell populations are crucial for bone physiology. Characterization of these cell types remains restricted to heterogenous bulk populations with limited information on whether they are unique or overlap with previously characterized cell types. Here we show, through comprehensive functional and single-cell transcriptomic analyses, that postnatal long bones of mice contain at least two types of bone progenitors with bona fide skeletal stem cell (SSC) characteristics. An early osteochondral SSC (ocSSC) facilitates long bone growth and repair, while a second type, a perivascular SSC (pvSSC), co-emerges with long bone marrow and contributes to shape the hematopoietic stem cell niche and regenerative demand. We establish that pvSSCs, but not ocSSCs, are the origin of bone marrow adipose tissue. Lastly, we also provide insight into residual SSC heterogeneity as well as potential crosstalk between the two spatially distinct cell populations. These findings comprehensively address previously unappreciated shortcomings of SSC research.
10.7554/eLife.66063
Spatial control of cell fate using synthetic surfaces to potentiate TGF-beta signaling.
Li Lingyin,Klim Joseph R,Derda Ratmir,Courtney Adam H,Kiessling Laura L
Proceedings of the National Academy of Sciences of the United States of America
In organisms, cell-fate decisions result from external cues presented by the extracellular microenvironment or the niche. In principle, synthetic niches can be engineered to give rise to patterned cell signaling, an advance that would transform the fields of tissue engineering and regenerative medicine. Biomaterials that display adhesive motifs are critical steps in this direction, but promoting localized signaling remains a major obstacle. We sought to exert precise spatial control over activation of TGF-β signaling. TGF-β signaling, which plays fundamental roles in development, tissue homeostasis, and cancer, is initiated by receptor oligomerization. We therefore hypothesized that preorganizing the transmembrane receptors would potentiate local TGF-β signaling. To generate surfaces that would nucleate the signaling complex, we employed defined self-assembled monolayers that present peptide ligands to TGF-β receptors. These displays of nondiffusible ligands do not compete with the growth factor but rather sensitize bound cells to subpicomolar concentrations of endogenous TGF-β. Cells adhering to the surfaces undergo TGF-β-mediated growth arrest and the epithelial to mesenchymal transition. Gene expression profiles reveal that the surfaces selectively regulate TGF-β responsive genes. This strategy provides access to tailored surfaces that can deliver signals with spatial control.
10.1073/pnas.1101454108
Transcription and Signaling Regulators in Developing Neuronal Subtypes of Mouse and Human Enteric Nervous System.
Memic Fatima,Knoflach Viktoria,Morarach Khomgrit,Sadler Rebecca,Laranjeira Catia,Hjerling-Leffler Jens,Sundström Erik,Pachnis Vassilis,Marklund Ulrika
Gastroenterology
BACKGROUND & AIMS:The enteric nervous system (ENS) regulates gastrointestinal function via different subtypes of neurons, organized into fine-tuned neural circuits. It is not clear how cell diversity is created within the embryonic ENS; information required for development of cell-based therapies and models of enteric neuropathies. We aimed to identify proteins that regulate ENS differentiation and network formation. METHODS:We generated and compared RNA expression profiles of the entire ENS, ENS progenitor cells, and non-ENS gut cells of mice, collected at embryonic days 11.5 and 15.5, when different subtypes of neurons are formed. Gastrointestinal tissues from R26ReYFP reporter mice crossed to Sox10-CreER or Wnt1-Cre mice were dissected and the 6 populations of cells were isolated by flow cytometry. We used histochemistry to map differentially expressed proteins in mouse and human gut tissues at different stages of development, in different regions. We examined enteric neuronal diversity and gastric function in Wnt1-Cre x Sox6 mice, which do not express the Sox6 gene in the ENS. RESULTS:We identified 147 transcription and signaling factors that varied in spatial and temporal expression during development of the mouse ENS. Of the factors also analyzed in human ENS, most were conserved. We uncovered 16 signaling pathways (such as fibroblast growth factor and Eph/ephrin pathways). Transcription factors were grouped according to their specific expression in enteric progenitor cells (such as MEF2C), enteric neurons (such as SOX4), or neuron subpopulations (such as SATB1 and SOX6). Lack of SOX6 in the ENS reduced the numbers of gastric dopamine neurons and delayed gastric emptying. CONCLUSIONS:Using transcriptome and histochemical analyses of the developing mouse and human ENS, we mapped expression patterns of transcription and signaling factors. Further studies of these candidate determinants might elucidate the mechanisms by which enteric stem cells differentiate into neuronal subtypes and form distinct connectivity patterns during ENS development. We found expression of SOX6 to be required for development of gastric dopamine neurons.
10.1053/j.gastro.2017.10.005
Interpretable dimensionality reduction of single cell transcriptome data with deep generative models.
Nature communications
Single-cell RNA-sequencing has great potential to discover cell types, identify cell states, trace development lineages, and reconstruct the spatial organization of cells. However, dimension reduction to interpret structure in single-cell sequencing data remains a challenge. Existing algorithms are either not able to uncover the clustering structures in the data or lose global information such as groups of clusters that are close to each other. We present a robust statistical model, scvis, to capture and visualize the low-dimensional structures in single-cell gene expression data. Simulation results demonstrate that low-dimensional representations learned by scvis preserve both the local and global neighbor structures in the data. In addition, scvis is robust to the number of data points and learns a probabilistic parametric mapping function to add new data points to an existing embedding. We then use scvis to analyze four single-cell RNA-sequencing datasets, exemplifying interpretable two-dimensional representations of the high-dimensional single-cell RNA-sequencing data.
10.1038/s41467-018-04368-5
Small RNA profiling of Xenopus embryos reveals novel miRNAs and a new class of small RNAs derived from intronic transposable elements.
Genome research
Small RNA control of gene expression is critical for developmental processes in vertebrate embryos. To determine the dynamics of small RNA expression and to uncover novel small RNAs in the early vertebrate embryo, we performed high-throughput sequencing of all small RNAs in Xenopus tropicalis embryos at three developmental time points and in dissected halves of gastrula embryos. This analysis allowed us to identify novel microRNAs and we show that microRNA expression is highly dynamic and spatially localized in early embryos. In addition, we have developed a microRNA prediction pipeline and demonstrate that it has the power to predict new miRNAs that are experimentally detectable in frogs, mice, and humans. By combining the small RNA sequencing with mRNA profiling at the different developmental stages, we identify a new class of small noncoding RNAs that we name siteRNAs, which align in clusters to introns of protein-coding genes. We show that siteRNAs are derived from remnants of transposable elements present in the introns. We find that genes containing clusters of siteRNAs are transcriptionally repressed as compared with all genes. Furthermore, we show that this is true for individual genes containing siteRNA clusters, and that these genes are enriched in specific repressive histone modifications. Our data thus suggest a new mechanism of siteRNA-mediated gene silencing in vertebrates, and provide an example of how mobile elements can affect gene regulation.
10.1101/gr.144469.112
Mapping the temporal and spatial dynamics of the human endometrium in vivo and in vitro.
Nature genetics
The endometrium, the mucosal lining of the uterus, undergoes dynamic changes throughout the menstrual cycle in response to ovarian hormones. We have generated dense single-cell and spatial reference maps of the human uterus and three-dimensional endometrial organoid cultures. We dissect the signaling pathways that determine cell fate of the epithelial lineages in the lumenal and glandular microenvironments. Our benchmark of the endometrial organoids reveals the pathways and cell states regulating differentiation of the secretory and ciliated lineages both in vivo and in vitro. In vitro downregulation of WNT or NOTCH pathways increases the differentiation efficiency along the secretory and ciliated lineages, respectively. We utilize our cellular maps to deconvolute bulk data from endometrial cancers and endometriotic lesions, illuminating the cell types dominating in each of these disorders. These mechanistic insights provide a platform for future development of treatments for common conditions including endometriosis and endometrial carcinoma.
10.1038/s41588-021-00972-2
A novel Drosophila antisense scaRNA with a predicted guide function.
Tortoriello Giuseppe,Accardo Maria Carmela,Scialò Filippo,Angrisani Alberto,Turano Mimmo,Furia Maria
Gene
A significant portion of eukaryotic small ncRNA transcriptome is composed by small nucleolar RNAs. From archaeal to mammalian cells, these molecules act as guides in the site-specific pseudouridylation or methylation of target RNAs. We used a bioinformatics search program to detect Drosophila putative orthologues of U79, one out of ten snoRNAs produced by GAS5, a human ncRNA involved in apoptosis, susceptibility to cancer and autoimmune diseases. This search led to the definition of a list of U79-related fly snoRNAs whose genomic organization, evolution and expression strategy are discussed here. We report that an intriguing novel specimen, named Dm46E3, is transcribed as a longer, unspliced precursor from the reverse strand of eiger, a fly regulatory gene that plays a key role in cell differentiation, apoptosis and immune response. Expression of Dm46E3 was found significantly up-regulated in a mutant strain in which eiger transcription is greatly reduced, suggesting that these two sense-antisense genes may be mutually regulated. Relevant to its function, Dm46E3 concentrated specifically in the Cajal bodies, followed a dynamic spatial expression profile during embryogenesis and displayed a degenerate antisense element that enables it to target U1b, a developmentally regulated isoform of the U1 spliceosomal snRNA that is particularly abundant in embryos.
10.1016/j.gene.2009.02.005
OX40 and 4-1BB delineate distinct immune profiles in sarcoma.
Oncoimmunology
Systemic relapse after radiotherapy and surgery is the major cause of disease-related mortality in sarcoma patients. Combining radiotherapy and immunotherapy is under investigation as a means to improve response rates. However, the immune contexture of sarcoma is understudied. Here, we use a retrospective cohort of sarcoma patients, treated with neoadjuvant radiotherapy, and TCGA data. We explore therapeutic targets of relevance to sarcoma, using genomics and multispectral immunohistochemistry to provide insights into the tumor immune microenvironment across sarcoma subtypes. Differential gene expression between radioresponsive myxoid liposarcoma (MLPS) and more radioresistant undifferentiated pleomorphic sarcoma (UPS) indicated UPS contained higher transcript levels of a number of immunotherapy targets (CD73/, CD39/, CD25/, and 4-1BB/). We focused on 4-1BB/ and other costimulatory molecules. In TCGA data, 4-1BB correlated to an inflamed and exhausted phenotype. OX40/ and 4-1BB/ were highly expressed in sarcoma subtypes versus other cancers. Despite OX40 and 4-1BB being described as Treg markers, we identified that they delineate distinct tumor immune profiles. This was true for sarcoma and other cancers. While only a limited number of samples could be analyzed, spatial analysis of OX40 expression identified two diverse phenotypes of OX40+ Tregs, one associated with and one independent of tertiary lymphoid structures (TLSs). Patient stratification is of intense interest for immunotherapies. We provide data supporting the viewpoint that a cohort of sarcoma patients, appropriately selected, are promising candidates for immunotherapies. Spatial profiling of OX40+ Tregs, in relation to TLSs, could be an additional metric to improve future patient stratification.
10.1080/2162402X.2022.2066050
A spatially resolved timeline of the human maternal-fetal interface.
Nature
Beginning in the first trimester, fetally derived extravillous trophoblasts (EVTs) invade the uterus and remodel its spiral arteries, transforming them into large, dilated blood vessels. Several mechanisms have been proposed to explain how EVTs coordinate with the maternal decidua to promote a tissue microenvironment conducive to spiral artery remodelling (SAR). However, it remains a matter of debate regarding which immune and stromal cells participate in these interactions and how this evolves with respect to gestational age. Here we used a multiomics approach, combining the strengths of spatial proteomics and transcriptomics, to construct a spatiotemporal atlas of the human maternal-fetal interface in the first half of pregnancy. We used multiplexed ion beam imaging by time-of-flight and a 37-plex antibody panel to analyse around 500,000 cells and 588 arteries within intact decidua from 66 individuals between 6 and 20 weeks of gestation, integrating this dataset with co-registered transcriptomics profiles. Gestational age substantially influenced the frequency of maternal immune and stromal cells, with tolerogenic subsets expressing CD206, CD163, TIM-3, galectin-9 and IDO-1 becoming increasingly enriched and colocalized at later time points. By contrast, SAR progression preferentially correlated with EVT invasion and was transcriptionally defined by 78 gene ontology pathways exhibiting distinct monotonic and biphasic trends. Last, we developed an integrated model of SAR whereby invasion is accompanied by the upregulation of pro-angiogenic, immunoregulatory EVT programmes that promote interactions with the vascular endothelium while avoiding the activation of maternal immune cells.
10.1038/s41586-023-06298-9
Decoding myofibroblast origins in human kidney fibrosis.
Nature
Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.
10.1038/s41586-020-2941-1
An Integrated Global Analysis of Compartmentalized HRAS Signaling.
Santra Tapesh,Herrero Ana,Rodriguez Javier,von Kriegsheim Alex,Iglesias-Martinez Luis F,Schwarzl Thomas,Higgins Des,Aye Thin-Thin,Heck Albert J R,Calvo Fernando,Agudo-Ibáñez Lorena,Crespo Piero,Matallanas David,Kolch Walter
Cell reports
Modern omics technologies allow us to obtain global information on different types of biological networks. However, integrating these different types of analyses into a coherent framework for a comprehensive biological interpretation remains challenging. Here, we present a conceptual framework that integrates protein interaction, phosphoproteomics, and transcriptomics data. Applying this method to analyze HRAS signaling from different subcellular compartments shows that spatially defined networks contribute specific functions to HRAS signaling. Changes in HRAS protein interactions at different sites lead to different kinase activation patterns that differentially regulate gene transcription. HRAS-mediated signaling is the strongest from the cell membrane, but it regulates the largest number of genes from the endoplasmic reticulum. The integrated networks provide a topologically and functionally resolved view of HRAS signaling. They reveal distinct HRAS functions including the control of cell migration from the endoplasmic reticulum and TP53-dependent cell survival when signaling from the Golgi apparatus.
10.1016/j.celrep.2019.02.038
Prognostic, therapeutic and diagnostic potential of microRNAs in non-small cell lung cancer.
Markou Athina,Liang Yu,Lianidou Evi
Clinical chemistry and laboratory medicine
Non-small cell lung carcinomas (NSCLC) account for about 80% of lung cancers and their remarkable heterogeneity manifests in histology, pathogenesis, prognosis, and response to treatments. Recent advances in molecular characterization help stratifying NSCLC patients for their potential benefit from targeting therapies. However, the fundamental mechanisms underlying the tumoral heterogeneity remain poorly understood. Expression profiling of many microRNAs (miRNAs) in various normal and disease tissues demonstrated unique spatial and temporal expression patterns and some miRNAs have been functionally characterized as oncogenes or tumor suppressor genes. Genome-wide screening identified specific miRNA expression signatures associated with clinical outcome of NSCLC patients. A group of miRNAs that has enriched expression in normal lung was found down regulated in NSCLC and may function as tumor suppressor genes. In this review we: a) summarize the current understanding of the critical role that miRNAs play in normal cell functions and in disease biology especially in lung cancer tumorigenesis, b) highlight their potential as biomarkers for lung cancer risk stratification, outcome prediction and classification of histologic subtypes, c) critically assess current knowledge on lung-enriched miRNAs and expression of their predicted target genes in NSCLC and d) evaluate their potential as circulating biomarkers and therapeutic targets in lung cancer.
10.1515/CCLM.2011.661
CYP79F1 and CYP79F2 have distinct functions in the biosynthesis of aliphatic glucosinolates in Arabidopsis.
Chen Sixue,Glawischnig Erich,Jørgensen Kirsten,Naur Peter,Jørgensen Bodil,Olsen Carl-Erik,Hansen Carsten H,Rasmussen Hasse,Pickett John A,Halkier Barbara A
The Plant journal : for cell and molecular biology
Cytochromes P450 of the CYP79 family catalyze the conversion of amino acids to oximes in the biosynthesis of glucosinolates, a group of natural plant products known to be involved in plant defense and as a source of flavor compounds, cancer-preventing agents and bioherbicides. We report a detailed biochemical analysis of the substrate specificity and kinetics of CYP79F1 and CYP79F2, two cytochromes P450 involved in the biosynthesis of aliphatic glucosinolates in Arabidopsis thaliana. Using recombinant CYP79F1 and CYP79F2 expressed in Escherichia coli and Saccharomyces cerevisiae, respectively, we show that CYP79F1 metabolizes mono- to hexahomomethionine, resulting in both short- and long-chain aliphatic glucosinolates. In contrast, CYP79F2 exclusively metabolizes long-chain elongated penta- and hexahomomethionines. CYP79F1 and CYP79F2 are spatially and developmentally regulated, with different gene expression patterns. CYP79F2 is highly expressed in hypocotyl and roots, whereas CYP79F1 is strongly expressed in cotyledons, rosette leaves, stems, and siliques. A transposon-tagged CYP79F1 knockout mutant completely lacks short-chain aliphatic glucosinolates, but has an increased level of long-chain aliphatic glucosinolates, especially in leaves and seeds. The level of long-chain aliphatic glucosinolates in a transposon-tagged CYP79F2 knockout mutant is substantially reduced, whereas the level of short-chain aliphatic glucosinolates is not affected. Biochemical characterization of CYP79F1 and CYP79F2, and gene expression analysis, combined with glucosinolate profiling of knockout mutants demonstrate the functional role of these enzymes. This provides valuable insights into the metabolic network leading to the biosynthesis of aliphatic glucosinolates, and into metabolic engineering of altered aliphatic glucosinolate profiles to improve nutritional value and pest resistance.
Disruption of NIPBL/Scc2 in Cornelia de Lange Syndrome provokes cohesin genome-wide redistribution with an impact in the transcriptome.
Nature communications
Cornelia de Lange syndrome (CdLS) is a rare disease affecting multiple organs and systems during development. Mutations in the cohesin loader, NIPBL/Scc2, were first described and are the most frequent in clinically diagnosed CdLS patients. The molecular mechanisms driving CdLS phenotypes are not understood. In addition to its canonical role in sister chromatid cohesion, cohesin is implicated in the spatial organization of the genome. Here, we investigate the transcriptome of CdLS patient-derived primary fibroblasts and observe the downregulation of genes involved in development and system skeletal organization, providing a link to the developmental alterations and limb abnormalities characteristic of CdLS patients. Genome-wide distribution studies demonstrate a global reduction of NIPBL at the NIPBL-associated high GC content regions in CdLS-derived cells. In addition, cohesin accumulates at NIPBL-occupied sites at CpG islands potentially due to reduced cohesin translocation along chromosomes, and fewer cohesin peaks colocalize with CTCF.
10.1038/s41467-021-24808-z
Molecular, spatial, and functional single-cell profiling of the hypothalamic preoptic region.
Science (New York, N.Y.)
The hypothalamus controls essential social behaviors and homeostatic functions. However, the cellular architecture of hypothalamic nuclei-including the molecular identity, spatial organization, and function of distinct cell types-is poorly understood. Here, we developed an imaging-based in situ cell-type identification and mapping method and combined it with single-cell RNA-sequencing to create a molecularly annotated and spatially resolved cell atlas of the mouse hypothalamic preoptic region. We profiled ~1 million cells, identified ~70 neuronal populations characterized by distinct neuromodulatory signatures and spatial organizations, and defined specific neuronal populations activated during social behaviors in male and female mice, providing a high-resolution framework for mechanistic investigation of behavior circuits. The approach described opens a new avenue for the construction of cell atlases in diverse tissues and organisms.
10.1126/science.aau5324
The Wnt3a/β-catenin target gene Mesogenin1 controls the segmentation clock by activating a Notch signalling program.
Nature communications
Segmentation is an organizing principle of body plans. The segmentation clock, a molecular oscillator best illustrated by the cyclic expression of Notch signalling genes, controls the periodic cleavage of somites from unsegmented presomitic mesoderm during vertebrate segmentation. Wnt3a controls the spatiotemporal expression of cyclic Notch genes; however, the underlying mechanisms remain obscure. Here we show by transcriptional profiling of Wnt3a (-/-) embryos that the bHLH transcription factor, Mesogenin1 (Msgn1), is a direct target gene of Wnt3a. To identify Msgn1 targets, we conducted genome-wide studies of Msgn1 activity in embryonic stem cells. We show that Msgn1 is a major transcriptional activator of a Notch signalling program and synergizes with Notch to trigger clock gene expression. Msgn1 also indirectly regulates cyclic genes in the Fgf and Wnt pathways. Thus, Msgn1 is a central component of a transcriptional cascade that translates a spatial Wnt3a gradient into a temporal pattern of clock gene expression.
10.1038/ncomms1381
Spatiotemporal MicroRNA-Gene Expression Network Related to Orofacial Clefts.
Journal of dental research
Craniofacial structures change dynamically in morphology during development through the coordinated regulation of various cellular molecules. However, it remains unclear how these complex mechanisms are regulated in a spatiotemporal manner. Here we applied natural cubic splines to model gene and microRNA (miRNA) expression from embryonic day (E) 10.5 to E14.5 in the proximal and distal regions of the maxillary processes to identify spatiotemporal patterns of gene and miRNA expression, followed by constructing corresponding regulatory networks. Three major groups of differentially expressed genes (DEGs) were identified, including 3,927 temporal, 314 spatial, and 494 spatiotemporal DEGs. Unsupervised clustering further resolved these spatiotemporal DEGs into 8 clusters with distinct expression patterns. Interestingly, we found 2 clusters of differentially expressed miRNAs: 1 had 80 miRNAs monotonically decreasing and the other had 97 increasing across developmental stages. To evaluate the phenotypic relevance of these DEGs during craniofacial development, we integrated data from the CleftGeneDB database and constructed the regulatory networks of genes related to orofacial clefts. Our analysis revealed 2 hub miRNAs, mmu-miR-325-3p and mmu-miR-384-5p, that repressed cleft-related genes , , , , and , while their expression increased over time. On the contrary, 2 hub miRNAs, mmu-miR-218-5p and mmu-miR-338-5p, repressed cleft-related genes , , , , and , while their expression decreased over time. Our experiments indicated that these miRNA mimics significantly inhibited cell proliferation in mouse embryonic palatal mesenchymal (MEPM) cells and O9-1 cells through the regulation of genes associated with cleft palate and validated the role of our regulatory networks in orofacial clefts. To facilitate interactive exploration of these data, we developed a user-friendly web tool to visualize the gene and miRNA expression patterns across developmental stages, as well as the regulatory networks (https://fyan.shinyapps.io/facebase_shiny/). Taken together, our results provide a valuable resource that serves as a reference map for future research in craniofacial development.
10.1177/00220345221105816
Single-cell transcriptome analyses reveal microglia types associated with proliferative retinopathy.
JCI insight
Pathological angiogenesis is a major cause of irreversible blindness in individuals of all age groups with proliferative retinopathy (PR). Mononuclear phagocytes (MPs) within neovascular areas contribute to aberrant retinal angiogenesis. Due to their cellular heterogeneity, defining the roles of MP subsets in PR onset and progression has been challenging. Here, we aimed to investigate the heterogeneity of microglia associated with neovascularization and to characterize the transcriptional profiles and metabolic pathways of proangiogenic microglia in a mouse model of oxygen-induced PR (OIR). Using transcriptional single-cell sorting, we comprehensively mapped all microglia populations in retinas of room air (RA) and OIR mice. We have unveiled several unique types of PR-associated microglia (PRAM) and identified markers, signaling pathways, and regulons associated with these cells. Among these microglia subpopulations, we found a highly proliferative microglia subset with high self-renewal capacity and a hypermetabolic microglia subset that expresses high levels of activating microglia markers, glycolytic enzymes, and proangiogenic Igf1. IHC staining shows that these PRAM were spatially located within or around neovascular tufts. These unique types of microglia have the potential to promote retinal angiogenesis, which may have important implications for future treatment of PR and other pathological ocular angiogenesis-related diseases.
10.1172/jci.insight.160940
Brg1 coordinates multiple processes during retinogenesis and is a tumor suppressor in retinoblastoma.
Aldiri Issam,Ajioka Itsuki,Xu Beisi,Zhang Jiakun,Chen Xiang,Benavente Claudia,Finkelstein David,Johnson Dianna,Akiyama Jennifer,Pennacchio Len A,Dyer Michael A
Development (Cambridge, England)
Retinal development requires precise temporal and spatial coordination of cell cycle exit, cell fate specification, cell migration and differentiation. When this process is disrupted, retinoblastoma, a developmental tumor of the retina, can form. Epigenetic modulators are central to precisely coordinating developmental events, and many epigenetic processes have been implicated in cancer. Studying epigenetic mechanisms in development is challenging because they often regulate multiple cellular processes; therefore, elucidating the primary molecular mechanisms involved can be difficult. Here we explore the role of Brg1 (Smarca4) in retinal development and retinoblastoma in mice using molecular and cellular approaches. Brg1 was found to regulate retinal size by controlling cell cycle length, cell cycle exit and cell survival during development. Brg1 was not required for cell fate specification but was required for photoreceptor differentiation and cell adhesion/polarity programs that contribute to proper retinal lamination during development. The combination of defective cell differentiation and lamination led to retinal degeneration in Brg1-deficient retinae. Despite the hypocellularity, premature cell cycle exit, increased cell death and extended cell cycle length, retinal progenitor cells persisted in Brg1-deficient retinae, making them more susceptible to retinoblastoma. ChIP-Seq analysis suggests that Brg1 might regulate gene expression through multiple mechanisms.
10.1242/dev.124800
Single-cell multi-omics analysis of human testicular germ cell tumor reveals its molecular features and microenvironment.
Nature communications
Seminoma is the most common malignant solid tumor in 14 to 44 year-old men. However, its molecular features and tumor microenvironment (TME) is largely unexplored. Here, we perform a series of studies via genomics profiling (single cell multi-omics and spatial transcriptomics) and functional examination using seminoma samples and a seminoma cell line. We identify key gene expression programs share between seminoma and primordial germ cells, and further characterize the functions of TFAP2C in promoting tumor invasion and migration. We also identify 15 immune cell subtypes in TME, and find that subtypes with exhaustion features were located closer to the tumor region through combined spatial transcriptome analysis. Furthermore, we identify key pathways and genes that may facilitate seminoma disseminating beyond the seminiferous tubules. These findings advance our knowledge of seminoma tumorigenesis and produce a multi-omics atlas of in situ human seminoma microenvironment, which could help discover potential therapy targets for seminoma.
10.1038/s41467-023-44305-9
Hierarchical patterning modes orchestrate hair follicle morphogenesis.
PLoS biology
Two theories address the origin of repeating patterns, such as hair follicles, limb digits, and intestinal villi, during development. The Turing reaction-diffusion system posits that interacting diffusible signals produced by static cells first define a prepattern that then induces cell rearrangements to produce an anatomical structure. The second theory, that of mesenchymal self-organisation, proposes that mobile cells can form periodic patterns of cell aggregates directly, without reference to any prepattern. Early hair follicle development is characterised by the rapid appearance of periodic arrangements of altered gene expression in the epidermis and prominent clustering of the adjacent dermal mesenchymal cells. We assess the contributions and interplay between reaction-diffusion and mesenchymal self-organisation processes in hair follicle patterning, identifying a network of fibroblast growth factor (FGF), wingless-related integration site (WNT), and bone morphogenetic protein (BMP) signalling interactions capable of spontaneously producing a periodic pattern. Using time-lapse imaging, we find that mesenchymal cell condensation at hair follicles is locally directed by an epidermal prepattern. However, imposing this prepattern's condition of high FGF and low BMP activity across the entire skin reveals a latent dermal capacity to undergo spatially patterned self-organisation in the absence of epithelial direction. This mesenchymal self-organisation relies on restricted transforming growth factor (TGF) β signalling, which serves to drive chemotactic mesenchymal patterning when reaction-diffusion patterning is suppressed, but, in normal conditions, facilitates cell movement to locally prepatterned sources of FGF. This work illustrates a hierarchy of periodic patterning modes operating in organogenesis.
10.1371/journal.pbio.2002117
Multi-axial self-organization properties of mouse embryonic stem cells into gastruloids.
Nature
The emergence of multiple axes is an essential element in the establishment of the mammalian body plan. This process takes place shortly after implantation of the embryo within the uterus and relies on the activity of gene regulatory networks that coordinate transcription in space and time. Whereas genetic approaches have revealed important aspects of these processes, a mechanistic understanding is hampered by the poor experimental accessibility of early post-implantation stages. Here we show that small aggregates of mouse embryonic stem cells (ESCs), when stimulated to undergo gastrulation-like events and elongation in vitro, can organize a post-occipital pattern of neural, mesodermal and endodermal derivatives that mimic embryonic spatial and temporal gene expression. The establishment of the three major body axes in these 'gastruloids' suggests that the mechanisms involved are interdependent. Specifically, gastruloids display the hallmarks of axial gene regulatory systems as exemplified by the implementation of collinear Hox transcriptional patterns along an extending antero-posterior axis. These results reveal an unanticipated self-organizing capacity of aggregated ESCs and suggest that gastruloids could be used as a complementary system to study early developmental events in the mammalian embryo.
10.1038/s41586-018-0578-0
LncRNAs in vertebrates: advances and challenges.
Mallory Allison C,Shkumatava Alena
Biochimie
Beyond the handful of classic and well-characterized long noncoding RNAs (lncRNAs), more recently, hundreds of thousands of lncRNAs have been identified in multiple species including bacteria, plants and vertebrates, and the number of newly annotated lncRNAs continues to increase as more transcriptomes are analyzed. In vertebrates, the expression of many lncRNAs is highly regulated, displaying discrete temporal and spatial expression patterns, suggesting roles in a wide range of developmental processes and setting them apart from classic housekeeping ncRNAs. In addition, the deregulation of a subset of these lncRNAs has been linked to the development of several diseases, including cancers, as well as developmental anomalies. However, the majority of vertebrate lncRNA functions remain enigmatic. As such, a major task at hand is to decipher the biological roles of lncRNAs and uncover the regulatory networks upon which they impinge. This review focuses on our emerging understanding of lncRNAs in vertebrate animals, highlighting some recent advances in their functional analyses across several species and emphasizing the current challenges researchers face to characterize lncRNAs and identify their in vivo functions.
10.1016/j.biochi.2015.03.014
Pharmacogenomics of phenolic antioxidant butylated hydroxyanisole (BHA) in the small intestine and liver of Nrf2 knockout and C57BL/6J mice.
Nair Sujit,Xu Changjiang,Shen Guoxiang,Hebbar Vidya,Gopalakrishnan Avantika,Hu Rong,Jain Mohit Raja,Lin Wen,Keum Young-Sam,Liew Celine,Chan Jefferson Y,Kong Ah-Ng Tony
Pharmaceutical research
PURPOSE:The objective of this study was to investigate the pharmacogenomics and the spatial regulation of global gene expression profiles elicited by cancer chemopreventive agent butylated hydroxyanisole (BHA) in mouse small intestine and liver as well as to identify BHA-modulated nuclear factor-E2-related factor 2 (Nrf2)-dependent genes. METHODS:C57BL/6J (+/+; wildtype) and C57BL/6J/Nrf2(-/-; knockout) mice were administered a single 200 mg/kg oral dose of BHA or only vehicle. Both small intestine and liver were collected at 3 h after treatment and total RNA was extracted. Gene expression profiles were analyzed using 45,000 Affymetrix mouse genome 430 2.0 array and GeneSpring 7.2 software. Microarray results were validated by quantitative real-time reverse transcription-PCR analyses. RESULTS:Clusters of genes that were either induced or suppressed more than two fold by BHA treatment compared with vehicle in C57BL/6J/Nrf2(-/-; knockout) and C57BL/6J Nrf2 (+/+; wildtype) mice genotypes were identified. Amongst these, in small intestine and liver, 1,490 and 493 genes respectively were identified as Nrf2-dependent and upregulated, and 1,090 and 824 genes respectively as Nrf2-dependent and downregulated. Based on their biological functions, these genes can be categorized into ubiquitination/proteolysis, apoptosis/cell cycle, electron transport, detoxification, cell growth/differentiation, transcription factors/interacting partners, kinases and phosphatases, transport, biosynthesis/metabolism, RNA/protein processing and nuclear assembly, and DNA replication genes. Phase II detoxification/antioxidant genes as well as novel molecular target genes, including putative interacting partners of Nrf2 such as nuclear corepressors and coactivators, were also identified as Nrf2-dependent genes. CONCLUSIONS:The identification of BHA-regulated and Nrf2-dependent genes not only provides potential novel insights into the gestalt biological effects of BHA on the pharmacogenomics and spatial regulation of global gene expression profiles in cancer chemoprevention, but also points to the pivotal role of Nrf2 in these biological processes.
10.1007/s11095-006-9099-x
Discriminative identification of transcriptional responses of promoters and enhancers after stimulus.
Kleftogiannis Dimitrios,Kalnis Panos,Arner Erik,Bajic Vladimir B
Nucleic acids research
Promoters and enhancers regulate the initiation of gene expression and maintenance of expression levels in spatial and temporal manner. Recent findings stemming from the Cap Analysis of Gene Expression (CAGE) demonstrate that promoters and enhancers, based on their expression profiles after stimulus, belong to different transcription response subclasses. One of the most promising biological features that might explain the difference in transcriptional response between subclasses is the local chromatin environment. We introduce a novel computational framework, PEDAL, for distinguishing effectively transcriptional profiles of promoters and enhancers using solely histone modification marks, chromatin accessibility and binding sites of transcription factors and co-activators. A case study on data from MCF-7 cell-line reveals that PEDAL can identify successfully the transcription response subclasses of promoters and enhancers from two different stimulations. Moreover, we report subsets of input markers that discriminate with minimized classification error MCF-7 promoter and enhancer transcription response subclasses. Our work provides a general computational approach for identifying effectively cell-specific and stimulation-specific promoter and enhancer transcriptional profiles, and thus, contributes to improve our understanding of transcriptional activation in human.
10.1093/nar/gkw1015
Spatially Resolved Bioenergetic and Genetic Reprogramming Through the Brain of Rats Bearing Implanted C6 Gliomas As Detected by Multinuclear High-Resolution Magic Angle Spinning and Genomic Analysis.
Righi Valeria,García-Martín María-Luisa,Mucci Adele,Schenetti Luisa,Tugnoli Vitaliano,Lopez-Larrubia Pilar,Cerdán Sebastián
Journal of proteome research
We used H, C HRMAS and genomic analysis to investigate regionally the transition from oxidative to glycolytic phenotype and its relationship with altered gene expression in adjacent biopsies through the brain of rats bearing C6 gliomas. Tumor-bearing animals were anesthetized and infused with a solution of [1-C]-glucose, and small adjacent biopsies were obtained spanning transversally from the contralateral hemisphere (regions I and II), the right and left peritumoral areas (regions III and V, respectively), and the tumor core (region IV). These biopsies were analyzed by H, C HRMAS and by quantitative gene expression techniques. Glycolytic metabolism, as reflected by the [3-C]-lactate content, increased clearly from regions I to IV, recovering partially to physiological levels in region V. In contrast, oxidative metabolism, as reflected by the [4-C]-glutamate labeling, decreased in regions I-IV, recovering partially in region V. This metabolic shift from normal to malignant metabolic phenotype paralleled changes in the expression of HIF1α, HIF2α, HIF3α genes, downstream transporters, and regulatory glycolytic, oxidative, and anaplerotic genes in the same regions. Together, our results indicate that genetic and metabolic alterations occurring in the brain of rats bearing C6 gliomas colocalize in situ and the profile of genetic alterations in every region can be inferred from the metabolomic profiles observed in situ by multinuclear HRMAS.
10.1021/acs.jproteome.8b00130
Comparative transcriptomic analysis reveals the regulatory mechanisms of catechins synthesis in different cultivars of Camellia sinensis.
Food research international (Ottawa, Ont.)
Camellia sinensis (L.) O. Kuntze is used to produce tea, a beverage consumed worldwide. Catechins are major medically active components of C. sinensis and can be used clinically to treat hyperglycaemia, hypertension, and cancer. In this study, we aimed to identify the genes involved in catechins biosynthesis. To this end, we analysed transcriptome data from two different cultivars of C. sinensis using DNBSEQ technology. In total,47,717 unigenes were obtained from two cultivars of C. sinensis, of which 9429 were predicted as new unigenes. In our analyses of the Kyoto Encyclopedia of Genes and Genomes database, 212 unigenes encoding 13 key enzymes involved in catechins biosynthesis were identified; the structures of leucoanthocyanidin reductase and anthocyanidin reductase were spatially modelled. Some of these key enzymes were verified by real-time quantitative polymerase chain reaction, and multiple genes encoding plant resistance proteins or transcription factors were identified and analysed. Furthermore, two microRNAs involved in the regulation of catechins biosynthesis were explored. Differentially expressed genes involved in the flavonoid biosynthesis pathway were identified from pairwise comparisons of genes from different cultivars of tea plants. Overall, our findings expanded the number of publicly available transcript datasets for this valuable plant species and identified candidate genes related to the biosynthesis of C. sinensis catechins, thereby establishing a foundation for further in-depth studies of catechins biosynthesis in varieties or cultivars of C. sinensis.
10.1016/j.foodres.2022.111375
Single cell transcriptomics of human epidermis identifies basal stem cell transition states.
Nature communications
How stem cells give rise to epidermis is unclear despite the crucial role the epidermis plays in barrier and appendage formation. Here we use single cell-RNA sequencing to interrogate basal stem cell heterogeneity of human interfollicular epidermis and find four spatially distinct stem cell populations at the top and bottom of rete ridges and transitional positions between the basal and suprabasal epidermal layers. Cell-cell communication modeling suggests that basal cell populations serve as crucial signaling hubs to maintain epidermal communication. Combining pseudotime, RNA velocity, and cellular entropy analyses point to a hierarchical differentiation lineage supporting multi-stem cell interfollicular epidermal homeostasis models and suggest that transitional basal stem cells are stable states essential for proper stratification. Finally, alterations in differentially expressed transitional basal stem cell genes result in severe thinning of human skin equivalents, validating their essential role in epidermal homeostasis and reinforcing the critical nature of basal stem cell heterogeneity.
10.1038/s41467-020-18075-7
Transcriptional profiling of macrophages in situ in metastatic melanoma reveals localization-dependent phenotypes and function.
Cell reports. Medicine
Modulation of immune function at the tumor site could improve patient outcomes. Here, we analyze patient samples of metastatic melanoma, a tumor responsive to T cell-based therapies, and find that tumor-infiltrating T cells are primarily juxtaposed to CD14 monocytes/macrophages rather than melanoma cells. Using immunofluorescence-guided laser capture microdissection, we analyze transcriptomes of CD3 T cells, CD14 monocytes/macrophages, and melanoma cells in non-dissociated tissue. Stromal CD14 cells display a specific transcriptional signature distinct from CD14 cells within tumor nests. This signature contains LY75, a gene linked with antigen capture and regulation of tolerance and immunity in dendritic cells (DCs). When applied to TCGA cohorts, this gene set can distinguish patients with significantly prolonged survival in metastatic cutaneous melanoma and other cancers. Thus, the stromal CD14 cell signature represents a candidate biomarker and suggests that reprogramming of stromal macrophages to acquire DC function may offer a therapeutic opportunity for metastatic cancers.
10.1016/j.xcrm.2022.100621
Dynamic profiling of immune microenvironment during anti-PD-1 immunotherapy for head and neck squamous cell carcinoma: the IPRICE study.
BMC cancer
BACKGROUND:Immune checkpoint inhibitors of programmed cell death protein 1 (PD-1) represent a significant breakthrough in treating head and neck squamous cell carcinoma (HNSCC), with long-lasting responses and prolonged survival observed in first- and second-line therapy. However, this is observed in < 20% of patients and high primary/secondary resistance may occur. The primary objective of the identification of predictive factors for the response to anti-PD-1 immunotherapy in head and neck squamous cell carcinoma (IPRICE) study is to identify predictive factors of response to anti-PD-1 immunotherapy. METHODS:The IPRICE study is a single-center, prospective, non-randomized, open-label, and interventional clinical trial. Liquid and tumor biopsies will be performed in 54 patients with recurrent/metastatic (R/M) HNSCC undergoing anti-PD-1 immunotherapy alone to compare the evolution of gene expression and immunological profile between responders and non-responders. We will use a multidisciplinary approach including spatial transcriptomics, single seq-RNA analysis, clinical data, and medical images. Genes, pathways, and transcription factors potentially involved in the immune response will also be analyzed, including genes involved in the interferon-gamma (IFN-γ) pathway, immunogenic cell death and mitophagy, hypoxia, circulating miRNA-mediated immunomodulation, cytokines, and immune repertoire within the tumor microenvironment (TME). With a follow-up period of 3-years, these data will help generate effective biomarkers to define optimal therapeutic strategy and new immunomodulatory agents based on a better understanding of primary/secondary resistance mechanisms. Tumor biopsy will be performed initially before the start of immunotherapy at the first tumor assessment and is only proposed at tumor progression. Clinical data will be collected using a dedicated Case Report Form (CRF). DISCUSSION:Identifying predictive factors of the response to anti-PD-1 immunotherapy and optimizing long-term immune response require a thorough understanding of the intrinsic and acquired resistance to immunotherapy. To achieve this, dynamic profiling of TME during anti-PD-1 immunotherapy based on analysis of tumor biopsy samples is critical. This will be accomplished through the anatomical localization of HNSCC, which will allow for the analysis of multiple biopsies during treatment and the emergence of breakthrough technologies including single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics. TRIAL REGISTRATION:Clinicaltrial.gov. Registered April 14, 2022, https://www. CLINICALTRIALS:gov/study/NCT05328024 .
10.1186/s12885-023-11672-x
ETS-family genes in pancreatic development.
Kobberup Sune,Nyeng Pia,Juhl Kirstine,Hutton John,Jensen Jan
Developmental dynamics : an official publication of the American Association of Anatomists
ETS-family factors play major roles in development and cancer, notably as critical targets for extra-cellular signaling pathways, including MAPK-signaling. Given the presently limited knowledge on the role of ETS-factors in pancreatic development, we here sought to characterize all 26 individual members of the ETS-family in relation to pancreatic development using a combination of genomics, RT-PCR, and histological techniques. This analysis uncovers 22 ETS family genes displaying select spatial and temporal expression patterns in the developing pancreas. Highly specific expression of ETS-family components is observed in pancreatic progenitor cells or the associated embryonic mesenchyme. Other members are linked to the differentiation of more mature pancreatic cells, including exocrine and endocrine cell types. We find that two members of the Etv subfamily, Etv4 and Etv5, are expressed in cells proximal to pancreatic mesenchyme, and, furthermore, induced in FGF10-arrested pancreatic progenitors suggesting that these factors mediate mesenchymal-to-epithelial signaling.
10.1002/dvdy.21292
Stroma AReactive Invasion Front Areas (SARIFA) proves prognostic relevance in gastric carcinoma and is based on a tumor-adipocyte interaction indicating an altered immune response.
Gastric cancer : official journal of the International Gastric Cancer Association and the Japanese Gastric Cancer Association
BACKGROUND:Recently, we presented Stroma AReactive Invasion Front Areas (SARIFA) as a new histomorphologic negative prognostic biomarker in gastric cancer. It is defined as direct contact between tumor cells and fat cells. The aim of this study was to further elucidate the underlying genomic, transcriptional, and immunological mechanisms of the SARIFA phenomenon. METHODS:To address these questions, SARIFA was classified on H&E-stained tissue sections of three cohorts: an external cohort (n = 489, prognostic validation), the TCGA-STAD cohort (n = 194, genomic and transcriptomic analysis), and a local cohort (n = 60, digital spatial profiling (whole transcriptome) and double RNA in situ hybridization/immunostaining of cytokines). RESULTS:SARIFA status proved to be an independent negative prognostic factor for overall survival in an external cohort of gastric carcinomas. In TCGA-STAD cohort, SARIFA is not driven by distinct genomic alterations, whereas the gene expression analyses showed an upregulation of FABP4 in SARIFA-positive tumors. In addition, the transcriptional regulations of white adipocyte differentiation, triglyceride metabolism, and catabolism were upregulated in pathway analyses. In the DSP analysis of SARIFA-positive tumors, FABP4 and the transcriptional regulation of white adipocyte differentiation were upregulated in macrophages. Additionally, a significantly lower expression of the cytokines IL6 and TNFα was observed at the invasion front. CONCLUSIONS:SARIFA proves to be a strong negative prognostic biomarker in advanced gastric cancer, implicating an interaction of tumor cells with tumor-promoting adipocytes with crucial changes in tumor cell metabolism. SARIFA is not driven by tumor genetics but is very likely driven by an altered immune response as a causative mechanism.
10.1007/s10120-023-01436-8
Distal Alternative Last Exons Localize mRNAs to Neural Projections.
Taliaferro J Matthew,Vidaki Marina,Oliveira Ruan,Olson Sara,Zhan Lijun,Saxena Tanvi,Wang Eric T,Graveley Brenton R,Gertler Frank B,Swanson Maurice S,Burge Christopher B
Molecular cell
Spatial restriction of mRNA to distinct subcellular locations enables local regulation and synthesis of proteins. However, the organizing principles of mRNA localization remain poorly understood. Here we analyzed subcellular transcriptomes of neural projections and soma of primary mouse cortical neurons and two neuronal cell lines and found that alternative last exons (ALEs) often confer isoform-specific localization. Surprisingly, gene-distal ALE isoforms were four times more often localized to neurites than gene-proximal isoforms. Localized isoforms were induced during neuronal differentiation and enriched for motifs associated with muscleblind-like (Mbnl) family RNA-binding proteins. Depletion of Mbnl1 and/or Mbnl2 reduced localization of hundreds of transcripts, implicating Mbnls in localization of mRNAs to neurites. We provide evidence supporting a model in which the linkage between genomic position of ALEs and subcellular localization enables coordinated induction of localization-competent mRNA isoforms through a post-transcriptional regulatory program that is induced during differentiation and reversed in cellular reprogramming and cancer.
10.1016/j.molcel.2016.01.020
Spatially distinct neutrophil responses within the inflammatory lesions of pneumonic plague.
Stasulli Nikolas M,Eichelberger Kara R,Price Paul A,Pechous Roger D,Montgomery Stephanie A,Parker Joel S,Goldman William E
mBio
UNLABELLED:During pneumonic plague, the bacterium Yersinia pestis elicits the development of inflammatory lung lesions that continue to expand throughout infection. This lesion development and persistence are poorly understood. Here, we examine spatially distinct regions of lung lesions using laser capture microdissection and transcriptome sequencing (RNA-seq) analysis to identify transcriptional differences between lesion microenvironments. We show that cellular pathways involved in leukocyte migration and apoptosis are downregulated in the center of lung lesions compared to the periphery. Probing for the bacterial factor(s) important for the alteration in neutrophil survival, we show both in vitro and in vivo that Y. pestis increases neutrophil survival in a manner that is dependent on the type III secretion system effector YopM. This research explores the complexity of spatially distinct host-microbe interactions and emphasizes the importance of cell relevance in assays in order to fully understand Y. pestis virulence. IMPORTANCE:Yersinia pestis is a high-priority pathogen and continues to cause outbreaks worldwide. The ability of Y. pestis to be transmitted via respiratory droplets and its history of weaponization has led to its classification as a select agent most likely to be used as a biological weapon. Unrestricted bacterial growth during the initial preinflammatory phase primes patients to be infectious once disease symptoms begin in the proinflammatory phase, and the rapid disease progression can lead to death before Y. pestis infection can be diagnosed and treated. Using in vivo analyses and focusing on relevant cell types during pneumonic plague infection, we can identify host pathways that may be manipulated to extend the treatment window for pneumonic plague patients.
10.1128/mBio.01530-15
Increased intron retention is a post-transcriptional signature associated with progressive aging and Alzheimer's disease.
Adusumalli Swarnaseetha,Ngian Zhen-Kai,Lin Wei-Qi,Benoukraf Touati,Ong Chin-Tong
Aging cell
Intron retention (IR) by alternative splicing is a conserved regulatory mechanism that can affect gene expression and protein function during adult development and age-onset diseases. However, it remains unclear whether IR undergoes spatial or temporal changes during different stages of aging or neurodegeneration like Alzheimer's disease (AD). By profiling the transcriptome of Drosophila head cells at different ages, we observed a significant increase in IR events for many genes during aging. Differential IR affects distinct biological functions at different ages and occurs at several AD-associated genes in older adults. The increased nucleosome occupancy at the differentially retained introns in young animals suggests that it may regulate the level of IR during aging. Notably, an increase in the number of IR events was also observed in healthy older mouse and human brain tissues, as well as in the cerebellum and frontal cortex from independent AD cohorts. Genes with differential IR shared many common features, including shorter intron length, no perturbation in their mRNA level, and enrichment for biological functions that are associated with mRNA processing and proteostasis. The differentially retained introns identified in AD frontal cortex have higher GC content, with many of their mRNA transcripts showing an altered level of protein expression compared to control samples. Taken together, our results suggest that an increased IR is an conserved signature that is associated with aging. By affecting pathways involved in mRNA and protein homeostasis, changes of IR pattern during aging may regulate the transition from healthy to pathological state in late-onset sporadic AD.
10.1111/acel.12928
Spatial regulation of gene expression during growth of articular cartilage in juvenile mice.
Lui Julian C,Chau Michael,Chen Weiping,Cheung Crystal S F,Hanson Jeffrey,Rodriguez-Canales Jaime,Nilsson Ola,Baron Jeffrey
Pediatric research
BACKGROUND:In juvenile mammals, the epiphyses of long bones grow by chondrogenesis within the articular cartilage. A better understanding of the molecular mechanisms that regulate the growth of articular cartilage may give insight into the antecedents of joint disease, such as osteoarthritis. METHODS:We used laser capture microdissection to isolate chondrocytes from the superficial, middle, and deep zones of growing tibial articular cartilage in the 1-wk-old mouse and then investigated expression patterns by microarray. To identify molecular markers for each zone of the growing articular cartilage, we found genes showing zone-specific expression and confirmed by real-time PCR and in situ hybridization. RESULTS:Bioinformatic analyses implicated ephrin receptor signaling, Wnt signaling, and bone morphogenetic protein signaling in the spatial regulation of chondrocyte differentiation during growth. Molecular markers were identified for superficial (e.g., Cilp, Prg4), middle (Cxcl14, Tnn), and deep zones (Sfrp5, Frzb). Comparison between juvenile articular and growth plate cartilage revealed that the superficial-to-deep zone transition showed similarity with the hypertrophic-to-resting zone transition. CONCLUSION:Laser capture microdissection combined with microarray analysis identified novel signaling pathways that are spatially regulated in growing mouse articular cartilage and revealed similarities between the molecular architecture of the growing articular cartilage and that of growth plate cartilage.
10.1038/pr.2014.208
TBR2 coordinates neurogenesis expansion and precise microcircuit organization via Protocadherin 19 in the mammalian cortex.
Lv Xiaohui,Ren Si-Qiang,Zhang Xin-Jun,Shen Zhongfu,Ghosh Tanay,Xianyu Anjin,Gao Peng,Li Zhizhong,Lin Susan,Yu Yang,Zhang Qiangqiang,Groszer Matthias,Shi Song-Hai
Nature communications
Cerebral cortex expansion is a hallmark of mammalian brain evolution; yet, how increased neurogenesis is coordinated with structural and functional development remains largely unclear. The T-box protein TBR2/EOMES is preferentially enriched in intermediate progenitors and supports cortical neurogenesis expansion. Here we show that TBR2 regulates fine-scale spatial and circuit organization of excitatory neurons in addition to enhancing neurogenesis in the mouse cortex. TBR2 removal leads to a significant reduction in neuronal, but not glial, output of individual radial glial progenitors as revealed by mosaic analysis with double markers. Moreover, in the absence of TBR2, clonally related excitatory neurons become more laterally dispersed and their preferential synapse development is impaired. Interestingly, TBR2 directly regulates the expression of Protocadherin 19 (PCDH19), and simultaneous PCDH19 expression rescues neurogenesis and neuronal organization defects caused by TBR2 removal. Together, these results suggest that TBR2 coordinates neurogenesis expansion and precise microcircuit assembly via PCDH19 in the mammalian cortex.
10.1038/s41467-019-11854-x
Single-cell RNA-seq reveals actinic keratosis-specific keratinocyte subgroups and their crosstalk with secretory-papillary fibroblasts.
Journal of the European Academy of Dermatology and Venereology : JEADV
BACKGROUND AND AIM:Actinic keratosis (AK) represents an intraepidermal malignant neoplasm with the proliferation of atypical keratinocytes. AK lesions are regarded as early in situ squamous cell carcinomas (SCCs) having the potential to progress into invasive SCC (iSCC) and metastasize, causing death. This study aimed to investigate the heterogeneity of keratinocytes and how this heterogeneity promoted AK development and progression. METHODS:We employed single-cell RNA sequencing (scRNA-seq) to examine the heterogeneity of keratinocytes and dermal fibroblast clusters in AKs and adjacent normal skins. Cell clustering, pseudotime trajectory construction, gene ontology enrichment analysis, transcription factor network analysis, and cell-cell communication were used to investigate the heterogeneity of keratinocytes in AK. The cellular identity and function were verified by immunohistochemical and immunofluorescence staining. RESULTS:Using scRNA-seq, we revealed 13 keratinocyte subgroups (clusters 0-12) in AK tissues and characterized 2 AK-specific clusters. Cluster 9 displayed high levels of IL1R2 and WFDC2, and cluster 11 showed high levels of FADS2 and FASN. The percentages of cells in these two clusters significantly increased in AK compared with normal tissues. The existence and spatial localization of AK-specific IL1R2+WFDC2+ cluster were verified by immunohistochemical and immunofluorescence staining. Functional studies indicated that the genes identified in the IL1R2+WFDC2+ cluster were crucial for epithelial cell proliferation, migration, and angiogenesis. Further immunofluorescent staining revealed the interactions between AK-specific keratinocytes and secretory-papillary fibroblasts mainly through ANGPTL4-ITGA5 signalling pathway rarely seen in normal tissues. CONCLUSION:The findings of this study might help better understand AK pathogenesis.
10.1111/jdv.19289
Pan-cancer computational histopathology reveals mutations, tumor composition and prognosis.
Nature cancer
We use deep transfer learning to quantify histopathological patterns across 17,355 hematoxylin and eosin-stained histopathology slide images from 28 cancer types and correlate these with matched genomic, transcriptomic and survival data. This approach accurately classifies cancer types and provides spatially resolved tumor and normal tissue distinction. Automatically learned computational histopathological features correlate with a large range of recurrent genetic aberrations across cancer types. This includes whole-genome duplications, which display universal features across cancer types, individual chromosomal aneuploidies, focal amplifications and deletions, as well as driver gene mutations. There are widespread associations between bulk gene expression levels and histopathology, which reflect tumor composition and enable the localization of transcriptomically defined tumor-infiltrating lymphocytes. Computational histopathology augments prognosis based on histopathological subtyping and grading, and highlights prognostically relevant areas such as necrosis or lymphocytic aggregates. These findings show the remarkable potential of computer vision in characterizing the molecular basis of tumor histopathology.
10.1038/s43018-020-0085-8
Molecular characterization of onychomatricoma: Spatial profiling reveals the role of onychofibroblasts in its pathogenesis.
Experimental dermatology
Onychomatricoma (OM) is a rare nail unit tumour with a characteristic presentation of finger-like projections arising from the nail matrix. Due to the lack of transcriptome information, the mechanisms underlying its development are largely unknown. To characterize molecular features involved in the disease pathogenesis, we used digital spatial profiling (DSP) in 2 cases of OM and normal control nail units. Based on the histological evaluation, we selectively profiled 69 regions of interest covering epithelial and stromal compartments of each tissue section. Dermoscopic and histopathologic findings were reviewed in 6 cases. Single-cell RNA sequencing of nail units and DSP were combined to define cell type contributions of OM. We identified 173 genes upregulated in stromal compartments of OM compared to onychodermis, specialized nail mesenchyme. Gene ontology analysis of the upregulated genes suggested the role of Wnt pathway activation in OM pathogenesis. We also found PLA2G2A, a known modulator of Wnt signalling, is strongly and specifically expressed in the OM stroma. The potential role of Wnt pathway was further supported by strong nuclear localization of β-catenin in OM. Compared to the nail matrix epithelium, only a few genes were increased in OM epithelium. Deconvolution of nail unit cell types showed that onychofibroblasts are the dominant cell type in OM stroma. Altogether, integrated spatial and single-cell multi-omics concluded that OM is a tumour that derives a significant proportion of its origin from onychofibroblasts and is associated with upregulation of Wnt signals, which play a key role in the disease pathogenesis.
10.1111/exd.14736
The molecular genetic landscape of human brain size variation.
Cell reports
Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.
10.1016/j.celrep.2023.113439
Radiomics Response Signature for Identification of Metastatic Colorectal Cancer Sensitive to Therapies Targeting EGFR Pathway.
Dercle Laurent,Lu Lin,Schwartz Lawrence H,Qian Min,Tejpar Sabine,Eggleton Peter,Zhao Binsheng,Piessevaux Hubert
Journal of the National Cancer Institute
BACKGROUND:The authors sought to forecast survival and enhance treatment decisions for patients with liver metastatic colorectal cancer by using on-treatment radiomics signature to predict tumor sensitiveness to irinotecan, 5-fluorouracil, and leucovorin (FOLFIRI) alone (F) or in combination with cetuximab (FC). METHODS:We retrospectively analyzed 667 metastatic colorectal cancer patients treated with F or FC. Computed tomography quality was classified as high (HQ) or standard (SD). Four datasets were created using the nomenclature (treatment) - (quality). Patients were randomly assigned (2:1) to training or validation sets: FCHQ: 78:38, FCSD: 124:62, FHQ: 78:51, FSD: 158:78. Four tumor-imaging biomarkers measured quantitative radiomics changes between standard of care computed tomography scans at baseline and 8 weeks. Using machine learning, the performance of the signature to classify tumors as treatment sensitive or treatment insensitive was trained and validated using receiver operating characteristic (ROC) curves. Hazard ratio and Cox regression models evaluated association with overall survival (OS). RESULTS:The signature (area under the ROC curve [95% confidence interval (CI)]) used temporal decrease in tumor spatial heterogeneity plus boundary infiltration to successfully predict sensitivity to antiepidermal growth factor receptor therapy (FCHQ: 0.80 [95% CI = 0.69 to 0.94], FCSD: 0.72 [95% CI = 0.59 to 0.83]) but failed with chemotherapy (FHQ: 0.59 [95% CI = 0.44 to 0.72], FSD: 0.55 [95% CI = 0.43 to 0.66]). In cetuximab-containing sets, radiomics signature outperformed existing biomarkers (KRAS-mutational status, and tumor shrinkage by RECIST 1.1) for detection of treatment sensitivity and was strongly associated with OS (two-sided P < .005). CONCLUSIONS:Radiomics response signature can serve as an intermediate surrogate marker of OS. The signature outperformed known biomarkers in providing an early prediction of treatment sensitivity and could be used to guide cetuximab treatment continuation decisions.
10.1093/jnci/djaa017
Paracrine WNT5A Signaling Inhibits Expansion of Tumor-Initiating Cells.
Borcherding Nicholas,Kusner David,Kolb Ryan,Xie Qing,Li Wei,Yuan Fang,Velez Gabriel,Askeland Ryan,Weigel Ronald J,Zhang Weizhou
Cancer research
It is not well understood how paracrine communication between basal and luminal cell populations in the mammary gland affects tumorigenesis. During ErbB2-induced mammary tumorigenesis, enriched mammary stem cells that represent a subpopulation of basal cells exhibit enhanced tumorigenic capacity compared with the corresponding luminal progenitors. Transcript profiling of tumors derived from basal and luminal tumor-initiating cells (TIC) revealed preferential loss of the noncanonical Wnt ligand WNT5A in basal TIC-derived tumors. Heterozygous loss of WNT5A was correlated with shorter survival of breast cancer patients. In a mouse model of ErbB2-induced breast cancer, Wnt5a heterozygosity promoted tumor multiplicity and pulmonary metastasis. As a TGFβ substrate, luminal cell-produced WNT5A induced a feed-forward loop to activate SMAD2 in a RYK and TGFβR1-dependent manner to limit the expansion of basal TIC in a paracrine fashion, a potential explanation for the suppressive effect of WNT5A in mammary tumorigenesis. Our results identify the WNT5A/RYK module as a spatial regulator of the TGFβ-SMAD signaling pathway in the context of mammary gland development and carcinogenesis, offering a new perspective on tumor suppression provided by basal-luminal cross-talk in normal mammary tissue.
10.1158/0008-5472.CAN-14-2761
Tumor-on-a-chip: a microfluidic model to study cell response to environmental gradients.
Lab on a chip
Limited blood supply and rapid tumor metabolism within solid tumors leads to nutrient starvation, waste product accumulation and the generation of pH gradients across the tumor mass. These environmental conditions modify multiple cellular functions, including metabolism, proliferation, and drug response. However, capturing the spatial metabolic and phenotypic heterogeneity of the tumor with classic in vitro models remains challenging. Thus, in this work a microfluidic tumor slice model was developed to study cell behavior under metabolic starvation gradients. The presented microdevice comprises a central chamber where tumor cells were cultured in a 3D collagen hydrogel. A lumen on the flank of the chamber was used to perfuse media, mimicking the vasculature. Under these circumstances, tumor cell metabolism led to the generation of viability, proliferation and pH gradients. The model decoupled the influence of oxygen from other nutrients, revealing that cell necrosis at the core of the model could be explained by nutrient starvation. The microdevice can be disassembled to retrieve the cells from the desired locations to study molecular adaptions due to nutrient starvation. When exposed to these pH gradients and low nutrient conditions, cancer cells showed multiple changes in their gene expression profile depending on their distance from the lumen. Those cells located further from the lumen upregulated several genes related to stress and survival response, whereas genes related to proliferation and DNA repair were downregulated. This model may help to identify new therapeutic opportunities to target the metabolic heterogeneity observed in solid tumors.
10.1039/c9lc00270g
Simultaneous Hypoxia and Low Extracellular pH Suppress Overall Metabolic Rate and Protein Synthesis In Vitro.
Sørensen Brita Singers,Busk Morten,Overgaard Jens,Horsman Michael R,Alsner Jan
PloS one
BACKGROUND:The tumor microenvironment is characterized by regions of hypoxia and acidosis which are linked to poor prognosis. This occurs due to an aberrant vasculature as well as high rates of glycolysis and lactate production in tumor cells even in the presence of oxygen (the Warburg effect), which weakens the spatial linkage between hypoxia and acidosis. METHODS:Five different human squamous cell carcinoma cell lines (SiHa, FaDuDD, UTSCC5, UTSCC14 and UTSCC15) were treated with hypoxia, acidosis (pH 6.3), or a combination, and gene expression analyzed using microarray. SiHa and FaDuDD were chosen for further characterization of cell energetics and protein synthesis. Total cellular ATP turnover and relative glycolytic dependency was determined by simultaneous measurements of oxygen consumption and lactate synthesis rates and total protein synthesis was determined by autoradiographic quantification of the incorporation of 35S-labelled methionine and cysteine into protein. RESULTS:Microarray analysis allowed differentiation between genes induced at low oxygen only at normal extracellular pH (pHe), genes induced at low oxygen at both normal and low pHe, and genes induced at low pHe independent of oxygen concentration. Several genes were found to be upregulated by acidosis independent of oxygenation. Acidosis resulted in a more wide-scale change in gene expression profiles than hypoxia including upregulation of genes involved in the translation process, for example Eukaryotic translation initiation factor 4A, isoform 2 (EIF4A2), and Ribosomal protein L37 (RPL37). Acidosis suppressed overall ATP turnover and protein synthesis by 50%. Protein synthesis, but not total ATP production, was also suppressed under hypoxic conditions. A dramatic decrease in ATP turnover (SiHa) and protein synthesis (both cell lines) was observed when hypoxia and low pHe were combined. CONCLUSIONS:We demonstrate here that the influence of hypoxia and acidosis causes different responses, both in gene expression and in de novo protein synthesis, depending on whether the two factors induced alone or overlapping, and as such it is important for in vivo studies to take this into account.
10.1371/journal.pone.0134955
Metabolic markers in relation to hypoxia; staining patterns and colocalization of pimonidazole, HIF-1α, CAIX, LDH-5, GLUT-1, MCT1 and MCT4.
Rademakers Saskia E,Lok Jasper,van der Kogel Albert J,Bussink Johan,Kaanders Johannes H A M
BMC cancer
BACKGROUND:The cellular response of malignant tumors to hypoxia is diverse. Several important endogenous metabolic markers are upregulated under hypoxic conditions. We examined the staining patterns and co-expression of HIF-1α, CAIX, LDH-5, GLUT-1, MCT1 and MCT4 with the exogenous hypoxic cell marker pimonidazole and the association of marker expression with clinicopathological characteristics. METHODS:20 biopsies of advanced head and neck carcinomas were immunohistochemically stained and analyzed. All patients were given the hypoxia marker pimonidazole intravenously 2 h prior to biopsy taking. The tumor area positive for each marker, the colocalization of the different markers and the distribution of the markers in relation to the blood vessels were assessed by semiautomatic quantitative analysis. RESULTS:MCT1 staining was present in hypoxic (pimonidazole stained) as well as non-hypoxic areas in almost equal amounts. MCT1 expression showed a significant overall correlation (r = 0.75, p < 0.001) and strong spatial relationship with CAIX. LDH-5 showed the strongest correlation with pimonidazole (r = 0.66, p = 0.002). MCT4 and GLUT-1 demonstrated a typical diffusion-limited hypoxic pattern and showed a high degree of colocalization. Both MCT4 and CAIX showed a higher expression in the primary tumor in node positive patients (p = 0.09 both). CONCLUSIONS:Colocalization and staining patterns of metabolic and hypoxia-related proteins provides valuable additional information over single protein analyses and can improve the understanding of their functions and environmental influences.
10.1186/1471-2407-11-167
A novel technique to combine and analyse spatial and temporal expression datasets: A case study with the sea anemone Nematostella vectensis to identify potential gene interactions.
Abdol Amir M,Röttinger Eric,Jansson Fredrik,Kaandorp Jaap A
Developmental biology
Understanding genetic interactions during early development of a given organism, is the first step toward unveiling gene regulatory networks (GRNs) that govern a biological process of interest. Predicting such interactions from large expression datasets by performing targeted knock-down/knock-out approaches is a challenging task. We use the currently available expression datasets (in situ hybridization images & qPCR time series) for a basal anthozoan the sea anemone N. vectensis to construct continuous spatiotemporal gene expression patterns during its early development. Moreover, by combining cluster results from each dataset we develop a method that provides testable hypotheses about potential genetic interactions. We show that the analysis of spatial gene expression patterns reveals functional regions of the embryo during the gastrulation. The clustering results from qPCR time series unveils significant temporal events and highlights genes potentially involved in N. vectensis gastrulation. Furthermore, we introduce a method for merging the clustering results from spatial and temporal datasets by which we can group genes that are expressed in the same region and at the time. We demonstrate that the merged clusters can be used to identify GRN interactions involved in various processes and to predict possible activators or repressors of any gene in the dataset. Finally, we validate our methods and results by predicting the repressor effect of NvErg on NvBra in the central domain during the gastrulation that has recently been confirmed by functional analysis.
10.1016/j.ydbio.2017.06.004
Single-nucleus RNA-seq2 reveals functional crosstalk between liver zonation and ploidy.
Nature communications
Single-cell RNA-seq reveals the role of pathogenic cell populations in development and progression of chronic diseases. In order to expand our knowledge on cellular heterogeneity, we have developed a single-nucleus RNA-seq2 method tailored for the comprehensive analysis of the nuclear transcriptome from frozen tissues, allowing the dissection of all cell types present in the liver, regardless of cell size or cellular fragility. We use this approach to characterize the transcriptional profile of individual hepatocytes with different levels of ploidy, and have discovered that ploidy states are associated with different metabolic potential, and gene expression in tetraploid mononucleated hepatocytes is conditioned by their position within the hepatic lobule. Our work reveals a remarkable crosstalk between gene dosage and spatial distribution of hepatocytes.
10.1038/s41467-021-24543-5
Single-cell Spatial Proteomic Revelations on the Multiparametric MRI Heterogeneity of Clinically Significant Prostate Cancer.
Pachynski Russell K,Kim Eric H,Miheecheva Natalia,Kotlov Nikita,Ramachandran Akshaya,Postovalova Ekaterina,Galkin Ilia,Svekolkin Viktor,Lyu Yang,Zou Qiong,Cao Dengfeng,Gaut Joseph,Ippolito Joseph E,Bagaev Alexander,Bruttan Maria,Gancharova Olga,Nomie Krystle,Tsiper Maria,Andriole Gerald L,Ataullakhanov Ravshan,Hsieh James J
Clinical cancer research : an official journal of the American Association for Cancer Research
PURPOSE:Multiparametric MRI (mpMRI) has become an indispensable radiographic tool in diagnosing prostate cancer. However, mpMRI fails to visualize approximately 15% of clinically significant prostate cancer (csPCa). The molecular, cellular, and spatial underpinnings of such radiographic heterogeneity in csPCa are unclear. EXPERIMENTAL DESIGN:We examined tumor tissues from clinically matched patients with mpMRI-invisible and mpMRI-visible csPCa who underwent radical prostatectomy. Multiplex immunofluorescence single-cell spatial imaging and gene expression profiling were performed. Artificial intelligence-based analytic algorithms were developed to examine the tumor ecosystem and integrate with corresponding transcriptomics. RESULTS:More complex and compact epithelial tumor architectures were found in mpMRI-visible than in mpMRI-invisible prostate cancer tumors. In contrast, similar stromal patterns were detected between mpMRI-invisible prostate cancer and normal prostate tissues. Furthermore, quantification of immune cell composition and tumor-immune interactions demonstrated a lack of immune cell infiltration in the malignant but not in the adjacent nonmalignant tissue compartments, irrespective of mpMRI visibility. No significant difference in immune profiles was detected between mpMRI-visible and mpMRI-invisible prostate cancer within our patient cohort, whereas expression profiling identified a 24-gene stromal signature enriched in mpMRI-invisible prostate cancer. Prostate cancer with strong stromal signature exhibited a favorable survival outcome within The Cancer Genome Atlas prostate cancer cohort. Notably, five recurrences in the 8 mpMRI-visible patients with csPCa and no recurrence in the 8 clinically matched patients with mpMRI-invisible csPCa occurred during the 5-year follow-up post-prostatectomy. CONCLUSIONS:Our study identified distinct molecular, cellular, and structural characteristics associated with mpMRI-visible csPCa, whereas mpMRI-invisible tumors were similar to normal prostate tissue, likely contributing to mpMRI invisibility.
10.1158/1078-0432.CCR-20-4217
Transcriptomics and RNA-Based Therapeutics as Potential Approaches to Manage SARS-CoV-2 Infection.
International journal of molecular sciences
SARS-CoV-2 is a coronavirus family member that appeared in China in December 2019 and caused the disease called COVID-19, which was declared a pandemic in 2020 by the World Health Organization. In recent months, great efforts have been made in the field of basic and clinical research to understand the biology and infection processes of SARS-CoV-2. In particular, transcriptome analysis has contributed to generating new knowledge of the viral sequences and intracellular signaling pathways that regulate the infection and pathogenesis of SARS-CoV-2, generating new information about its biology. Furthermore, transcriptomics approaches including spatial transcriptomics, single-cell transcriptomics and direct RNA sequencing have been used for clinical applications in monitoring, detection, diagnosis, and treatment to generate new clinical predictive models for SARS-CoV-2. Consequently, RNA-based therapeutics and their relationship with SARS-CoV-2 have emerged as promising strategies to battle the SARS-CoV-2 pandemic with the assistance of novel approaches such as CRISPR-CAS, ASOs, and siRNA systems. Lastly, we discuss the importance of precision public health in the management of patients infected with SARS-CoV-2 and establish that the fusion of transcriptomics, RNA-based therapeutics, and precision public health will allow a linkage for developing health systems that facilitate the acquisition of relevant clinical strategies for rapid decision making to assist in the management and treatment of the SARS-CoV-2-infected population to combat this global public health problem.
10.3390/ijms231911058
Reactive astrocytes potentiate tumor aggressiveness in a murine glioma resection and recurrence model.
Okolie Onyinyechukwu,Bago Juli R,Schmid Ralf S,Irvin David M,Bash Ryan E,Miller C Ryan,Hingtgen Shawn D
Neuro-oncology
BACKGROUND:Surgical resection is a universal component of glioma therapy. Little is known about the postoperative microenvironment due to limited preclinical models. Thus, we sought to develop a glioma resection and recurrence model in syngeneic immune-competent mice to understand how surgical resection influences tumor biology and the local microenvironment. METHODS:We genetically engineered cells from a murine glioma mouse model to express fluorescent and bioluminescent reporters. Established allografts were resected using image-guided microsurgery. Postoperative tumor recurrence was monitored by serial imaging, and the peritumoral microenvironment was characterized by histopathology and immunohistochemistry. Coculture techniques were used to explore how astrocyte injury influences tumor aggressiveness in vitro. Transcriptome and secretome alterations in injured astrocytes was examined by RNA-seq and Luminex. RESULTS:We found that image-guided resection achieved >90% reduction in tumor volume but failed to prevent both local and distant tumor recurrence. Immunostaining for glial fibrillary acidic protein and nestin showed that resection-induced injury led to temporal and spatial alterations in reactive astrocytes within the peritumoral microenvironment. In vitro, we found that astrocyte injury induced transcriptome and secretome alterations and promoted tumor proliferation, as well as migration. CONCLUSIONS:This study demonstrates a unique syngeneic model of glioma resection and recurrence in immune-competent mice. Furthermore, this model provided insights into the pattern of postsurgical tumor recurrence and changes in the peritumoral microenvironment, as well as the impact of injured astrocytes on glioma growth and invasion. A better understanding of the postsurgical tumor microenvironment will allow development of targeted anticancer agents that improve surgery-mediated effects on tumor biology.
10.1093/neuonc/now117
Spatiotemporal dissection of the cell cycle with single-cell proteogenomics.
Nature
The cell cycle, over which cells grow and divide, is a fundamental process of life. Its dysregulation has devastating consequences, including cancer. The cell cycle is driven by precise regulation of proteins in time and space, which creates variability between individual proliferating cells. To our knowledge, no systematic investigations of such cell-to-cell proteomic variability exist. Here we present a comprehensive, spatiotemporal map of human proteomic heterogeneity by integrating proteomics at subcellular resolution with single-cell transcriptomics and precise temporal measurements of individual cells in the cell cycle. We show that around one-fifth of the human proteome displays cell-to-cell variability, identify hundreds of proteins with previously unknown associations with mitosis and the cell cycle, and provide evidence that several of these proteins have oncogenic functions. Our results show that cell cycle progression explains less than half of all cell-to-cell variability, and that most cycling proteins are regulated post-translationally, rather than by transcriptomic cycling. These proteins are disproportionately phosphorylated by kinases that regulate cell fate, whereas non-cycling proteins that vary between cells are more likely to be modified by kinases that regulate metabolism. This spatially resolved proteomic map of the cell cycle is integrated into the Human Protein Atlas and will serve as a resource for accelerating molecular studies of the human cell cycle and cell proliferation.
10.1038/s41586-021-03232-9
Measuring cell identity in noisy biological systems.
Birnbaum Kenneth D,Kussell Edo
Nucleic acids research
Global gene expression measurements are increasingly obtained as a function of cell type, spatial position within a tissue and other biologically meaningful coordinates. Such data should enable quantitative analysis of the cell-type specificity of gene expression, but such analyses can often be confounded by the presence of noise. We introduce a specificity measure Spec that quantifies the information in a gene's complete expression profile regarding any given cell type, and an uncertainty measure dSpec, which measures the effect of noise on specificity. Using global gene expression data from the mouse brain, plant root and human white blood cells, we show that Spec identifies genes with variable expression levels that are nonetheless highly specific of particular cell types. When samples from different individuals are used, dSpec measures genes' transcriptional plasticity in each cell type. Our approach is broadly applicable to mapped gene expression measurements in stem cell biology, developmental biology, cancer biology and biomarker identification. As an example of such applications, we show that Spec identifies a new class of biomarkers, which exhibit variable expression without compromising specificity. The approach provides a unifying theoretical framework for quantifying specificity in the presence of noise, which is widely applicable across diverse biological systems.
10.1093/nar/gkr591
Mapping human development at single-cell resolution.
Behjati Sam,Lindsay Susan,Teichmann Sarah A,Haniffa Muzlifah
Development (Cambridge, England)
Human development is regulated by spatiotemporally restricted molecular programmes and is pertinent to many areas of basic biology and human medicine, such as stem cell biology, reproductive medicine and childhood cancer. Mapping human development has presented significant technological, logistical and ethical challenges. The availability of established human developmental biorepositories and the advent of cutting-edge single-cell technologies provide new opportunities to study human development. Here, we present a working framework for the establishment of a human developmental cell atlas exploiting single-cell genomics and spatial analysis. We discuss how the development atlas will benefit the scientific and clinical communities to advance our understanding of basic biology, health and disease.
10.1242/dev.152561
A computational approach towards a gene regulatory network for the developing Nematostella vectensis gut.
Botman Daniel,Röttinger Eric,Martindale Mark Q,de Jong Johann,Kaandorp Jaap A
PloS one
BACKGROUND:The starlet sea anemone Nematostella vectensis is a diploblastic cnidarian that expresses a set of conserved genes for gut formation during its early development. During the last decade, the spatial distribution of many of these genes has been visualized with RNA hybridization or protein immunolocalization techniques. However, due to N. vectensis' curved and changing morphology, quantification of these spatial data is problematic. A method is developed for two-dimensional gene expression quantification, which enables a numerical analysis and dynamic modeling of these spatial patterns. METHODS/RESULT:In this work, first standardized gene expression profiles are generated from publicly available N. vectensis embryo images that display mRNA and/or protein distributions. Then, genes expressed during gut formation are clustered based on their expression profiles, and further grouped based on temporal appearance of their gene products in embryonic development. Representative expression profiles are manually selected from these clusters, and used as input for a simulation-based optimization scheme. This scheme iteratively fits simulated profiles to the selected profiles, leading to an optimized estimation of the model parameters. Finally, a preliminary gene regulatory network is derived from the optimized model parameters. OUTLOOK:While the focus of this study is N. vectensis, the approach outlined here is suitable for inferring gene regulatory networks in the embryonic development of any animal, thus allowing to comparatively study gene regulation of gut formation in silico across various species.
10.1371/journal.pone.0103341
Correlating global gene regulation to angiogenesis in the developing chick extra-embryonic vascular system.
Javerzat Sophie,Franco Mélanie,Herbert John,Platonova Natalia,Peille Anne-Lise,Pantesco Véronique,De Vos John,Assou Said,Bicknell Roy,Bikfalvi Andreas,Hagedorn Martin
PloS one
BACKGROUND:Formation of blood vessels requires the concerted regulation of an unknown number of genes in a spatial-, time- and dosage-dependent manner. Determining genes, which drive vascular maturation is crucial for the identification of new therapeutic targets against pathological angiogenesis. METHODOLOGY/PRINCIPAL FINDINGS:[corrected] We accessed global gene regulation throughout maturation of the chick chorio-allantoic membrane (CAM), a highly vascularized tissue, using pan genomic microarrays. Seven percent of analyzed genes showed a significant change in expression (>2-fold, FDR<5%) with a peak occurring from E7 to E10, when key morphogenetic and angiogenic genes such as BMP4, SMO, HOXA3, EPAS1 and FGFR2 were upregulated, reflecting the state of an activated endothelium. At later stages, a general decrease in gene expression occurs, including genes encoding mitotic factors or angiogenic mediators such as CYR61, EPAS1, MDK and MYC. We identified putative human orthologs for 77% of significantly regulated genes and determined endothelial cell enrichment for 20% of the orthologs in silico. Vascular expression of several genes including ENC1, FSTL1, JAM2, LDB2, LIMS1, PARVB, PDE3A, PRCP, PTRF and ST6GAL1 was demonstrated by in situ hybridization. Up to 9% of the CAM genes were also overexpressed in human organs with related functions, such as placenta and lung or the thyroid. 21-66% of CAM genes enriched in endothelial cells were deregulated in several human cancer types (P<.0001). Interfering with PARVB (encoding parvin, beta) function profoundly changed human endothelial cell shape, motility and tubulogenesis, suggesting an important role of this gene in the angiogenic process. CONCLUSIONS/SIGNIFICANCE:Our study underlines the complexity of gene regulation in a highly vascularized organ during development. We identified a restricted number of novel genes enriched in the endothelium of different species and tissues, which may play crucial roles in normal and pathological angiogenesis.
10.1371/journal.pone.0007856
Detection and removal of spatial bias in multiwell assays.
Lachmann Alexander,Giorgi Federico M,Alvarez Mariano J,Califano Andrea
Bioinformatics (Oxford, England)
MOTIVATION:Multiplex readout assays are now increasingly being performed using microfluidic automation in multiwell format. For instance, the Library of Integrated Network-based Cellular Signatures (LINCS) has produced gene expression measurements for tens of thousands of distinct cell perturbations using a 384-well plate format. This dataset is by far the largest 384-well gene expression measurement assay ever performed. We investigated the gene expression profiles of a million samples from the LINCS dataset and found that the vast majority (96%) of the tested plates were affected by a significant 2D spatial bias. RESULTS:Using a novel algorithm combining spatial autocorrelation detection and principal component analysis, we could remove most of the spatial bias from the LINCS dataset and show in parallel a dramatic improvement of similarity between biological replicates assayed in different plates. The proposed methodology is fully general and can be applied to any highly multiplexed assay performed in multiwell format. CONTACT:ac2248@columbia.edu SUPPLEMENTARY INFORMATION:Supplementary data are available at Bioinformatics online.
10.1093/bioinformatics/btw092
Human distal lung maps and lineage hierarchies reveal a bipotent progenitor.
Nature
Mapping the spatial distribution and molecular identity of constituent cells is essential for understanding tissue dynamics in health and disease. We lack a comprehensive map of human distal airways, including the terminal and respiratory bronchioles (TRBs), which are implicated in respiratory diseases. Here, using spatial transcriptomics and single-cell profiling of microdissected distal airways, we identify molecularly distinct TRB cell types that have not-to our knowledge-been previously characterized. These include airway-associated LGR5 fibroblasts and TRB-specific alveolar type-0 (AT0) cells and TRB secretory cells (TRB-SCs). Connectome maps and organoid-based co-cultures reveal that LGR5 fibroblasts form a signalling hub in the airway niche. AT0 cells and TRB-SCs are conserved in primates and emerge dynamically during human lung development. Using a non-human primate model of lung injury, together with human organoids and tissue specimens, we show that alveolar type-2 cells in regenerating lungs transiently acquire an AT0 state from which they can differentiate into either alveolar type-1 cells or TRB-SCs. This differentiation programme is distinct from that identified in the mouse lung. Our study also reveals mechanisms that drive the differentiation of the bipotent AT0 cell state into normal or pathological states. In sum, our findings revise human lung cell maps and lineage trajectories, and implicate an epithelial transitional state in primate lung regeneration and disease.
10.1038/s41586-022-04541-3
Establishing neuronal diversity in the spinal cord: a time and a place.
Sagner Andreas,Briscoe James
Development (Cambridge, England)
The vertebrate spinal cord comprises multiple functionally distinct neuronal cell types arranged in characteristic positions. During development, these different types of neurons differentiate from transcriptionally distinct neural progenitors that are arrayed in discrete domains along the dorsal-ventral and anterior-posterior axes of the embryonic spinal cord. This organization arises in response to morphogen gradients acting upstream of a gene regulatory network, the architecture of which determines the spatial and temporal pattern of gene expression. In recent years, substantial progress has been made in deciphering the regulatory network that underlies the specification of distinct progenitor and neuronal cell identities. In this Review, we outline how distinct neuronal cell identities are established in response to spatial and temporal patterning systems, and outline novel experimental approaches to study the emergence and function of neuronal diversity in the spinal cord.
10.1242/dev.182154
Chromatin remodeling enzyme Snf2h regulates embryonic lens differentiation and denucleation.
Development (Cambridge, England)
Ocular lens morphogenesis is a model for investigating mechanisms of cellular differentiation, spatial and temporal gene expression control, and chromatin regulation. Brg1 (Smarca4) and Snf2h (Smarca5) are catalytic subunits of distinct ATP-dependent chromatin remodeling complexes implicated in transcriptional regulation. Previous studies have shown that Brg1 regulates both lens fiber cell differentiation and organized degradation of their nuclei (denucleation). Here, we employed a conditional Snf2h(flox) mouse model to probe the cellular and molecular mechanisms of lens formation. Depletion of Snf2h induces premature and expanded differentiation of lens precursor cells forming the lens vesicle, implicating Snf2h as a key regulator of lens vesicle polarity through spatial control of Prox1, Jag1, p27(Kip1) (Cdkn1b) and p57(Kip2) (Cdkn1c) gene expression. The abnormal Snf2h(-/-) fiber cells also retain their nuclei. RNA profiling of Snf2h(-/) (-) and Brg1(-/-) eyes revealed differences in multiple transcripts, including prominent downregulation of those encoding Hsf4 and DNase IIβ, which are implicated in the denucleation process. In summary, our data suggest that Snf2h is essential for the establishment of lens vesicle polarity, partitioning of prospective lens epithelial and fiber cell compartments, lens fiber cell differentiation, and lens fiber cell nuclear degradation.
10.1242/dev.135285
Human corneal epithelial subpopulations: oxygen dependent ex vivo expansion and transcriptional profiling.
Bath Chris
Acta ophthalmologica
Corneal epithelium is being regenerated throughout life by limbal epithelial stem cells (LESCs) believed to be located in histologically defined stem cell niches in corneal limbus. Defective or dysfunctional LESCs result in limbal stem cell deficiency (LSCD) causing pain and decreased visual acuity. Since the first successful treatment of LSCD by transplantation of ex vivo expanded LESCs in 1997, many attempts have been carried out to optimize culture conditions to improve the outcome of surgery. To date, progress in this field of bioengineering is substantially hindered by both the lack of specific biomarkers of LESCs and the lack of a precise molecular characterization of in situ epithelial subpopulations. The aim of this dissertation was to optimize culture systems with regard to the environmental oxygen concentration for selective ex vivo expansion of LESCs and to analyse in situ subpopulations in human corneal epithelium using a combination of laser capture microdissection and RNA sequencing for global transcriptomic profiling. We compared dissociation cultures, using either expansion on γ-irradiated NIH/3T3 feeder cells in serum-rich medium or expansion directly on plastic in serum-free EpiLife medium, using a range of physiologically relevant oxygen concentrations (2%, 5%, 10%, 15% and 20%). Using immunocytochemistry and advanced fluorescence microscopy, cells were characterized regarding growth, cell cycle distribution, colony-forming efficiency (CFE), phenotypes and cytomorphometry. Limbal epithelial cells expanded in 2% O2 exhibited slow growth, low fraction of cells in S/G2 , high CFE, high expression of stem cell markers ABCG2 and p63α, and low fraction of differentiation marker CK3 resembling a LESC phenotype. The effect of hypoxia to maintain LESCs in culture was not dependent on the system used for propagation (Bath et al. 2013a). Laser capture microdissection was used to isolate cellular subpopulations in situ from the spatially defined differentiation pathway in human corneal epithelium according to an optimized protocol for maintenance of expression profiles. Isolated total RNA from basal limbal crypts (BLCs), superficial limbal crypts (SLCs), paracentral/central cornea and limbal stroma was amplified and converted to fragmented cDNA libraries for use in deep paired-end next-generation sequencing. Global transcriptional profiling was carried out using bioinformatics. The location of primitive cells in BLCs, migratory and activated cells in SLCs and differentiated cells in paracentral/central cornea was evident from mapping of significantly upregulated genes in each compartment to the gene ontology (GO). Interestingly, many GO terms in BLCs were also involved in neurogenic processes, whereas many GO terms in SLCs were related to vasculature. Mapping upregulated genes in BLCs to pathway annotations in Kyoto Encyclopedia of Genes and Genomes described many active pathways as signalling and cancer-associated pathways. We supply extensive information on possible novel biomarkers, reveal insight into both active pathways and novel regulators of LESCs such as Lrig1 and SOX9 and provide an immense amount of data for future exploration (Bath et al. 2013b). Selective ex vivo expansion of LESCs in hypoxia and the comprehensive molecular characterization of corneal epithelial subpopulations in situ are expected to be beneficial for the future treatment of LSCD by cultured limbal epithelial transplantation.
10.1111/aos.12157
Pixelated spatial gene expression analysis from tissue.
Ganguli A,Ornob A,Spegazzini N,Liu Y,Damhorst G,Ghonge T,Thornton B,Konopka C J,Dobrucki W,Clare S E,Bhargava R,Smith A M,Kosari F,Bashir R
Nature communications
Here, we present a technique that performs on-chip picoliter real-time reverse transcriptase loop mediated isothermal amplification (RT-LAMP) reactions on a histological tissue section without any analyte purification while preserving the native spatial location of the nucleic acid molecules. We demonstrate this method by amplifying TOP2A messenger RNA (mRNA) in a prostate cancer xenograft with 100 µm spatial resolution and by visualizing the variation in threshold time of amplification across the tissue. The on-chip reaction was validated by mRNA fluorescence in situ hybridization (mFISH) from cells in the tissue section. The entire process, from tissue loading on microchip to results from RT-LAMP can be carried out in less than 2 h. We anticipate that this technique, with its ease of use, fast turnaround, and quantitative molecular outputs, would become an invaluable tissue analysis tool for researchers and clinicians in the biomedical arena.
10.1038/s41467-017-02623-9
Global analysis of gene expression patterns in developing mouse neocortex using serial analysis of gene expression.
Gunnersen J M,Augustine C,Spirkoska V,Kim M,Brown M,Tan S-S
Molecular and cellular neurosciences
Molecular inventories of the developing mouse neocortex before and after birth were generated using the global gene expression profiling tool serial analysis of gene expression (SAGE). Libraries were generated from embryonic day 15 and postnatal day 1 mouse neocortex and more than 40,000 tags were collected (20,211 and 22,001 tags, representing 11,706 and 12,402 transcripts, respectively). Comparison of the two libraries resulted in the identification of 321 transcripts that were differentially expressed (P < 0.05). Differential expression was independently verified for selected genes by Northern blotting, and in situ hybridization revealed spatial expression patterns in the neocortex. Differentially expressed transcripts included genes known to be important in neocortical development (e.g., brain factor 1, neuroD2, and Id2), genes not previously associated with neocortical development (such as brahma-related gene 1, receptor for activated C-kinase I, hypermethylated in cancer 2, and Evi9), and genes of unknown identity or function.
10.1006/mcne.2001.1098
Spatial profiling of chromatin accessibility in mouse and human tissues.
Nature
Cellular function in tissue is dependent on the local environment, requiring new methods for spatial mapping of biomolecules and cells in the tissue context. The emergence of spatial transcriptomics has enabled genome-scale gene expression mapping, but the ability to capture spatial epigenetic information of tissue at the cellular level and genome scale is lacking. Here we describe a method for spatially resolved chromatin accessibility profiling of tissue sections using next-generation sequencing (spatial-ATAC-seq) by combining in situ Tn5 transposition chemistry and microfluidic deterministic barcoding. Profiling mouse embryos using spatial-ATAC-seq delineated tissue-region-specific epigenetic landscapes and identified gene regulators involved in the development of the central nervous system. Mapping the accessible genome in the mouse and human brain revealed the intricate arealization of brain regions. Applying spatial-ATAC-seq to tonsil tissue resolved the spatially distinct organization of immune cell types and states in lymphoid follicles and extrafollicular zones. This technology progresses spatial biology by enabling spatially resolved chromatin accessibility profiling to improve our understanding of cell identity, cell state and cell fate decision in relation to epigenetic underpinnings in development and disease.
10.1038/s41586-022-05094-1
Spatial transcriptomics stratifies psoriatic disease severity by emergent cellular ecosystems.
Science immunology
Whereas the cellular and molecular features of human inflammatory skin diseases are well characterized, their tissue context and systemic impact remain poorly understood. We thus profiled human psoriasis (PsO) as a prototypic immune-mediated condition with a high predilection for extracutaneous involvement. Spatial transcriptomics (ST) analyses of 25 healthy, active lesion, and clinically uninvolved skin biopsies and integration with public single-cell transcriptomics data revealed marked differences in immune microniches between healthy and inflamed skin. Tissue-scale cartography further identified core disease features across all active lesions, including the emergence of an inflamed suprabasal epidermal state and the presence of B lymphocytes in lesional skin. Both lesional and distal nonlesional samples were stratified by skin disease severity and not by the presence of systemic disease. This segregation was driven by macrophage-, fibroblast-, and lymphatic-enriched spatial regions with gene signatures associated with metabolic dysfunction. Together, these findings suggest that mild and severe forms of PsO have distinct molecular features and that severe PsO may profoundly alter the cellular and metabolic composition of distal unaffected skin sites. In addition, our study provides a valuable resource for the research community to study spatial gene organization of healthy and inflamed human skin.
10.1126/sciimmunol.abq7991
Trastuzumab and pertuzumab produce changes in morphology and estrogen receptor signaling in ovarian cancer xenografts revealing new treatment strategies.
Faratian Dana,Zweemer Annelien J M,Nagumo Yoko,Sims Andrew H,Muir Morwenna,Dodds Michael,Mullen Peter,Um Inhwa,Kay Charlene,Hasmann Max,Harrison David J,Langdon Simon P
Clinical cancer research : an official journal of the American Association for Cancer Research
PURPOSE:The aim of this study was to investigate the antitumor effects of HER2-directed combination therapy in ovarian cancer xenograft models to evaluate their potential. The combinations of trastuzumab and pertuzumab, and trastuzumab and aromatase inhibitor therapy were investigated. EXPERIMENTAL DESIGN:The effects of trastuzumab, pertuzumab, and letrozole on growth response, apoptosis, morphology, and gene and protein expression were evaluated in the SKOV3 ovarian cancer cell line xenograft and a panel of five human ovarian xenografts derived directly from clinical specimens. RESULTS:The combination of HER2-directed antibodies showed enhanced antitumor activity compared with single antibody therapy in the SKOV3 xenograft model. Apoptosis, morphology, and estrogen-regulated gene expression were modulated by these antibodies in both spatial and temporal manners. A panel of ovarian cancer xenografts showed differential growth responses to the combination of trastuzumab and pertuzumab. High HER2 expression and increasing HER3 protein expression on treatment were associated with growth response. In trastuzumab-treated SKOV3 tumors, there was a change in tumor morphology, with a reduction in frequency of estrogen receptor alpha (ERα)-negative clear cell areas. Trastuzumab, but not pertuzumab, increased expression of ERα in SKOV3 xenografts when analyzed by quantitative immunofluorescence. ERα and downstream signaling targets were modulated by trastuzumab alone and in combination. Trastuzumab enhanced the responsiveness of SKOV3 xenografts to letrozole when given in combination. CONCLUSIONS:These data suggest that trastuzumab in combination with pertuzumab could be an effective approach in high HER2-expressing ovarian cancers and could also enhance sensitivity to endocrine therapy in ERα-positive ovarian cancer.
10.1158/1078-0432.CCR-10-2461
A hidden Ising model for ChIP-chip data analysis.
Mo Qianxing,Liang Faming
Bioinformatics (Oxford, England)
MOTIVATION:Chromatin immunoprecipitation (ChIP) coupled with tiling microarray (chip) experiments have been used in a wide range of biological studies such as identification of transcription factor binding sites and investigation of DNA methylation and histone modification. Hidden Markov models are widely used to model the spatial dependency of ChIP-chip data. However, parameter estimation for these models is typically either heuristic or suboptimal, leading to inconsistencies in their applications. To overcome this limitation and to develop an efficient software, we propose a hidden ferromagnetic Ising model for ChIP-chip data analysis. RESULTS:We have developed a simple, but powerful Bayesian hierarchical model for ChIP-chip data via a hidden Ising model. Metropolis within Gibbs sampling algorithm is used to simulate from the posterior distribution of the model parameters. The proposed model naturally incorporates the spatial dependency of the data, and can be used to analyze data with various genomic resolutions and sample sizes. We illustrate the method using three publicly available datasets and various simulated datasets, and compare it with three closely related methods, namely TileMap HMM, tileHMM and BAC. We find that our method performs as well as TileMap HMM and BAC for the high-resolution data from Affymetrix platform, but significantly outperforms the other three methods for the low-resolution data from Agilent platform. Compared with the BAC method which also involves MCMC simulations, our method is computationally much more efficient. AVAILABILITY:A software called iChip is freely available at http://www.bioconductor.org/. CONTACT:moq@mskcc.org.
10.1093/bioinformatics/btq032
Spatial and temporal intratumour heterogeneity has potential consequences for single biopsy-based neuroblastoma treatment decisions.
Nature communications
Intratumour heterogeneity is a major cause of treatment failure in cancer. We present in-depth analyses combining transcriptomic and genomic profiling with ultra-deep targeted sequencing of multiregional biopsies in 10 patients with neuroblastoma, a devastating childhood tumour. We observe high spatial and temporal heterogeneity in somatic mutations and somatic copy-number alterations which are reflected on the transcriptomic level. Mutations in some druggable target genes including ALK and FGFR1 are heterogeneous at diagnosis and/or relapse, raising the issue whether current target prioritization and molecular risk stratification procedures in single biopsies are sufficiently reliable for therapy decisions. The genetic heterogeneity in gene mutations and chromosome aberrations observed in deep analyses from patient courses suggest clonal evolution before treatment and under treatment pressure, and support early emergence of metastatic clones and ongoing chromosomal instability during disease evolution. We report continuous clonal evolution on mutational and copy number levels in neuroblastoma, and detail its implications for therapy selection, risk stratification and therapy resistance.
10.1038/s41467-021-26870-z
Mutational patterns in chemotherapy resistant muscle-invasive bladder cancer.
Liu David,Abbosh Philip,Keliher Daniel,Reardon Brendan,Miao Diana,Mouw Kent,Weiner-Taylor Amaro,Wankowicz Stephanie,Han Garam,Teo Min Yuen,Cipolla Catharine,Kim Jaegil,Iyer Gopa,Al-Ahmadie Hikmat,Dulaimi Essel,Chen David Y T,Alpaugh R Katherine,Hoffman-Censits Jean,Garraway Levi A,Getz Gad,Carter Scott L,Bellmunt Joaquim,Plimack Elizabeth R,Rosenberg Jonathan E,Van Allen Eliezer M
Nature communications
Despite continued widespread use, the genomic effects of cisplatin-based chemotherapy and implications for subsequent treatment are incompletely characterized. Here, we analyze whole exome sequencing of matched pre- and post-neoadjuvant cisplatin-based chemotherapy primary bladder tumor samples from 30 muscle-invasive bladder cancer patients. We observe no overall increase in tumor mutational burden post-chemotherapy, though a significant proportion of subclonal mutations are unique to the matched pre- or post-treatment tumor, suggesting chemotherapy-induced and/or spatial heterogeneity. We subsequently identify and validate a novel mutational signature in post-treatment tumors consistent with known characteristics of cisplatin damage and repair. We find that post-treatment tumor heterogeneity predicts worse overall survival, and further observe alterations in cell-cycle and immune checkpoint regulation genes in post-treatment tumors. These results provide insight into the clinical and genomic dynamics of tumor evolution with cisplatin-based chemotherapy, suggest mechanisms of clinical resistance, and inform development of clinically relevant biomarkers and trials of combination therapies.
10.1038/s41467-017-02320-7
Candidate genes for panhypopituitarism identified by gene expression profiling.
Physiological genomics
Mutations in the transcription factors PROP1 and PIT1 (POU1F1) lead to pituitary hormone deficiency and hypopituitarism in mice and humans. The dysmorphology of developing Prop1 mutant pituitaries readily distinguishes them from those of Pit1 mutants and normal mice. This and other features suggest that Prop1 controls the expression of genes besides Pit1 that are important for pituitary cell migration, survival, and differentiation. To identify genes involved in these processes we used microarray analysis of gene expression to compare pituitary RNA from newborn Prop1 and Pit1 mutants and wild-type littermates. Significant differences in gene expression were noted between each mutant and their normal littermates, as well as between Prop1 and Pit1 mutants. Otx2, a gene critical for normal eye and pituitary development in humans and mice, exhibited elevated expression specifically in Prop1 mutant pituitaries. We report the spatial and temporal regulation of Otx2 in normal mice and Prop1 mutants, and the results suggest Otx2 could influence pituitary development by affecting signaling from the ventral diencephalon and regulation of gene expression in Rathke's pouch. The discovery that Otx2 expression is affected by Prop1 deficiency provides support for our hypothesis that identifying molecular differences in mutants will contribute to understanding the molecular mechanisms that control pituitary organogenesis and lead to human pituitary disease.
10.1152/physiolgenomics.00080.2011
Spatial and transcriptional heterogeneity of pancreatic beta cell neogenesis revealed by a time-resolved reporter system.
Diabetologia
AIMS/HYPOTHESIS:While pancreatic beta cells have been shown to originate from endocrine progenitors in ductal regions, it remains unclear precisely where beta cells emerge from and which transcripts define newborn beta cells. We therefore investigated characteristics of newborn beta cells extracted by a time-resolved reporter system. METHODS:We established a mouse model, 'Ins1-GFP; Timer', which provides spatial information during beta cell neogenesis with high temporal resolution. Single-cell RNA-sequencing (scRNA-seq) was performed on mouse beta cells sorted by fluorescent reporter to uncover transcriptomic profiles of newborn beta cells. scRNA-seq of human embryonic stem cell (hESC)-derived beta-like cells was also performed to compare newborn beta cell features between mouse and human. RESULTS:Fluorescence imaging of Ins1-GFP; Timer mouse pancreas successfully dissected newly generated beta cells as green fluorescence-dominant cells. This reporter system revealed that, as expected, some newborn beta cells arise close to the ducts (β); unexpectedly, the others arise away from the ducts and adjacent to blood vessels (β). Single-cell transcriptomic analyses demonstrated five distinct populations among newborn beta cells, confirming spatial heterogeneity of beta cell neogenesis such as high probability of glucagon-positive β, musculoaponeurotic fibrosarcoma oncogene family B (MafB)-positive β and musculoaponeurotic fibrosarcoma oncogene family A (MafA)-positive β cells. Comparative analysis with scRNA-seq data of mouse newborn beta cells and hESC-derived beta-like cells uncovered transcriptional similarity between mouse and human beta cell neogenesis including microsomal glutathione S-transferase 1 (MGST1)- and synaptotagmin 13 (SYT13)-highly-expressing state. CONCLUSIONS/INTERPRETATION:The combination of time-resolved histological imaging with single-cell transcriptional mapping demonstrated novel features of spatial and transcriptional heterogeneity in beta cell neogenesis, which will lead to a better understanding of beta cell differentiation for future cell therapy. DATA AVAILABILITY:Raw and processed single-cell RNA-sequencing data for this study has been deposited in the Gene Expression Omnibus under accession number GSE155742.
10.1007/s00125-022-05662-0
Cell growth in aggregates determines gene expression, proliferation, survival, chemoresistance, and sensitivity to immune effectors in follicular lymphoma.
Gravelle Pauline,Jean Christine,Familiades Julien,Decaup Emilie,Blanc Amandine,Bezombes-Cagnac Christine,Laurent Camille,Savina Ariel,Fournié Jean-Jacques,Laurent Guy
The American journal of pathology
Lymphomas grow as dense aggregates in patients, but whether this spatial organization affects lymphoma cell biology is unknown. We grew follicular lymphoma (FL) cells in vitro as multicellular aggregates of lymphoma cells to investigate this question. Gene expression analysis revealed that 612 genes were differentially expressed when cells grew in multicellular aggregates of lymphoma cells rather than in suspension. These genes correspond to several GO biological processes, such as hypoxia, activation of NF-κB pathway, and negative regulation of cell cycle, a gene signature also found in the transcriptomes from FL biopsies. Pimonidazole staining, HIF-1A accumulation, and VEGFA release confirmed that cells in multicellular aggregates of lymphoma cells actually respond to hypoxia. In adaptation to such conditions, they also displayed an activated NF-κB pathway and a quiescent status far more frequently than in suspension. When cultured in three dimensions, FL cells display resistance to doxorubicin and bendamustine, two drugs largely used in FL therapy, compared to FL cultured in suspension. Finally, multicellular aggregates of lymphoma cells were also found to be less sensitive to purified natural killer cells. To conclude, our study shows that in FL, spatial organization results in dramatic changes in FL biology, including gene expression, proliferation, drug resistance, and immune escape.
10.1016/j.ajpath.2013.09.018
Local network component analysis for quantifying transcription factor activities.
Shi Qianqian,Zhang Chuanchao,Guo Weifeng,Zeng Tao,Lu Lina,Jiang Zhonglin,Wang Ziming,Liu Juan,Chen Luonan
Methods (San Diego, Calif.)
Transcription factors (TFs) could regulate physiological transitions or determine stable phenotypic diversity. The accurate estimation on TF regulatory signals or functional activities is of great significance to guide biological experiments or elucidate molecular mechanisms, but still remains challenging. Traditional methods identify TF regulatory signals at the population level, which masks heterogeneous regulation mechanisms in individuals or subgroups, thus resulting in inaccurate analyses. Here, we propose a novel computational framework, namely local network component analysis (LNCA), to exploit data heterogeneity and automatically quantify accurate transcription factor activity (TFA) in practical terms, through integrating the partitioned expression sets (i.e., local information) and prior TF-gene regulatory knowledge. Specifically, LNCA adopts an adaptive optimization strategy, which evaluates the local similarities of regulation controls and corrects biases during data integration, to construct the TFA landscape. In particular, we first numerically demonstrate the effectiveness of LNCA for the simulated data sets, compared with traditional methods, such as FastNCA, ROBNCA and NINCA. Then, we apply our model to two real data sets with implicit temporal or spatial regulation variations. The results show that LNCA not only recognizes the periodic mode along the S. cerevisiae cell cycle process, but also substantially outperforms over other methods in terms of accuracy and consistency. In addition, the cross-validation study for glioblastomas multiforme (GBM) indicates that the TFAs, identified by LNCA, can better distinguish clinically distinct tumor groups than the expression values of the corresponding TFs, thus opening a new way to classify tumor subtypes and also providing a novel insight into cancer heterogeneity. AVAILABILITY:LNCA was implemented as a Matlab package, which is available at http://sysbio.sibcb.ac.cn/cb/chenlab/software.htm/LNCApackage_0.1.rar.
10.1016/j.ymeth.2017.06.018
Molecular portraits of intratumoral heterogeneity in human ovarian cancer.
Choi Yoon Pyo,Shim Hyo Sup,Gao Ming-Qing,Kang Suki,Cho Nam Hoon
Cancer letters
One of the most common characteristic profiles of cancer is intratumoral heterogeneity (ITH). We aimed to clarify the molecular profiles and biological significance of ITH with relation to cancer stem cell (CSC). We analyzed five primary cultured clones generated from different spatial zones, front and rear zone, of a fresh-frozen ovarian tumor tissue, performing ATP-CRA, conventional RT-PCR, side population (SP) analysis, flow cytometry immunophenotyping, and cell proliferation assays. We also carried out array CGH and Ingenuity Pathways Analysis (IPA) between SP and non-SP (NSP) cells. Clones from tumor front zone showed phenotypically and genetically distinct subpopulations with relatively higher SP proportions, CD24(+) and CD117(+) expression, and chemotherapeutic resistance. We demonstrate that phenotype of SP cells in heterogeneous clones of human ovarian cancer was closely related to CD24(+), CD117(+), and combined CD117(+)/CD24(+) fractions. Chromosomal alterations in SP cells relative to NSP cells were closely related to the novel core networks of cancer stem cell-like cells (CSCs), such as cycle checkpoint regulation, notch, PTEN, wnt/β-catenin, PI3K/AKT, integrin, and cytokine and chemokine signaling. ITH could arise from clonal diversity closely related to CSC-like molecules, as evidenced by accumulated genetic, transcriptional and gene products alterations in SP.
10.1016/j.canlet.2011.03.018