logo logo
Differential expression of the embryo/cancer gene ECSA(DPPA2), the cancer/testis gene BORIS and the pluripotency structural gene OCT4, in human preimplantation development. Monk Marilyn,Hitchins Megan,Hawes Susan Molecular human reproduction In this paper, we examine the expression profiles of two new putative pluripotent stem cell genes, the embryo/cancer sequence A gene (ECSA) and the cancer/testis gene Brother Of the Regulator of Imprinted Sites (BORIS), in human oocytes, preimplantation embryos, primordial germ cells (PGCs) and embryo stem (ES) cells. Their expression profiles are compared with that of the well-known pluripotency gene, OCT4, using a primer design that avoids amplification of the multiple OCT4 pseudogenes. As expected, OCT4 is high in human oocytes, down-regulated in early cleavage stages and then expressed de novo in human blastocysts and PGCs. BORIS and ECSA show distinct profiles of expression in that BORIS is predominantly expressed in the early stages of preimplantation development, in oocytes and 4-cell embryos, whereas ECSA is predominantly expressed in the later stages, blastocysts and PGCs. BORIS is not detected in blastocysts, PGCs or other fetal and adult somatic tissue tested. Thus, BORIS and ECSA may be involved in two different aspects of reprogramming in development, viz., in late gametogenesis, and at the time of formation of the ES cells (inner cell mass (ICM) and PGC), respectively. However, in human ES cells, where a deprogrammed stem cell state is stably established in culture, an immunofluoresence study shows that all three genes are co-expressed at the protein level. Thus, following their derivation from ICM cells, ES cells may undergo further transformation in culture to express a number of embryo and germ line stem cell functions, which, in normal development, show different temporal and spatial specificity of expression. 10.1093/molehr/gan025
The Diagnostic, Prognostic, and Therapeutic Potential of Adaptive Immune Receptor Repertoire Profiling in Cancer. Cowell Lindsay G Cancer research Lymphocytes play a critical role in antitumor immune responses. They are directly targeted by some therapies, and the composition and spatial organization of intratumor T-cell populations is prognostic in some cancer types. A better understanding of lymphocyte population dynamics over the course of disease and in response to therapy is urgently needed to guide therapy decisions and to develop new therapy targets. Deep sequencing of the repertoire of antigen receptor-encoding genes expressed in a lymphocyte population has become a widely used approach for profiling the population's immune status. Lymphocyte antigen receptor repertoire deep sequencing data can be used to assess the clonal richness and diversity of lymphocyte populations; to track clone members over time, between tissues, and across lymphocyte subsets; to detect clonal expansion; and to detect the recruitment of new clones into a tissue. Repertoire sequencing is thus a critical complement to other methods of lymphocyte and immune profiling in cancer. This review describes the current state of knowledge based on repertoire sequencing studies conducted on human cancer patients, with a focus on studies of the T-cell receptor beta chain locus. The review then outlines important questions left unanswered and suggests future directions for the field. 10.1158/0008-5472.CAN-19-1457
Multiparametric analysis of colorectal cancer immune responses. World journal of gastroenterology Colorectal cancer (CRC) is a heterogeneous disease, with a diverse and plastic immune cell infiltrate. These immune cells play an important role in regulating tumour growth - progression or elimination. Some populations of cells have a strong correlation with disease-free survival, making them useful prognostic markers. In particular, the infiltrate of CD3 and CD8 T cells into CRC tumours has been validated worldwide as a valuable indicator of patient prognosis. However, the heterogeneity of the immune response, both between patients with tumours of different molecular subtypes, and within the tumour itself, necessitates the use of multiparametric analysis in the investigation of tumour-specific immune responses. This review will outline the multiparametric analysis techniques that have been developed and applied to studying the role of immune cells in the tumour, with a focus on colorectal cancer. Because much of the data in this disease relates to T cell subsets and heterogeneity, we have used T cell populations as examples throughout. Flow and mass cytometry give a detailed representation of the cells within the tumour in a single-cell suspension on a per-cell basis. Imaging technologies, such as imaging mass cytometry, are used to investigate increasing numbers of markers whilst retaining the spatial and structural information of the tumour section and the infiltrating immune cells. Together, the analyses of multiple immune parameters can provide valuable information to guide clinical decision-making in CRC. 10.3748/wjg.v24.i27.2995
Notch2-mediated plasticity between marginal zone and follicular B cells. Nature communications Follicular B (FoB) and marginal zone B (MZB) cells are functionally and spatially distinct mature B cell populations in the spleen, originating from a Notch2-dependent fate decision after splenic influx of immature transitional B cells. In the B cell follicle, a Notch2-signal is provided by DLL-1-expressing fibroblasts. However, it is unclear whether FoB cells, which are in close contact with these DLL-1 expressing fibroblasts, can also differentiate to MZB cells if they receive a Notch2-signal. Here, we show induced Notch2IC-expression in FoB cells re-programs mature FoB cells into bona fide MZB cells as is evident from the surface phenotype, localization, immunological function and transcriptome of these cells. Furthermore, the lineage conversion from FoB to MZB cells occurs in immunocompetent wildtype mice. These findings demonstrate plasticity between mature FoB and MZB cells that can be driven by a singular signaling event, the activation of Notch2. 10.1038/s41467-021-21359-1
Distinct and shared gene expression for human innate versus adaptive helper lymphoid cells. Ercolano Giuseppe,Wyss Tania,Salomé Bérengère,Romero Pedro,Trabanelli Sara,Jandus Camilla Journal of leukocyte biology Innate lymphoid cells (ILCs) are the latest identified innate immune cell family. Given their similarity in transcription factor expression and cytokine secretion profiles, ILCs have been considered as the innate phenocopy of CD4 Th cells. Here, we explored the transcriptome of circulating human ILC subsets as opposed to CD4 Th cell subsets. We describe transcriptomic differences between total ILCs and total CD4 Th cells, as well as between paired innate and adaptive cell subsets (ILC1 vs. Th1; ILC2 vs. Th2; and ILC3 vs. Th17 cells). In particular, we observed differences in expression of genes involved in cell trafficking such as CCR1, CCR6 and CXCR3, innate activation and inhibitory functions, including CD119, 2B4, TIGIT, and CTLA-4, and neuropeptide receptors, such as VIPR2. Moreover, we report for the first time on distinct expression of long noncoding RNAs (lncRNAs) in innate vs. adaptive cells, arguing for a potential role of lncRNA in shaping human ILC biology. Altogether, our results point for unique, rather than redundant gene organization in ILCs compared to CD4 Th cells, in regard to kinetics, fine-tuning and spatial organization of the immune response. 10.1002/JLB.5MA0120-209R
HPV+ oropharyngeal squamous cell carcinomas from patients with two tumors display synchrony of viral genomes yet discordant mutational profiles and signatures. Carcinogenesis Human papillomavirus (HPV) positive oropharyngeal squamous cell carcinoma (HPV + OPSCC) is increasing in prevalence in the USA, as are cases of patients with multiple HPV + OPSCCs (mHPV + OPSCC). mHPV + OPSCCs present a unique opportunity to examine HPV + OPSCC mutation acquisition and evolution. We performed sequencing of the viral genome, somatic exome and somatic transcriptome from 8 patients each with 2 spatially distinct HPV + OPSCCs, and 37 'traditional' HPV + OPSCCs to first address if paired tumors are caused by the same viral isolate and next, if acquired alterations, and the underlying processes driving mutagenesis, are shared within pairs. All tumor pairs contained viral genomes from the same HPV type 16 sublineage and differed by 0-2 clonal single nucleotide polymorphisms (SNPs), suggesting infection with the same viral isolate. Despite this, there was significant discordance in expression profiles, mutational burden and mutational profiles between tumors in a pair, with only two pairs sharing any overlapping mutations (3/3343 variants). Within tumor pairs there was a striking discrepancy of mutational signatures, exemplified by no paired tumors sharing high APOBEC mutational burden. Here, leveraging mHPV + OPSCCs as a model system to study mutation acquisition in virally mediated tumors, in which the germline, environmental exposures, immune surveillance and tissue/organ type were internally controlled, we demonstrate that despite infection by the same viral isolate, paired mHPV + OPSCCs develop drastically different somatic alterations and even more strikingly, appear to be driven by disparate underlying mutational processes. Thus, despite a common starting point, HPV + OPSCCs evolve through variable mutational processes with resultant stochastic mutational profiles. 10.1093/carcin/bgaa111
Genomic profiling in Down syndrome acute lymphoblastic leukemia identifies histone gene deletions associated with altered methylation profiles. Loudin M G,Wang J,Leung H-C Eastwood,Gurusiddappa S,Meyer J,Condos G,Morrison D,Tsimelzon A,Devidas M,Heerema N A,Carroll A J,Plon S E,Hunger S P,Basso G,Pession A,Bhojwani D,Carroll W L,Rabin K R Leukemia Patients with Down syndrome (DS) and acute lymphoblastic leukemia (ALL) have distinct clinical and biological features. Whereas most DS-ALL cases lack the sentinel cytogenetic lesions that guide risk assignment in childhood ALL, JAK2 mutations and CRLF2 overexpression are highly enriched. To further characterize the unique biology of DS-ALL, we performed genome-wide profiling of 58 DS-ALL and 68 non-DS (NDS) ALL cases by DNA copy number, loss of heterozygosity, gene expression and methylation analyses. We report a novel deletion within the 6p22 histone gene cluster as significantly more frequent in DS-ALL, occurring in 11 DS (22%) and only 2 NDS cases (3.1%) (Fisher's exact P=0.002). Homozygous deletions yielded significantly lower histone expression levels, and were associated with higher methylation levels, distinct spatial localization of methylated promoters and enrichment of highly methylated genes for specific pathways and transcription factor-binding motifs. Gene expression profiling demonstrated heterogeneity of DS-ALL cases overall, with supervised analysis defining a 45-transcript signature associated with CRLF2 overexpression. Further characterization of pathways associated with histone deletions may identify opportunities for novel targeted interventions. 10.1038/leu.2011.128
Ventral prostate predominant l, a novel mouse gene expressed exclusively in the prostate. Wubah Judith A,Fischer Carolyn M,Rolfzen Laura N,Khalili May,Kang Jason,Green Jeffrey E,Bieberich Charles J The Prostate BACKGROUND:Despite the region-specific nature of human prostate disease, there is a paucity of information regarding the molecular basis of prostate regionalization and patterning. To elucidate genetic mechanisms that underlie prostate growth and development, we investigated differential gene expression in mouse prostate lobes. METHODS:mRNA differential display analysis was used to identify differentially expressed genes during development of ventral, anterior, and dorsolateral prostate lobes. Differential gene expression was confirmed by Northern blot analysis and RT-PCR. RESULTS:A novel gene, Ventral prostate predominant1 (Vpp1) was identified. Vpp1 mRNA was evident in all lobes but accumulated predominantly in the ventral prostate, and was detected on postnatal day 7 through adulthood exclusively in the prostate gland. The steady-state level of Vpp1 mRNA decreased markedly in response to castration, suggesting androgen regulation of Vpp1 expression. Analysis of TRAMP tumors demonstrated a dramatic decrease in the level of Vpp1 mRNA. CONCLUSIONS:The spatial distribution and early postnatal onset of Vpp1 expression is consistent with a role for this gene in prostate regionalization. The absolute prostate specificity of Vpp1 expression may allow this gene to serve as a paradigm to study the molecular basis of gene expression that is restricted exclusively to the prostate gland.
Geometric compensation applied to image analysis of cell populations with morphological variability: a new role for a classical concept. Figueiredo Joana,Rodrigues Isabel,Ribeiro João,Fernandes Maria Sofia,Melo Soraia,Sousa Bárbara,Paredes Joana,Seruca Raquel,Sanches João M Scientific reports Immunofluorescence is the gold standard technique to determine the level and spatial distribution of fluorescent-tagged molecules. However, quantitative analysis of fluorescence microscopy images faces crucial challenges such as morphologic variability within cells. In this work, we developed an analytical strategy to deal with cell shape and size variability that is based on an elastic geometric alignment algorithm. Firstly, synthetic images mimicking cell populations with morphological variability were used to test and optimize the algorithm, under controlled conditions. We have computed expression profiles specifically assessing cell-cell interactions (IN profiles) and profiles focusing on the distribution of a marker throughout the intracellular space of single cells (RD profiles). To experimentally validate our analytical pipeline, we have used real images of cell cultures stained for E-cadherin, tubulin and a mitochondria dye, selected as prototypes of membrane, cytoplasmic and organelle-specific markers. The results demonstrated that our algorithm is able to generate a detailed quantitative report and a faithful representation of a large panel of molecules, distributed in distinct cellular compartments, independently of cell's morphological features. This is a simple end-user method that can be widely explored in research and diagnostic labs to unravel protein regulation mechanisms or identify protein expression patterns associated with disease. 10.1038/s41598-018-28570-z
Gene-expression and immunohistochemical study of specific T-cell subsets and accessory cell types in the transformation and prognosis of follicular lymphoma. Glas Annuska M,Knoops Laurent,Delahaye Leonie,Kersten Marie José,Kibbelaar Robby E,Wessels Lodewyk A,van Laar Ryan,van Krieken J Han J M,Baars Joke W,Raemaekers John,Kluin Philip M,van't Veer Laura J,de Jong Daphne Journal of clinical oncology : official journal of the American Society of Clinical Oncology PURPOSE:Despite the generally favorable clinical course in follicular lymphoma (FL), a minority of patients have a poor prognosis-with death within 3 years of diagnosis-most often due to transformation to aggressive disease. PATIENTS AND METHODS:In this study, we analyzed the potential of predicting early transformation on the basis of gene expression and immunologic parameters in FL biopsy samples taken at diagnosis. RESULTS:At the gene-expression level, FL is a highly uniform disease at the time of diagnosis, precluding the detection of sufficiently validated prognostic gene-expression profiles suitable for a clinical setting. Combinations of differentially expressed genes indicate that immunologic mechanisms play a differential role in the risk of early transformation. Using immunohistochemistry for specific cell populations, the spatial distribution to neoplastic follicles and the activation of CD4-positive T-helper cells (P = .002) and specifically T-helper 1 (P = .004) were shown to be highly discriminatory to predict early transformation. A role for functional modulation of follicular dendritic cells could also be supported (P = .04). Other cell populations, including CD68-positive macrophages and regulatory T cells, were not differentially present. CONCLUSION:These results support the identification of FL as an immunologically functional disease in which an interaction of the tumor cells and the functional composition of the microenvironment determines the clinical behavior. 10.1200/JCO.2006.06.1648
Molecular Dissection of CD8 T-Cell Dysfunction. Wang Chao,Singer Meromit,Anderson Ana C Trends in immunology Chronic viral infections and cancer often lead to the emergence of dysfunctional or 'exhausted' CD8 T cells, and the restoration of their functions is currently the focus of therapeutic interventions. In this review, we detail recent advances in the annotation of the gene modules and the epigenetic landscape associated with T-cell dysfunction. Together with analysis of single-cell transcriptomes, these findings have enabled a deeper and more precise understanding of the transcriptional mechanisms that induce and maintain the dysfunctional state and highlight the heterogeneity of CD8 T-cell phenotypes present in chronically inflamed tissue. We discuss the relevance of these findings for understanding the transcriptional and spatial regulation of dysfunctional T cells and for the design of therapeutics. 10.1016/j.it.2017.05.008
Cell-free DNA captures tumor heterogeneity and driver alterations in rapid autopsies with pre-treated metastatic cancer. Nature communications In patients with metastatic cancer, spatial heterogeneity of somatic alterations may lead to incomplete assessment of a cancer's mutational profile when analyzing a single tumor biopsy. In this study, we perform sequencing of cell-free DNA (cfDNA) and distinct metastatic tissue samples from ten rapid autopsy cases with pre-treated metastatic cancer. We show that levels of heterogeneity in genetic biomarkers vary between patients but that gene expression signatures representative of the tumor microenvironment are more consistent. Across nine patients with plasma samples available, we are able to detect 62/62 truncal and 47/121 non-truncal point mutations in cfDNA. We observe that mutation clonality in cfDNA is correlated with the number of metastatic lesions in which the mutation is detected and use this result to derive a clonality threshold to classify truncal and non-truncal driver alterations with reasonable specificity. In contrast, mutation truncality is more often incorrectly assigned when studying single tissue samples. Our results demonstrate the utility of a single cfDNA sample relative to that of single tissue samples when treating patients with metastatic cancer. 10.1038/s41467-021-23394-4
Multiplatform Molecular Profiling Reveals Epigenomic Intratumor Heterogeneity in Ependymoma. Liu S John,Magill Stephen T,Vasudevan Harish N,Hilz Stephanie,Villanueva-Meyer Javier E,Lastella Sydney,Daggubati Vikas,Spatz Jordan,Choudhury Abrar,Orr Brent A,Demaree Benjamin,Seo Kyounghee,Ferris Sean P,Abate Adam R,Oberheim Bush Nancy Ann,Bollen Andrew W,McDermott Michael W,Costello Joseph F,Raleigh David R Cell reports Ependymomas exist within distinct genetic subgroups, but the molecular diversity within individual ependymomas is unknown. We perform multiplatform molecular profiling of 6 spatially distinct samples from an ependymoma with C11orf95-RELA fusion. DNA methylation and RNA sequencing distinguish clusters of samples according to neuronal development gene expression programs that could also be delineated by differences in magnetic resonance blood perfusion. Exome sequencing and phylogenetic analysis reveal epigenomic intratumor heterogeneity and suggest that chromosomal structural alterations may precede accumulation of single-nucleotide variants during ependymoma tumorigenesis. In sum, these findings shed light on the oncogenesis and intratumor heterogeneity of ependymoma. 10.1016/j.celrep.2020.01.018
Inferring causal molecular networks: empirical assessment through a community-based effort. Hill Steven M,Heiser Laura M,Cokelaer Thomas,Unger Michael,Nesser Nicole K,Carlin Daniel E,Zhang Yang,Sokolov Artem,Paull Evan O,Wong Chris K,Graim Kiley,Bivol Adrian,Wang Haizhou,Zhu Fan,Afsari Bahman,Danilova Ludmila V,Favorov Alexander V,Lee Wai Shing,Taylor Dane,Hu Chenyue W,Long Byron L,Noren David P,Bisberg Alexander J, ,Mills Gordon B,Gray Joe W,Kellen Michael,Norman Thea,Friend Stephen,Qutub Amina A,Fertig Elana J,Guan Yuanfang,Song Mingzhou,Stuart Joshua M,Spellman Paul T,Koeppl Heinz,Stolovitzky Gustavo,Saez-Rodriguez Julio,Mukherjee Sach Nature methods It remains unclear whether causal, rather than merely correlational, relationships in molecular networks can be inferred in complex biological settings. Here we describe the HPN-DREAM network inference challenge, which focused on learning causal influences in signaling networks. We used phosphoprotein data from cancer cell lines as well as in silico data from a nonlinear dynamical model. Using the phosphoprotein data, we scored more than 2,000 networks submitted by challenge participants. The networks spanned 32 biological contexts and were scored in terms of causal validity with respect to unseen interventional data. A number of approaches were effective, and incorporating known biology was generally advantageous. Additional sub-challenges considered time-course prediction and visualization. Our results suggest that learning causal relationships may be feasible in complex settings such as disease states. Furthermore, our scoring approach provides a practical way to empirically assess inferred molecular networks in a causal sense. 10.1038/nmeth.3773
Distinct Patterns of Myeloid Cell Infiltration in Patients With hrHPV-Positive and hrHPV-Negative Penile Squamous Cell Carcinoma: The Importance of Assessing Myeloid Cell Densities Within the Spatial Context of the Tumor. Rafael Tynisha S,de Vries Hielke M,Ottenhof Sarah R,Hofland Ingrid,Broeks Annegien,de Jong Jeroen,Bekers Elise,Horenblas Simon,de Menezes Renée X,Jordanova Ekaterina S,Brouwer Oscar R Frontiers in immunology Comprehensive analysis of tumor infiltrating myeloid cells in the tumor microenvironment of penile squamous cell carcinoma (PSCC) is lacking. In this retrospective study, for the first time, PSCC resection specimens (N = 103) were annotated into the following compartments: intratumoral tumor (IT Tumor), intratumoral stroma (IT Stroma), peritumoral tumor (PT Tumor) and peritumoral stroma (PT Stroma) compartments. We then quantified CD14+, CD68+ and CD163+ myeloid cells within these compartments using an image analysis software and assessed their association with various clinical parameters, including high-risk human papillomavirus (hrHPV) status. In the total cohort, hrHPV status, grade of differentiation, age and tumor size were associated with myeloid cell densities. hrHPV tumors had higher infiltration rates of CD14+, CD68+ and CD163+ myeloid cells in the IT tumor compartment (p < 0.001, for all) compared to hrHPV tumors. Furthermore, when examining the association between compartment-specific infiltration and differentiation grade, increased myeloid cell densities in the IT tumor compartment were associated with a more advanced histological grade (p < 0.001, for all). This association remained significant when the hrHPV cohort (N = 60) was analyzed (CD14+ p = 0.001; CD68+ p < 0.001; CD163+ p = 0.004). Subgroup analysis in the hrHPV group (N = 43) showed that high infiltration rates of CD68+ and CD163+ cells in the PT tumor compartment were associated with lymph node (LN) metastasis (p = 0.031 and p = 0.026, respectively). Regarding the association between myeloid cell densities and disease-specific survival, the risk of death was found to decrease slightly as the number of myeloid cells in the IT tumor compartment increased (CD14+ p 0.04; CD68+ p = 0.05; CD163+ p = 0.02). However, after adjusting for hrHPV, no independent association between myeloid densities and disease-specific survival were found. Altogether, these findings demonstrate the importance of assessing myeloid cell densities within the spatial context of the tumor. Further studies are needed to unravel the specific phenotype of myeloid cells residing in the different compartments, their effect on clinical parameters and the impact of hrHPV on the recruitment of myeloid cell populations in PSCC. 10.3389/fimmu.2021.682030
Expression profiles of 39 HOX genes in normal human adult organs and anaplastic thyroid cancer cell lines by quantitative real-time RT-PCR system. Takahashi Yoko,Hamada Jun-ichi,Murakawa Katsuhiko,Takada Minoru,Tada Mitsuhiro,Nogami Ikuko,Hayashi Nobuyasu,Nakamori Shoji,Monden Morito,Miyamoto Masaki,Katoh Hiroyuki,Moriuchi Tetsuya Experimental cell research HOX genes are well known as master control genes in embryonic morphogenesis. We hypothesized that HOX genes give cells spatial information to maintain tissue- or organ-specificity in adult body and that the deregulated expression of HOX genes results in tumor development. We established a comprehensive analysis system to quantify expression of 39 human HOX genes based on the real-time reverse transcription PCR (RT-PCR) method. Analysis of 39 HOX genes of 20 normal adult organs by this system revealed that 5' HOX genes were expressed in organs in the caudal parts of the body, and that the more caudal regions the more numbers of HOX genes were expressed. It was also found that the expression patterns of HOX genes were more similar in the adjacent genes on the same cluster rather than in those belonging to the same paralogs. Compared with normal thyroid tissues, thyroid cancer cell lines showed the altered expression of some HOX genes, especially Abd-B homeobox family genes. Our results showed that HOX genes were organ-specifically expressed in adult body and that the deregulated expressions of Abd-B family genes were implicated in thyroid tumor development. 10.1016/j.yexcr.2003.09.024
Spatial clustering and common regulatory elements correlate with coordinated gene expression. Zhang Jingyu,Chen Hengyu,Li Ruoyan,Taft David A,Yao Guang,Bai Fan,Xing Jianhua PLoS computational biology Many cellular responses to surrounding cues require temporally concerted transcriptional regulation of multiple genes. In prokaryotic cells, a single-input-module motif with one transcription factor regulating multiple target genes can generate coordinated gene expression. In eukaryotic cells, transcriptional activity of a gene is affected by not only transcription factors but also the epigenetic modifications and three-dimensional chromosome structure of the gene. To examine how local gene environment and transcription factor regulation are coupled, we performed a combined analysis of time-course RNA-seq data of TGF-β treated MCF10A cells and related epigenomic and Hi-C data. Using Dynamic Regulatory Events Miner (DREM), we clustered differentially expressed genes based on gene expression profiles and associated transcription factors. Genes in each class have similar temporal gene expression patterns and share common transcription factors. Next, we defined a set of linear and radial distribution functions, as used in statistical physics, to measure the distributions of genes within a class both spatially and linearly along the genomic sequence. Remarkably, genes within the same class despite sometimes being separated by tens of million bases (Mb) along genomic sequence show a significantly higher tendency to be spatially close despite sometimes being separated by tens of Mb along the genomic sequence than those belonging to different classes do. Analyses extended to the process of mouse nervous system development arrived at similar conclusions. Future studies will be able to test whether this spatial organization of chromosomes contributes to concerted gene expression. 10.1371/journal.pcbi.1006786
Characterizing the molecular spatial and temporal field of injury in early-stage smoker non-small cell lung cancer patients after definitive surgery by expression profiling. Kadara Humam,Shen Li,Fujimoto Junya,Saintigny Pierre,Chow Chi-Wan,Lang Wenhua,Chu Zuoming,Garcia Melinda,Kabbout Mohamed,Fan You-Hong,Behrens Carmen,Liu Diane A,Mao Li,Lee J Jack,Gold Kathryn A,Wang Jing,Coombes Kevin R,Kim Edward S,Hong Waun Ki,Wistuba Ignacio I Cancer prevention research (Philadelphia, Pa.) Gene expression alterations in response to cigarette smoke have been characterized in normal-appearing bronchial epithelium of healthy smokers, and it has been suggested that adjacent histologically normal tissue displays tumor-associated molecular abnormalities. We sought to delineate the spatial and temporal molecular lung field of injury in smoker patients with early-stage non-small cell lung cancer (NSCLC; n = 19) who were accrued into a surveillance clinical trial for annual follow-up and bronchoscopies within 1 year after definitive surgery. Bronchial brushings and biopsies were obtained from six different sites in the lung at the time of inclusion in the study and at 12, 24, and 36 months after the first time point. Affymetrix Human Gene 1.0 ST arrays were used for whole-transcript expression profiling of airways (n = 391). Microarray analysis identified gene features (n = 1,165) that were nonuniform by site and differentially expressed between airways adjacent to tumors relative to more distant samples as well as those (n = 1,395) that were significantly altered with time up to 3 years. In addition, gene interaction networks mediated by phosphoinositide 3-kinase (PI3K) and extracellular signal-regulated kinase (ERK)1/2 were modulated in adjacent compared with contralateral airways and the latter network with time. Furthermore, phosphorylated AKT and ERK1/2 immunohistochemical expression were significantly increased with time (nuclear pAKT, P = 0.03; cytoplasmic pAKT, P < 0.0001; pERK1/2, P = 0.02) and elevated in adjacent compared with more distant airways (nuclear pAKT, P = 0.04; pERK1/2, P = 0.03). This study highlights spatial and temporal cancer-associated expression alterations in the molecular field of injury of patients with early-stage NSCLCs after definitive surgery that warrant further validation in independent studies. 10.1158/1940-6207.CAPR-12-0290
Whole-Body Single-Cell Sequencing Reveals Transcriptional Domains in the Annelid Larval Body. Molecular biology and evolution Animal bodies comprise diverse arrays of cells. To characterize cellular identities across an entire body, we have compared the transcriptomes of single cells randomly picked from dissociated whole larvae of the marine annelid Platynereis dumerilii. We identify five transcriptionally distinct groups of differentiated cells, each expressing a unique set of transcription factors and effector genes that implement cellular phenotypes. Spatial mapping of cells into a cellular expression atlas, and wholemount in situ hybridization of group-specific genes reveals spatially coherent transcriptional domains in the larval body, comprising, for example, apical sensory-neurosecretory cells versus neural/epidermal surface cells. These domains represent new, basic subdivisions of the annelid body based entirely on differential gene expression, and are composed of multiple, transcriptionally similar cell types. They do not represent clonal domains, as revealed by developmental lineage analysis. We propose that the transcriptional domains that subdivide the annelid larval body represent families of related cell types that have arisen by evolutionary diversification. Their possible evolutionary conservation makes them a promising tool for evo-devo research. 10.1093/molbev/msx336
Primer on the human genome. Tsai Kenneth Y,Tsao Hensin Journal of the American Academy of Dermatology UNLABELLED:The study of the expression patterns of many genes, or even the entire genome, is now routinely possible. Such powerful tools have enabled hypothesis-generating research at a scale never before possible. Moreover, spatially or temporally linked gene and protein expression, implying co-regulation and functional relatedness, has led to the identification of particular clusters of genes important for fundamental biologic processes, such as development and cancer. Not only is this expected to yield further mechanistic insights into disease processes, but perhaps most exciting, it will likely establish the foundation of predictive medicine, in which understanding of individual genomic signatures leads to the use of appropriately targeted therapy. LEARNING OBJECTIVE:At the conclusion of this learning activity, participants should be able to understand the fundamental tenets of molecular biology as they relate to the field of genomics. 10.1016/j.jaad.2006.10.985
Pro-inflammatory chemokine-chemokine receptor interactions within the Ewing sarcoma microenvironment determine CD8(+) T-lymphocyte infiltration and affect tumour progression. Berghuis Dagmar,Santos Susy J,Baelde Hans J,Taminiau Antonie Hm,Egeler R Maarten,Schilham Marco W,Hogendoorn Pancras Cw,Lankester Arjan C The Journal of pathology Ewing sarcoma is an aggressive round cell sarcoma with poor patient prognosis, particularly in cases of advanced-stage disease. Dynamic tumor-host immune interations within the tumor microenvironment may polarize in situ immune responses and shape tumor development and/or progression. To gain insight into the nature of tumour-host immune interactions within the Ewing sarcoma microenvironment, the presence and spatial distribution of infiltrating CD8(+) /CD4(+) T-lymphocytes were evaluated in therapy-naive Ewing sarcoma. Expression profiling of 40 different chemokines and several chemokine receptors was performed in therapy-naive tumours and cell lines by qPCR, immunohistochemistry, and flow cytometry. Considerable inter-tumour variation was observed regarding density, type, and distribution of infiltrating T-lymphocytes. Tumour-infiltrating T-cells contained significantly higher percentages of CD8(+) T-lymphocytes as compared to stroma-infiltrating cells, suggesting preferential migration of this T-cell type into tumour areas. Gene expression levels of several type 1-associated, pro-inflammatory chemokines (CXCR3- and CCR5-ligands CXCL9, CXCL10, and CCL5) correlated positively with infiltrating (CD8(+) ) T-lymphocyte numbers expressing corresponding chemokine receptors. Survival analyses demonstrated an impact of tumour-infiltrating, and not stroma-infiltrating, CD8(+) T-lymphocytes on tumour progression. At protein level, both tumour and stromal cells expressed the IFNγ-inducible chemokines CXCL9 and CXCL10. CCR5-ligand CCL5 was exclusively expressed by non-tumoural stromal/infiltrating cells. Together, our results indicate that an inflammatory immune microenvironment with high expression of type 1-associated chemokines may be critical for the recruitment of (CD8(+) ) T-lymphocytes expressing corresponding chemokine receptors. The observed impact of tumour-infiltrating (CD8(+) ) T-lymphocytes is consistent with a role for adaptive anti-tumour immunity in the prevention of Ewing sarcoma progression. Recognition of the merits and exploitation/induction of an inflammatory microenvironment may improve the efficacy of natural immune responses against, and (adoptive) immunotherapeutic approaches for, Ewing sarcoma. 10.1002/path.2819
A cellular hierarchy in melanoma uncouples growth and metastasis. Nature Although melanoma is notorious for its high degree of heterogeneity and plasticity, the origin and magnitude of cell-state diversity remains poorly understood. Equally, it is unclear whether growth and metastatic dissemination are supported by overlapping or distinct melanoma subpopulations. Here, by combining mouse genetics, single-cell and spatial transcriptomics, lineage tracing and quantitative modelling, we provide evidence of a hierarchical model of tumour growth that mirrors the cellular and molecular logic underlying the cell-fate specification and differentiation of the embryonic neural crest. We show that tumorigenic competence is associated with a spatially localized perivascular niche, a phenotype acquired through an intercellular communication pathway established by endothelial cells. Consistent with a model in which only a fraction of cells are fated to fuel growth, temporal single-cell tracing of a population of melanoma cells with a mesenchymal-like state revealed that these cells do not contribute to primary tumour growth but, instead, constitute a pool of metastatic initiating cells that switch cell identity while disseminating to secondary organs. Our data provide a spatially and temporally resolved map of the diversity and trajectories of melanoma cell states and suggest that the ability to support growth and metastasis are limited to distinct pools of cells. The observation that these phenotypic competencies can be dynamically acquired after exposure to specific niche signals warrant the development of therapeutic strategies that interfere with the cancer cell reprogramming activity of such microenvironmental cues. 10.1038/s41586-022-05242-7
Integrated Genomic Analysis Identifies Driver Genes and Cisplatin-Resistant Progenitor Phenotype in Pediatric Liver Cancer. Cancer discovery Pediatric liver cancers (PLC) comprise diverse diseases affecting infants, children, and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of patients with hepatoblastoma, all before three years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between "hepatocytic," "liver progenitor," and "mesenchymal" molecular subgroups of hepatoblastoma. We showed that during chemotherapy, "liver progenitor" cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies. SIGNIFICANCE: PLCs are deadly when they resist chemotherapy, with limited alternative treatment options. Using a multiomics approach, we identified PLC driver genes and the cellular phenotype at the origin of cisplatin resistance. We validated new treatments targeting these molecular features in cell lines and xenografts.. 10.1158/2159-8290.CD-20-1809
Molecular imaging of hepatocellular carcinoma. Gharib Ahmed M,Thomasson David,Li King C P Gastroenterology Molecular imaging may be defined as spatially localized and/or temporally resolved sensing of molecular and cellular processes in vivo. An imageable molecular event may be the result of the overexpression of a gene that produces a specific messenger RNA. The overexpressed protein could be an enzyme or could be incorporated into cell-surface transporters or receptors. Any step of this process is a potential target for molecular imaging. Current molecular imaging modalities include magnetic resonance imaging, nuclear imaging, ultrasound, and optical imaging. Nuclear medicine has been at the forefront of molecular imaging because of the relatively high sensitivity to detect nanomolar or picomolar quantities of the radiolabeled imaging probe. Imaging has had a central role in the diagnosis of hepatocellular carcinoma (HCC), which is considered the fifth most frequent malignancy worldwide. Nuclear imaging was one of the earlier modalities used for liver imaging. Traditional tracers included technetium 99m ( 99m Tc) sulfur colloid, gallium 67, and 99m Tc iminodiacetate acid analogues. Other less traditional probes include 99m Tc diethylenetriamine pentaacetic acid-galactosyl-human serum albumin for evaluation of functional liver volume and 99m Tc-labeled tetrofosmin and methoxisobutylisonitrile for detecting drug resistance. Fluorodeoxyglucose is the most widely used probe for positron emission tomography (PET) tumor imaging; however, carbon 11-labeled acetate appears to show improved sensitivity and specificity for HCC. Oxygen 15 PET imaging allows for the measurement of hepatic and tumor blood flow. Difficulties developing specific imaging methods for HCC are caused by the lack of obvious specific molecular targets, problems with drug delivery, and poor contrast-to-noise. No magic molecular imaging method exists today to accurately detect, characterize, and monitor HCC in vivo. 10.1053/j.gastro.2004.09.029
How should clinicians address intratumour heterogeneity in clear cell renal cell carcinoma? Current opinion in urology PURPOSE OF REVIEW:Despite the availability of multiple targeted therapies, the 5-year survival rate of patients with metastatic clear cell renal cell carcinoma (ccRCC) rarely exceeds 10%. Recent insights into the mutational landscape and evolutionary dynamics of ccRCC have offered up a plausible explanation for these outcomes. The purpose of this review is to link the research findings to potential changes in clinical practice. RECENT FINDINGS:Intratumour heterogeneity (ITH) dominates the evolutionary landscape in ccRCC at the genetic, transcriptomic and proteomic level. Spatial and temporal separation of tumour subclones within the primary tumour as well as between primary and metastatic sites has been demonstrated at single nucleotide resolution. In the cases analysed to date, approximately two-thirds of somatic mutations are not shared between multiple biopsies from the same primary tumour. Very few of the key disease-driving events are shared across all primary tumour regions (with the exception of VHL and loss of chromosome 3p), whereas the majority are restricted to one or more tumour regions (TP53, SETD2, BAP1, PTEN, mTOR, PIK3CA and KDM5C). SUMMARY:ITH must be considered in the management of ccRCC with respect to diagnostic procedures, prognostic and predictive biomarkers and drug development. 10.1097/MOU.0000000000000204
Benzo[a]pyrene Exposure Reduces Cell-Type Diversity and Stimulates Sex-Biased Damage Pathways in End Organs of Lupus-Prone Mice. International journal of molecular sciences Studies indicate that genetic factors only account for approximately thirty percent of all autoimmune diseases, while the rest of autoimmune pathogenesis is attributed to environmental factors including toxic chemicals. To understand if and how environmental pollutants trigger autoimmunity, we investigated the effect of benzo[a]pyrene (BaP) exposure on the development of autoimmune phenotypes in the lupus-prone MRL strain. The exposure of MRL mice to BaP over the course of 8 weeks before lupus onset resulted in total body weight loss in males, while marginal changes in anti-dsDNA levels occurred. Multi-organ analyses of BaP-treated and control MRL mice suggested that the kidney is a major organ directly affected by the metabolism of benzene-containing compounds, with increased expression of BaP-target genes including and . Intriguingly, spatial transcriptomic data showed that BaP caused a drastic reduction in cell-type diversity in both the kidneys and spleen of MRL mice. Further analysis of the molecular pathways affected suggested a sex-biased effect of BaP treatment, with the upregulated expression of angiogenesis genes in the lungs and an increased deposition of C3 in the kidneys of male mice. While SLE is more common in women, the disease is more severe in male patients, with an increased risk of disease progression to renal failure and lung cancer. Our results reveal sex-biased molecular pathways stimulated by BaP which may help explain the increased likelihood of end organ damage in males with lupus. 10.3390/ijms24076163
Mapping dynamic histone acetylation patterns to gene expression in nanog-depleted murine embryonic stem cells. Markowetz Florian,Mulder Klaas W,Airoldi Edoardo M,Lemischka Ihor R,Troyanskaya Olga G PLoS computational biology Embryonic stem cells (ESC) have the potential to self-renew indefinitely and to differentiate into any of the three germ layers. The molecular mechanisms for self-renewal, maintenance of pluripotency and lineage specification are poorly understood, but recent results point to a key role for epigenetic mechanisms. In this study, we focus on quantifying the impact of histone 3 acetylation (H3K9,14ac) on gene expression in murine embryonic stem cells. We analyze genome-wide histone acetylation patterns and gene expression profiles measured over the first five days of cell differentiation triggered by silencing Nanog, a key transcription factor in ESC regulation. We explore the temporal and spatial dynamics of histone acetylation data and its correlation with gene expression using supervised and unsupervised statistical models. On a genome-wide scale, changes in acetylation are significantly correlated to changes in mRNA expression and, surprisingly, this coherence increases over time. We quantify the predictive power of histone acetylation for gene expression changes in a balanced cross-validation procedure. In an in-depth study we focus on genes central to the regulatory network of Mouse ESC, including those identified in a recent genome-wide RNAi screen and in the PluriNet, a computationally derived stem cell signature. We find that compared to the rest of the genome, ESC-specific genes show significantly more acetylation signal and a much stronger decrease in acetylation over time, which is often not reflected in a concordant expression change. These results shed light on the complexity of the relationship between histone acetylation and gene expression and are a step forward to dissect the multilayer regulatory mechanisms that determine stem cell fate. 10.1371/journal.pcbi.1001034
The Tumor Immune Microenvironment Drives a Prognostic Relevance That Correlates with Bladder Cancer Subtypes. Pfannstiel Carolin,Strissel Pamela L,Chiappinelli Katherine B,Sikic Danijel,Wach Sven,Wirtz Ralph M,Wullweber Adrian,Taubert Helge,Breyer Johannes,Otto Wolfgang,Worst Thomas,Burger Maximilian,Wullich Bernd,Bolenz Christian,Fuhrich Nicole,Geppert Carol I,Weyerer Veronika,Stoehr Robert,Bertz Simone,Keck Bastian,Erlmeier Franziska,Erben Philipp,Hartmann Arndt,Strick Reiner,Eckstein Markus, , , , Cancer immunology research Muscle-invasive bladder cancer (MIBC) represents approximately two-thirds of invasive urothelial bladder cancers (UBC) and has high morbidity and mortality. Men are over 3-fold more frequently affected by UBC than women. Despite intensive efforts to improve patient treatment and outcome, two-thirds of patients with UBC will have a recurrence or disease progression within 5 years. We demonstrated that the quantity and spatial distribution of stromal tumor-infiltrating lymphocytes (sTIL) within the tumor immune microenvironment (TIME) predict stages of tumor inflammation, subtypes, and patient survival and correlate with expression of immune checkpoints in an analysis of 542 patients with MIBC. High sTILs indicated an inflamed subtype with an 80% 5-year DSS, and a lack of immune infiltrates identified an uninflamed subtype with a survival rate of less than 25%. A separate immune evading phenotype with upregulated immune checkpoints associated with poor survival. Within the TIME are tertiary lymphoid structures (TLS), which can mediate antitumor activity via immune cells. High TLS amounts and close tumor distance correlated significantly with an inflamed phenotype and favorable survival. The uninflamed and evasion phenotypes showed lowest TLS numbers, farthest tumor distances, and shortest survival. High inflammation also correlated with increased neoantigen load and mutational burden. Patients treated with adjuvant chemotherapy showed a favorable prognosis, which was dependent on high sTILs. Determination of sTILs and tumor subtypes may stratify therapy success and patient survival, and considering sTILs can easily be quantified using simple morphologic parameters, like hematoxylin and eosin, sTILs can be implemented for predicting patient survival in a routine manner. 10.1158/2326-6066.CIR-18-0758
Raman-guided subcellular pharmaco-metabolomics for metastatic melanoma cells. Du Jiajun,Su Yapeng,Qian Chenxi,Yuan Dan,Miao Kun,Lee Dongkwan,Ng Alphonsus H C,Wijker Reto S,Ribas Antoni,Levine Raphael D,Heath James R,Wei Lu Nature communications Non-invasively probing metabolites within single live cells is highly desired but challenging. Here we utilize Raman spectro-microscopy for spatial mapping of metabolites within single cells, with the specific goal of identifying druggable metabolic susceptibilities from a series of patient-derived melanoma cell lines. Each cell line represents a different characteristic level of cancer cell de-differentiation. First, with Raman spectroscopy, followed by stimulated Raman scattering (SRS) microscopy and transcriptomics analysis, we identify the fatty acid synthesis pathway as a druggable susceptibility for differentiated melanocytic cells. We then utilize hyperspectral-SRS imaging of intracellular lipid droplets to identify a previously unknown susceptibility of lipid mono-unsaturation within de-differentiated mesenchymal cells with innate resistance to BRAF inhibition. Drugging this target leads to cellular apoptosis accompanied by the formation of phase-separated intracellular membrane domains. The integration of subcellular Raman spectro-microscopy with lipidomics and transcriptomics suggests possible lipid regulatory mechanisms underlying this pharmacological treatment. Our method should provide a general approach in spatially-resolved single cell metabolomics studies. 10.1038/s41467-020-18376-x
Long Non-Coding RNAs, Cell Cycle, and Human Breast Cancer. Human gene therapy The long non-coding RNAs (lncRNAs) constitute an important class of the human transcriptome. The discovery of lncRNAs provided one of many unexpected results of the post-genomic era and uncovered a huge number of previously ignored transcriptional events. In recent years, lncRNAs are known to be linked with human diseases, with particular focus on cancer. Growing evidence has indicated that dysregulation of lncRNAs in breast cancer (BC) is strongly associated with the occurrence, development, and progress. Increasing numbers of lncRNAs have been found to interact with cell cycle progression and tumorigenesis in BC. The lncRNAs can exert their effect as a tumor suppressor or oncogene and regulate tumor development through direct or indirect regulation of cancer-related modulators and signaling pathways. What is more, lncRNAs are excellent candidates for promising therapeutic targets in BC due to the features of high tissue and cell-type specific expression. However, the underlying mechanisms of lncRNAs in BC still remain largely undefined. Here, we concisely summarize and sort out the current understanding of research progress in relationships of the roles for lncRNA in regulating the cell cycle. We also summarize the evidence for aberrant lncRNA expression in BC, and the potential for lncRNA to improve BC therapy is also discussed. Together, lncRNAs can be considered as exciting therapeutic candidates whose expression can be altered to impede BC progression. 10.1089/hum.2023.074
MicroRNA-182 plays an onco-miRNA role in cervical cancer. Tang Tao,Wong Hoi Kin,Gu Wenyi,Yu Mei-yung,To Ka-Fai,Wang Chi Chiu,Wong Yick Fu,Cheung Tak Hong,Chung Tony Kwok Hung,Choy Kwong Wai Gynecologic oncology OBJECTIVES:The purposes of this study were to identify aberrantly expressed miRNAs and investigate their pathogenic roles in cervical cancer. METHODS:miRNA expression was assessed in cervical cancer cell lines, micro-dissected normal cervical epithelium cells and primary cervical carcinoma by TaqMan RT-PCR. Spatial expression of miR-182 in cervical carcinoma and normal cervix was explored by in situ hybridization. HeLa xenograft mice model was used for evaluation of the effect on tumor growth of miR-182 inhibitor. Western blot, flow cytometry and gene expression analysis were used for identification of the functional role of miR-182 in HeLa cells. RESULTS:Two up-regulated (miR-182 and -183) and nine down-regulated (miR-211, 145, 223, 150, 142-5p, 328, 195, 199b, 142-3p) microRNAs were consistently identified in cervical cancer cell lines. Further investigation confirmed the most up-regulated miRNA (miR-182) was significantly elevated in primary cervical carcinoma and discovered a significant correlation between the increased expression of miR-182 and advanced stages of cervical cancer. In HeLa xenograft mouse model, we demonstrated that inhibition of the miR-182 could exert the effect of tumor growth regression. Western blot, flow cytometry and pathway analysis for the HeLa cells with miR-182 over/down-expression in vitro showed that miR-182 was involved in apoptosis and cell cycle pathways, it also associated with the regulation of FOXO1. CONCLUSIONS:Our findings indicated that miR-182 plays an onco-miRNA role in cervical cancer and its alteration is associated with cervical cancer pathogenesis by disrupting cell proliferation. 10.1016/j.ygyno.2012.12.043
Multiregion ultra-deep sequencing reveals early intermixing and variable levels of intratumoral heterogeneity in colorectal cancer. Suzuki Yuka,Ng Sarah Boonhsi,Chua Clarinda,Leow Wei Qiang,Chng Jermain,Liu Shi Yang,Ramnarayanan Kalpana,Gan Anna,Ho Dan Liang,Ten Rachel,Su Yan,Lezhava Alexandar,Lai Jiunn Herng,Koh Dennis,Lim Kiat Hon,Tan Patrick,Rozen Steven G,Tan Iain Beehuat Molecular oncology Intratumor heterogeneity (ITH) contributes to cancer progression and chemoresistance. We sought to comprehensively describe ITH of somatic mutations, copy number, and transcriptomic alterations involving clinically and biologically relevant gene pathways in colorectal cancer (CRC). We performed multiregion, high-depth (384× on average) sequencing of 799 cancer-associated genes in 24 spatially separated primary tumor and nonmalignant tissues from four treatment-naïve CRC patients. We then used ultra-deep sequencing (17 075× on average) to accurately verify the presence or absence of identified somatic mutations in each sector. We also digitally measured gene expression and copy number alterations using NanoString assays. We identified the subclonal point mutations and determined the mutational timing and phylogenetic relationships among spatially separated sectors of each tumor. Truncal mutations, those shared by all sectors in the tumor, affected the well-described driver genes such as APC, TP53, and KRAS. With sequencing at 17 075×, we found that mutations first detected at a sequencing depth of 384× were in fact more widely shared among sectors than originally assessed. Interestingly, ultra-deep sequencing also revealed some mutations that were present in all spatially dispersed sectors, but at subclonal levels. Ultra-high-depth validation sequencing, copy number analysis, and gene expression profiling provided a comprehensive and accurate genomic landscape of spatial heterogeneity in CRC. Ultra-deep sequencing allowed more sensitive detection of somatic mutations and a more accurate assessment of ITH. By detecting the subclonal mutations with ultra-deep sequencing, we traced the genomic histories of each tumor and the relative timing of mutational events. We found evidence of early mixing, in which the subclonal ancestral mutations intermixed across the sectors before the acquisition of subsequent nontruncal mutations. Our findings also indicate that different CRC patients display markedly variable ITH, suggesting that each patient's tumor possesses a unique genomic history and spatial organization. 10.1002/1878-0261.12012
Early gene expression analysis in 9L orthotopic tumor-bearing rats identifies immune modulation in molecular response to synchrotron microbeam radiation therapy. Bouchet Audrey,Sakakini Nathalie,El Atifi Michèle,Le Clec'h Céline,Brauer Elke,Moisan Anaïck,Deman Pierre,Rihet Pascal,Le Duc Géraldine,Pelletier Laurent PloS one Synchrotron Microbeam Radiation Therapy (MRT) relies on the spatial fractionation of the synchrotron photon beam into parallel micro-beams applying several hundred of grays in their paths. Several works have reported the therapeutic interest of the radiotherapy modality at preclinical level, but biological mechanisms responsible for the described efficacy are not fully understood to date. The aim of this study was to identify the early transcriptomic responses of normal brain and glioma tissue in rats after MRT irradiation (400Gy). The transcriptomic analysis of similarly irradiated normal brain and tumor tissues was performed 6 hours after irradiation of 9 L orthotopically tumor-bearing rats. Pangenomic analysis revealed 1012 overexpressed and 497 repressed genes in the irradiated contralateral normal tissue and 344 induced and 210 repressed genes in tumor tissue. These genes were grouped in a total of 135 canonical pathways. More than half were common to both tissues with a predominance for immunity or inflammation (64 and 67% of genes for normal and tumor tissues, respectively). Several pathways involving HMGB1, toll-like receptors, C-type lectins and CD36 may serve as a link between biochemical changes triggered by irradiation and inflammation and immunological challenge. Most immune cell populations were involved: macrophages, dendritic cells, natural killer, T and B lymphocytes. Among them, our results highlighted the involvement of Th17 cell population, recently described in tumor. The immune response was regulated by a large network of mediators comprising growth factors, cytokines, lymphokines. In conclusion, early response to MRT is mainly based on inflammation and immunity which appear therefore as major contributors to MRT efficacy. 10.1371/journal.pone.0081874
A human fetal lung cell atlas uncovers proximal-distal gradients of differentiation and key regulators of epithelial fates. Cell We present a multiomic cell atlas of human lung development that combines single-cell RNA and ATAC sequencing, high-throughput spatial transcriptomics, and single-cell imaging. Coupling single-cell methods with spatial analysis has allowed a comprehensive cellular survey of the epithelial, mesenchymal, endothelial, and erythrocyte/leukocyte compartments from 5-22 post-conception weeks. We identify previously uncharacterized cell states in all compartments. These include developmental-specific secretory progenitors and a subtype of neuroendocrine cell related to human small cell lung cancer. Our datasets are available through our web interface (https://lungcellatlas.org). To illustrate its general utility, we use our cell atlas to generate predictions about cell-cell signaling and transcription factor hierarchies which we rigorously test using organoid models. 10.1016/j.cell.2022.11.005
RNAi microarray analysis in cultured mammalian cells. Mousses Spyro,Caplen Natasha J,Cornelison Robert,Weaver Don,Basik Mark,Hautaniemi Sampsa,Elkahloun Abdel G,Lotufo Roberto A,Choudary Ashish,Dougherty Edward R,Suh Ed,Kallioniemi Olli Genome research RNA interference (RNAi) mediated by small interfering RNAs (siRNAs) is a powerful new tool for analyzing gene knockdown phenotypes in living mammalian cells. To facilitate large-scale, high-throughput functional genomics studies using RNAi, we have developed a microarray-based technology for highly parallel analysis. Specifically, siRNAs in a transfection matrix were first arrayed on glass slides, overlaid with a monolayer of adherent cells, incubated to allow reverse transfection, and assessed for the effects of gene silencing by digital image analysis at a single cell level. Validation experiments with HeLa cells stably expressing GFP showed spatially confined, sequence-specific, time- and dose-dependent inhibition of green fluorescence for those cells growing directly on microspots containing siRNA targeting the GFP sequence. Microarray-based siRNA transfections analyzed with a custom-made quantitative image analysis system produced results that were identical to those from traditional well-based transfection, quantified by flow cytometry. Finally, to integrate experimental details, image analysis, data display, and data archiving, we developed a prototype information management system for high-throughput cell-based analyses. In summary, this RNAi microarray platform, together with ongoing efforts to develop large-scale human siRNA libraries, should facilitate genomic-scale cell-based analyses of gene function. 10.1101/gr.1478703
Spatial transcriptomic profiling reveals the pathophysiology of early-stage hidradenitis suppurativa. The British journal of dermatology 10.1093/bjd/ljad274
A hypergraph-based learning algorithm for classifying gene expression and arrayCGH data with prior knowledge. Tian Ze,Hwang TaeHyun,Kuang Rui Bioinformatics (Oxford, England) MOTIVATION:Incorporating biological prior knowledge into predictive models is a challenging data integration problem in analyzing high-dimensional genomic data. We introduce a hypergraph-based semi-supervised learning algorithm called HyperPrior to classify gene expression and array-based comparative genomic hybridization (arrayCGH) data using biological knowledge as constraints on graph-based learning. HyperPrior is a robust two-step iterative method that alternatively finds the optimal labeling of the samples and the optimal weighting of the features, guided by constraints encoding prior knowledge. The prior knowledge for analyzing gene expression data is that cancer-related genes tend to interact with each other in a protein-protein interaction network. Similarly, the prior knowledge for analyzing arrayCGH data is that probes that are spatially nearby in their layout along the chromosomes tend to be involved in the same amplification or deletion event. Based on the prior knowledge, HyperPrior imposes a consistent weighting of the correlated genomic features in graph-based learning. RESULTS:We applied HyperPrior to test two arrayCGH datasets and two gene expression datasets for both cancer classification and biomarker identification. On all the datasets, HyperPrior achieved competitive classification performance, compared with SVMs and the other baselines utilizing the same prior knowledge. HyperPrior also identified several discriminative regions on chromosomes and discriminative subnetworks in the PPI, both of which contain cancer-related genomic elements. Our results suggest that HyperPrior is promising in utilizing biological prior knowledge to achieve better classification performance and more biologically interpretable findings in gene expression and arrayCGH data. AVAILABILITY:http://compbio.cs.umn.edu/HyperPrior CONTACT:kuang@cs.umn.edu SUPPLEMENTARY INFORMATION:Supplementary data are available at bioinformatics online. 10.1093/bioinformatics/btp467
Cancer-Associated Fibroblast-Like Tumor Cells Remodel the Ewing Sarcoma Tumor Microenvironment. Clinical cancer research : an official journal of the American Association for Cancer Research PURPOSE:Despite limited genetic and histologic heterogeneity, Ewing sarcoma (EwS) tumor cells are transcriptionally heterogeneous and display varying degrees of mesenchymal lineage specification in vitro. In this study, we investigated if and how transcriptional heterogeneity of EwS cells contributes to heterogeneity of tumor phenotypes in vivo. EXPERIMENTAL DESIGN:Single-cell proteogenomic-sequencing of EwS cell lines was performed and integrated with patient tumor transcriptomic data. Cell subpopulations were isolated by FACS for assessment of gene expression and phenotype. Digital spatial profiling and human whole transcriptome analysis interrogated transcriptomic heterogeneity in EwS xenografts. Tumor cell subpopulations and matrix protein deposition were evaluated in xenografts and patient tumors using multiplex immunofluorescence staining. RESULTS:We identified CD73 as a biomarker of highly mesenchymal EwS cell subpopulations in tumor models and patient biopsies. CD73+ tumor cells displayed distinct transcriptional and phenotypic properties, including selective upregulation of genes that are repressed by EWS::FLI1, and increased migratory potential. CD73+ cells were distinguished in vitro and in vivo by increased expression of matrisomal genes and abundant deposition of extracellular matrix (ECM) proteins. In epithelial-derived malignancies, ECM is largely deposited by cancer-associated fibroblasts (CAF), and we thus labeled CD73+ EwS cells, CAF-like tumor cells. Marked heterogeneity of CD73+ EwS cell frequency and distribution was detected in tumors in situ, and CAF-like tumor cells and associated ECM were observed in peri-necrotic regions and invasive foci. CONCLUSIONS:EwS tumor cells can adopt CAF-like properties, and these distinct cell subpopulations contribute to tumor heterogeneity by remodeling the tumor microenvironment. See related commentary by Kuo and Amatruda, p. 5002. 10.1158/1078-0432.CCR-23-1111
Somatic pairing of chromosome 19 in renal oncocytoma is associated with deregulated EGLN2-mediated [corrected] oxygen-sensing response. Koeman Julie M,Russell Ryan C,Tan Min-Han,Petillo David,Westphal Michael,Koelzer Katherine,Metcalf Julie L,Zhang Zhongfa,Matsuda Daisuke,Dykema Karl J,Houseman Heather L,Kort Eric J,Furge Laura L,Kahnoski Richard J,Richard Stéphane,Vieillefond Annick,Swiatek Pamela J,Teh Bin Tean,Ohh Michael,Furge Kyle A PLoS genetics Chromosomal abnormalities, such as structural and numerical abnormalities, are a common occurrence in cancer. The close association of homologous chromosomes during interphase, a phenomenon termed somatic chromosome pairing, has been observed in cancerous cells, but the functional consequences of somatic pairing have not been established. Gene expression profiling studies revealed that somatic pairing of chromosome 19 is a recurrent chromosomal abnormality in renal oncocytoma, a neoplasia of the adult kidney. Somatic pairing was associated with significant disruption of gene expression within the paired regions and resulted in the deregulation of the prolyl-hydroxylase EGLN2 [corrected] a key protein that regulates the oxygen-dependent degradation of hypoxia-inducible factor (HIF). Overexpression of EGLN2 [corrected] in renal oncocytoma increased ubiquitin-mediated destruction of HIF and concomitantly suppressed the expression of several HIF-target genes, including the pro-death BNIP3L gene. The transcriptional changes that are associated with somatic pairing of chromosome 19 mimic the transcriptional changes that occur following DNA amplification. Therefore, in addition to numerical and structural chromosomal abnormalities, alterations in chromosomal spatial dynamics should be considered as genomic events that are associated with tumorigenesis. The identification of EGLN2 as a significantly deregulated gene that maps within the paired chromosome region directly implicates defects in the oxygen-sensing network to the biology of renal oncocytoma. 10.1371/journal.pgen.1000176
Immunogenomics of Colorectal Cancer Response to Checkpoint Blockade: Analysis of the KEYNOTE 177 Trial and Validation Cohorts. Gastroenterology BACKGROUND & AIMS:Colorectal cancer (CRC) shows variable response to immune checkpoint blockade, which can only partially be explained by high tumor mutational burden (TMB). We conducted an integrated study of the cancer tissue and associated tumor microenvironment (TME) from patients treated with pembrolizumab (KEYNOTE 177 clinical trial) or nivolumab to dissect the cellular and molecular determinants of response to anti- programmed cell death 1 (PD1) immunotherapy. METHODS:We selected multiple regions per tumor showing variable T-cell infiltration for a total of 738 regions from 29 patients, divided into discovery and validation cohorts. We performed multiregional whole-exome and RNA sequencing of the tumor cells and integrated these with T-cell receptor sequencing, high-dimensional imaging mass cytometry, detection of programmed death-ligand 1 (PDL1) interaction in situ, multiplexed immunofluorescence, and computational spatial analysis of the TME. RESULTS:In hypermutated CRCs, response to anti-PD1 immunotherapy was not associated with TMB but with high clonality of immunogenic mutations, clonally expanded T cells, low activation of Wnt signaling, deregulation of the interferon gamma pathway, and active immune escape mechanisms. Responsive hypermutated CRCs were also rich in cytotoxic and proliferating PD1CD8 T cells interacting with PDL1 antigen-presenting macrophages. CONCLUSIONS:Our study clarified the limits of TMB as a predictor of response of CRC to anti-PD1 immunotherapy. It identified a population of antigen-presenting macrophages interacting with CD8 T cells that consistently segregate with response. We therefore concluded that anti-PD1 agents release the PD1-PDL1 interaction between CD8 T cells and macrophages to promote cytotoxic antitumor activity. 10.1053/j.gastro.2021.06.064
Three-dimensional spatial transcriptomics uncovers cell type localizations in the human rheumatoid arthritis synovium. Communications biology The inflamed rheumatic joint is a highly heterogeneous and complex tissue with dynamic recruitment and expansion of multiple cell types that interact in multifaceted ways within a localized area. Rheumatoid arthritis synovium has primarily been studied either by immunostaining or by molecular profiling after tissue homogenization. Here, we use Spatial Transcriptomics, where tissue-resident RNA is spatially labeled in situ with barcodes in a transcriptome-wide fashion, to study local tissue interactions at the site of chronic synovial inflammation. We report comprehensive spatial RNA-Seq data coupled to cell type-specific localization patterns at and around organized structures of infiltrating leukocyte cells in the synovium. Combining morphological features and high-throughput spatially resolved transcriptomics may be able to provide higher statistical power and more insights into monitoring disease severity and treatment-specific responses in seropositive and seronegative rheumatoid arthritis. 10.1038/s42003-022-03050-3
Heterogeneity of tumor-infiltrating lymphocytes ascribed to local immune status rather than neoantigens by multi-omics analysis of glioblastoma multiforme. Feng Lin,Qian Haipeng,Yu Xuexin,Liu Kan,Xiao Ting,Zhang Chengli,Kuang Manchao,Cheng Shujun,Li Xueji,Wan Jinghai,Zhang Kaitai Scientific reports Hypothetically, intratumoral genomic heterogeneity has the potential to foster tumor-infiltrating lymphocyte (TIL) diversity; however, no study has directly tested this hypothesis by simultaneously investigating somatic mutations, TIL diversity, and immune response activity. Thus, we performed whole-exome sequencing, immune repertoire sequencing and gene expression on ten spatially separated tumor samples obtained from two tumor masses excised from a glioblastoma multiforme (GBM) patient, and we included peripheral blood as control. We found that although the multi-region samples from one tumor shared more common mutations than those from different tumors, the TIL populations did not. TIL repertoire diversity did not significantly correlate with the number of non-synonymous mutations; however, TIL diversity was highly correlated with local immune activity, as the pathways were all immune-related pathways that highly positive correlated with local TIL diversity. Twenty-three genes with expression largely unaffected by the intratumor heterogeneity were extracted from these pathways. Fifty GBM patients were stratified into two clusters by the expression of these genes with significant difference in prognosis. This finding was validated by The Cancer Genome Atlas (TCGA) GBM dataset, which indicated that despite the heterogeneity of intra-tumor immune status, the overall level of the immune response in GBM could be connected with prognosis. 10.1038/s41598-017-05538-z
In situ protein expression in tumour spheres: development of an immunostaining protocol for confocal microscopy. Weiswald Louis-Bastien,Guinebretière Jean-Marc,Richon Sophie,Bellet Dominique,Saubaméa Bruno,Dangles-Marie Virginie BMC cancer BACKGROUND:Multicellular tumour sphere models have been shown to closely mimic phenotype characteristics of in vivo solid tumours, or to allow in vitro propagation of cancer stem cells (CSCs). CSCs are usually characterized by the expression of specific membrane markers using flow cytometry (FC) after enzymatic dissociation. Consequently, the spatial location of positive cells within spheres is not documented. Confocal microscopy is the best technique for the imaging of thick biological specimens after multi-labelling but suffers from poor antibody penetration. Thus, we describe here a new protocol for in situ confocal imaging of protein expression in intact spheroids. METHODS:Protein expression in whole spheroids (150 mum in diameter) from two human colon cancer cell lines, HT29 and CT320X6, has been investigated with confocal immunostaining, then compared with profiles obtained through paraffin immunohistochemistry (pIHC) and FC. Target antigens, relevant for colon cancer and with different expression patterns, have been studied. RESULTS:We first demonstrate that our procedure overcomes the well-known problem of antibody penetration in compact structures by performing immunostaining of EpCAM, a membrane protein expressed by all cells within our spheroids. EpCAM expression is detected in all cells, even the deepest ones. Likewise, antibody access is confirmed with CK20 and CD44 immunostaining. Confocal imaging shows that 100% of cells express beta-catenin, mainly present in the plasma membrane with also cytoplasmic and nuclear staining, in agreement with FC and pIHC data. pIHC and confocal imaging show similar CA 19-9 cytoplasmic and membranar expression profile in a cell subpopulation. CA 19-9+ cell count confirms confocal imaging as a highly sensitive method (75%, 62% and 51%, for FC, confocal imaging and pIHC, respectively). Finally, confocal imaging reveals that the weak expression of CD133, a putative colon CSC marker, is restricted to the luminal cell surface of colorectal cancer acini, with CD133+ cellular debris into glandular lumina. CONCLUSION:The present protocol enables in situ visualization of protein expression in compact three-dimensional models by whole mount confocal imaging, allowing the accurate localization and quantification of cells expressing specific markers. It should prove useful to study rare events like CSCs within tumour spheres. 10.1186/1471-2407-10-106
Integrative Genomic Analyses Identify LncRNA Regulatory Networks across Pediatric Leukemias and Solid Tumors. Cancer research Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies. SIGNIFICANCE:Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks. 10.1158/0008-5472.CAN-22-3186
DNA Methylation Landscape Reflects the Spatial Organization of Chromatin in Different Cells. Zhang Ling,Xie Wen Jun,Liu Sirui,Meng Luming,Gu Chan,Gao Yi Qin Biophysical journal The relationship between DNA methylation and chromatin structure is still largely unknown. By analyzing a large set of published sequencing data, we observed a long-range power law correlation of DNA methylation with cell class-specific scaling exponents in the range of tens of kilobases. We showed that such cell class-specific scaling exponents are caused by different patchiness of DNA methylation in different cells. By modeling the chromatin structure using high-resolution chromosome conformation capture data and mapping the methylation level onto the modeled structure, we demonstrated that the patchiness of DNA methylation is related to chromatin structure. The scaling exponents of the power law correlation are thus a display of the spatial organization of chromatin. Besides the long-range correlation, we also showed that the local correlation of DNA methylation is associated with nucleosome positioning. The local correlation of partially methylated domains is different from that of nonpartially methylated domains, suggesting that their chromatin structures differ at the scale of several hundred base pairs (covering a few nucleosomes). Our study provides a novel, to our knowledge, view of the spatial organization of chromatin structure from a perspective of DNA methylation, in which both long-range and local correlations of DNA methylation along the genome reflect the spatial organization of chromatin. 10.1016/j.bpj.2017.08.019
Spatial genomic heterogeneity in diffuse intrinsic pontine and midline high-grade glioma: implications for diagnostic biopsy and targeted therapeutics. Hoffman Lindsey M,DeWire Mariko,Ryall Scott,Buczkowicz Pawel,Leach James,Miles Lili,Ramani Arun,Brudno Michael,Kumar Shiva Senthil,Drissi Rachid,Dexheimer Phillip,Salloum Ralph,Chow Lionel,Hummel Trent,Stevenson Charles,Lu Q Richard,Jones Blaise,Witte David,Aronow Bruce,Hawkins Cynthia E,Fouladi Maryam Acta neuropathologica communications INTRODUCTION:Diffuse intrinsic pontine glioma (DIPG) and midline high-grade glioma (mHGG) are lethal childhood brain tumors. Spatial genomic heterogeneity has been well-described in adult HGG but has not been comprehensively characterized in pediatric HGG. We performed whole exome sequencing on 38-matched primary, contiguous, and metastatic tumor sites from eight children with DIPG (n = 7) or mHGG (n = 1) collected using a unique MRI-guided autopsy protocol. Validation was performed using Sanger sequencing, Droplet Digital polymerase-chain reaction, immunohistochemistry, and fluorescent in-situ hybridization. RESULTS:Median age at diagnosis was 6.1 years (range: 2.9-23.3 years). Median overall survival was 13.2 months (range: 11.2-32.2 months). Contiguous tumor infiltration and distant metastases were observed in seven and six patients, respectively, including leptomeningeal dissemination in three DIPGs. Histopathological heterogeneity was evident in seven patients, including intra-pontine heterogeneity in two DIPGs, ranging from World Health Organization grade II to IV astrocytoma. We found conservation of heterozygous K27M mutations in H3F3A (n = 4) or HIST1H3B (n = 3) across all primary, contiguous, and metastatic tumor sites in all DIPGs. ACVR1 (n = 2), PIK3CA (n = 2), FGFR1 (n = 2), and MET (n = 1) were also intra-tumorally conserved. ACVR1 was co-mutated with HIST1H3B (n = 2). In contrast, PDGFRA amplification and mutation were spatially heterogeneous, as were mutations in BCOR (n = 1), ATRX (n = 2), and MYC (n = 1). TP53 aberrations (n = 3 patients) varied by type and location between primary and metastatic tumors sites but were intra-tumorally conserved. CONCLUSION:Spatial conservation of prognostically-relevant and therapeutically-targetable somatic mutations in DIPG and mHGG contrasts the significant heterogeneity of driver mutations seen in adult HGG and supports uniform implementation of diagnostic biopsy in DIPG and mHGG to classify molecular risk groups and guide therapeutic strategy. 10.1186/s40478-015-0269-0
Radiogenomics of hepatocellular carcinoma: multiregion analysis-based identification of prognostic imaging biomarkers by integrating gene data-a preliminary study. Xia Wei,Chen Ying,Zhang Rui,Yan Zhuangzhi,Zhou Xiaobo,Zhang Bo,Gao Xin Physics in medicine and biology Our objective was to identify prognostic imaging biomarkers for hepatocellular carcinoma in contrast-enhanced computed tomography (CECT) with biological interpretations by associating imaging features and gene modules. We retrospectively analyzed 371 patients who had gene expression profiles. For the 38 patients with CECT imaging data, automatic intra-tumor partitioning was performed, resulting in three spatially distinct subregions. We extracted a total of 37 quantitative imaging features describing intensity, geometry, and texture from each subregion. Imaging features were selected after robustness and redundancy analysis. Gene modules acquired from clustering were chosen for their prognostic significance. By constructing an association map between imaging features and gene modules with Spearman rank correlations, the imaging features that significantly correlated with gene modules were obtained. These features were evaluated with Cox's proportional hazard models and Kaplan-Meier estimates to determine their prognostic capabilities for overall survival (OS). Eight imaging features were significantly correlated with prognostic gene modules, and two of them were associated with OS. Among these, the geometry feature volume fraction of the subregion, which was significantly correlated with all prognostic gene modules representing cancer-related interpretation, was predictive of OS (Cox p  =  0.022, hazard ratio  =  0.24). The texture feature cluster prominence in the subregion, which was correlated with the prognostic gene module representing lipid metabolism and complement activation, also had the ability to predict OS (Cox p  =  0.021, hazard ratio  =  0.17). Imaging features depicting the volume fraction and textural heterogeneity in subregions have the potential to be predictors of OS with interpretable biological meaning. 10.1088/1361-6560/aaa609
Comparison of the Agilent, ROMA/NimbleGen and Illumina platforms for classification of copy number alterations in human breast tumors. Baumbusch L O,Aarøe J,Johansen F E,Hicks J,Sun H,Bruhn L,Gunderson K,Naume B,Kristensen V N,Liestøl K,Børresen-Dale A-L,Lingjaerde O C BMC genomics BACKGROUND:Microarray Comparative Genomic Hybridization (array CGH) provides a means to examine DNA copy number aberrations. Various platforms, brands and underlying technologies are available, facing the user with many choices regarding platform sensitivity and number, localization, and density distribution of probes. RESULTS:We evaluate three different platforms presenting different nature and arrangement of the probes: The Agilent Human Genome CGH Microarray 44 k, the ROMA/NimbleGen Representational Oligonucleotide Microarray 82 k, and the Illumina Human-1 Genotyping 109 k BeadChip, with Agilent being gene oriented, ROMA/NimbleGen being genome oriented, and Illumina being genotyping oriented. We investigated copy number changes in 20 human breast tumor samples representing different gene expression subclasses, using a suite of graphical and statistical methods designed to work across platforms. Despite substantial differences in the composition and spatial distribution of probes, the comparison revealed high overall concordance. Notably however, some short amplifications and deletions of potential biological importance were not detected by all platforms. Both correlation and cluster analysis indicate a somewhat higher similarity between ROMA/NimbleGen and Illumina than between Agilent and the other two platforms. The programs developed for the analysis are available from http://www.ifi.uio.no/bioinf/Projects/. CONCLUSION:We conclude that platforms based on different technology principles reveal similar aberration patterns, although we observed some unique amplification or deletion peaks at various locations, only detected by one of the platforms. The correct platform choice for a particular study is dependent on whether the appointed research intention is gene, genome, or genotype oriented. 10.1186/1471-2164-9-379
Whole-genome microarrays of fission yeast: characteristics, accuracy, reproducibility, and processing of array data. Lyne Rachel,Burns Gavin,Mata Juan,Penkett Chris J,Rustici Gabriella,Chen Dongrong,Langford Cordelia,Vetrie David,Bähler Jürg BMC genomics BACKGROUND:The genome of the fission yeast Schizosaccharomyces pombe has recently been sequenced, setting the stage for the post-genomic era of this increasingly popular model organism. We have built fission yeast microarrays, optimised protocols to improve array performance, and carried out experiments to assess various characteristics of microarrays. RESULTS:We designed PCR primers to amplify specific probes (180-500 bp) for all known and predicted fission yeast genes, which are printed in duplicate onto separate regions of glass slides together with control elements (approximately 13,000 spots/slide). Fluorescence signal intensities depended on the size and intragenic position of the array elements, whereas the signal ratios were largely independent of element properties. Only the coding strand is covalently linked to the slides, and our array elements can discriminate transcriptional direction. The microarrays can distinguish sequences with up to 70% identity, above which cross-hybridisation contributes to the signal intensity. We tested the accuracy of signal ratios and measured the reproducibility of array data caused by biological and technical factors. Because the technical variability is lower, it is best to use samples prepared from independent biological experiments to obtain repeated measurements with swapping of fluorochromes to prevent dye bias. We also developed a script that discards unreliable data and performs a normalization to correct spatial artefacts. CONCLUSIONS:This paper provides data for several microarray properties that are rarely measured. The results define critical parameters for microarray design and experiments and provide a framework to optimise and interpret array data. Our arrays give reproducible and accurate expression ratios with high sensitivity. The scripts for primer design and initial data processing as well as primer sequences and detailed protocols are available from our website. 10.1186/1471-2164-4-27
Functional genomics guided with MR imaging: mouse tumor model study. Guccione Samira,Yang Yi-Shan,Shi Gongyi,Lee Daniel Y,Li King C P,Bednarski Mark D Radiology To gain a better understanding of gene expression patterns in tumors, the authors used contrast material-enhanced magnetic resonance (MR) imaging to noninvasively characterize regions within the same tumor to provide a correlate for genomic analysis. Gene expression profiles of samples from a mouse tumor model obtained from contrast-enhanced and nonenhanced regions within the same tumor were compared with MR imaging and functional genomics. From these samples, 11000 genes were analyzed: 10 genes were up-regulated in the contrast-enhanced areas, and one gene was up-regulated in the nonenhanced regions. Several of these genes encode extracellular matrix proteins. Findings in this study demonstrate that MR imaging can serve as a powerful noninvasive tool for characterizing different regions of tumors to guide genomic analysis with high spatial and temporal resolution. 10.1148/radiol.2282020907
Single-Cell and Spatial Transcriptomic Analysis of Human Skin Delineates Intercellular Communication and Pathogenic Cells. The Journal of investigative dermatology Epidermal homeostasis is governed by a balance between keratinocyte proliferation and differentiation with contributions from cell-cell interactions, but conserved or divergent mechanisms governing this equilibrium across species and how an imbalance contributes to skin disease are largely undefined. To address these questions, human skin single-cell RNA sequencing and spatial transcriptomics data were integrated and compared with mouse skin data. Human skin cell-type annotation was improved using matched spatial transcriptomics data, highlighting the importance of spatial context in cell-type identity, and spatial transcriptomics refined cellular communication inference. In cross-species analyses, we identified a human spinous keratinocyte subpopulation that exhibited proliferative capacity and a heavy metal processing signature, which was absent in mouse and may account for species differences in epidermal thickness. This human subpopulation was expanded in psoriasis and zinc-deficiency dermatitis, attesting to disease relevance and suggesting a paradigm of subpopulation dysfunction as a hallmark of the disease. To assess additional potential subpopulation drivers of skin diseases, we performed cell-of-origin enrichment analysis within genodermatoses, nominating pathogenic cell subpopulations and their communication pathways, which highlighted multiple potential therapeutic targets. This integrated dataset is encompassed in a publicly available web resource to aid mechanistic and translational studies of normal and diseased skin. 10.1016/j.jid.2023.02.040
SCAN: Spatiotemporal Cloud Atlas for Neural cells. Nucleic acids research The nervous system is one of the most complicated and enigmatic systems within the animal kingdom. Recently, the emergence and development of spatial transcriptomics (ST) and single-cell RNA sequencing (scRNA-seq) technologies have provided an unprecedented ability to systematically decipher the cellular heterogeneity and spatial locations of the nervous system from multiple unbiased aspects. However, efficiently integrating, presenting and analyzing massive multiomic data remains a huge challenge. Here, we manually collected and comprehensively analyzed high-quality scRNA-seq and ST data from the nervous system, covering 10 679 684 cells. In addition, multi-omic datasets from more than 900 species were included for extensive data mining from an evolutionary perspective. Furthermore, over 100 neurological diseases (e.g. Alzheimer's disease, Parkinson's disease, Down syndrome) were systematically analyzed for high-throughput screening of putative biomarkers. Differential expression patterns across developmental time points, cell types and ST spots were discerned and subsequently subjected to extensive interpretation. To provide researchers with efficient data exploration, we created a new database with interactive interfaces and integrated functions called the Spatiotemporal Cloud Atlas for Neural cells (SCAN), freely accessible at http://47.98.139.124:8799 or http://scanatlas.net. SCAN will benefit the neuroscience research community to better exploit the spatiotemporal atlas of the neural system and promote the development of diagnostic strategies for various neurological disorders. 10.1093/nar/gkad895
Single-cell transcriptome atlas of the human corpus cavernosum. Nature communications The corpus cavernosum is the most important structure for penile erection, and its dysfunction causes many physiological and psychological problems. However, its cellular heterogeneity and signalling networks at the molecular level are poorly understood because of limited access to samples. Here, we profile 64,993 human cavernosal single-cell transcriptomes from three males with normal erection and five organic erectile dysfunction patients. Cell communication analysis reveals that cavernosal fibroblasts are central to the paracrine signalling network and regulate microenvironmental homeostasis. Combining with immunohistochemical staining, we reveal the cellular heterogeneity and describe a detailed spatial distribution map for each fibroblast, smooth muscle and endothelial subcluster in the corpus cavernosum. Furthermore, comparative analysis and related functional experiments identify candidate regulatory signalling pathways in the pathological process. Our study provides an insight into the human corpus cavernosum microenvironment and a reference for potential erectile dysfunction therapies. 10.1038/s41467-022-31950-9
Single-cell transcriptomics of human-skin-equivalent organoids. Cell reports Several methods for generating human-skin-equivalent (HSE) organoid cultures are in use to study skin biology; however, few studies thoroughly characterize these systems. To fill this gap, we use single-cell transcriptomics to compare in vitro HSEs, xenograft HSEs, and in vivo epidermis. By combining differential gene expression, pseudotime analyses, and spatial localization, we reconstruct HSE keratinocyte differentiation trajectories that recapitulate known in vivo epidermal differentiation pathways and show that HSEs contain major in vivo cellular states. However, HSEs also develop unique keratinocyte states, an expanded basal stem cell program, and disrupted terminal differentiation. Cell-cell communication modeling shows aberrant epithelial-to-mesenchymal transition (EMT)-associated signaling pathways that alter upon epidermal growth factor (EGF) supplementation. Last, xenograft HSEs at early time points post transplantation significantly rescue many in vitro deficits while undergoing a hypoxic response that drives an alternative differentiation lineage. This study highlights the strengths and limitations of organoid cultures and identifies areas for potential innovation. 10.1016/j.celrep.2023.112511
Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics. Sottoriva Andrea,Spiteri Inmaculada,Piccirillo Sara G M,Touloumis Anestis,Collins V Peter,Marioni John C,Curtis Christina,Watts Colin,Tavaré Simon Proceedings of the National Academy of Sciences of the United States of America Glioblastoma (GB) is the most common and aggressive primary brain malignancy, with poor prognosis and a lack of effective therapeutic options. Accumulating evidence suggests that intratumor heterogeneity likely is the key to understanding treatment failure. However, the extent of intratumor heterogeneity as a result of tumor evolution is still poorly understood. To address this, we developed a unique surgical multisampling scheme to collect spatially distinct tumor fragments from 11 GB patients. We present an integrated genomic analysis that uncovers extensive intratumor heterogeneity, with most patients displaying different GB subtypes within the same tumor. Moreover, we reconstructed the phylogeny of the fragments for each patient, identifying copy number alterations in EGFR and CDKN2A/B/p14ARF as early events, and aberrations in PDGFRA and PTEN as later events during cancer progression. We also characterized the clonal organization of each tumor fragment at the single-molecule level, detecting multiple coexisting cell lineages. Our results reveal the genome-wide architecture of intratumor variability in GB across multiple spatial scales and patient-specific patterns of cancer evolution, with consequences for treatment design. 10.1073/pnas.1219747110
MiR-486-5p Downregulation Marks an Early Event in Colorectal Carcinogenesis. Diseases of the colon and rectum BACKGROUND:MicroRNAs are dysregulated in colorectal cancer and subsets correlated with advanced tumor stage and metastasis. Data are lacking on microRNA dysregulation from early to late-stage disease. OBJECTIVE:The purpose of this study was to identify a microRNA signature associated with the primary tumor and metastatic site in stage IV disease and to examine whether the signature is evident in earlier stages. DESIGN:A microRNA profile was generated and then explored in normal colon tissue (n = 5), early stage (stage I and II; n = 10), and late-stage (stage III and IV; n = 14) colorectal primary tumors via polymerase chain reaction to delineate molecular events that may promote colorectal carcinogenesis. SETTING:Genome-wide microRNA expression profiling was performed. PATIENTS:A total of 14 patient-matched stage IV primary colorectal cancer tumors and corresponding liver metastases were included. MAIN OUTCOME MEASURES:MicroRNA array technology was used to identify microRNA expression-predictive metastatic potential in the primary tumor. RESULTS:A distinct 9-member signature group of microRNAs was concurrent in stage IV primary colorectal cancer and their corresponding liver metastases, when compared with surrounding unaffected colon and liver tissue (microRNA-18b, microRNA-93, microRNA-182, microRNA-183, microRNA21, microRNA-486-5p, microRNA-500a, microRNA-552, and microRNA-941). Of the microRNA panel, only microRNA486-5p was differentially expressed in early stage colorectal cancer samples compared with normal tissue (p = 0.001) and additionally differentially expressed between late-stage colorectal cancer samples and normal tissue (p < 0.01). LIMITATIONS:Our microRNA profile was generated in a small subset of patients and will require validation in more samples. CONCLUSIONS:We identified a distinct microRNA signature in primary colon and matched metastatic disease. On additional investigation, 1 microRNA was differentially expressed in both early and late-stage cancer patient samples, and it may herald an early event in colorectal carcinogenesis. This study warrants additional investigation with a larger patient cohort to better understand the effect of microRNAs in carcinogenesis. See Video Abstract at http://links.lww.com/DCR/A723. 10.1097/DCR.0000000000001192
Development and Validation of a Combined Hypoxia and Immune Prognostic Classifier for Head and Neck Cancer. Clinical cancer research : an official journal of the American Association for Cancer Research PURPOSE:Intratumoral hypoxia and immunity have been correlated with patient outcome in various tumor settings. However, these factors are not currently considered for treatment selection in head and neck cancer (HNC) due to lack of validated biomarkers. Here we sought to develop a hypoxia-immune classifier with potential application in patient prognostication and prediction of response to targeted therapy. EXPERIMENTAL DESIGN:A 54-gene hypoxia-immune signature was constructed on the basis of literature review. Gene expression was analyzed using the The Cancer Genome Atlas (TCGA) HNC dataset ( = 275) and validated using two independent cohorts ( = 130 and 123). IHC was used to investigate the utility of a simplified protein signature. The spatial distribution of hypoxia and immune markers was examined using multiplex immunofluorescence staining. RESULTS:Unsupervised hierarchical clustering of TCGA dataset (development cohort) identified three patient subgroups with distinct hypoxia-immune phenotypes and survival profiles: hypoxia/immune, hypoxia/immune, and mixed, with 5-year overall survival (OS) rates of 71%, 51%, and 49%, respectively ( = 0.0015). The prognostic relevance of the hypoxia-immune gene signature was replicated in two independent validation cohorts. Only PD-L1 and intratumoral CD3 protein expression were associated with improved OS on multivariate analysis. Hypoxia/immune and hypoxia/immune tumors were overrepresented in "inflamed" and "immune-desert" microenvironmental profiles, respectively. Multiplex staining demonstrated an inverse correlation between CA-IX expression and prevalence of intratumoral CD3 T cells ( = -0.5464; = 0.0377), further corroborating the transcription-based classification. CONCLUSIONS:We developed and validated a hypoxia-immune prognostic transcriptional classifier, which may have clinical application to guide the use of hypoxia modification and targeted immunotherapies for the treatment of HNC. 10.1158/1078-0432.CCR-18-3314
Myeloid-specific KDM6B inhibition sensitizes glioblastoma to PD1 blockade. Nature cancer Glioblastoma (GBM) tumors are enriched in immune-suppressive myeloid cells and are refractory to immune checkpoint therapy (ICT). Targeting epigenetic pathways to reprogram the functional phenotype of immune-suppressive myeloid cells to overcome resistance to ICT remains unexplored. Single-cell and spatial transcriptomic analyses of human GBM tumors demonstrated high expression of an epigenetic enzyme-histone 3 lysine 27 demethylase (KDM6B)-in intratumoral immune-suppressive myeloid cell subsets. Importantly, myeloid cell-specific Kdm6b deletion enhanced proinflammatory pathways and improved survival in GBM tumor-bearing mice. Mechanistic studies showed that the absence of Kdm6b enhances antigen presentation, interferon response and phagocytosis in myeloid cells by inhibition of mediators of immune suppression including Mafb, Socs3 and Sirpa. Further, pharmacological inhibition of KDM6B mirrored the functional phenotype of Kdm6b-deleted myeloid cells and enhanced anti-PD1 efficacy. This study thus identified KDM6B as an epigenetic regulator of the functional phenotype of myeloid cell subsets and a potential therapeutic target for enhanced response to ICT. 10.1038/s43018-023-00620-0
Elevation of WNT5A expression in polyp formation in Lkb1+/- mice and Peutz-Jeghers syndrome. Lai Cecilia,Robinson James,Clark Sue,Stamp Gordon,Poulsom Richard,Silver Andrew The Journal of pathology Peutz-Jeghers syndrome (PJS) is a rare, inherited disease caused by germline mutation of the LKB1 gene. Patients with PJS develop characteristic polyps in the digestive tract and carry an elevated risk of cancers in multiple organs, including the intestinal tract. While LKB1 is capable of phosphorylating AMPK and regulates the mTOR pathway, it is also known to be a multitasking protein that can influence other cellular processes, including cell polarity. We hypothesized that there may be other biological pathways directly or indirectly affected by the loss of LKB1 in PJS and aimed to investigate this possibility through transcriptional profiling of polyps harvested from an Lkb1(+/-) mouse model of PJS and from PJS patients. We identified alterations in the mRNA level of a wide range of genes, including some that are involved in Wnt signalling (Wnt5a, Wif1, Dixdc1, Wnt11, Ccnd1, and Ccnd2), although we did not observe nuclear localization of β-catenin in over 93 human PJS intestinal polyps or in 24 gastric polyps from Lkb1(+/-) mice. Among these genes, WNT5A, a non-canonical and non-transforming Wnt, is consistently up-regulated in both Lkb1(+/-) mice and human PJS polyps at a high level. We performed in situ hybridization to further define the spatial expression pattern of WNT5A and observed a strong signal in the stroma of mouse and human polyps compared to no or very low expression in the mucosa. Our findings indicate that WNT5A plays an important role in PJS polyposis. 10.1002/path.2835
Fast alignment and preprocessing of chromatin profiles with Chromap. Nature communications As sequencing depth of chromatin studies continually grows deeper for sensitive profiling of regulatory elements or chromatin spatial structures, aligning and preprocessing of these sequencing data have become the bottleneck for analysis. Here we present Chromap, an ultrafast method for aligning and preprocessing high throughput chromatin profiles. Chromap is comparable to BWA-MEM and Bowtie2 in alignment accuracy and is over 10 times faster than traditional workflows on bulk ChIP-seq/Hi-C profiles and than 10x Genomics' CellRanger v2.0.0 pipeline on single-cell ATAC-seq profiles. 10.1038/s41467-021-26865-w
Methylation profiling defines an extensive field defect in histologically normal prostate tissues associated with prostate cancer. Yang Bing,Bhusari Sachin,Kueck Jessica,Weeratunga Pushpa,Wagner Jennifer,Leverson Glen,Huang Wei,Jarrard David F Neoplasia (New York, N.Y.) Prostate cancer (PCa) is typically found as a multifocal disease suggesting the potential for molecular defects within the morphologically normal tissue. The frequency and spatial extent of DNA methylation changes encompassing a potential field defect are unknown. A comparison of non-tumor-associated (NTA) prostate to histologically indistinguishable tumor-associated (TA) prostate tissues detected a distinct profile of DNA methylation alterations (0.2%) using genome-wide DNA arrays based on the Encyclopedia of DNA Elements 18 sequence that tile both gene-rich and poor regions. Hypomethylation (87%) occurred more frequently than hypermethylation (13%). Several of the most significantly altered loci (CAV1, EVX1, MCF2L, and FGF1) were then used as probes to map the extent of these DNA methylation changes in normal tissues from prostates containing cancer. In TA tissues, the extent of methylation was similar both adjacent (2 mm) and at a distance (>1 cm) from tumor foci. These loci were also able to distinguish NTA from TA tissues in a validation set of patient samples. These mapping studies indicate that a spatially widespread epigenetic defect occurs in the peripheral prostate tissues of men who have PCa that may be useful in the detection of this disease. 10.1593/neo.13280
Synaptic and peptidergic connectome of a neurosecretory center in the annelid brain. Williams Elizabeth A,Verasztó Csaba,Jasek Sanja,Conzelmann Markus,Shahidi Réza,Bauknecht Philipp,Mirabeau Olivier,Jékely Gáspár eLife Neurosecretory centers in animal brains use peptidergic signaling to influence physiology and behavior. Understanding neurosecretory center function requires mapping cell types, synapses, and peptidergic networks. Here we use transmission electron microscopy and gene expression mapping to analyze the synaptic and peptidergic connectome of an entire neurosecretory center. We reconstructed 78 neurosecretory neurons and mapped their synaptic connectivity in the brain of larval , a marine annelid. These neurons form an anterior neurosecretory center expressing many neuropeptides, including hypothalamic peptide orthologs and their receptors. Analysis of peptide-receptor pairs in spatially mapped single-cell transcriptome data revealed sparsely connected networks linking specific neuronal subsets. We experimentally analyzed one peptide-receptor pair and found that a neuropeptide can couple neurosecretory and synaptic brain signaling. Our study uncovered extensive networks of peptidergic signaling within a neurosecretory center and its connection to the synaptic brain. 10.7554/eLife.26349
Osteoprotegerin in bone metastases: mathematical solution to the puzzle. Ryser Marc D,Qu Yiding,Komarova Svetlana V PLoS computational biology Bone is a common site for cancer metastasis. To create space for their growth, cancer cells stimulate bone resorbing osteoclasts. Cytokine RANKL is a key osteoclast activator, while osteoprotegerin (OPG) is a RANKL decoy receptor and an inhibitor of osteoclastogenesis. Consistently, systemic application of OPG decreases metastatic tumor burden in bone. However, OPG produced locally by cancer cells was shown to enhance osteolysis and tumor growth. We propose that OPG produced by cancer cells causes a local reduction in RANKL levels, inducing a steeper RANKL gradient away from the tumor and towards the bone tissue, resulting in faster resorption and tumor expansion. We tested this hypothesis using a mathematical model of nonlinear partial differential equations describing the spatial dynamics of OPG, RANKL, PTHrP, osteoclasts, tumor and bone mass. We demonstrate that at lower expression rates, tumor-derived OPG enhances the chemotactic RANKL gradient and osteolysis, whereas at higher expression rates OPG broadly inhibits RANKL and decreases osteolysis and tumor burden. Moreover, tumor expression of a soluble mediator inducing RANKL in the host tissue, such as PTHrP, is important for correct orientation of the RANKL gradient. A meta-analysis of OPG, RANKL and PTHrP expression in normal prostate, carcinoma and metastatic tissues demonstrated an increase in expression of OPG, but not RANKL, in metastatic prostate cancer, and positive correlation between OPG and PTHrP in metastatic prostate cancer. The proposed mechanism highlights the importance of the spatial distribution of receptors, decoys and ligands, and can be applied to other systems involving regulation of spatially anisotropic processes. 10.1371/journal.pcbi.1002703
High-content image analysis to study phenotypic heterogeneity in endothelial cell monolayers. Chesnais Francois,Hue Jonas,Roy Errin,Branco Marco,Stokes Ruby,Pellon Aize,Le Caillec Juliette,Elbahtety Eyad,Battilocchi Matteo,Danovi Davide,Veschini Lorenzo Journal of cell science Endothelial cells (ECs) are heterogeneous across and within tissues, reflecting distinct, specialised functions. EC heterogeneity has been proposed to underpin EC plasticity independently from vessel microenvironments. However, heterogeneity driven by contact-dependent or short-range cell-cell crosstalk cannot be evaluated with single cell transcriptomic approaches, as spatial and contextual information is lost. Nonetheless, quantification of EC heterogeneity and understanding of its molecular drivers is key to developing novel therapeutics for cancer, cardiovascular diseases and for revascularisation in regenerative medicine. Here, we developed an EC profiling tool (ECPT) to examine individual cells within intact monolayers. We used ECPT to characterise different phenotypes in arterial, venous and microvascular EC populations. In line with other studies, we measured heterogeneity in terms of cell cycle, proliferation, and junction organisation. ECPT uncovered a previously under-appreciated single-cell heterogeneity in NOTCH activation. We correlated cell proliferation with different NOTCH activation states at the single-cell and population levels. The positional and relational information extracted with our novel approach is key to elucidating the molecular mechanisms underpinning EC heterogeneity. 10.1242/jcs.259104
The spatial association of gene expression evolves from synchrony to asynchrony and stochasticity with age. Wang Qi,Huang Jianhua,Zhang Xinmin,Wu Bin,Liu Xiaoyu,Shen Ziyin PloS one For multicellular organisms, different tissues coordinate to integrate physiological functions, although this systematically and gradually declines in the aging process. Therefore, an association exists between tissue coordination and aging, and investigating the evolution of tissue coordination with age is of interest. In the past decade, both common and heterogeneous aging processes among tissues were extensively investigated. The results on spatial association of gene changes that determine lifespan appear complex and paradoxical. To reconcile observed commonality and heterogeneity of gene changes among tissues and to address evolution feature of tissue coordination with age, we introduced a new analytical strategy to systematically analyze genome-wide spatio-temporal gene expression profiles. We first applied the approach to natural aging process in three species (Rat, Mouse and Drosophila) and then to anti-aging process in Mouse. The results demonstrated that temporal gene expression alteration in different tissues experiences a progressive association evolution from spatial synchrony to asynchrony and stochasticity with age. This implies that tissue coordination gradually declines with age. Male mice showed earlier spatial asynchrony in gene expression than females, suggesting that male animals are more prone to aging than females. The confirmed anti-aging interventions (resveratrol and caloric restriction) enhanced tissue coordination, indicating their underlying anti-aging mechanism on multiple tissue levels. Further, functional analysis suggested asynchronous DNA/protein damage accumulation as well as asynchronous repair, modification and degradation of DNA/protein in tissues possibly contributes to asynchronous and stochastic changes of tissue microenvironment. This increased risk for a variety of age-related diseases such as neurodegeneration and cancer that eventually accelerate organismal aging and death. Our study suggests a novel molecular event occurring in aging process of multicellular species that may represent an intrinsic molecular mechanism of aging. 10.1371/journal.pone.0024076
Resolving cancer-stroma interfacial signalling and interventions with micropatterned tumour-stromal assays. Shen Keyue,Luk Samantha,Hicks Daniel F,Elman Jessica S,Bohr Stefan,Iwamoto Yoshiko,Murray Ryan,Pena Kristen,Wang Fangjing,Seker Erkin,Weissleder Ralph,Yarmush Martin L,Toner Mehmet,Sgroi Dennis,Parekkadan Biju Nature communications Tumour-stromal interactions are a determining factor in cancer progression. In vivo, the interaction interface is associated with spatially resolved distributions of cancer and stromal phenotypes. Here, we establish a micropatterned tumour-stromal assay (μTSA) with laser capture microdissection to control the location of co-cultured cells and analyse bulk and interfacial tumour-stromal signalling in driving cancer progression. μTSA reveals a spatial distribution of phenotypes in concordance with human oestrogen receptor-positive (ER+) breast cancer samples, and heterogeneous drug activity relative to the tumour-stroma interface. Specifically, an unknown mechanism of reversine is shown in targeting tumour-stromal interfacial interactions using ER+ MCF-7 breast cancer and bone marrow-derived stromal cells. Reversine suppresses MCF-7 tumour growth and bone metastasis in vivo by reducing tumour stromalization including collagen deposition and recruitment of activated stromal cells. This study advocates μTSA as a platform for studying tumour microenvironmental interactions and cancer field effects with applications in drug discovery and development. 10.1038/ncomms6662
Proteomic profiling of the oncogenic septin 9 reveals isoform-specific interactions in breast cancer cells. Proteomics Septins are a family of multimeric GTP-binding proteins, which are abnormally expressed in cancer. Septin 9 (SEPT9) is an essential and ubiquitously expressed septin with multiple isoforms, which have differential expression patterns and effects in breast cancer cells. It is unknown, however, if SEPT9 isoforms associate with different molecular networks and functions. Here, we performed a proteomic screen in MCF-7 breast cancer cells to identify the interactome of GFP-SEPT9 isoforms 1, 4 and 5, which vary significantly in their N-terminal extensions. While all three isoforms associated with SEPT2 and SEPT7, the truncated SEPT9_i4 and SEPT9_i5 interacted with septins of the SEPT6 group more promiscuously than SEPT9_i1, which bound predominately SEPT8. Spatial mapping and functional clustering of non-septin partners showed isoform-specific differences in interactions with proteins of distinct subcellular organelles (e.g., nuclei, centrosomes, cilia) and functions such as cell signalling and ubiquitination. The interactome of the full length SEPT9_i1 was more enriched in cytoskeletal regulators, while the truncated SEPT9_i4 and SEPT9_i5 exhibited preferential and isoform-specific interactions with nuclear, signalling, and ubiquitinating proteins. These data provide evidence for isoform-specific interactions, which arise from truncations in the N-terminal extensions of SEPT9, and point to novel roles in the pathogenesis of breast cancer. 10.1002/pmic.202100155
A single-cell atlas of mouse lung development. Development (Cambridge, England) Lung organogenesis requires precise timing and coordination to effect spatial organization and function of the parenchymal cells. To provide a systematic broad-based view of the mechanisms governing the dynamic alterations in parenchymal cells over crucial periods of development, we performed a single-cell RNA-sequencing time-series yielding 102,571 epithelial, endothelial and mesenchymal cells across nine time points from embryonic day 12 to postnatal day 14 in mice. Combining computational fate-likelihood prediction with RNA in situ hybridization and immunofluorescence, we explore lineage relationships during the saccular to alveolar stage transition. The utility of this publicly searchable atlas resource (www.sucrelab.org/lungcells) is exemplified by discoveries of the complexity of type 1 pneumocyte function and characterization of mesenchymal Wnt expression patterns during the saccular and alveolar stages - wherein major expansion of the gas-exchange surface occurs. We provide an integrated view of cellular dynamics in epithelial, endothelial and mesenchymal cell populations during lung organogenesis. 10.1242/dev.199512
Integration of Single-Cell Genomics Datasets. Adey Andrew C Cell Welch et al. and Stuart et al. present novel techniques for the integration of single-cell RNA-seq datasets across multiple platforms, individuals, and species. They both extend these strategies to map cell types between RNA-seq datasets with epigenetic properties and in situ transcript profiling. The ability to transfer information between datasets and spatial methods will enable more comprehensive profiling and comparisons of cell populations in complex biological systems. 10.1016/j.cell.2019.05.034
Disease-gene discovery by integration of 3D gene expression and transcription factor binding affinities. Piro Rosario M,Molineris Ivan,Di Cunto Ferdinando,Eils Roland,König Rainer Bioinformatics (Oxford, England) MOTIVATION:The computational evaluation of candidate genes for hereditary disorders is a non-trivial task. Several excellent methods for disease-gene prediction have been developed in the past 2 decades, exploiting widely differing data sources to infer disease-relevant functional relationships between candidate genes and disorders. We have shown recently that spatially mapped, i.e. 3D, gene expression data from the mouse brain can be successfully used to prioritize candidate genes for human Mendelian disorders of the central nervous system. RESULTS:We improved our previous work 2-fold: (i) we demonstrate that condition-independent transcription factor binding affinities of the candidate genes' promoters are relevant for disease-gene prediction and can be integrated with our previous approach to significantly enhance its predictive power; and (ii) we define a novel similarity measure-termed Relative Intensity Overlap-for both 3D gene expression patterns and binding affinity profiles that better exploits their disease-relevant information content. Finally, we present novel disease-gene predictions for eight loci associated with different syndromes of unknown molecular basis that are characterized by mental retardation. 10.1093/bioinformatics/bts720
Quantitative evidence for early metastatic seeding in colorectal cancer. Nature genetics Both the timing and molecular determinants of metastasis are unknown, hindering treatment and prevention efforts. Here we characterize the evolutionary dynamics of this lethal process by analyzing exome-sequencing data from 118 biopsies from 23 patients with colorectal cancer with metastases to the liver or brain. The data show that the genomic divergence between the primary tumor and metastasis is low and that canonical driver genes were acquired early. Analysis within a spatial tumor growth model and statistical inference framework indicates that early disseminated cells commonly (81%, 17 out of 21 evaluable patients) seed metastases while the carcinoma is clinically undetectable (typically, less than 0.01 cm). We validated the association between early drivers and metastasis in an independent cohort of 2,751 colorectal cancers, demonstrating their utility as biomarkers of metastasis. This conceptual and analytical framework provides quantitative in vivo evidence that systemic spread can occur early in colorectal cancer and illuminates strategies for patient stratification and therapeutic targeting of the canonical drivers of tumorigenesis. 10.1038/s41588-019-0423-x
Gene expression profiles of NO- and HNO-donor treated breast cancer cells: insights into tumor response and resistance pathways. Cheng Robert Y S,Basudhar Debashree,Ridnour Lisa A,Heinecke Julie L,Kesarwala Aparna H,Glynn Sharon,Switzer Christopher H,Ambs Stefan,Miranda Katrina M,Wink David A Nitric oxide : biology and chemistry Nitric oxide (NO) synthase 2 (NOS2), a major inflammatory protein, modulates disease progression via NO in a number of pathologies, including cancer. The role of NOS2-derived NO is not only flux-dependent, which is higher in mouse vs human cells, but also varies based on spatial and temporal distribution both within tumor cells and in the tumor microenvironment. NO donors have been utilized to mimic NO flux conditions and to investigate the effects of varied NO concentrations. As a wide range of effects mediated by NO and other nitrogen oxides such as nitroxyl (HNO) have been elucidated, multiple NO- and HNO-releasing compounds have been developed as potential therapeutics, including as tumor modulators. One of the challenges is to determine differences in biomarker expression from extracellular vs intracellular generation of NO or HNO. Taking advantage of new NO and HNO releasing agents, we have characterized the gene expression profile of estrogen receptor-negative human breast cancer (MDA-MB-231) cells following exposure to aspirin, the NO donor DEA/NO, the HNO donor IPA/NO andtheir intracellularly-activated prodrug conjugates DEA/NO-aspirin and IPA/NO-aspirin. Comparison of the gene expression profiles demonstrated that several genes were uniquely expressed with respect to NO or HNO, such as miR-21, HSP70, cystathionine γ-lyase and IL24. These findings provide insight into targets and pathways that could be therapeutically exploited by the redox related species NO and HNO. 10.1016/j.niox.2014.08.003
Evaluation of targeted therapies in advanced breast cancer: the need for large-scale molecular screening and transformative clinical trial designs. Fadoukhair Z,Zardavas D,Chad M A,Goulioti T,Aftimos P,Piccart M Oncogene Breast cancer (BC) has been classified into four intrinsic subtypes through seminal studies employing gene expression profiling analysis of primary tumours, namely the luminal A and B subtypes, the human epidermal growth factor receptor 2-like subtype and the basal-like subtype. More recently, the emergence of high-throughput genomic sequencing techniques, such as next-generation or massive parallel sequencing has expanded our understanding of the complex genomic landscapes of BC, with marked intertumour heterogeneity seen among different patients. In addition, increasing evidence indicates intratumour heterogeneity, with molecular differences observed within one patient, both spatially and longitudinally. These phenomena have an impact on the clinical development of molecularly targeted agents, with the classical paradigm of population-based clinical trials being no longer efficient. In the era of genomically driven oncology, three complementary tools can accelerate the clinical development of targeted agents for advanced BC as follows: (i) the implementation of molecular profiling of metastatic tumour lesions, as exemplified by the AURORA (Aiming to Understand the Molecular Aberrations in Metastatic Breast Cancer) programme; (ii) serial assessments of circulating tumour DNA, allowing a more thorough molecular interrogation of metastatic tumour burden; and (iii) new innovative clinical trial designs able to address the challenges of the increasing molecular fragmentation of BC. 10.1038/onc.2015.249
Integrated miRNA-mRNA spatial signature for oral squamous cell carcinoma: a prospective profiling study of Narrow Band Imaging guided resection. Farah Camile S,Fox Simon A,Dalley Andrew J Scientific reports Oral squamous cell carcinoma (OSCC) is a common malignancy for which there is poor prognosis and limited therapeutic options. The objective was to identify mRNA targets of dysregulated miRNAs in OSCC using integrated analysis and understand molecular abnormality in surgical margins. We used biopsies along the spatial axis from normal tissue defined by narrow band imaging (NBI) through conventional white light (WL) margins to tumour from 18 patients undergoing surgical resection for OSCC. Overall 119 miRNA and 4794 mRNA were differentially expressed along the adjacent normal tissue to tumour axis. Analysis of miRNA profiles demonstrated the NBI margins were molecularly distinct from both the tumour and WL margin. Integrated analysis identified 193 miRNA-mRNA interactions correlated to the spatial axis of NBI-WL-T. We used cross-validation analysis to derive a spatial interactome signature of OSCC comprising 100 putative miRNA-mRNA interactions between 40 miRNA and 96 mRNA. Bioinformatic analysis suggests that miRNA dysregulation in OSCC may contribute to activation of the oncostatin M, BDNF and TGF-β pathways. Our data demonstrates that surgical margins defined by NBI leave less potentially malignant residual tissue. The miRNA-mRNA interactome provides insight into dysregulated miRNA signalling in OSCC and supports molecular definition of tumour margins. 10.1038/s41598-018-19341-x
Overlapping gene expression profiles of cell migration and tumor invasion in human bladder cancer identify metallothionein 1E and nicotinamide N-methyltransferase as novel regulators of cell migration. Wu Y,Siadaty M S,Berens M E,Hampton G M,Theodorescu D Oncogene Cell migration is essential to cancer invasion and metastasis and is spatially and temporally integrated through transcriptionally dependent and independent mechanisms. As cell migration is studied in vitro, it is important to identify genes that both drive cell migration and are biologically relevant in promoting invasion and metastasis in patients with cancer. Here, gene expression profiling and a high-throughput cell migration system answers this question in human bladder cancer. In vitro migration rates of 40 microarray-profiled human bladder cancer cell lines were measured by radial migration assay. Genes whose expression was either directly or inversely associated with cell migration rate were identified and subsequently evaluated for their association with cancer stage in 61 patients. This analysis identified genes known to be associated with cell invasion such as versican, and novel ones, including metallothionein 1E (MT1E) and nicotinamide N-methyltransferase (NNMT), whose expression correlated positively with cancer cell migration and tumor stage. Using loss of function analysis, we show that MT1E and NNMT are necessary for cancer cell migration. These studies provide a general approach to identify the clinically relevant genes in cancer cell migration and mechanistically implicate two novel genes in this process in human bladder cancer. 10.1038/onc.2008.264
Spatial transcriptomic analysis reveals inflammatory foci defined by senescent cells in the white matter, hippocampi and cortical grey matter in the aged mouse brain. GeroScience There is strong evidence that aging is associated with an increased presence of senescent cells in the brain. The finding that treatment with senolytic drugs improves cognitive performance of aged laboratory mice suggests that increased cellular senescence is causally linked to age-related cognitive decline. The relationship between senescent cells and their relative locations within the brain is critical to understanding the pathology of age-related cognitive decline and dementia. To assess spatial distribution of cellular senescence in the aged mouse brain, spatially resolved whole transcriptome mRNA expression was analyzed in sections of brains derived from young (3 months old) and aged (28 months old) C57BL/6 mice while capturing histological information in the same tissue section. Using this spatial transcriptomics (ST)-based method, microdomains containing senescent cells were identified on the basis of their senescence-related gene expression profiles (i.e., expression of the senescence marker cyclin-dependent kinase inhibitor p16 encoded by the Cdkn2a gene) and were mapped to different anatomical brain regions. We confirmed that brain aging is associated with increased cellular senescence in the white matter, the hippocampi and the cortical grey matter. Transcriptional analysis of the senescent cell-containing ST spots shows that presence of senescent cells is associated with a gene expression signature suggestive of neuroinflammation. GO enrichment analysis of differentially expressed genes in the outer region of senescent cell-containing ST spots ("neighboring ST spots") also identified functions related to microglia activation and neuroinflammation. In conclusion, senescent cells accumulate with age in the white matter, the hippocampi and cortical grey matter and likely contribute to the genesis of inflammatory foci in a paracrine manner. 10.1007/s11357-022-00521-7
Spatial transcriptomics reveals distinct tissue niches linked with steroid responsiveness in acute gastrointestinal GVHD. Blood Severe acute graft-versus-host disease (aGVHD) is associated with significant mortality and morbidity, especially in steroid-resistant (SR) cases. Spatial transcriptomic technology can elucidate tissue-based interactions in vivo and possibly identify predictors of treatment response. Tissue sections from 32 treatment-naïve patients with biopsy-confirmed lower gastrointestinal (GI) aGVHD were obtained. The GeoMx digital spatial profiler was used to capture transcriptome profiles of >18 000 genes from different foci of immune infiltrates, colonic epithelium, and vascular endothelium. Each tissue compartment sampled showed 2 distinct clusters that were analyzed for differential expression and spatially resolved correlation of gene signatures. Classic cell-mediated immunity signatures, normal differentiated epithelial cells, and inflamed vasculature dominated foci sampled from steroid-sensitive cases. In contrast, a neutrophil predominant noncanonical inflammation with regenerative epithelial cells and some indication of angiogenic endothelial response was overrepresented in areas from SR cases. Evaluation of potential prognostic biomarkers identified ubiquitin specific peptidase 17-like (USP17L) family of genes as being differentially expressed in immune cells from patients with worsened survival. In summary, we demonstrate distinct tissue niches with unique gene expression signatures within lower GI tissue from patients with aGVHD and provide evidence of a potential prognostic biomarker. 10.1182/blood.2023020644
Integration of spatial and single-cell transcriptomic data elucidates mouse organogenesis. Nature biotechnology Molecular profiling of single cells has advanced our knowledge of the molecular basis of development. However, current approaches mostly rely on dissociating cells from tissues, thereby losing the crucial spatial context of regulatory processes. Here, we apply an image-based single-cell transcriptomics method, sequential fluorescence in situ hybridization (seqFISH), to detect mRNAs for 387 target genes in tissue sections of mouse embryos at the 8-12 somite stage. By integrating spatial context and multiplexed transcriptional measurements with two single-cell transcriptome atlases, we characterize cell types across the embryo and demonstrate that spatially resolved expression of genes not profiled by seqFISH can be imputed. We use this high-resolution spatial map to characterize fundamental steps in the patterning of the midbrain-hindbrain boundary (MHB) and the developing gut tube. We uncover axes of cell differentiation that are not apparent from single-cell RNA-sequencing (scRNA-seq) data, such as early dorsal-ventral separation of esophageal and tracheal progenitor populations in the gut tube. Our method provides an approach for studying cell fate decisions in complex tissues and development. 10.1038/s41587-021-01006-2
Dynamic MAPK signaling activity underlies a transition from growth arrest to proliferation in mutant tumors. Ji Tiantian,Zhang Lina,Deng Mingxi,Huang Shengshuo,Wang Ying,Pham Tri Thanh,Smith Andrew Alan,Sridhar Varun,Cabernard Clemens,Wang Jiguang,Yan Yan Disease models & mechanisms Human tumors exhibit plasticity and evolving capacity over time. It is difficult to study the mechanisms of how tumors change over time in human patients, in particular during the early stages when a few oncogenic cells are barely detectable. Here, we used a tumor model caused by loss of (), a highly conserved apicobasal cell polarity gene, to investigate the spatial-temporal dynamics of early tumorigenesis events. The fly mutant tumors have been successfully used to model many aspects of tumorigenesis processes. However, it is still unknown whether mutant tumors exhibit plasticity and evolvability along the temporal axis. We found that mutant tumors displayed different growth rates and cell cycle profiles over time, indicative of a growth arrest-to-proliferation transition as the mutant tumors progress. Longitudinal bulk and single-cell transcriptomic analysis of mutant tumors revealed that the MAPK pathway, including JNK and ERK signaling activities, showed quantitative changes over time. We found that high JNK signaling activity caused G2/M cell cycle arrest in early mutant tumors. In addition, JNK signaling activity displayed a radial polarity with the JNK cells located at the periphery of mutant tumors, providing an inherent mechanism that leads to an overall decrease in JNK signaling activity over time. We also found that ERK signaling activity, in contrast to JNK activity, increased over time and promoted growth in late-stage mutant tumors. Furthermore, high JNK signaling activity repressed ERK signaling activity in early mutant tumors. Together, these data demonstrate that dynamic MAPK signaling activity, fueled by intratumor heterogeneity derived from tissue topological differences, drives a growth arrest-to-proliferation transition in mutant tumors.This article has an associated First Person interview with the joint first authors of the paper. 10.1242/dmm.040147
Gene expression profiling of human diseases by serial analysis of gene expression. Ye Shui Q,Usher David C,Zhang Li Q Journal of biomedical science Until recently, the approach to understanding the molecular basis of complex syndromes such as cancer, coronary artery disease, and diabetes was to study the behavior of individual genes. However, it is generally recognized that expression of a number of genes is coordinated both spatially and temporally and that this coordination changes during the development and progression of diseases. Newly developed functional genomic approaches, such as serial analysis of gene expression (SAGE) and DNA microarrays have enabled researchers to determine the expression pattern of thousands of genes simultaneously. One attractive feature of SAGE compared to microarrays is its ability to quantify gene expression without prior sequence information or information about genes that are thought to be expressed. SAGE has been successfully applied to the gene expression profiling of a number of human diseases. In this review, we will first discuss SAGE technique and contrast it to microarray. We will then highlight new biological insights that have emerged from its application to the study of human diseases. 10.1007/bf02256531
Spatial profiling of early primate gastrulation in utero. Nature Gastrulation controls the emergence of cellular diversity and axis patterning in the early embryo. In mammals, this transformation is orchestrated by dynamic signalling centres at the interface of embryonic and extraembryonic tissues. Elucidating the molecular framework of axis formation in vivo is fundamental for our understanding of human development and to advance stem-cell-based regenerative approaches. Here we illuminate early gastrulation of marmoset embryos in utero using spatial transcriptomics and stem-cell-based embryo models. Gaussian process regression-based 3D transcriptomes delineate the emergence of the anterior visceral endoderm, which is hallmarked by conserved (HHEX, LEFTY2, LHX1) and primate-specific (POSTN, SDC4, FZD5) factors. WNT signalling spatially coordinates the formation of the primitive streak in the embryonic disc and is counteracted by SFRP1 and SFRP2 to sustain pluripotency in the anterior domain. Amnion specification occurs at the boundaries of the embryonic disc through ID1, ID2 and ID3 in response to BMP signalling, providing a developmental rationale for amnion differentiation of primate pluripotent stem cells (PSCs). Spatial identity mapping demonstrates that primed marmoset PSCs exhibit the highest similarity to the anterior embryonic disc, whereas naive PSCs resemble the preimplantation epiblast. Our 3D transcriptome models reveal the molecular code of lineage specification in the primate embryo and provide an in vivo reference to decipher human development. 10.1038/s41586-022-04953-1
Heterogeneous expression of zinc-finger E-box-binding homeobox 1 plays a pivotal role in metastasis via regulation of miR-200c in epithelial-mesenchymal transition. Muto Yuta,Suzuki Koichi,Kato Takaharu,Tsujinaka Shingo,Ichida Kosuke,Takayama Yuji,Fukui Taro,Kakizawa Nao,Watanabe Fumiaki,Saito Masaaki,Futsuhara Kazushige,Noda Hiroshi,Miyakura Yasuyuki,Konishi Fumio,Rikiyama Toshiki International journal of oncology Although epithelial-mesenchymal transition (EMT) has been implicated as the pivotal event in metastasis, there is insufficient evidence related to EMT in clinical settings. Intratumor heterogeneity may lead to underestimation of gene expression representing EMT. In the present study, we investigated the expression of EMT-associated genes and microRNAs in primary colorectal cancer while considering intratumor heterogeneity. One-hundred and thirty-three multiple spatially separated samples were obtained from 8 patients with metastatic colorectal cancers and 8 with non-metastatic colorectal cancers, from the tumor center (TC), invasive front (IF) and metastasis. Differences in gene and microRNA expression were investigated by microarray and quantitative reverse-transcription PCR. Gene expression microarray analysis detected 7920 sites showing differing levels of gene expression among the TC, IF and metastasis. Expression of the EMT-associated gene zinc-finger E-box-binding homeobox 1 (ZEB1) significantly increased in the IF (p<0.01). To exclude individual differences, the expression ratio between TC and IF in each tumor was applied to analysis. This approach enabled recognition of the activation of the VEGF and Wnt signaling pathways, which were involved in metastasis via promotion of EMT. While no activation of these pathways was seen at the TC, regardless of whether tumors were metastatic or non-metastatic, they were preferentially activated at the IF in metastatic tumors, where high ZEB1 expression was seen in connection with decreased miR-200c expression. Multiple sampling in a tumor revealed that heterogeneous ZEB1 expression induced by EMT-associated signaling pathways played a pivotal role in metastasis via regulation of miR-200c. 10.3892/ijo.2016.3583
Evaluation of genetic patterns in different tumor areas of intermediate-grade prostatic adenocarcinomas by high-resolution genomic array analysis. van Dekken Herman,Paris Pamela L,Albertson Donna G,Alers Janneke C,Andaya Armann,Kowbel David,van der Kwast Theodorus H,Pinkel Daniel,Schröder Fritz H,Vissers Kees J,Wildhagen Mark F,Collins Colin Genes, chromosomes & cancer Prostate cancer is known for its highly heterogeneous histological appearance. Data concerning the cytogenetic content of areas with different histology are sparse. We have genetically evaluated 10 prostatic adenocarcinomas with intermediate histopathological grades (Gleason score 7) that showed two distinctive growth patterns with different pathologies, that is, Gleason grades 3 and 4 (G3 and G4). The G3 and G4 tumor specimens were taken from spatially separated regions within the cancer mass. Array-based comparative genomic hybridization (aCGH) was performed to obtain genotypes from the 10 pairs of G3 and G4 cancer areas. The cancer DNAs were retrieved from formalin-fixed and paraffin-embedded tissues allowing optimal recognition and selection of target cells. A genome-wide 2,400-element BAC array that provided high-resolution detection of both deletions and amplifications was used. In the 20 G3 and G4 areas, 252 genomic aberrations (88 gains, 164 deletions) were noted, of which 86 were concurrent in G3 and G4 areas (34% overlap). Ninety-five of the 252 alterations were defined by a single BAC clone (54 gains, 41 deletions). Overlapping changes were more frequent for deletions (46%) than for gains (13%). Frequent coinciding deletions (> or = 20% of tumors) were seen on 8p (60%), 6q (30%), 1p (20%), 2q (20%), proximal 8q (20%), 10q (20%), 13q (20%), 16q (20%), and 18q (20%). A frequent overlapping gain (> or = 20% of tumors) was detected on distal 13q (20%). The patterns of imbalance could be found to coincide in the G3 and G4 areas of the majority of cancers. Array-based CGH can be used as a tool for the evaluation of genetic patterns in prostate cancer. 10.1002/gcc.20001
A holistic view of the malignant organism we call glioblastoma. Cell Tumors are not simply a chaotic mass of mutated cells but can follow complex organizational principles, including in space. In this issue of Cell, Mathur and colleagues reconstruct a 3D genomic, epigenomic, and transcriptomic spatial cartograph of glioblastoma, offering a "whole-tumor" perspective with patterns of clonal expansion that are embedded in neurodevelopmental hierarchy. 10.1016/j.cell.2023.12.021
DNMAD: web-based diagnosis and normalization for microarray data. Vaquerizas Juan M,Dopazo Joaquín,Díaz-Uriarte Ramón Bioinformatics (Oxford, England) SUMMARY:We present a web server for Diagnosis and Normalization of MicroArray Data (DNMAD). DNMAD includes several common data transformations such as spatial and global robust local regression or multiple slide normalization, and allows for detecting several kinds of errors that result from the manipulation and the image analysis of the arrays. This tool offers a user-friendly interface, and is completely integrated within the Gene Expression Pattern Analysis Suite (GEPAS). AVAILABILITY:The tool is accessible on-line at http://dnmad.bioinfo.cnio.es. 10.1093/bioinformatics/bth401
In-channel printing-device opening assay for micropatterning multiple cells and gene analysis. Zhou Hao,Zhao Liang,Zhang Xueji Analytical chemistry Herein we report an easy but versatile method for patterning different cells on a single substrate by using a microfluidic approach that allows not only spatial and temporal control of multiple microenvironments but also retrieval of specific treated cells to profile their expressed genetic information at around 10-cell resolution. By taking advantages of increased surface area of gold nanoparticles on a poly(dimethylsiloxane) (PDMS) coated substrate, cell adhesive-promotive protein, human fibronectin (hFN) can be significantly accumulated on designed regions where cells can recognize the protein and spread out. Moreover, the whole device can be easily opened by hand without any loss of patterned cells which could be retrieved by mouth-pipet. Consequently, we demonstrate the possibility of analyzing the difference of gene expression patterns between wild type MCF-7 cell and MCF/Adr (drug-resistant cell line) from less than 400 cells in total for a single comprehensive assay, including parallel experiments, controls, and multiple dose treatments. Certain genes, especially the P-glycoprotein coding gene (ABCB1), show high expression level in resistant cells compared with the wild type, suggesting a possible pathway that may contribute to the antidrug mechanism. 10.1021/ac504823s
Clustering gene expression data using a diffraction-inspired framework. Dinger Steven C,Van Wyk Michael A,Carmona Sergio,Rubin David M Biomedical engineering online BACKGROUND:The recent developments in microarray technology has allowed for the simultaneous measurement of gene expression levels. The large amount of captured data challenges conventional statistical tools for analysing and finding inherent correlations between genes and samples. The unsupervised clustering approach is often used, resulting in the development of a wide variety of algorithms. Typical clustering algorithms require selecting certain parameters to operate, for instance the number of expected clusters, as well as defining a similarity measure to quantify the distance between data points. The diffraction-based clustering algorithm however is designed to overcome this necessity for user-defined parameters, as it is able to automatically search the data for any underlying structure. METHODS:The diffraction-based clustering algorithm presented in this paper is tested using five well-known expression datasets pertaining to cancerous tissue samples. The clustering results are then compared to those results obtained from conventional algorithms such as the k-means, fuzzy c-means, self-organising map, hierarchical clustering algorithm, Gaussian mixture model and density-based spatial clustering of applications with noise (DBSCAN). The performance of each algorithm is measured using an average external criterion and an average validity index. RESULTS:The diffraction-based clustering algorithm is shown to be independent of the number of clusters as the algorithm searches the feature space and requires no form of parameter selection. The results show that the diffraction-based clustering algorithm performs significantly better on the real biological datasets compared to the other existing algorithms. CONCLUSION:The results of the diffraction-based clustering algorithm presented in this paper suggest that the method can provide researchers with a new tool for successfully analysing microarray data. 10.1186/1475-925X-11-85
How many clusters: a validation index for arbitrary-shaped clusters. Bayá Ariel E,Granitto Pablo M IEEE/ACM transactions on computational biology and bioinformatics Clustering validation indexes are intended to assess the goodness of clustering results. Many methods used to estimate the number of clusters rely on a validation index as a key element to find the correct answer. This paper presents a new validation index based on graph concepts, which has been designed to find arbitrary shaped clusters by exploiting the spatial layout of the patterns and their clustering label. This new clustering index is combined with a solid statistical detection framework, the gap statistic. The resulting method is able to find the right number of arbitrary-shaped clusters in diverse situations, as we show with examples where this information is available. A comparison with several relevant validation methods is carried out using artificial and gene expression data sets. The results are very encouraging, showing that the underlying structure in the data can be more accurately detected with the new clustering index. Our gene expression data results also indicate that this new index is stable under perturbation of the input data. 10.1109/TCBB.2013.32
Tumor slices as a model to evaluate doxorubicin in vitro treatment and expression of trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 in canine mammary gland cancer. Sobral Renata A,Honda Suzana T,Katayama Maria Lucia H,Brentani Helena,Brentani M Mitzi,Patrão Diogo F C,Folgueira Maria Aparecida A K Acta veterinaria Scandinavica BACKGROUND:In women with breast cancer submitted to neoadjuvant chemotherapy based in doxorubicin, tumor expression of groups of three genes (PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2) have classified them as responsive or resistant. We have investigated whether expression of these trios of genes could predict mammary carcinoma response in dogs and whether tumor slices, which maintain epithelial-mesenchymal interactions, could be used to evaluate drug response in vitro. METHODS:Tumors from 38 dogs were sliced and cultured with or without doxorubicin 1 muM for 24 h. Tumor cells were counted by two observers to establish a percentage variation in cell number, between slices. Based on these results, a reduction in cell number between treated and control samples > or = 21.7%, arbitrarily classified samples, as drug responsive. Tumor expression of PRSS11, MTSS1, CLPTM1 and SMYD2, was evaluated by real time PCR. Relative expression results were then transformed to their natural logarithm values, which were spatially disposed according to the expression of trios of genes, comprising PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. Fisher linear discrimination test was used to generate a separation plane between responsive and non-responsive tumors. RESULTS:Culture of tumor slices for 24 h was feasible. Nine samples were considered responsive and 29 non-responsive to doxorubicin, considering the pre-established cut-off value of cell number reduction > or = 21.7%, between doxorubicin treated and control samples. Relative gene expression was evaluated and tumor samples were then spatially distributed according to the expression of the trios of genes: PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2. A separation plane was generated. However, no clear separation between responsive and non-responsive samples could be observed. CONCLUSION:Three-dimensional distribution of samples according to the expression of the trios of genes PRSS11, MTSS1, CLPTM1 and PRSS11, MTSS1, SMYD2 could not predict doxorubicin in vitro responsiveness. Short term culture of mammary gland cancer slices may be an interesting model to evaluate chemotherapy activity. 10.1186/1751-0147-50-27
Head and Neck Squamous Cell Carcinomas Are Characterized by a Stable Immune Signature Within the Primary Tumor Over Time and Space. Clinical cancer research : an official journal of the American Association for Cancer Research Genetic and morphologic heterogeneity is well-documented in solid cancers. Immune cells are also variably distributed within the tumor; this heterogeneity is difficult to assess in small biopsies, and may confound our understanding of the determinants of successful immunotherapy. We examined the transcriptomic variability of the immunologic signature in head and neck squamous cell carcinoma (HNSCC) within individual tumors using transcriptomic and IHC assessments. Forty-four tumor biopsies from 16 HNSCC patients, taken at diagnosis and later at resection, were analyzed using RNA-sequencing. Variance filtering was used to identify the top 4,000 most variable genes. Principal component analysis, hierarchical clustering, and correlation analysis were performed. Gene expression of was correlated to IHC analysis. Analysis of immunologic gene expression was highly consistent in replicates from the same cancer. Across the cohort, samples from the same patient were most similar to each other, both spatially (at diagnosis) and, notably, over time (diagnostic biopsy compared with resection); comparison of global gene expression by hierarchical clustering ( ≤ 0.0001) and correlation analysis [median intrapatient = 0.82; median interpatient = 0.63]. gene transcript counts were highly correlated with CD8 T-cell counts by IHC ( = 0.82). Our data demonstrate that in HNSCC the global tumor and adaptive immune signatures are stable between discrete parts of the same tumor and also at different timepoints. This suggests that immunologic heterogeneity may not be a key reason for failure of immunotherapy and underpins the use of transcriptomics for immunologic evaluation of novel agents in HNSCC patients. . 10.1158/1078-0432.CCR-17-0373
Customized oligonucleotide array-based comparative genomic hybridization as a clinical assay for genomic profiling of chronic lymphocytic leukemia. The Journal of molecular diagnostics : JMD Chromosome gains and losses used for risk stratification in chronic lymphocytic leukemia (CLL) are commonly assessed by multiprobe fluorescence in situ hybridization (FISH) studies. We designed and validated a customized array-comparative genomic hybridization (aCGH) platform as a clinical assay for CLL genomic profiling. A 60-mer, 44,000-probe oligonucleotide array with a 50-kb average spatial resolution was augmented with high-density probe tiling at loci that are frequently aberrant in CLL. Aberrations identified by aCGH were compared with those identified by a FISH panel, including locus-specific probes to ATM (11q22.3), the centromeric region of chromosome 12 (12p11.1-q11), D13S319 (13q14.3), LAMP1 (13q34), and TP53 (17p13.1). In 100 CLL samples, aCGH/FISH concordance was seen for 89% of FISH-called aberrations at the ATM (n=18), D13S319 (n=42), LAMP (n=12), and TP53 (n=22) loci and for chromosome 12 (n=14). Eighty-four percentage of FISH/aCGH discordant calls were in samples either at or below the limit of aCGH sensitivity (10% to 25% FISH aberration-containing cells). Therefore, aCGH profiling is a feasible routine clinical test with comparable results to multiprobe FISH studies; however, it may be less sensitive than FISH in cases with low-level aberrations. Further, a customized array design can provide comprehensive genomic profiling with additional accuracy in both identifying and defining the extent of small aberrations at target loci. 10.2353/jmoldx.2009.080037
Microarray analysis of cutaneous squamous cell carcinomas reveals enhanced expression of epidermal differentiation complex genes. Hudson Laurie G,Gale James M,Padilla R Steven,Pickett Gavin,Alexander Bryan E,Wang Jing,Kusewitt Donna F Molecular carcinogenesis Gene expression profiles were determined for 12 cutaneous squamous cell carcinomas (SCC) removed from sun-exposed sites on nonimmunosuppressed patients. Gene expression in each SCC was compared to that in sun-exposed skin from the same patient using the Affymetrix HGU133 2.0 PlusGeneChip. We identified 440 genes with increased expression in SCC and 738 with decreased expression; overall we identified a large number of small changes in gene expression rather than a few marked changes that distinguished SCC from sun-exposed skin. Analyzing this robust data set according to biofunctional pathways using DAVID, transcriptional control elements using oPOSSUM, and chromosomal location using GSEA suggested genetic and epigenetic mechanisms of gene expression regulation in SCC. Some altered patterns of gene expression in SCC were consistent with regulation of spatially separated genes by a number of developmentally important transcription factors (forkhead, HMG, and homeo factors) that negatively regulated gene expression and to a few factors that positively regulated expression (Creb-1, NFkappaB, RelA, and Sp-1). We also found that coordinately enhanced expression of epidermal differentiation complex genes on chromosome 1q21 was a hallmark of SCC. A novel finding in our study was enhanced expression of keratin 13 in SCC, a result validated by immunohistochemical staining of an SCC tumor tissue array. 10.1002/mc.20636
Expression imbalance map: a new visualization method for detection of mRNA expression imbalance regions. Kano Makoto,Nishimura Kunihiro,Ishikawa Shumpei,Tsutsumi Shuichi,Hirota Koichi,Hirose Michitaka,Aburatani Hiroyuki Physiological genomics We describe the development of a new visualization method, called the expression imbalance map (EIM), for detecting mRNA expression imbalance regions, reflecting genomic losses and gains at a much higher resolution than conventional technologies such as comparative genomic hybridization (CGH). Simple spatial mapping of the microarray expression profiles on chromosomal location provides little information about genomic structure, because mRNA expression levels do not completely reflect genomic copy number and some microarray probes would be of low quality. The EIM, which does not employ arbitrary selection of thresholds in conjunction with hypergeometric distribution-based algorithm, has a high tolerance of these complex factors. The EIM could detect regionally underexpressed or overexpressed genes (called, here, an expression imbalance region) in lung cancer specimens from their gene expression data of oligonucleotide microarray. Many known as well as potential loci with frequent genomic losses or gains were detected as expression imbalance regions by the EIM. Therefore, the EIM should provide the user with further insight into genomic structure through mRNA expression. 10.1152/physiolgenomics.00116.2002
Variation in MicroRNA Expression Profile of Uterine Leiomyoma with Endometrial Cavity Distortion and Endometrial Cavity Non-Distortion. International journal of molecular sciences The expression profile of microRNA (miRNA) in uterine leiomyoma (UL) cells is different from that in normal uterine myometrial (UM) cells. The effect of UL cells on uterine receptivity might vary according to their ability to distort the uterine endometrial cavity. However, the variation in miRNA expression profiles between endometrial cavity-distorting leiomyoma (ECDL) and endometrial cavity non-distorting leiomyoma (ECNDL) cells remains unknown. This study aimed to elucidate whether the expression profile of miRNAs in ECDL cells is dissimilar to that of ECNDL cells in uterus. Pelviscopic myomectomy was performed to obtain tissue samples of UL and their corresponding normal UM tissues (matched) from patients with UL ( = 26), among whom women with ECNDL and ECDL numbered 15 and 11, respectively. The relative expression of hsa-miR-15b, -29a, -29b, -29c, -197, and -200c as well as the candidate target genes in UL cells was compared to those in the matched UM cells using qRT-PCR to assess their ability to cause ECD. The spatial expression of miRNAs and target genes in the UL tissues was analyzed using in situ hybridization. Target gene expression was analyzed using qPCR after transfection with the mimics and inhibitors of miRNAs in UL cells. The relative expression level of miR-15b was upregulated, and the relative expression levels of miR-29a, -29b, -29c, -197, and -200c were downregulated in UL cells compared to those in UM cells. The relative expression levels of progesterone receptor, estrogen receptor, and matrix metalloproteinases (MMPs) were upregulated in UL cells compared to those in UM cells. The relative expression levels of miR-29c and -200c were downregulated, and the relative expression levels of estrogen receptor, MMPs and tissue inhibitors of metalloproteinases (TIMPs) were upregulated in ECDL cells compared to those in ECNDL cells. The expression profile of miRNAs in UL cells varied with respect to the occurrence or absence of endometrial cavity distortion. The biochemical properties of UL might be regulated by miRNAs in order to alter their effect on structural homeostasis of the uterus. 10.3390/ijms19092524
Enriching the molecular definition of the airway "field of cancerization:" establishing new paradigms for the patient at risk for lung cancer. Gomperts Brigitte N,Walser Tonya C,Spira Avrum,Dubinett Steven M Cancer prevention research (Philadelphia, Pa.) The "field of cancerization" refers to histologically normal-appearing tissue adjacent to neoplastic tissue that displays molecular abnormalities, some of which are the same as those of the tumor. Improving our understanding of these molecular events is likely to increase our understanding of carcinogenesis. Kadara and colleagues attempt to characterize the molecular events occurring temporally and spatially within the field of cancerization of patients with early-stage non-small cell lung cancer (NSCLC) following definitive surgery. They followed patients with bronchoscopies annually after tumor resection and extracted RNA from the serial brushings from different endobronchial sites. They then conducted microarray analysis to identify gene expression differences over time and in different sites in the airway. Candidate genes were found that may have biologic relevance to the field of cancerization. For example, expression of phosphorylated AKT and ERK1/2 was found to increase in the airway epithelium with time. Although there are limitations in the study design, this investigation demonstrates the utility of identifying molecular changes in histologically normal airway epithelium in lung cancer. In addition to increasing our understanding of lung cancer biology, studying the field of cancerization has the potential to identify biomarkers from samples obtained in a minimally invasive manner. 10.1158/1940-6207.CAPR-12-0470
Probing the transcriptome of Boehmeria nivea reveals candidate genes associated with the biosynthesis of chlorogenic acid. Gene Boehmeria nivea (L.) Gaudich is used in traditional Chinese medicine. Chlorogenic acids are major medically active components of Boehmeria nivea, which can be used clinically to treat hyperglycemia, pneumonia, and cancer. To identify the genes involved in chlorogenic acid biosynthesis, we analyzed transcriptome data from leaf, root, and stem tissues of Boehmeria nivea using the Illumina Hi-Seq 4000 platform. A total of 146,790 unigenes were obtained from Boehmeria nivea, of which 106,786 were annotated in public databases. In analyses of the KEGG (Kyoto Encyclopedia of Genes and Genome) database, 484 unigenes that encode the five key enzymes involved in chlorogenic acid biosynthesis were identified, and shikimate O-hydroxycinnamoyl transferase was spatially simulated. Some of these key enzyme unigenes expression levels were verified by RT-qPCR (real-time quantitative Polymerase Chain Reaction). Furthermore, multiple genes encoding plant resistance proteins or transcription factors were identified and analyzed. Differentially expressed genes were identified by performing pairwise comparison of genes between tissues. This study increases the number of public transcript datasets of this species and identifies candidate genes related to the biosynthesis of chlorogenic acid, laying a foundation for the further exploration of this pathway in Boehmeria nivea. 10.1016/j.gene.2022.146579
Use of gene expression profiling to direct in vivo molecular imaging of lung cancer. Grimm Jan,Kirsch David G,Windsor Stephen D,Kim Carla F Bender,Santiago Philip M,Ntziachristos Vasilis,Jacks Tyler,Weissleder Ralph Proceedings of the National Academy of Sciences of the United States of America Using gene expression profiling, we identified cathepsin cysteine proteases as highly up-regulated genes in a mouse model of human lung adenocarcinoma. Overexpression of cathepsin proteases in these lung tumors was confirmed by immunohistochemistry and Western blotting. Therefore, an optical probe activated by cathepsin proteases was selected to detect murine lung tumors in vivo as small as 1 mm in diameter and spatially separated. We generated 3D maps of the fluorescence signal and fused them with anatomical computed tomography images to show a close correlation between fluorescence signal and tumor burden. By serially imaging the same mouse, optical imaging was used to follow tumor progression. This study demonstrates the capability for molecular imaging of a primary lung tumor by using endogenous proteases expressed by a tumor. It also highlights the feasibility of using gene expression profiling to identify molecular targets for imaging lung cancer. 10.1073/pnas.0503920102
Combined gene expression and protein interaction analysis of dynamic modularity in glioma prognosis. Zhang Xiaoyu,Yang Hongbin,Gong Binsheng,Jiang Chuanlu,Yang Lizhuang Journal of neuro-oncology Because of the variety of factors affecting glioma prognosis, prediction of patient survival is particularly difficult. Protein-protein interaction (PPI) networks have been considered with regard to how their spatial characteristics relate to glioma. However, the dynamic nature of PPIs in vivo makes them temporally and spatially complex events. Integration of prognosis-specific co-expression information adds further dynamic features to these networks. Although some biomarkers for glioma prognosis have been identified, none is sufficient for accurate prediction of either prognosis or improved survival. We have established co-expressed protein-interaction networks that integrate protein-protein interactions with glioma gene-expression profiles related to different survival times. Biomarkers related to glioma prognosis were identified by comparative analysis of the dynamic features of the glioma prognosis network, particularly subnetworks. Four significantly differently expressed genes (SDEGs) are upregulated and ten SDEGs downregulated as lifetime is extended. In addition, 97 enhanced differently co-expressed protein interactions (DCPIs) and 99 weakened DCPIs were associated with glioma patient lifetime extension. We propose a method for estimating glioma prognosis on the basis of the construction of a dynamic modular network. We have used this method to identify dynamic genes and interactions related to glioma prognosis. Among these, enhanced MYC expression was related to lifetime extension, as were interactions between E2F1 and RB1 and between EGFR and p38. This method is a novel means of studying the molecular mechanisms determining prognosis in glioma. 10.1007/s11060-011-0757-4
Distinct evolutionary trajectories of primary high-grade serous ovarian cancers revealed through spatial mutational profiling. Bashashati Ali,Ha Gavin,Tone Alicia,Ding Jiarui,Prentice Leah M,Roth Andrew,Rosner Jamie,Shumansky Karey,Kalloger Steve,Senz Janine,Yang Winnie,McConechy Melissa,Melnyk Nataliya,Anglesio Michael,Luk Margaret T Y,Tse Kane,Zeng Thomas,Moore Richard,Zhao Yongjun,Marra Marco A,Gilks Blake,Yip Stephen,Huntsman David G,McAlpine Jessica N,Shah Sohrab P The Journal of pathology High-grade serous ovarian cancer (HGSC) is characterized by poor outcome, often attributed to the emergence of treatment-resistant subclones. We sought to measure the degree of genomic diversity within primary, untreated HGSCs to examine the natural state of tumour evolution prior to therapy. We performed exome sequencing, copy number analysis, targeted amplicon deep sequencing and gene expression profiling on 31 spatially and temporally separated HGSC tumour specimens (six patients), including ovarian masses, distant metastases and fallopian tube lesions. We found widespread intratumoural variation in mutation, copy number and gene expression profiles, with key driver alterations in genes present in only a subset of samples (eg PIK3CA, CTNNB1, NF1). On average, only 51.5% of mutations were present in every sample of a given case (range 10.2-91.4%), with TP53 as the only somatic mutation consistently present in all samples. Complex segmental aneuploidies, such as whole-genome doubling, were present in a subset of samples from the same individual, with divergent copy number changes segregating independently of point mutation acquisition. Reconstruction of evolutionary histories showed one patient with mixed HGSC and endometrioid histology, with common aetiologic origin in the fallopian tube and subsequent selection of different driver mutations in the histologically distinct samples. In this patient, we observed mixed cell populations in the early fallopian tube lesion, indicating that diversity arises at early stages of tumourigenesis. Our results revealed that HGSCs exhibit highly individual evolutionary trajectories and diverse genomic tapestries prior to therapy, exposing an essential biological characteristic to inform future design of personalized therapeutic solutions and investigation of drug-resistance mechanisms. 10.1002/path.4230
Quantitative nuclear histomorphometric features are predictive of Oncotype DX risk categories in ductal carcinoma in situ: preliminary findings. Li Haojia,Whitney Jon,Bera Kaustav,Gilmore Hannah,Thorat Mangesh A,Badve Sunil,Madabhushi Anant Breast cancer research : BCR BACKGROUND:Oncotype DX (ODx) is a 12-gene assay assessing the recurrence risk (high, intermediate, and low) of ductal carcinoma in situ (pre-invasive breast cancer), which guides clinicians regarding prescription of radiotherapy. However, ODx is expensive, time-consuming, and tissue-destructive. In addition, the actual prognostic meaning for the intermediate ODx risk category remains unclear. METHODS:In this work, we evaluated the ability of quantitative nuclear histomorphometric features extracted from hematoxylin and eosin-stained slide images of 62 ductal carcinoma in situ (DCIS) patients to distinguish between the corresponding ODx risk categories. The prognostic value of the identified image signature was further evaluated on an independent validation set of 30 DCIS patients in its ability to distinguish those DCIS patients who progressed to invasive carcinoma versus those who did not. Following nuclear segmentation and feature extraction, feature ranking strategies were employed to identify the most discriminating features between individual ODx risk categories. The selected features were then combined with machine learning classifiers to establish models to predict ODx risk categories. The model performance was evaluated using the average area under the receiver operating characteristic curve (AUC) using cross validation. In addition, an unsupervised clustering approach was also implemented to evaluate the ability of nuclear histomorphometric features to discriminate between the ODx risk categories. RESULTS:Features relating to spatial distribution, orientation disorder, and texture of nuclei were identified as most discriminating between the high ODx and the intermediate, low ODx risk categories. Additionally, the AUC of the most discriminating set of features for the different classification tasks was as follows: (1) high vs low ODx (0.68), (2) high vs. intermediate ODx (0.67), (3) intermediate vs. low ODx (0.57), (4) high and intermediate vs. low ODx (0.63), (5) high vs. low and intermediate ODx (0.66). Additionally, the unsupervised clustering resulted in intermediate ODx risk category patients being co-clustered with low ODx patients compared to high ODx. CONCLUSION:Our results appear to suggest that nuclear histomorphometric features can distinguish high from low and intermediate ODx risk category patients. Additionally, our findings suggest that histomorphometric features for intermediate ODx were more similar to low ODx compared to high ODx risk category. 10.1186/s13058-019-1200-6
In silico dissection of cell-type-associated patterns of gene expression in prostate cancer. Stuart Robert O,Wachsman William,Berry Charles C,Wang-Rodriguez Jessica,Wasserman Linda,Klacansky Igor,Masys Dan,Arden Karen,Goodison Steven,McClelland Michael,Wang Yipeng,Sawyers Anne,Kalcheva Iveta,Tarin David,Mercola Dan Proceedings of the National Academy of Sciences of the United States of America Prostate tumors are complex entities composed of malignant cells mixed and interacting with nonmalignant cells. However, molecular analyses by standard gene expression profiling are limited because spatial information and nontumor cell types are lost in sample preparation. We scored 88 prostate specimens for relative content of tumor, benign hyperplastic epithelium, stroma, and dilated cystic glands. The proportions of these cell types were then linked in silico to gene expression levels determined by microarray analysis, revealing unique cell-specific profiles. Gene expression differences for malignant and nonmalignant epithelial cells (tumor versus benign hyperplastic epithelium) could be identified without being confounded by contributions from stroma that dominate many samples or sacrificing possible paracrine influences. Cell-specific expression of selected genes was validated by immunohistochemistry and quantitative PCR. The results provide patterns of gene expression for these three lineages with relevance to pathogenetic, diagnostic, and therapeutic considerations. 10.1073/pnas.2536479100
Quantification of differential ErbB1 and ErbB2 cell surface expression and spatial nanoclustering through plasmon coupling. Nano letters Cell surface receptors play ubiquitous roles in cell signaling and communication and their expression levels are important biomarkers for many diseases. Expression levels are, however, only one factor that determines the physiological activity of a receptor. For some surface receptors, their distribution on the cell surface, especially their clustering, provides additional mechanisms for regulation. To access this spatial information robust assays are required that provide detailed insight into the organization of cell surface receptors on nanometer length scales. In this manuscript, we demonstrate through combination of scattering spectroscopy, electron microscopy, and generalized multiple particle Mie theory (GMT) simulations that the density- and morphology-dependent spectral response of Au nanoparticle (NP) immunolabels bound to the epidermal growth factor receptors ErbB1 and ErbB2 encodes quantitative information of both the cell surface expression and spatial clustering of the two receptors in different unliganded in vitro cancer cell lines (SKBR3, MCF7, A431). A systematic characterization of the collective spectral responses of NPs targeted at ErbB1 and ErbB2 at various NP concentrations indicates differences in the large-scale organization of ErbB1 and ErbB2 in cell lines that overexpress these receptors. Validation experiments in the scanning electron microscope (SEM) confirm that NPs targeted at ErbB1 on A431 are more strongly clustered than NPs bound to ErbB2 on SKBR3 or MCF7 at overall comparable NP surface densities. This finding is consistent with the existence of larger receptor clusters for ErbB1 than for ErbB2 in the plasma membranes of the respective cells. 10.1021/nl3012227
Preclinical applications of imaging for cancer gene therapy. Briat Arnaud,Vassaux Georges Expert reviews in molecular medicine Gene therapy is a very attractive strategy in experimental cancer therapy. Ideally, the approach aims to deliver therapeutic genes selectively to cancer cells. However, progress in the improvement of gene therapy formulations has been hampered by difficulties in measuring transgene delivery and in quantifying transgene expression in vivo. In clinical trials, endpoints rely almost exclusively on the analysis of biopsies by molecular and histopathological methods, which provide limited information. Therefore, to ensure the rational development of gene therapy, a crucial issue is the utilisation of technologies for the non-invasive monitoring of spatial and temporal gene expression in vivo upon administration of a gene delivery vector. Such imaging technologies would allow the generation of quantitative information about gene expression and the assessment of cancer gene therapy efficacy. In the past decade, progress has been made in the field of in vivo molecular imaging. This review highlights the various methods currently being developed in preclinical models. 10.1017/S1462399406000044
Global gene expression profiling of healthy human brain and its application in studying neurological disorders. Negi Simarjeet K,Guda Chittibabu Scientific reports Brain function is governed by precise regulation of gene expression across its anatomically distinct structures; however, the expression patterns of genes across hundreds of brain structures are not clearly understood. Here, we describe a gene expression model, which is representative of the healthy human brain transcriptome by using data from the Allen Brain Atlas. Our in-depth gene expression profiling revealed that 84% of genes are expressed in at least one of the 190 brain structures studied. Hierarchical clustering based on gene expression profiles delineated brain regions into structurally tiered spatial groups and we observed striking enrichment for region-specific processes. Further, weighted co-expression network analysis identified 19 robust modules of highly correlated genes enriched with functional associations for neurogenesis, dopamine signaling, immune regulation and behavior. Also, structural distribution maps of major neurotransmission systems in the brain were generated. Finally, we developed a supervised classification model, which achieved 84% and 81% accuracies for predicting autism- and Parkinson's-implicated genes, respectively, using our expression model as a baseline. This study represents the first use of global gene expression profiling from healthy human brain to develop a disease gene prediction model and this generic methodology can be applied to study any neurological disorder. 10.1038/s41598-017-00952-9
Quantification of intrinsic subtype ambiguity in Luminal A breast cancer and its relationship to clinical outcomes. Kumar Neeraj,Zhao Dan,Bhaumik Dulal,Sethi Amit,Gann Peter H BMC cancer BACKGROUND:PAM50 gene profiling assigns each cancer to a single intrinsic subtype. However, individual cancers vary in their adherence to a prototype, and due to bulk tissue sampling, some may exhibit expression patterns that indicate intra-tumor admixture of multiple subtypes. Our objective was to develop admixture metrics from PAM50 gene expression profiles in order to stratify Luminal A (LumA) cases according to their degree of subtype admixture, and then relate such admixture to clinical and molecular variables. METHODS:We re-constructed scaled, normalized PAM50 profiles for 1980 cases (674 LumA) in the METABRIC cohort and for each case computed its Mahalanobis (M-) distance from its assigned centroid and M-distance from all other centroids. We used t-SNE plots to visualize overlaps in subtype clustering. With Normal-like cases excluded, we developed two metrics: Median Distance Criteria (MDC) classified pure cases as those located within the 50th percentile of the LumA centroid and > =50th percentile from any other centroid. Distance Ratio Criteria (DRC) was computed as the ratio of M-distances from the LumA centroid to the nearest non-assigned centroid. Pure and admixed LumA cases were compared on clinical/molecular traits. TCGA LumA cases (n = 509) provided independent validation. RESULTS:Compared to pure cases in METABRIC, admixed ones had older age at diagnosis, larger tumor size, and higher grade and stage. These associations were stronger for the DRC metric compared to MDC. Admixed cases were associated with HER2 gain, high proliferation, higher PAM50 recurrence scores, more frequent TP53 mutation, and less frequent PIK3CA mutation. Similar results were observed in the TCGA validation cohort, which also showed a positive association between admixture and number of clonal populations estimated by PyClone. LumA-LumB confusion predominated, but other combinations were also present. Degree of admixture was associated with overall survival in both cohorts, as was disease-free survival in TCGA, independent of age, grade and stage (HR = 2.85, Tertile 3 vs.1). CONCLUSIONS:Luminal A breast cancers subgrouped based on PAM50 subtype purity support the hypothesis that admixed cases have worse clinical features and survival. Future analyses will explore more extensive genomic metrics for admixture and their spatial significance within a single tumor. 10.1186/s12885-019-5392-z
Mapping genomic and transcriptomic alterations spatially in epithelial cells adjacent to human breast carcinoma. Abdalla Moustafa,Tran-Thanh Danh,Moreno Juan,Iakovlev Vladimir,Nair Ranju,Kanwar Nisha,Abdalla Mohamed,Lee Jennifer P Y,Kwan Jennifer Yin Yee,Cawthorn Thomas R,Warren Keisha,Arneson Nona,Wang Dong-Yu,Fox Natalie S,Youngson Bruce J,Miller Naomi A,Easson Alexandra M,McCready David,Leong Wey L,Boutros Paul C,Done Susan J Nature communications Almost all genomic studies of breast cancer have focused on well-established tumours because it is technically challenging to study the earliest mutational events occurring in human breast epithelial cells. To address this we created a unique dataset of epithelial samples ductoscopically obtained from ducts leading to breast carcinomas and matched samples from ducts on the opposite side of the nipple. Here, we demonstrate that perturbations in mRNA abundance, with increasing proximity to tumour, cannot be explained by copy number aberrations. Rather, we find a possibility of field cancerization surrounding the primary tumour by constructing a classifier that evaluates where epithelial samples were obtained relative to a tumour (cross-validated micro-averaged AUC = 0.74). We implement a spectral co-clustering algorithm to define biclusters. Relating to over-represented bicluster pathways, we further validate two genes with tissue microarrays and in vitro experiments. We highlight evidence suggesting that bicluster perturbation occurs early in tumour development. 10.1038/s41467-017-01357-y
A novel machine learning approach reveals latent vascular phenotypes predictive of renal cancer outcome. Scientific reports Gene expression signatures are commonly used as predictive biomarkers, but do not capture structural features within the tissue architecture. Here we apply a 2-step machine learning framework for quantitative imaging of tumor vasculature to derive a spatially informed, prognostic gene signature. The trained algorithms classify endothelial cells and generate a vascular area mask (VAM) in H&E micrographs of clear cell renal cell carcinoma (ccRCC) cases from The Cancer Genome Atlas (TCGA). Quantification of VAMs led to the discovery of 9 vascular features (9VF) that predicted disease-free-survival in a discovery cohort (n = 64, HR = 2.3). Correlation analysis and information gain identified a 14 gene expression signature related to the 9VF's. Two generalized linear models with elastic net regularization (14VF and 14GT), based on the 14 genes, separated independent cohorts of up to 301 cases into good and poor disease-free survival groups (14VF HR = 2.4, 14GT HR = 3.33). For the first time, we successfully applied digital image analysis and targeted machine learning to develop prognostic, morphology-based, gene expression signatures from the vascular architecture. This novel morphogenomic approach has the potential to improve previous methods for biomarker development. 10.1038/s41598-017-13196-4
Transcriptomics and Spatial Proteomics for Discovery and Validation of Missing Proteins in the Human Ovary. Journal of proteome research Efforts to understand the complexities of human biology encompass multidimensional aspects, with proteins emerging as crucial components. However, studying the human ovary introduces unique challenges due to its complex dynamics and changes over a lifetime, varied cellular composition, and limited sample access. Here, four new RNA-seq samples of ovarian cortex spanning ages of 7 to 32 were sequenced and added to the existing data in the Human Protein Atlas (HPA) database www.proteinatlas.org, opening the doors to unique possibilities for exploration of oocyte-specific proteins. Based on transcriptomics analysis of the four new tissue samples representing both prepubertal girls and women of fertile age, we selected 20 protein candidates that lacked previous evidence at the protein level, so-called "missing proteins" (MPs). The proteins were validated using high-resolution antibody-based profiling and single-cell transcriptomics. Fourteen proteins exhibited consistent single-cell expression patterns in oocytes and granulosa cells, confirming their presence in the ovary and suggesting that these proteins play important roles in ovarian function, thus proposing that these 14 proteins should no longer be classified as MPs. This research significantly advances the understanding of MPs, unearthing fresh avenues for prospective exploration. By integrating innovative methodologies and leveraging the wealth of data in the HPA database, these insights contribute to refining our understanding of protein roles within the human ovary and opening the doors for further investigations into missing proteins and human reproduction. 10.1021/acs.jproteome.3c00545
Single-cell transcriptomics identifies Keap1-Nrf2 regulated collective invasion in a tumor model. eLife Apicobasal cell polarity loss is a founding event in epithelial-mesenchymal transition and epithelial tumorigenesis, yet how pathological polarity loss links to plasticity remains largely unknown. To understand the mechanisms and mediators regulating plasticity upon polarity loss, we performed single-cell RNA sequencing of ovaries, where inducing polarity-gene -knockdown (Lgl-KD) causes invasive multilayering of the follicular epithelia. Analyzing the integrated Lgl-KD and transcriptomes, we discovered the cells specific to the various discernible phenotypes and characterized the underlying gene expression. A genetic requirement of Keap1-Nrf2 signaling in promoting multilayer formation of Lgl-KD cells was further identified. Ectopic expression of Keap1 increased the volume of delaminated follicle cells that showed enhanced invasive behavior with significant changes to the cytoskeleton. Overall, our findings describe the comprehensive transcriptome of cells within the follicle cell tumor model at the single-cell resolution and identify a previously unappreciated link between Keap1-Nrf2 signaling and cell plasticity at early tumorigenesis. 10.7554/eLife.80956
Quantitative image analysis of cellular heterogeneity in breast tumors complements genomic profiling. Yuan Yinyin,Failmezger Henrik,Rueda Oscar M,Ali H Raza,Gräf Stefan,Chin Suet-Feung,Schwarz Roland F,Curtis Christina,Dunning Mark J,Bardwell Helen,Johnson Nicola,Doyle Sarah,Turashvili Gulisa,Provenzano Elena,Aparicio Sam,Caldas Carlos,Markowetz Florian Science translational medicine Solid tumors are heterogeneous tissues composed of a mixture of cancer and normal cells, which complicates the interpretation of their molecular profiles. Furthermore, tissue architecture is generally not reflected in molecular assays, rendering this rich information underused. To address these challenges, we developed a computational approach based on standard hematoxylin and eosin-stained tissue sections and demonstrated its power in a discovery and validation cohort of 323 and 241 breast tumors, respectively. To deconvolute cellular heterogeneity and detect subtle genomic aberrations, we introduced an algorithm based on tumor cellularity to increase the comparability of copy number profiles between samples. We next devised a predictor for survival in estrogen receptor-negative breast cancer that integrated both image-based and gene expression analyses and significantly outperformed classifiers that use single data types, such as microarray expression signatures. Image processing also allowed us to describe and validate an independent prognostic factor based on quantitative analysis of spatial patterns between stromal cells, which are not detectable by molecular assays. Our quantitative, image-based method could benefit any large-scale cancer study by refining and complementing molecular assays of tumor samples. 10.1126/scitranslmed.3004330
Intra-tumour heterogeneity of diffuse large B-cell lymphoma involves the induction of diversified stroma-tumour interfaces. Sangaletti Sabina,Iannelli Fabio,Zanardi Federica,Cancila Valeria,Portararo Paola,Botti Laura,Vacca Davide,Chiodoni Claudia,Di Napoli Arianna,Valenti Cesare,Rizzello Celeste,Vegliante Maria Carmela,Pisati Federica,Gulino Alessandro,Ponzoni Maurilio,Colombo Mario Paolo,Tripodo Claudio EBioMedicine BACKGROUND:Intra-tumour heterogeneity in lymphoid malignancies encompasses selection of genetic events and epigenetic regulation of transcriptional programs. Clonal-related neoplastic cell populations are unsteadily subjected to immune editing and metabolic adaptations within different tissue microenvironments. How tissue-specific mesenchymal cells impact on the diversification of aggressive lymphoma clones is still unknown. METHODS:Combining in situ quantitative immunophenotypical analyses and RNA sequencing we investigated the intra-tumour heterogeneity and the specific mesenchymal modifications that are associated with A20 diffuse large B-cell lymphoma (DLBCL) cells seeding of different tissue microenvironments. Furthermore, we characterized features of lymphoma-associated stromatogenesis in human DLBCL samples using Digital Spatial Profiling, and established their relationship with prognostically relevant variables, such as MYC. FINDINGS:We found that the tissue microenvironment casts a relevant influence over A20 transcriptional landscape also impacting on Myc and DNA damage response programs. Extending the investigation to mice deficient for the matricellular protein SPARC, a stromal prognostic factor in human DLBCL, we demonstrated a different immune imprint on A20 cells according to stromal Sparc proficiency. Through Digital Spatial Profiling of 87 immune and stromal genes on human nodal DLBCL regions characterized by different mesenchymal composition, we demonstrate intra-lesional heterogeneity arising from diversified mesenchymal contextures and impacting on the stromal and immune milieu. INTERPRETATION:Our study provides experimental evidence that stromal microenvironment generates topological determinants of intra-tumour heterogeneity in DLBCL involving key transcriptional pathways such as Myc expression, damage response programs and immune checkpoints. FUNDING:This study has been supported by the Italian Foundation for Cancer Research (AIRC) (grants 15999 and 22145 to C. Tripodo) and by the University of Palermo. 10.1016/j.ebiom.2020.103055
Hyper-Graph Regularized Constrained NMF for Selecting Differentially Expressed Genes and Tumor Classification. Jiao Cui-Na,Gao Ying-Lian,Yu Na,Liu Jin-Xing,Qi Lian-Yong IEEE journal of biomedical and health informatics Non-negative Matrix Factorization (NMF) is a dimensionality reduction approach for learning a parts-based and linear representation of non-negative data. It has attracted more attention because of that. In practice, NMF not only neglects the manifold structure of data samples, but also overlooks the priori label information of different classes. In this paper, a novel matrix decomposition method called Hyper-graph regularized Constrained Non-negative Matrix Factorization (HCNMF) is proposed for selecting differentially expressed genes and tumor sample classification. The advantage of hyper-graph learning is to capture local spatial information in high dimensional data. This method incorporates a hyper-graph regularization constraint to consider the higher order data sample relationships. The application of hyper-graph theory can effectively find pathogenic genes in cancer datasets. Besides, the label information is further incorporated in the objective function to improve the discriminative ability of the decomposition matrix. Supervised learning with label information greatly improves the classification effect. We also provide the iterative update rules and convergence proofs for the optimization problems of HCNMF. Experiments under The Cancer Genome Atlas (TCGA) datasets confirm the superiority of HCNMF algorithm compared with other representative algorithms through a set of evaluations. 10.1109/JBHI.2020.2975199
Cell type specificity of signaling: view from membrane receptors distribution and their downstream transduction networks. He Ying,Yu Zhonghao,Ge Dongya,Wang-Sattler Rui,Thiesen Hans-Jürgen,Xie Lu,Li Yixue Protein & cell Studies on cell signaling pay more attention to spatial dynamics and how such diverse organization can relate to high order of cellular capabilities. To overview the specificity of cell signaling, we integrated human receptome data with proteome spatial expression profiles to systematically investigate the specificity of receptors and receptor-triggered transduction networks across 62 normal cell types and 14 cancer types. Six percent receptors showed cell-type-specific expression, and 4% signaling networks presented enriched cell-specific proteins induced by the receptors. We introduced a concept of "response context" to annotate the cell-type dependent signaling networks. We found that most cells respond similarly to the same stimulus, as the "response contexts" presented high functional similarity. Despite this, the subtle spatial diversity can be observed from the difference in network architectures. The architecture of the signaling networks in nerve cells displayed less completeness than that in glandular cells, which indicated cellular-context dependent signaling patterns are elaborately spatially organized. Likewise, in cancer cells most signaling networks were generally dysfunctional and less complete than that in normal cells. However, glioma emerged hyper-activated transduction mechanism in malignant state. Receptor ATP6AP2 and TNFRSF21 induced rennin-angiotensin and apoptosis signaling were found likely to explain the glioma-specific mechanism. This work represents an effort to decipher context-specific signaling network from spatial dimension. Our results indicated that although a majority of cells engage general signaling response with subtle differences, the spatial dynamics of cell signaling can not only deepen our insights into different signaling mechanisms, but also help understand cell signaling in disease. 10.1007/s13238-012-2049-y
Activity of distinct growth factor receptor network components in breast tumors uncovers two biologically relevant subtypes. Genome medicine BACKGROUND:The growth factor receptor network (GFRN) plays a significant role in driving key oncogenic processes. However, assessment of global GFRN activity is challenging due to complex crosstalk among GFRN components, or pathways, and the inability to study complex signaling networks in patient tumors. Here, pathway-specific genomic signatures were used to interrogate GFRN activity in breast tumors and the consequent phenotypic impact of GRFN activity patterns. METHODS:Novel pathway signatures were generated in human primary mammary epithelial cells by overexpressing key genes from GFRN pathways (HER2, IGF1R, AKT1, EGFR, KRAS (G12V), RAF1, BAD). The pathway analysis toolkit Adaptive Signature Selection and InteGratioN (ASSIGN) was used to estimate pathway activity for GFRN components in 1119 breast tumors from The Cancer Genome Atlas (TCGA) and across 55 breast cancer cell lines from the Integrative Cancer Biology Program (ICBP43). These signatures were investigated for their relationship to pro- and anti-apoptotic protein expression and drug response in breast cancer cell lines. RESULTS:Application of these signatures to breast tumor gene expression data identified two novel discrete phenotypes characterized by concordant, aberrant activation of either the HER2, IGF1R, and AKT pathways ("the survival phenotype") or the EGFR, KRAS (G12V), RAF1, and BAD pathways ("the growth phenotype"). These phenotypes described a significant amount of the variability in the total expression data across breast cancer tumors and characterized distinctive patterns in apoptosis evasion and drug response. The growth phenotype expressed lower levels of BIM and higher levels of MCL-1 proteins. Further, the growth phenotype was more sensitive to common chemotherapies and targeted therapies directed at EGFR and MEK. Alternatively, the survival phenotype was more sensitive to drugs inhibiting HER2, PI3K, AKT, and mTOR, but more resistant to chemotherapies. CONCLUSIONS:Gene expression profiling revealed a bifurcation pattern in GFRN activity represented by two discrete phenotypes. These phenotypes correlate to unique mechanisms of apoptosis and drug response and have the potential of pinpointing targetable aberration(s) for more effective breast cancer treatments. 10.1186/s13073-017-0429-x
Identification of transcriptome alterations in the prefrontal cortex, hippocampus, amygdala and hippocampus of suicide victims. Glavan Daniela,Gheorman Victor,Gresita Andrei,Hermann Dirk M,Udristoiu Ion,Popa-Wagner Aurel Scientific reports Suicide is one of the leading causes of death globally for all ages, and as such presents a very serious problem for clinicians worldwide. However, the underlying neurobiological pathology remains to a large extent unknown. In order to address this gap, we have carried out a genome-wide investigation of the gene expression in the amygdala, hippocampus, prefrontal cortex and thalamus in post-mortem brain samples obtained from 20 suicide completers and 7 control subjects. By KEGG enrichment analysis indicated we identified novel clusters of downregulated pathways involved in antigen neutralization and autoimmune thyroid disease (amygdala, thalamus), decreased axonal plasticity in the hippocampus. Two upregulated pathways were involved in neuronal death in the hippocampus and olfactory transduction in the thalamus and the prefrontal cortex. Autoimmune thyroid disease pathway was downregulated only in females. Metabolic pathways involved in Notch signaling amino acid metabolism and unsaturated lipid synthesis were thalamus-specific. Suicide-associated changes in the expression of several genes and pseudogenes that point to various functional mechanisms possibly implicated in the pathology of suicide. Two genes (SNORA13 and RNU4-2) involved in RNA processing were common to all brain regions analyzed. Most of the identified gene expression changes were related to region-specific dysregulated manifestation of genetic and epigenetic mechanisms underlying neurodevelopmental disorders (SNORD114-10, SUSd1), motivation, addiction and motor disorders (CHRNA6), long-term depression (RAB3B), stress response, major depression and schizophrenia (GFAP), signal transduction at the neurovascular unit (NEXN) and inhibitory neurotransmission in spatial learning, neural plasticity (CALB2; CLIC6, ENPP1). Some of the differentially expressed genes were brain specific non-coding RNAs involved in the regulation of translation (SNORA13). One, (PARM1) is a potential oncogene and prognostic biomarker for colorectal cancer with no known function in the brain. Disturbed gene expression involved in antigen neutralization, autoimmunity, neural plasticity, stress response, signal transduction at the neurovascular unit, dysregulated nuclear RNA processing and translation and epigenetic imprinting signatures is associated with suicide and point to regulatory non-coding RNAs as potential targets of new drugs development. 10.1038/s41598-021-98210-6
Multi-region sequencing unveils novel actionable targets and spatial heterogeneity in esophageal squamous cell carcinoma. Yan Ting,Cui Heyang,Zhou Yong,Yang Bin,Kong Pengzhou,Zhang Yingchun,Liu Yiqian,Wang Bin,Cheng Yikun,Li Jiayi,Guo Shixing,Xu Enwei,Liu Huijuan,Cheng Caixia,Zhang Ling,Chen Ling,Zhuang Xiaofei,Qian Yu,Yang Jian,Ma Yanchun,Li Hongyi,Wang Fang,Liu Jing,Liu Xuefeng,Su Dan,Wang Yan,Sun Ruifang,Guo Shiping,Li Yaoping,Cheng Xiaolong,Liu Zhihua,Zhan Qimin,Cui Yongping Nature communications Esophageal squamous cell carcinoma (ESCC) ranks fourth among cancer-related deaths in China due to the lack of actionable molecules. We performed whole-exome and T-cell receptor (TCR) repertoire sequencing on multi-regional tumors, normal tissues and blood samples from 39 ESCC patients. The data revealed 12.8% of ERBB4 mutations at patient level and functional study supported its oncogenic role. 18% of patients with early BRCA1/2 variants were associated with high-level contribution of signature 3, which was validated in an independent large cohort (n = 508). Furthermore, knockdown of BRCA1/2 dramatically increased sensitivity to cisplatin in ESCC cells. 5% of patients harbored focal high-level amplification of CD274 that led to massive expression of PD-L1, and might be more sensitive to immune checkpoint blockade. Finally, we found a tight correlation between genomic and TCR repertoire intra-tumor heterogeneity (ITH). Collectively, we reveal high-level ITH in ESCC, identify several potential actionable targets and may provide novel insight into ESCC treatment. 10.1038/s41467-019-09255-1
Multiplatform genomic profiling and magnetic resonance imaging identify mechanisms underlying intratumor heterogeneity in meningioma. Nature communications Meningiomas are the most common primary intracranial tumors, but the molecular drivers of meningioma tumorigenesis are poorly understood. We hypothesized that investigating intratumor heterogeneity in meningiomas would elucidate biologic drivers and reveal new targets for molecular therapy. To test this hypothesis, here we perform multiplatform molecular profiling of 86 spatially-distinct samples from 13 human meningiomas. Our data reveal that regional alterations in chromosome structure underlie clonal transcriptomic, epigenomic, and histopathologic signatures in meningioma. Stereotactic co-registration of sample coordinates to preoperative magnetic resonance images further suggest that high apparent diffusion coefficient (ADC) distinguishes meningioma regions with proliferating cells enriched for developmental gene expression programs. To understand the function of these genes in meningioma, we develop a human cerebral organoid model of meningioma and validate the high ADC marker genes CDH2 and PTPRZ1 as potential targets for meningioma therapy using live imaging, single cell RNA sequencing, CRISPR interference, and pharmacology. 10.1038/s41467-020-18582-7
Gene expression profiling of periodontitis-affected gingival tissue by spatial transcriptomics. Lundmark Anna,Gerasimcik Natalija,Båge Tove,Jemt Anders,Mollbrink Annelie,Salmén Fredrik,Lundeberg Joakim,Yucel-Lindberg Tülay Scientific reports Periodontitis is a highly prevalent chronic inflammatory disease of the periodontium, leading ultimately to tooth loss. In order to characterize the gene expression of periodontitis-affected gingival tissue, we have here simultaneously quantified and localized gene expression in periodontal tissue using spatial transcriptomics, combining RNA sequencing with histological analysis. Our analyses revealed distinct clusters of gene expression, which were identified to correspond to epithelium, inflamed areas of connective tissue, and non-inflamed areas of connective tissue. Moreover, 92 genes were identified as significantly up-regulated in inflamed areas of the gingival connective tissue compared to non-inflamed tissue. Among these, immunoglobulin lambda-like polypeptide 5 (IGLL5), signal sequence receptor subunit 4 (SSR4), marginal zone B and B1 cell specific protein (MZB1), and X-box binding protein 1 (XBP1) were the four most highly up-regulated genes. These genes were also verified as significantly higher expressed in gingival tissue of patients with periodontitis compared to healthy controls, using reverse transcription quantitative polymerase chain reaction. Moreover, the protein expressions of up-regulated genes were verified in gingival biopsies by immunohistochemistry. In summary, in this study, we report distinct gene expression signatures within periodontitis-affected gingival tissue, as well as specific genes that are up-regulated in inflamed areas compared to non-inflamed areas of gingival tissue. The results obtained from this study may add novel information on the genes and cell types contributing to pathogenesis of the chronic inflammatory disease periodontitis. 10.1038/s41598-018-27627-3
Topographic mapping of the glioblastoma proteome reveals a triple-axis model of intra-tumoral heterogeneity. Nature communications Glioblastoma is an aggressive form of brain cancer with well-established patterns of intra-tumoral heterogeneity implicated in treatment resistance and progression. While regional and single cell transcriptomic variations of glioblastoma have been recently resolved, downstream phenotype-level proteomic programs have yet to be assigned across glioblastoma's hallmark histomorphologic niches. Here, we leverage mass spectrometry to spatially align abundance levels of 4,794 proteins to distinct histologic patterns across 20 patients and propose diverse molecular programs operational within these regional tumor compartments. Using machine learning, we overlay concordant transcriptional information, and define two distinct proteogenomic programs, MYC- and KRAS-axis hereon, that cooperate with hypoxia to produce a tri-dimensional model of intra-tumoral heterogeneity. Moreover, we highlight differential drug sensitivities and relative chemoresistance in glioblastoma cell lines with enhanced KRAS programs. Importantly, these pharmacological differences are less pronounced in transcriptional glioblastoma subgroups suggesting that this model may provide insights for targeting heterogeneity and overcoming therapy resistance. 10.1038/s41467-021-27667-w
De novo transcriptome assembly of Angelica dahurica and characterization of coumarin biosynthesis pathway genes. Zhao Liqiang,Zhang Shengxiang,Shan Chunmiao,Shi Yuanyuan,Wu Huan,Wu Jiawen,Peng Daiyin Gene Angelica dahurica (Hoffm.) Benth. & Hook.f. ex Franch. & Sav (A. dahurica) is a famous Chinese herb known for the production of coumarins, important secondary metabolites with wide-ranging pharmacological activities. In particular, the methoxylated coumarins like those produced by A. dahurica are known for their anti-inflammatory, anti-cancer, and anti-oxidant pharmacological effects. However, the molecular mechanism of coumarin biosynthesis in A. dahurica has not been studied. Such investigation could help scientists harness the biosynthesis potential of methoxylated coumarins. Here we present, three transcriptomes corresponding to leaf, root, and stem tissues of A. dahurica. A total of 114,310 unigenes with an average length of 1118 bp were de novo assembled, and 81,404 (71.21%) of those unigenes were annotated. Then, 181 unigenes encoding the seven key enzymes involved were identified, for which COMT (Caffeic acid 3-O-methyltransferase) was spatially used in a phylogenetic analysis, and some of these key enzyme genes were verified by qRT-PCR. Differentially expressed genes and root-specific-expressed genes were identified, by comparing genes' profile activity between roots and other tissues. Furthermore, multiple genes encoding key enzymes or transcription factors related to coumarin biosynthesis were identified and analyzed. This study is the first to report comprehensive gene information of A. dahurica at the transcriptional level, and to distinguish candidate genes related to its biosynthesis of coumarin, thus laying a foundation for this pathway's further exploration in A. dahurica. 10.1016/j.gene.2021.145713
Context Specificity in Causal Signaling Networks Revealed by Phosphoprotein Profiling. Hill Steven M,Nesser Nicole K,Johnson-Camacho Katie,Jeffress Mara,Johnson Aimee,Boniface Chris,Spencer Simon E F,Lu Yiling,Heiser Laura M,Lawrence Yancey,Pande Nupur T,Korkola James E,Gray Joe W,Mills Gordon B,Mukherjee Sach,Spellman Paul T Cell systems Signaling networks downstream of receptor tyrosine kinases are among the most extensively studied biological networks, but new approaches are needed to elucidate causal relationships between network components and understand how such relationships are influenced by biological context and disease. Here, we investigate the context specificity of signaling networks within a causal conceptual framework using reverse-phase protein array time-course assays and network analysis approaches. We focus on a well-defined set of signaling proteins profiled under inhibition with five kinase inhibitors in 32 contexts: four breast cancer cell lines (MCF7, UACC812, BT20, and BT549) under eight stimulus conditions. The data, spanning multiple pathways and comprising ∼70,000 phosphoprotein and ∼260,000 protein measurements, provide a wealth of testable, context-specific hypotheses, several of which we experimentally validate. Furthermore, the data provide a unique resource for computational methods development, permitting empirical assessment of causal network learning in a complex, mammalian setting. 10.1016/j.cels.2016.11.013
Intratumoural immune heterogeneity as a hallmark of tumour evolution and progression in hepatocellular carcinoma. Nature communications The clinical relevance of immune landscape intratumoural heterogeneity (immune-ITH) and its role in tumour evolution remain largely unexplored. Here, we uncover significant spatial and phenotypic immune-ITH from multiple tumour sectors and decipher its relationship with tumour evolution and disease progression in hepatocellular carcinomas (HCC). Immune-ITH is associated with tumour transcriptomic-ITH, mutational burden and distinct immune microenvironments. Tumours with low immune-ITH experience higher immunoselective pressure and escape via loss of heterozygosity in human leukocyte antigens and immunoediting. Instead, the tumours with high immune-ITH evolve to a more immunosuppressive/exhausted microenvironment. This gradient of immune pressure along with immune-ITH represents a hallmark of tumour evolution, which is closely linked to the transcriptome-immune networks contributing to disease progression and immune inactivation. Remarkably, high immune-ITH and its transcriptomic signature are predictive for worse clinical outcome in HCC patients. This in-depth investigation of ITH provides evidence on tumour-immune co-evolution along HCC progression. 10.1038/s41467-020-20171-7
Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection. BMC genomics BACKGROUND:Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown. RESULTS:We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No mRNA sequence variants were associated with expression variation. CONCLUSIONS:Expression variation among localities may reflect morphological differences in tumors that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival. 10.1186/s12864-021-07994-4
Molecular asymmetry in the cephalochordate embryo revealed by single-blastomere transcriptome profiling. Lin Che-Yi,Lu Mei-Yeh Jade,Yue Jia-Xing,Li Kun-Lung,Le Pétillon Yann,Yong Luok Wen,Chen Yi-Hua,Tsai Fu-Yu,Lyu Yu-Feng,Chen Cheng-Yi,Hwang Sheng-Ping L,Su Yi-Hsien,Yu Jr-Kai PLoS genetics Studies in various animals have shown that asymmetrically localized maternal transcripts play important roles in axial patterning and cell fate specification in early embryos. However, comprehensive analyses of the maternal transcriptomes with spatial information are scarce and limited to a handful of model organisms. In cephalochordates (amphioxus), an early branching chordate group, maternal transcripts of germline determinants form a compact granule that is inherited by a single blastomere during cleavage stages. Further blastomere separation experiments suggest that other transcripts associated with the granule are likely responsible for organizing the posterior structure in amphioxus; however, the identities of these determinants remain unknown. In this study, we used high-throughput RNA sequencing of separated blastomeres to examine asymmetrically localized transcripts in two-cell and eight-cell stage embryos of the amphioxus Branchiostoma floridae. We identified 111 and 391 differentially enriched transcripts at the 2-cell stage and the 8-cell stage, respectively, and used in situ hybridization to validate the spatial distribution patterns for a subset of these transcripts. The identified transcripts could be categorized into two major groups: (1) vegetal tier/germ granule-enriched and (2) animal tier/anterior-enriched transcripts. Using zebrafish as a surrogate model system, we showed that overexpression of one animal tier/anterior-localized amphioxus transcript, zfp665, causes a dorsalization/anteriorization phenotype in zebrafish embryos by downregulating the expression of the ventral gene, eve1, suggesting a potential function of zfp665 in early axial patterning. Our results provide a global transcriptomic blueprint for early-stage amphioxus embryos. This dataset represents a rich platform to guide future characterization of molecular players in early amphioxus development and to elucidate conservation and divergence of developmental programs during chordate evolution. 10.1371/journal.pgen.1009294
Patient-Derived Xenograft Models Reveal Intratumor Heterogeneity and Temporal Stability in Neuroblastoma. Braekeveldt Noémie,von Stedingk Kristoffer,Fransson Susanne,Martinez-Monleon Angela,Lindgren David,Axelson Håkan,Levander Fredrik,Willforss Jakob,Hansson Karin,Øra Ingrid,Backman Torbjörn,Börjesson Anna,Beckman Siv,Esfandyari Javanshir,Berbegall Ana P,Noguera Rosa,Karlsson Jenny,Koster Jan,Martinsson Tommy,Gisselsson David,Påhlman Sven,Bexell Daniel Cancer research Patient-derived xenografts (PDX) and the Avatar, a single PDX mirroring an individual patient, are emerging tools in preclinical cancer research. However, the consequences of intratumor heterogeneity for PDX modeling of biomarkers, target identification, and treatment decisions remain underexplored. In this study, we undertook serial passaging and comprehensive molecular analysis of neuroblastoma orthotopic PDXs, which revealed strong intrinsic genetic, transcriptional, and phenotypic stability for more than 2 years. The PDXs showed preserved neuroblastoma-associated gene signatures that correlated with poor clinical outcome in a large cohort of patients with neuroblastoma. Furthermore, we captured spatial intratumor heterogeneity using ten PDXs from a single high-risk patient tumor. We observed diverse growth rates, transcriptional, proteomic, and phosphoproteomic profiles. PDX-derived transcriptional profiles were associated with diverse clinical characteristics in patients with high-risk neuroblastoma. These data suggest that high-risk neuroblastoma contains elements of both temporal stability and spatial intratumor heterogeneity, the latter of which complicates clinical translation of personalized PDX-Avatar studies into preclinical cancer research. These findings underpin the complexity of PDX modeling as a means to advance translational applications against neuroblastoma. . 10.1158/0008-5472.CAN-18-0527
Infrared Spectroscopic Imaging Visualizes a Prognostic Extracellular Matrix-Related Signature in Breast Cancer. Tiwari Saumya,Triulzi Tiziana,Holton Sarah,Regondi Viola,Paolini Biagio,Tagliabue Elda,Bhargava Rohit Scientific reports Molecular analysis techniques such as gene expression analysis and proteomics have contributed greatly to our understanding of cancer heterogeneity. In prior studies, gene expression analysis was shown to stratify patient outcome on the basis of tumor-microenvironment associated genes. A specific gene expression profile, referred to as ECM3 (Extracellular Matrix Cluster 3), indicated poorer survival in patients with grade III tumors. In this work, we aimed to visualize the downstream effects of this gene expression profile onto the tissue, thus providing a spatial context to altered gene expression profiles. Using infrared spectroscopic imaging, we identified spectral patterns specific to the ECM3 gene expression profile, achieving a high spectral classification performance of 0.87 as measured by the area under the curve of the receiver operating characteristic curve. On a patient level, we correctly identified 20 out of 22 ECM3 group patients and 19 out of 20 non-ECM3 group patients by using this spectroscopic imaging-based classifier. By comparing pixels that were identified as ECM3 or non-ECM3 with H&E and IHC images, we were also able to observe an association between tissue morphology and the gene expression clusters, showing the ability of our method to capture broad outcome associated features from infrared images. 10.1038/s41598-020-62403-2
Systematic analysis reveals long noncoding RNAs regulating neighboring transcription factors in human cancers. Liu Zhi,Dai Juncheng,Shen Hongbing Biochimica et biophysica acta. Molecular basis of disease Long noncoding RNAs (lncRNAs) are proposed to play essential roles in regulating gene transcription. Moreover, a subset has been implicated in modulating the expression of the nearby loci. Here we systematically evaluated the relationship between lncRNAs and their neighboring genes based on transcriptome expression profiles from 4900 samples across 12 cancer types. Our findings reveal that lncRNAs, especially those of high syntenic conservation across species, are spatially correlated with transcription factors across the genome. Combining the methods of conservation, co-expression, and causal inference test, we identified a list of 28 lncRNA/TF regulatory pairs across 12 TCGA cancer types, and 19 of which were further confirmed in additional cancer cell lines. Several of these pairs, including PTV1/MYC and GATA6-AS1/GATA6, show prior evidence of regulatory relationships. Other candidates such as LINC00261/FOXA2 and PITRM1-AS1/KLF6 were novel. Our study highlights the significant roles of lncRNAs in tumorigenesis and provides a comprehensive overview of lncRNA regulation on its neighboring TF genes in human cancers. 10.1016/j.bbadis.2018.05.006
Genome-wide RNA Tomography in the zebrafish embryo. Junker Jan Philipp,Noël Emily S,Guryev Victor,Peterson Kevin A,Shah Gopi,Huisken Jan,McMahon Andrew P,Berezikov Eugene,Bakkers Jeroen,van Oudenaarden Alexander Cell Advancing our understanding of embryonic development is heavily dependent on identification of novel pathways or regulators. Although genome-wide techniques such as RNA sequencing are ideally suited for discovering novel candidate genes, they are unable to yield spatially resolved information in embryos or tissues. Microscopy-based approaches, using in situ hybridization, for example, can provide spatial information about gene expression, but are limited to analyzing one or a few genes at a time. Here, we present a method where we combine traditional histological techniques with low-input RNA sequencing and mathematical image reconstruction to generate a high-resolution genome-wide 3D atlas of gene expression in the zebrafish embryo at three developmental stages. Importantly, our technique enables searching for genes that are expressed in specific spatial patterns without manual image annotation. We envision broad applicability of RNA tomography as an accurate and sensitive approach for spatially resolved transcriptomics in whole embryos and dissected organs. 10.1016/j.cell.2014.09.038
A molecular atlas of the developing ectoderm defines neural, neural crest, placode, and nonneural progenitor identity in vertebrates. Plouhinec Jean-Louis,Medina-Ruiz Sofía,Borday Caroline,Bernard Elsa,Vert Jean-Philippe,Eisen Michael B,Harland Richard M,Monsoro-Burq Anne H PLoS biology During vertebrate neurulation, the embryonic ectoderm is patterned into lineage progenitors for neural plate, neural crest, placodes and epidermis. Here, we use Xenopus laevis embryos to analyze the spatial and temporal transcriptome of distinct ectodermal domains in the course of neurulation, during the establishment of cell lineages. In order to define the transcriptome of small groups of cells from a single germ layer and to retain spatial information, dorsal and ventral ectoderm was subdivided along the anterior-posterior and medial-lateral axes by microdissections. Principal component analysis on the transcriptomes of these ectoderm fragments primarily identifies embryonic axes and temporal dynamics. This provides a genetic code to define positional information of any ectoderm sample along the anterior-posterior and dorsal-ventral axes directly from its transcriptome. In parallel, we use nonnegative matrix factorization to predict enhanced gene expression maps onto early and mid-neurula embryos, and specific signatures for each ectoderm area. The clustering of spatial and temporal datasets allowed detection of multiple biologically relevant groups (e.g., Wnt signaling, neural crest development, sensory placode specification, ciliogenesis, germ layer specification). We provide an interactive network interface, EctoMap, for exploring synexpression relationships among genes expressed in the neurula, and suggest several strategies to use this comprehensive dataset to address questions in developmental biology as well as stem cell or cancer research. 10.1371/journal.pbio.2004045
Developmentally regulated long non-coding RNAs in Xenopus tropicalis. Forouzmand Elmira,Owens Nick D L,Blitz Ira L,Paraiso Kitt D,Khokha Mustafa K,Gilchrist Michael J,Xie Xiaohui,Cho Ken W Y Developmental biology Advances in RNA sequencing technologies have led to the surprising discovery that a vast number of transcripts emanate from regions of the genome that are not part of coding genes. Although some of the smaller ncRNAs such as microRNAs have well-characterized functions, the majority of long ncRNA (lncRNA) functions remain poorly understood. Understanding the significance of lncRNAs is an important challenge facing biology today. A powerful approach to uncovering the function of lncRNAs is to explore temporal and spatial expression profiling. This may be particularly useful for classes of lncRNAs that have developmentally important roles as the expression of such lncRNAs will be expected to be both spatially and temporally regulated during development. Here, we take advantage of our ultra-high frequency (temporal) sampling of Xenopus embryos to analyze gene expression trajectories of lncRNA transcripts over the first 3 days of development. We computationally identify 5689 potential single- and multi-exon lncRNAs. These lncRNAs demonstrate clear dynamic expression patterns. A subset of them displays highly correlative temporal expression profiles with respect to those of the neighboring genes. We also identified spatially localized lncRNAs in the gastrula stage embryo. These results suggest that lncRNAs have regulatory roles during early embryonic development. 10.1016/j.ydbio.2016.06.016
Distinct tumor microenvironments of lytic and blastic bone metastases in prostate cancer patients. Journal for immunotherapy of cancer The most common metastatic lesions of prostate cancer are in bone and can be classified into three distinct pathology subtypes: lytic, blastic, and an indeterminate mixture of both. We investigated a cohort of decalcified formalin-fixed and paraffin-embedded (FFPE) patient specimens from the bone that contained metastatic prostate cancer with lytic or blastic features. These tissue sections were utilized for immunohistochemistry (IHC) staining, isolation of RNA for gene expression, and Digital Spatial Profiling (DSP) of changes in both the tumor and microenvironment. A diverse set of unique immune cell populations and signaling pathways to both lytic and blastic types of prostate cancer metastases were present. In blastic lesions immune cells were enriched for pSTAT3 and components of the JAK-STAT pathway. In lytic-type lesions, immune cells were enriched for pAKT activity and components of the PI3K-AKT pathway. Enrichment for immune checkpoints including PD-L1, B7-H4, OX40L, and IDO-1 were identified in blastic prostate cancer, providing new therapeutic targets for patients with bone metastases. Biopsies could guide selection of patients into appropriate therapeutic interventions based on protein levels and RNA expression of desired targets in metastatic disease. Molecular pathology has been an excellent complement to the diagnosis, treatment, and management of primary tumors and could be successfully extended to patients with metastatic lesions. 10.1186/s40425-019-0753-3
DNA copy number motifs are strong and independent predictors of survival in breast cancer. Communications biology Somatic copy number alterations are a frequent sign of genome instability in cancer. A precise characterization of the genome architecture would reveal underlying instability mechanisms and provide an instrument for outcome prediction and treatment guidance. Here we show that the local spatial behavior of copy number profiles conveys important information about this architecture. Six filters were defined to characterize regional traits in copy number profiles, and the resulting Copy Aberration Regional Mapping Analysis (CARMA) algorithm was applied to tumors in four breast cancer cohorts (n = 2919). The derived motifs represent a layer of information that complements established molecular classifications of breast cancer. A score reflecting presence or absence of motifs provided a highly significant independent prognostic predictor. Results were consistent between cohorts. The nonsite-specific occurrence of the detected patterns suggests that CARMA captures underlying replication and repair defects and could have a future potential in treatment stratification. 10.1038/s42003-020-0884-6
New developments in profiling and imaging of proteins from tissue sections by MALDI mass spectrometry. Chaurand Pierre,Norris Jeremy L,Cornett D Shannon,Mobley James A,Caprioli Richard M Journal of proteome research Molecular imaging of tissue by MALDI mass spectrometry is a powerful tool for visualizing the spatial distribution of constituent analytes with high molecular specificity. Although the technique is relatively young, it has already contributed to the understanding of many diverse areas of human health. In recent years, a great many advances in the practice of imaging mass spectrometry have taken place, making the technique more sensitive, robust, and ultimately useful. The purpose of this review is to highlight some of the more recent technological advances that have improved the efficiency of imaging mass spectrometry for clinical applications. Advances in the way MALDI mass spectrometry is integrated with histology, improved methods for automation, and better tools for data analysis are outlined in this review. Refined top-down strategies for the identification and validation of candidate biomarkers found in tissue sections are discussed. A clinical example highlighting the application of these methods to a cohort of clinical samples is described. 10.1021/pr060346u
Transcription Factor Profiling Identifies Spatially Heterogenous Mediators of Follicular Thyroid Cancer Invasion. Nicolson Norman G,Paulsson Johan O,Juhlin C Christofer,Carling Tobias,Korah Reju Endocrine pathology While minimally invasive follicular thyroid cancer (miFTC) generally has low risk of recurrence or death, encapsulated angioinvasive (eaFTC) or widely invasive (wiFTC) histological subtypes display significantly worse prognosis. Drivers of invasion are incompletely understood. Therefore, tissue samples including miFTC, eaFTC, and wiFTC tumors, as well as histologically normal thyroid adjacent to benign follicular adenomas, were selected from a cohort (n = 21) of thyroid tumor patients, and the gene expression of selected transcription factors was characterized with quantitative PCR. Invasion-relevant spatial expression patterns of selected transcription factors were subsequently characterized with immunohistochemistry. E2F1 was over-expressed in all 3 subtypes (p<0.01). SP1 was differentially expressed in eaFTC and wiFTC compared with normal (p=0.01 and 0.04, respectively). TCF7L2 was significantly upregulated in wiFTC specifically (p<0.05). While these findings were mRNA specific, immunohistochemistry of additional cancer-associated transcription factors revealed differential expression along the tumor invasive front relative to the central tumor, and histone acetylation modulators emerged as putative invasion markers. These findings may have significant implications for the interpretation of bulk gene expression analysis of thyroid tumor samples or for the development of targeted therapeutics for this rare but aggressive thyroid cancer variant. 10.1007/s12022-020-09651-0
Prospective patient stratification into robust cancer-cell intrinsic subtypes from colorectal cancer biopsies. The Journal of pathology Colorectal cancer (CRC) biopsies underpin accurate diagnosis, but are also relevant for patient stratification in molecularly-guided clinical trials. The consensus molecular subtypes (CMSs) and colorectal cancer intrinsic subtypes (CRISs) transcriptional signatures have potential clinical utility for improving prognostic/predictive patient assignment. However, their ability to provide robust classification, particularly in pretreatment biopsies from multiple regions or at different time points, remains untested. In this study, we undertook a comprehensive assessment of the robustness of CRC transcriptional signatures, including CRIS and CMS, using a range of tumour sampling methodologies currently employed in clinical and translational research. These include analyses using (i) laser-capture microdissected CRC tissue, (ii) eight publically available rectal cancer biopsy data sets (n = 543), (iii) serial biopsies (from AXEBeam trial, NCT00828672; n = 10), (iv) multi-regional biopsies from colon tumours (n = 29 biopsies, n = 7 tumours), and (v) pretreatment biopsies from the phase II rectal cancer trial COPERNCIUS (NCT01263171; n = 44). Compared to previous results obtained using CRC resection material, we demonstrate that CMS classification in biopsy tissue is significantly less capable of reliably classifying patient subtype (43% unknown in biopsy versus 13% unknown in resections, p = 0.0001). In contrast, there was no significant difference in classification rate between biopsies and resections when using the CRIS classifier. Additionally, we demonstrated that CRIS provides significantly better spatially- and temporally- robust classification of molecular subtypes in CRC primary tumour tissue compared to CMS (p = 0.003 and p = 0.02, respectively). These findings have potential to inform ongoing biopsy-based patient stratification in CRC, enabling robust and stable assignment of patients into clinically-informative arms of prospective multi-arm, multi-stage clinical trials. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland. 10.1002/path.5051
CTCF loss mediates unique DNA hypermethylation landscapes in human cancers. Damaschke Nathan A,Gawdzik Joseph,Avilla Mele,Yang Bing,Svaren John,Roopra Avtar,Luo Jian-Hua,Yu Yan P,Keles Sunduz,Jarrard David F Clinical epigenetics BACKGROUND:The chromatin insulator CCCTC-binding factor (CTCF) displays tissue-specific DNA binding sites that regulate transcription and chromatin organization. Despite evidence linking CTCF to the protection of epigenetic states through barrier insulation, the impact of CTCF loss on genome-wide DNA methylation sites in human cancer remains undefined. RESULTS:Here, we demonstrate that prostate and breast cancers within The Cancer Genome Atlas (TCGA) exhibit frequent copy number loss of CTCF and that this loss is associated with increased DNA methylation events that occur preferentially at CTCF binding sites. CTCF sites differ among tumor types and result in tissue-specific methylation patterns with little overlap between breast and prostate cancers. DNA methylation and transcriptome profiling in vitro establish that forced downregulation of CTCF leads to spatially distinct DNA hypermethylation surrounding CTCF binding sites, loss of CTCF binding, and decreased gene expression that is also seen in human tumors. DNA methylation inhibition reverses loss of expression at these CTCF-regulated genes. CONCLUSION:These findings establish CTCF loss as a major mediator in directing localized DNA hypermethylation events in a tissue-specific fashion and further support its role as a driver of the cancer phenotype. 10.1186/s13148-020-00869-7
A single cell transcriptome atlas of the developing zebrafish hindbrain. Development (Cambridge, England) Segmentation of the vertebrate hindbrain leads to the formation of rhombomeres, each with a distinct anteroposterior identity. Specialised boundary cells form at segment borders that act as a source or regulator of neuronal differentiation. In zebrafish, there is spatial patterning of neurogenesis in which non-neurogenic zones form at boundaries and segment centres, in part mediated by Fgf20 signalling. To further understand the control of neurogenesis, we have carried out single cell RNA sequencing of the zebrafish hindbrain at three different stages of patterning. Analyses of the data reveal known and novel markers of distinct hindbrain segments, of cell types along the dorsoventral axis, and of the transition of progenitors to neuronal differentiation. We find major shifts in the transcriptome of progenitors and of differentiating cells between the different stages analysed. Supervised clustering with markers of boundary cells and segment centres, together with RNA-seq analysis of Fgf-regulated genes, has revealed new candidate regulators of cell differentiation in the hindbrain. These data provide a valuable resource for functional investigations of the patterning of neurogenesis and the transition of progenitors to neuronal differentiation. 10.1242/dev.184143
A Graphical Model of Smoking-Induced Global Instability in Lung Cancer. Wang Yanbo,Qian Weikang,Yuan Bo IEEE/ACM transactions on computational biology and bioinformatics Smoking is the major cause of lung cancer and the leading cause of cancer-related death in the world. The most current view about lung cancer is no longer limited to individual genes being mutated by any carcinogenic insults from smoking. Instead, tumorigenesis is a phenotype conferred by many systematic and global alterations, leading to extensive heterogeneity and variation for both the genotypes and phenotypes of individual cancer cells. Thus, strategically it is foremost important to develop a methodology to capture any consistent and global alterations presumably shared by most of the cancerous cells for a given population. This is particularly true that almost all of the data collected from solid cancers (including lung cancers) are usually distant apart over a large span of temporal or even spatial contexts. Here, we report a multiple non-Gaussian graphical model to reconstruct the gene interaction network using two previously published gene expression datasets. Our graphical model aims to selectively detect gross structural changes at the level of gene interaction networks. Our methodology is extensively validated, demonstrating good robustness, as well as the selectivity and specificity expected based on our biological insights. In summary, gene regulatory networks are still relatively stable during presumably the early stage of neoplastic transformation. But drastic structural differences can be found between lung cancer and its normal control, including the gain of functional modules for cellular proliferations such as EGFR and PDGFRA, as well as the lost of the important IL6 module, supporting their roles as potential drug targets. Interestingly, our method can also detect early modular changes, with the ALDH3A1 and its associated interactions being strongly implicated as a potential early marker, whose activations appear to alter LCN2 module as well as its interactions with the important TP53-MDM2 circuitry. Our strategy using the graphical model to reconstruct gene interaction work with biologically-inspired constraints exemplifies the importance and beauty of biology in developing any bio-computational approach. 10.1109/TCBB.2016.2599867
Identification of two genes differentially expressed upon different spatial configuration of the MGH-U3 human bladder cancer cells. Champetier Serge,Fradet Yves,Bachvarov Dimcho Urologic oncology By using mRNA differential display technology, we have compared gene expression in neoconfluent (actively proliferating) and 14-day confluent (non-proliferative) MGH-U3 cells. Two clones (6.17.6 and 10.9.2) displayed significantly increased expression in the 14-day confluent MGH-U3 cells, as confirmed by Northern blot analysis. Nucleotide sequence showed that clone 6.17.6 was part of the human insulin-like growth factor binding protein-5 (IGFBP-5) cDNA, while clone 10.9.2 was part of a human gene with unknown function (clone HSI12602). The differential IGFBP-5 expression in this model system was also evidenced by Western analysis. The possible anti-proliferative role of IGFBP-5 in this model system is briefly discussed. 10.1016/s1078-1439(01)00130-2
Systems Biology Analysis of Heterocellular Signaling. Tape Christopher J Trends in biotechnology Tissues comprise multiple heterotypic cell types (e.g., epithelial, mesenchymal, and immune cells). Communication between heterotypic cell types is essential for biological cohesion and is frequently dysregulated in disease. Despite the importance of heterocellular communication, most systems biology techniques do not report cell-specific signaling data from mixtures of cells. As a result, our existing perspective of cellular behavior under-represents the influence of heterocellular signaling. Recent technical advances now permit the resolution of systems-level cell-specific signaling data. This review discusses how new physical, spatial, and isotopic resolving methods are facilitating unique systems biology studies of heterocellular communication. 10.1016/j.tibtech.2016.02.016
Long Noncoding RNAs: Insights from Biological Features and Functions to Diseases. Li Xiaolei,Wu Zhiqiang,Fu Xiaobing,Han Weidong Medicinal research reviews Over the past decade, genome-wide transcriptomic studies have shown that the mammalian genome is pervasively transcribed and produces many thousands of transcriptomes without bias from previous genome annotations. This finding, together with the discovery of a plethora of unexpected RNAs that have no obvious coding capacities, have challenged the traditional views that proteins are the main protagonists of cellular functions and that RNA is merely an intermediary between DNA sequence and its encoded protein. There are many different kinds of products that are generated by this pervasive transcription; this review focuses on long noncoding RNAs (lncRNAs) that have shown spatial and temporal specific patterns of expression and regulation in a wide variety of cells and tissues, adding significant complexity to the understanding of their biological roles. Recent research has shed new light onto the biological function significance of lncRNAs. Here, we review the rapidly advancing field of lncRNAs, describing their biological features and their roles in regulation of gene expression. Moreover, we highlight some recent advances in our understanding of ncRNA-mediated regulation of stem cell pluripotency, morphogenesis, and development, focusing mainly on the regulatory roles of lncRNAs. Finally, we consider the potential medical implications, and the potential use of lncRNAs in drug development and discovery and in the identification of molecular markers of diseases, including cancer. 10.1002/med.21254
Cell population characterization and discovery using single-cell technologies in endocrine systems. Cheung Leonard Y M,Rizzoti Karine Journal of molecular endocrinology In the last 15 years, single-cell technologies have become robust and indispensable tools to investigate cell heterogeneity. Beyond transcriptomic, genomic and epigenome analyses, technologies are constantly evolving, in particular toward multi-omics, where analyses of different source materials from a single cell are combined, and spatial transcriptomics, where resolution of cellular heterogeneity can be detected in situ. While some of these techniques are still being optimized, single-cell RNAseq has commonly been used because the examination of transcriptomes allows characterization of cell identity and, therefore, unravel previously uncharacterized diversity within cell populations. Most endocrine organs have now been investigated using this technique, and this has given new insights into organ embryonic development, characterization of rare cell types, and disease mechanisms. Here, we highlight recent studies, particularly on the hypothalamus and pituitary, and examine recent findings on the pancreas and reproductive organs where many single-cell experiments have been performed. 10.1530/JME-19-0276
Spatial transcriptomics of B cell and T cell receptors reveals lymphocyte clonal dynamics. Science (New York, N.Y.) The spatial distribution of lymphocyte clones within tissues is critical to their development, selection, and expansion. We have developed spatial transcriptomics of variable, diversity, and joining (VDJ) sequences (Spatial VDJ), a method that maps B cell and T cell receptor sequences in human tissue sections. Spatial VDJ captures lymphocyte clones that match canonical B and T cell distributions and amplifies clonal sequences confirmed by orthogonal methods. We found spatial congruency between paired receptor chains, developed a computational framework to predict receptor pairs, and linked the expansion of distinct B cell clones to different tumor-associated gene expression programs. Spatial VDJ delineates B cell clonal diversity and lineage trajectories within their anatomical niche. Thus, Spatial VDJ captures lymphocyte spatial clonal architecture across tissues, providing a platform to harness clonal sequences for therapy. 10.1126/science.adf8486
Meet the author: Troy McEachron. Molecular cell Troy McEachron, corresponding author of the technology review article "Applicability of spatial transcriptional profiling to cancer research," discusses his motivation for becoming a scientist and pursing pediatric cancer research, the challenges he faces as a Black scientist, and strategies for writing a review article. 10.1016/j.molcel.2021.03.040
Single-cell and spatial transcriptomics approaches of the bone marrow microenvironment. Al-Sabah Jude,Baccin Chiara,Haas Simon Current opinion in oncology PURPOSE OF REVIEW:The bone marrow is home to hematopoietic stem cells responsible for lifelong blood production, alongside mesenchymal stem cells required for skeletal regeneration. In the bone marrow, a unique combination of signals derived from a multitude of cell types results in the establishment of so-called niches that regulate stem-cell maintenance and differentiation. Recently, single-cell and spatially resolved transcriptomics technologies have been utilized to characterize the murine bone marrow microenvironment during homeostasis, stress and upon cancer-induced remodeling. In this review, we summarize the major findings of these studies. RECENT FINDINGS:Single-cell technologies applied to bone marrow provided the first systematic and label-free identification of bone marrow cell types, enabled their molecular and spatial characterization, and clarified the cellular sources of key prohematopoietic factors. Large transcriptional heterogeneity and novel subpopulations were observed in compartments previously thought to be homogenous. For example, Lepr Cxcl12-abundant reticular cells were shown to constitute the major source of prohematopoietic factors, but consist of subpopulations differing in their adipogenic versus osteogenic priming, morphology and localization. These subpopulations were suggested to act as professional cytokine secreting cells, thereby establishing distinct bone marrow niches. SUMMARY:Single-cell and spatially resolved transcriptomics approaches have clarified the molecular identity and localization of bone marrow-resident cell types, paving the road for a deeper exploration of bone marrow niches in the mouse and humans. 10.1097/CCO.0000000000000602
Applying single-cell technologies to clinical pathology: progress in nephropathology. The Journal of pathology Cells represent the basic building blocks of living organisms. Accurate characterisation of cellular phenotype, intercellular signalling networks, and the spatial organisation of cells within organs is crucial to deliver a better understanding of the processes underpinning physiology, and the perturbations that lead to disease. Single-cell methodologies have increased rapidly in scale and scope in recent years and are set to generate important insights into human disease. Here, we review current practices in nephropathology, which are dominated by relatively simple morphological descriptions of tissue biopsies based on their appearance using light microscopy. Bulk transcriptomics have more recently been used to explore glomerular and tubulointerstitial kidney disease, renal cancer, and the responses to injury and alloimmunity in kidney transplantation, generating novel disease insights and prognostic biomarkers. These studies set the stage for single-cell transcriptomic approaches that reveal cell-type-specific gene expression patterns in health and disease. These technologies allow genome-wide disease susceptibility genes to be interpreted with the knowledge of the specific cell populations within organs that express them, identifying candidate cell types for further study. Single-cell technologies are also moving beyond assaying individual cellular transcriptomes, to measuring the epigenetic landscape of single cells. Single-cell antigen-receptor gene sequencing also enables specific T- and B-cell clones to be tracked in different tissues and disease states. In the coming years these rich 'multi-omic' descriptions of kidney disease will enable histopathological descriptions to be comprehensively integrated with molecular phenotypes, enabling better disease classification and prognostication and the application of personalised treatment strategies. © 2020 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland. 10.1002/path.5417
Clinical management of breast cancer heterogeneity. Nature reviews. Clinical oncology Traditionally, intertumour heterogeneity in breast cancer has been documented in terms of different histological subtypes, treatment sensitivity profiles, and clinical outcomes among different patients. Results of high-throughput molecular profiling studies have subsequently revealed the true extent of this heterogeneity. Further complicating this scenario, the heterogeneous expression of the oestrogen receptor (ER), progesterone receptor (PR), and HER2 has been reported in different areas of the same tumour. Furthermore, discordance, in terms of ER, PR and HER2 expression, has also been reported between primary tumours and their matched metastatic lesions. High-throughput molecular profiling studies have confirmed that spatial and temporal intratumour heterogeneity of breast cancers exist at a level beyond common expectations. We describe the different levels of tumour heterogeneity, and discuss the strategies that can be adopted by clinicians to tackle treatment response and resistance issues associated with such heterogeneity, including a rationally selected combination of agents that target driver mutations, the targeting of deleterious passenger mutations, identifying and eradicating the 'lethal' clone, targeting the tumour microenvironment, or using adaptive treatments and immunotherapy. The identification of the most-appropriate strategies and their implementation in the clinic will prove highly challenging and necessitate the adoption of radically new practices for the optimal clinical management of breast malignancies. 10.1038/nrclinonc.2015.73
Hormone-regulated transcriptomes: lessons learned from estrogen signaling pathways in breast cancer cells. Hah Nasun,Kraus W Lee Molecular and cellular endocrinology Recent rapid advances in next generation sequencing technologies have expanded our understanding of steroid hormone signaling to a genome-wide level. In this review, we discuss the use of a novel genomic approach, global nuclear run-on coupled with massively parallel sequencing (GRO-seq), to explore new facets of the steroid hormone-regulated transcriptome, especially estrogen responses in breast cancer cells. GRO-seq is a high throughput sequencing method adapted from conventional nuclear run-on methodologies, which is used to obtain a map of the position and orientation of all transcriptionally engaged RNA polymerases across the genome with extremely high spatial resolution. GRO-seq, which is an excellent tool for examining transcriptional responses to extracellular stimuli, has been used to comprehensively assay the effects of estrogen signaling on the transcriptome of ERα-positive MCF-7 human breast cancer cells. These studies have revealed new details about estrogen-dependent transcriptional regulation, including effects on transcription by all three RNA polymerases, complex transcriptional dynamics in response to estrogen signaling, and identification novel, unannotated non-coding RNAs. Collectively, these studies have been useful in discerning the molecular logic of the estrogen-regulated mitogenic response. 10.1016/j.mce.2013.06.021
Spatial Transcriptomics Reveal Pitfalls and Opportunities for the Detection of Rare High-Plasticity Breast Cancer Subtypes. Laboratory investigation; a journal of technical methods and pathology Breast cancer is one of the most prominent types of cancers, in which therapeutic resistance is a major clinical concern. Specific subtypes, such as claudin-low and metaplastic breast carcinoma (MpBC), have been associated with high nongenetic plasticity, which can facilitate resistance. The similarities and differences between these orthogonal subtypes, identified by molecular and histopathological analyses, respectively, remain insufficiently characterized. Furthermore, adequate methods to identify high-plasticity tumors to better anticipate resistance are lacking. Here, we analyzed 11 triple-negative breast tumors, including 3 claudin-low and 4 MpBC, via high-resolution spatial transcriptomics. We combined pathological annotations and deconvolution approaches to precisely identify tumor spots, on which we performed signature enrichment, differential expression, and copy number analyses. We used The Cancer Genome Atlas and Cancer Cell Line Encyclopedia public databases for external validation of expression markers. By focusing our spatial transcriptomic analyses on tumor cells in MpBC samples, we bypassed the negative impact of stromal contamination and identified specific markers that are neither expressed in other breast cancer subtypes nor expressed in stromal cells. Three markers (BMPER, POPDC3, and SH3RF3) were validated in external expression databases encompassing bulk tumor material and stroma-free cell lines. We unveiled that existing bulk expression signatures of high-plasticity breast cancers are relevant in mesenchymal transdifferentiated compartments but can be hindered by abundant stromal cells in tumor samples, negatively impacting their clinical applicability. Spatial transcriptomic analyses constitute powerful tools to identify specific expression markers and could thus enhance diagnosis and clinical care of rare high-plasticity breast cancers. 10.1016/j.labinv.2023.100258
3D-cardiomics: A spatial transcriptional atlas of the mammalian heart. Journal of molecular and cellular cardiology Understanding the spatial gene expression and regulation in the heart is key to uncovering its developmental and physiological processes, during homeostasis and disease. Numerous techniques exist to gain gene expression and regulation information in organs such as the heart, but few utilize intuitive true-to-life three-dimensional representations to analyze and visualise results. Here we combined transcriptomics with 3D-modelling to interrogate spatial gene expression in the mammalian heart. For this, we microdissected and sequenced transcriptome-wide 18 anatomical sections of the adult mouse heart. Our study has unveiled known and novel genes that display complex spatial expression in the heart sub-compartments. We have also created 3D-cardiomics, an interface for spatial transcriptome analysis and visualization that allows the easy exploration of these data in a 3D model of the heart. 3D-cardiomics is accessible from http://3d-cardiomics.erc.monash.edu/. 10.1016/j.yjmcc.2021.09.011
Discovering hematoma-stimulated circuits for secondary brain injury after intraventricular hemorrhage by spatial transcriptome analysis. Frontiers in immunology Introduction:Central nervous system (CNS) diseases, such as neurodegenerative disorders and brain diseases caused by acute injuries, are important, yet challenging to study due to disease lesion locations and other complexities. Methods:Utilizing the powerful method of spatial transcriptome analysis together with novel algorithms we developed for the study, we report here for the first time a 3D trajectory map of gene expression changes in the brain following acute neural injury using a mouse model of intraventricular hemorrhage (IVH). IVH is a common and representative complication after various acute brain injuries with severe mortality and mobility implications. Results:Our data identified three main 3D global pseudospace-time trajectory bundles that represent the main neural circuits from the lateral ventricle to the hippocampus and primary cortex affected by experimental IVH stimulation. Further analysis indicated a rapid response in the primary cortex, as well as a direct and integrated effect on the hippocampus after IVH stimulation. Discussion:These results are informative for understanding the pathophysiological changes, including the spatial and temporal patterns of gene expression changes, in IVH patients after acute brain injury, strategizing more effective clinical management regimens, and developing novel bioinformatics strategies for the study of other CNS diseases. The algorithm strategies used in this study are searchable via a web service (www.combio-lezhang.online/3dstivh/home). 10.3389/fimmu.2023.1123652
A recent survey on colon cancer detection techniques. Rathore Saima,Hussain Mutawarra,Ali Ahmad,Khan Asifullah IEEE/ACM transactions on computational biology and bioinformatics Colon cancer causes deaths of about half a million people every year. Common method of its detection is histopathological tissue analysis, which, though leads to vital diagnosis, is significantly correlated to the tiredness, experience, and workload of the pathologist. Researchers have been working since decades to get rid of manual inspection, and to develop trustworthy systems for detecting colon cancer. Several techniques, based on spectral/spatial analysis of colon biopsy images, and serum and gene analysis of colon samples, have been proposed in this regard. Due to rapid evolution of colon cancer detection techniques, a latest review of recent research in this field is highly desirable. The aim of this paper is to discuss various colon cancer detection techniques. In this survey, we categorize the techniques on the basis of the adopted methodology and underlying data set, and provide detailed description of techniques in each category. Additionally, this study provides an extensive comparison of various colon cancer detection categories, and of multiple techniques within each category. Further, most of the techniques have been evaluated on similar data set to provide a fair performance comparison. Analysis reveals that neither of the techniques is perfect; however, research community is progressively inching toward the finest possible solution. 10.1109/TCBB.2013.84
Cancer meets neuroscience: the association between glioma occurrence and intrinsic brain features. Brain : a journal of neurology 10.1093/brain/awad034
Identification and validation of core genes in tumor-educated platelets for human gastrointestinal tumor diagnosis using network-based transcriptomic analysis. Platelets Gastrointestinal (GI) tumors have increasing incidence worldwide with their underlying mechanisms still not being fully understood. The use of tumor-educated platelets (TEPs) in liquid biopsy is a newly-emerged blood-based cancer diagnostic method. Herein, we aimed to investigate the genomic changes of TEPs in GI tumor development and their potential functions using network-based meta-analysis combined with bioinformatic methods. We used a total of three eligible RNA-seq datasets, which were integrated using multiple meta-analysis methods on the NetworkAnalyst website, and identified 775 DEGs (differentially expressed genes; 51 up-regulated and 724 down-regulated genes) in GI tumor relative to healthy control (HC) samples. These TEP DEGs were mostly enriched in bone marrow-derived cell types and associated with gene ontology (GO) of "carcinoma" and could affect pathways of "Integrated Cancer Pathway" and "Generic transcription pathway" respectively for highly and lowly expressed DEGs. Combined network-based meta-analysis and protein-protein interaction (PPI) analysis identified cyclin dependent kinase 1 (CDK1) and heat shock protein family A (Hsp70) member 5 (HSPA5) to be the hub genes with the highest degree centrality (DC), being up-regulated and down-regulated in TEPs, respectively. GO and Kyoto Encyclopedia of Genes and Genomes (KEGG) results showed that the hub genes were primarily related to cell cycle and division, nucleobase-containing compound and carbohydrate transport, and endoplasmic reticulum unfolded protein response. Additionally, the nomogram model suggested that the two-gene signature owns extraordinary predictive power for GI tumor diagnosis. Further, the two-gene signature was demonstrated to have potential value for metastatic GI tumor diagnosis. The expression levels of CDK1 and HSPA5 in clinical platelet samples were verified to be consistent with the results from bioinformatic analysis. This study identified a two-gene signature encompassing CDK1 and HSPA5 that can be used as a biomarker for GI tumor diagnosis and maybe even cancer-associated thrombosis (CAT)-related prognosis. 10.1080/09537104.2023.2212071
Molecular-genetic imaging based on reporter gene expression. Kang Joo Hyun,Chung June-Key Journal of nuclear medicine : official publication, Society of Nuclear Medicine Molecular imaging includes proteomic, metabolic, cellular biologic process, and genetic imaging. In a narrow sense, molecular imaging means genetic imaging and can be called molecular-genetic imaging. Imaging reporter genes play a leading role in molecular-genetic imaging. There are 3 major methods of molecular-genetic imaging, based on optical, MRI, and nuclear medicine modalities. For each of these modalities, various reporter genes and probes have been developed, and these have resulted in successful transitions from bench to bedside applications. Each of these imaging modalities has its unique advantages and disadvantages. Fluorescent and bioluminescent optical imaging modalities are simple, less expensive, more convenient, and more user friendly than other imaging modalities. Another advantage, especially of bioluminescence imaging, is its ability to detect low levels of gene expression. MRI has the advantage of high spatial resolution, whereas nuclear medicine methods are highly sensitive and allow data from small-animal imaging studies to be translated to clinical practice. Moreover, multimodality imaging reporter genes will allow us to choose the imaging technologies that are most appropriate for the biologic problem at hand and facilitate the clinical application of reporter gene technologies. Reporter genes can be used to visualize the levels of expression of particular exogenous and endogenous genes and several intracellular biologic phenomena, including specific signal transduction pathways, nuclear receptor activities, and protein-protein interactions. This technique provides a straightforward means of monitoring tumor mass and can visualize the in vivo distributions of target cells, such as immune cells and stem cells. Molecular imaging has gradually evolved into an important tool for drug discovery and development, and transgenic mice with an imaging reporter gene can be useful during drug and stem cell therapy development. Moreover, instrumentation improvements, the identification of novel targets and genes, and imaging probe developments suggest that molecular-genetic imaging is likely to play an increasingly important role in the diagnosis and therapy of cancer. 10.2967/jnumed.107.045955
Matrix-assisted laser desorption/ionization imaging mass spectrometry: a promising technique for reproductive research. Lagarrigue Mélanie,Lavigne Régis,Guével Blandine,Com Emmanuelle,Chaurand Pierre,Pineau Charles Biology of reproduction Matrix-assisted laser desorption/ionization (MALDI) tissue imaging mass spectrometry is particularly promising among the numerous applications of mass spectrometry. It is used for probing and analyzing the spatial arrangement of a wide range of molecules, including proteins, peptides, lipids, drugs, and metabolites, directly in thin slices of tissue. In the field of proteomics, the technology avoids tedious and time-consuming extraction and fractionation steps classically required for sample analysis. MALDI imaging mass spectrometry is increasingly recognized as a powerful method for clinical proteomics, particularly in cancer research. The technology has particular potential for the discovery of new tissue biomarker candidates, classification of tumors, early diagnosis or prognosis, elucidating pathogenesis pathways, and therapy monitoring. Over recent years, MALDI imaging mass spectrometry has been used for molecular profiling and imaging directly in male and female reproductive tissues. This review will consider some of the recent publications in the field, addressing a range of issues covering embryo development, gene expression product profiling during gametogenesis, and seeking and identifying biomarkers of reproductive cancers. The wealth of advances in mass spectrometry imaging will inevitably attract biologists and clinicians as the advantages and power of this technology become more widely known. This review will also discuss bottlenecks and the many technical issues that remain to be resolved before laboratories in the field can adopt the technology. We foresee that MALDI imaging mass spectrometry will have a major impact in reproductive research by opening new avenues to the understanding of various molecular mechanisms and the diagnosis of reproductive pathologies. 10.1095/biolreprod.111.094896
The Human Lung Cell Atlas: A High-Resolution Reference Map of the Human Lung in Health and Disease. American journal of respiratory cell and molecular biology Lung disease accounts for every sixth death globally. Profiling the molecular state of all lung cell types in health and disease is currently revolutionizing the identification of disease mechanisms and will aid the design of novel diagnostic and personalized therapeutic regimens. Recent progress in high-throughput techniques for single-cell genomic and transcriptomic analyses has opened up new possibilities to study individual cells within a tissue, classify these into cell types, and characterize variations in their molecular profiles as a function of genetics, environment, cell-cell interactions, developmental processes, aging, or disease. Integration of these cell state definitions with spatial information allows the in-depth molecular description of cellular neighborhoods and tissue microenvironments, including the tissue resident structural and immune cells, the tissue matrix, and the microbiome. The Human Cell Atlas consortium aims to characterize all cells in the healthy human body and has prioritized lung tissue as one of the flagship projects. Here, we present the rationale, the approach, and the expected impact of a Human Lung Cell Atlas. 10.1165/rcmb.2018-0416TR
Spatial transcriptomics reveals novel genes during the remodelling of the embryonic human arterial valves. PLoS genetics Abnormalities of the arterial valves, including bicuspid aortic valve (BAV) are amongst the most common congenital defects and are a significant cause of morbidity as well as predisposition to disease in later life. Despite this, and compounded by their small size and relative inaccessibility, there is still much to understand about how the arterial valves form and remodel during embryogenesis, both at the morphological and genetic level. Here we set out to address this in human embryos, using Spatial Transcriptomics (ST). We show that ST can be used to investigate the transcriptome of the developing arterial valves, circumventing the problems of accurately dissecting out these tiny structures from the developing embryo. We show that the transcriptome of CS16 and CS19 arterial valves overlap considerably, despite being several days apart in terms of human gestation, and that expression data confirm that the great majority of the most differentially expressed genes are valve-specific. Moreover, we show that the transcriptome of the human arterial valves overlaps with that of mouse atrioventricular valves from a range of gestations, validating our dataset but also highlighting novel genes, including four that are not found in the mouse genome and have not previously been linked to valve development. Importantly, our data suggests that valve transcriptomes are under-represented when using commonly used databases to filter for genes important in cardiac development; this means that causative variants in valve-related genes may be excluded during filtering for genomic data analyses for, for example, BAV. Finally, we highlight "novel" pathways that likely play important roles in arterial valve development, showing that mouse knockouts of RBP1 have arterial valve defects. Thus, this study has confirmed the utility of ST for studies of the developing heart valves and broadens our knowledge of the genes and signalling pathways important in human valve development. 10.1371/journal.pgen.1010777
System-wide transcriptome damage and tissue identity loss in COVID-19 patients. Cell reports. Medicine The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections. 10.1016/j.xcrm.2022.100522
Decoding functional cell-cell communication events by multi-view graph learning on spatial transcriptomics. Briefings in bioinformatics Cell-cell communication events (CEs) are mediated by multiple ligand-receptor (LR) pairs. Usually only a particular subset of CEs directly works for a specific downstream response in a particular microenvironment. We name them as functional communication events (FCEs) of the target responses. Decoding FCE-target gene relations is: important for understanding the mechanisms of many biological processes, but has been intractable due to the mixing of multiple factors and the lack of direct observations. We developed a method HoloNet for decoding FCEs using spatial transcriptomic data by integrating LR pairs, cell-type spatial distribution and downstream gene expression into a deep learning model. We modeled CEs as a multi-view network, developed an attention-based graph learning method to train the model for generating target gene expression with the CE networks, and decoded the FCEs for specific downstream genes by interpreting trained models. We applied HoloNet on three Visium datasets of breast cancer and liver cancer. The results detangled the multiple factors of FCEs by revealing how LR signals and cell types affect specific biological processes, and specified FCE-induced effects in each single cell. We conducted simulation experiments and showed that HoloNet is more reliable on LR prioritization in comparison with existing methods. HoloNet is a powerful tool to illustrate cell-cell communication landscapes and reveal vital FCEs that shape cellular phenotypes. HoloNet is available as a Python package at https://github.com/lhc17/HoloNet. 10.1093/bib/bbad359
Research Techniques Made Simple: Use of Imaging Mass Cytometry for Dermatological Research and Clinical Applications. The Journal of investigative dermatology Traditional immunohistochemistry (IHC) is inherently limited by its ability to analyze only several markers within a histological tissue section at a given time, which hinders in-depth characterization and phenotyping of tissues. Imaging mass cytometry (IMC), which combines IHC using metal-labeled antibodies with laser ablation and detection using mass cytometry by time-of-flight, overcomes this limitation with the capability to simultaneously analyze up to 40 protein markers to generate high-dimensional images from a single tissue section. IMC analysis preserves tissue architecture and spatial cellular relationships that would otherwise be lost or significantly altered in applications requiring tissue dissociation, such as flow cytometry or single-cell RNA sequencing. Resulting high-dimensional histological images permit spatially conserved analysis to identify unique cell populations, cellular interactions and avoidances, and insight into activation and behavioral status based on tissue location. IMC can be performed on both frozen and formalin-fixed paraffin-embedded tissue, allowing for previously banked samples to be analyzed and correlated with known clinical outcomes. Expectedly, IMC will change the landscape of investigative pathology, particularly when used in coordination with multiomic platforms to combine transcriptomic and proteomic data at a single-cell resolution. Here, we aim to highlight the potential utility of IMC within dermatologic research and clinical applications. 10.1016/j.jid.2020.12.008
Molecular biomarkers in pediatric glial tumors: a needed wind of change. Fontebasso Adam M,Bechet Denise,Jabado Nada Current opinion in oncology PURPOSE OF REVIEW:Glial tumors of the central nervous system (CNS) are the leading cause of cancer-related death and morbidity in children. Their diagnosis/prognosis relies mainly on clinical and histopathological factors. However, pathological grading is particularly challenging as there is substantial molecular heterogeneity in pediatric CNS tumors, which results in variable biological behavior in tumors with potentially identical histological diagnoses or limited reliable measures of classification for given subgroups. Novel molecular markers/pathways identified by integrated genomic/transcriptomic/epigenomic studies of cohorts of pediatric gliomas are revolutionizing this field and are summarized herein. RECENT FINDINGS:Studies of pediatric gliomas have identified unexpected oncogenic pathways implicated in gliomagenesis. These range from a single pathway/molecule defect such as abnormalities of the mitogen-activated-protein-kinase pathway considered to be a hallmark of pilocytic astrocytomas, to alterations in epigenomic modulators in higher-grade tumors. Importantly, the type, timing, and spatial clustering of these molecular alterations provide a better understanding of the pathogenesis of gliomas and critical markers for therapy that will help refine pathological grading. SUMMARY:Reappraisal of glioma classification using these novel biomarkers will likely change practice toward molecular pathology and their integration into clinical trials will enable personalized therapies based on the molecular fingerprint of individual tumors. 10.1097/CCO.0000000000000007
Murine breast cancers disorganize the liver transcriptome in a zonated manner. Communications biology The spatially organized gene expression program within the liver specifies hepatocyte functions according to their relative distances to the bloodstream (i.e., zonation), contributing to liver homeostasis. Despite the knowledge that solid cancers remotely disrupt liver homeostasis, it remains unexplored whether solid cancers affect liver zonation. Here, using spatial transcriptomics, we thoroughly investigate the abundance and zonation of hepatic genes in cancer-bearing mice. We find that breast cancers affect liver zonation in various distinct manners depending on biological pathways. Aspartate metabolism and triglyceride catabolic processes retain relatively intact zonation patterns, but the zonation of xenobiotic catabolic process genes exhibits a strong disruption. The acute phase response is induced in zonated manners. Furthermore, we demonstrate that breast cancers activate innate immune cells in particular neutrophils in distinct zonated manners, rather than in a uniform fashion within the liver. Collectively, breast cancers disorganize hepatic transcriptomes in zonated manners, thereby disrupting zonated functions of the liver. 10.1038/s42003-023-04479-w
Clinical molecular imaging. Miller Janet C,Thrall James H, Journal of the American College of Radiology : JACR This review summarizes the rapidly growing field of molecular imaging, the spatially localized and/or temporally resolved sensing of molecular and cellular processes in vivo. Molecular imaging is used to map the anatomic locations of specific molecules of interest within living tissue and has enormous potential as a powerful means to diagnose and monitor disease. Molecular imaging agents comprise a targeting component that confers localization and a component that enables external detectability with an imaging modality, such as PET, SPECT, MRI, optical, and ultrasound. The advantages and disadvantages of each of these modalities are discussed in regard to spatial resolution, temporal resolution, sensitivity, and cost. Molecular imaging agents can be divided into three categories, Type A, which bind directly to a target molecule, Type B, which are accumulated by molecular or cellular activity by the target, and Type C, which are undetectable when injected but can be imaged after they are activated by the target. The current status of clinical molecular imaging agents is presented as well as examples of some preclinical applications. The value of molecular imaging is illustrated by some examples for diseases such as cancer, neurological and psychiatric disorders, cardiovascular disease, infection and inflammation, and the monitoring of gene therapy and stem cell therapy. 10.1016/S1546-1440(03)00025-5
Spatial Transcriptomic Analysis of Focal and Normal Areas of Myocyte Disarray in Human Hypertrophic Cardiomyopathy. International journal of molecular sciences Hypertrophic Cardiomyopathy (HCM) is a common inherited disorder that can lead to heart failure and sudden cardiac death, characterized at the histological level by focal areas of myocyte disarray, hypertrophy and fibrosis, and only a few disease-targeted therapies exist. To identify the focal and spatially restricted alterations in the transcriptional pathways and reveal novel therapeutic targets, we performed a spatial transcriptomic analysis of the areas of focal myocyte disarray compared to areas of normal tissue using a commercially available platform (GeoMx, nanoString). We analyzed surgical myectomy tissue from four patients with HCM and the control interventricular septum tissue from two unused organ donor hearts that were free of cardiovascular disease. Histological sections were reviewed by an expert pathologist, and 72 focal areas with varying degrees of myocyte disarray (normal, mild, moderate, severe) were chosen for analysis. Areas of interest were interrogated with the Human Cancer Transcriptome Atlas designed to profile 1800 transcripts. Differential expression analysis revealed significant changes in gene expression between HCM and the control tissue, and functional enrichment analysis indicated that these genes were primarily involved in interferon production and mitochondrial energetics. Within the HCM tissue, differentially expressed genes between areas of normal and severe disarray were enriched for genes related to mitochondrial energetics and the extracellular matrix in severe disarray. An analysis of the gene expression of the ligand-receptor pair revealed that the HCM tissue exhibited downregulation of platelet-derived growth factor (PDGF), NOTCH, junctional adhesion molecule, and CD46 signaling while showing upregulation of fibronectin, CD99, cadherin, and amyloid precursor protein signaling. A deconvolution analysis utilizing the matched single nuclei RNA-sequencing (snRNA-seq) data to determine cell type composition in areas of interest revealed significant differences in fibroblast and vascular cell composition in areas of severe disarray when compared to normal areas in HCM samples. Cell composition in the normal areas of the control tissue was also divergent from the normal areas in HCM samples, which was consistent with the differential expression results. Overall, our data identify novel and potential disease-modifying targets for therapy in HCM. 10.3390/ijms241612625
Advanced Genomics-Based Approaches for Defining Allograft Rejection With Single Cell Resolution. Shi Tiffany,Roskin Krishna,Baker Brian M,Woodle E Steve,Hildeman David Frontiers in immunology Solid organ transplant recipients require long-term immunosuppression for prevention of rejection. Calcineurin inhibitor (CNI)-based immunosuppressive regimens have remained the primary means for immunosuppression for four decades now, yet little is known about their effects on graft resident and infiltrating immune cell populations. Similarly, the understanding of rejection biology under specific types of immunosuppression remains to be defined. Furthermore, development of innovative, rationally designed targeted therapeutics for mitigating or preventing rejection requires a fundamental understanding of the immunobiology that underlies the rejection process. The established use of microarray technologies in transplantation has provided great insight into gene transcripts associated with allograft rejection but does not characterize rejection on a single cell level. Therefore, the development of novel genomics tools, such as single cell sequencing techniques, combined with powerful bioinformatics approaches, has enabled characterization of immune processes at the single cell level. This can provide profound insights into the rejection process, including identification of resident and infiltrating cell transcriptomes, cell-cell interactions, and T cell receptor α/β repertoires. In this review, we discuss genomic analysis techniques, including microarray, bulk RNAseq (bulkSeq), single-cell RNAseq (scRNAseq), and spatial transcriptomic (ST) techniques, including considerations of their benefits and limitations. Further, other techniques, such as chromatin analysis assay for transposase-accessible chromatin sequencing (ATACseq), bioinformatic regulatory network analyses, and protein-based approaches are also examined. Application of these tools will play a crucial role in redefining transplant rejection with single cell resolution and likely aid in the development of future immunomodulatory therapies in solid organ transplantation. 10.3389/fimmu.2021.750754
Spatial sorting enables comprehensive characterization of liver zonation. Ben-Moshe Shani,Shapira Yonatan,Moor Andreas E,Manco Rita,Veg Tamar,Bahar Halpern Keren,Itzkovitz Shalev Nature metabolism The mammalian liver is composed of repeating hexagonal units termed lobules. Spatially resolved single-cell transcriptomics revealed that about half of hepatocyte genes are differentially expressed across the lobule, yet technical limitations impeded reconstructing similar global spatial maps of other hepatocyte features. Here, we show how zonated surface markers can be used to sort hepatocytes from defined lobule zones with high spatial resolution. We apply transcriptomics, miRNA array measurements and mass spectrometry proteomics to reconstruct spatial atlases of multiple zonated features. We demonstrate that protein zonation largely overlaps with mRNA zonation, with the periportal HNF4α as an exception. We identify zonation of miRNAs such as miR-122, and inverse zonation of miRNAs and their hepatocyte target genes, highlighting potential regulation of protein levels through zonated mRNA degradation. Among the targets we find the pericentral Wnt receptors Fzd7 and Fzd8 and the periportal Wnt inhibitors Tcf7l1 and Ctnnbip1. Our approach facilitates reconstructing spatial atlases of multiple cellular features in the liver and other structured tissues. 10.1038/s42255-019-0109-9
Transcriptome Wide Analysis of Magnetic Resonance Imaging-targeted Biopsy and Matching Surgical Specimens from High-risk Prostate Cancer Patients Treated with Radical Prostatectomy: The Target Must Be Hit. Radtke Jan Philipp,Takhar Mandeep,Bonekamp David,Kesch Claudia,Erho Nicholas,du Plessis Marguerite,Buerki Christine,Ong Kaye,Davicioni Elai,Hohenfellner Markus,Hadaschik Boris A European urology focus BACKGROUND:The most suspicious lesions on multiparametric magnetic resonance imaging (MRI) may be representative of final pathology. OBJECTIVE:We connect imaging with high-precision spatial annotation of biopsies and genomic cancer signatures to compare the genomic signals of the index lesion and biopsy cores of adjacent and far away locations. DESIGN, SETTING, AND PARTICIPANTS:Eleven patients diagnosed with high-risk prostate cancer on MRI/transrectal ultrasound-fusion biopsy (Bx) and treated with radical prostatectomy (RP). Five tissue specimens were collected from each patient. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS:Whole transcriptome RNA-expression was profiled for each sample. Genomic signatures were used to compare signals in MRI invisible versus visible foci using Pearson's correlation and to assess intratumoral heterogeneity using hierarchical clustering. RESULTS AND LIMITATIONS:Ten RP and 27 Bx-samples passed quality control. Gene expression between RP and index Bx, but not adjacent benign samples, was highly correlated. Genomic Gleason grade classifier features measured across the different samples showed concordant expression across Bx and RP tumor samples, while an inverse expression pattern was observed between tumor and benign samples indicating the lack of a strong field-effect. The distribution of low and high Prostate Imaging Reporting and Data System (PI-RADS) samples was 10 and 11, respectively. Genomics of all low PI-RADS samples resembled benign tissue and most high PI-RADS samples resembled cancer tissue. A strong association was observed between PI-RADS version 2 and Decipher as well as the genomic Gleason grade classifier score. Clustering analysis showed that most samples cluster tightly by patient. One patient showed unique tumor biology in index versus secondary lesion suggesting the presence of intrapatient heterogeneity and the utility in profiling multiple foci identified by MRI. CONCLUSIONS:MRI-targeted Bx-genomics show excellent correlation with RP-genomics and confirm the information captured by PI-RADS. Sampling of the target lesion must be precise as correlation between index and benign lesions was not seen. PATIENT SUMMARY:In this report, we tested if targeted prostate sampling using magnetic resonance imaging-fusion biopsy allows to genetically describe index tumors of prostate cancer. We found that imaging genomics correlated well with final prostatectomy provided that the target is hit precisely. 10.1016/j.euf.2017.01.005
Simultaneous profiling of host expression and microbial abundance by spatial metatranscriptome sequencing. Genome research We developed an analysis pipeline that can extract microbial sequences from spatial transcriptomic (ST) data and assign taxonomic labels, generating a spatial microbial abundance matrix in addition to the default host expression matrix, enabling simultaneous analysis of host expression and microbial distribution. We called the pipeline spatial metatranscriptome (SMT) and applied it on both human and murine intestinal sections and validated the spatial microbial abundance information with alternative assays. Biological insights were gained from these novel data that showed host-microbe interaction at various spatial scales. Finally, we tested experimental modification that can increase microbial capture while preserving host spatial expression quality and, by use of positive controls, quantitatively showed the capture efficiency and recall of our methods. This proof-of-concept work shows the feasibility of SMT analysis and paves the way for further experimental optimization and application. 10.1101/gr.277178.122
Maximizing the Utility of Transcriptomics Data in Inflammatory Skin Diseases. Frontiers in immunology Inflammatory skin diseases are induced by disorders of the host defense system of the skin, which is composed of a barrier, innate and acquired immunity, as well as the cutaneous microbiome. These disorders are characterized by recurrent cutaneous lesions and intense itch, which seriously affecting life quality of people across all ages and ethnicities. To elucidate molecular factors for typical inflammatory skin diseases (such as psoriasis and atopic dermatitis), transcriptomic profiling assays have been largely performed. Additionally, single-cell RNA sequencing (scRNA-seq) as well as spatial transcriptomic profiling have revealed multiple potential translational targets and offered guides to improve diagnosis and treatment strategies for inflammatory skin diseases. High-throughput transcriptomics data has shown unprecedented power to disclose the complex pathophysiology of inflammatory skin diseases. Here, we will summarize discoveries from transcriptomics data and discuss how to maximize the transcriptomics data to propel the development of diagnostic biomarkers and therapeutic targets in inflammatory skin diseases. 10.3389/fimmu.2021.761890
Microbial Pathogenesis in the Era of Spatial Omics. Infection and immunity The biology of a cell, whether it is a unicellular organism or part of a multicellular network, is influenced by cell type, temporal changes in cell state, and the cell's environment. Spatial cues play a critical role in the regulation of microbial pathogenesis strategies. Information about where the pathogen is-in a tissue or in proximity to a host cell-regulates gene expression and the compartmentalization of gene products in the microbe and the host. Our understanding of host and pathogen identity has bloomed with the accessibility of transcriptomics and proteomics techniques. A missing piece of the puzzle has been our ability to evaluate global transcript and protein expression in the context of the subcellular niche, primary cell, or native tissue environment during infection. This barrier is now lower with the advent of new spatial omics techniques to understand how location regulates cellular functions. This review will discuss how recent advances in spatial proteomics and transcriptomics approaches can address outstanding questions in microbial pathogenesis. 10.1128/iai.00442-22
Exercise Training Reduces the Inflammatory Response and Promotes Intestinal Mucosa-Associated Immunity in Lynch Syndrome. Clinical cancer research : an official journal of the American Association for Cancer Research PURPOSE:Lynch syndrome (LS) is a hereditary condition with a high lifetime risk of colorectal and endometrial cancers. Exercise is a non-pharmacologic intervention to reduce cancer risk, though its impact on patients with LS has not been prospectively studied. Here, we evaluated the impact of a 12-month aerobic exercise cycling intervention in the biology of the immune system in LS carriers. PATIENTS AND METHODS:To address this, we enrolled 21 patients with LS onto a non-randomized, sequential intervention assignation, clinical trial to assess the effect of a 12-month exercise program that included cycling classes 3 times weekly for 45 minutes versus usual care with a one-time exercise counseling session as control. We analyzed the effects of exercise on cardiorespiratory fitness, circulating, and colorectal-tissue biomarkers using metabolomics, gene expression by bulk mRNA sequencing, and spatial transcriptomics by NanoString GeoMx. RESULTS:We observed a significant increase in oxygen consumption (VO2peak) as a primary outcome of the exercise and a decrease in inflammatory markers (prostaglandin E) in colon and blood as the secondary outcomes in the exercise versus usual care group. Gene expression profiling and spatial transcriptomics on available colon biopsies revealed an increase in the colonic mucosa levels of natural killer and CD8+ T cells in the exercise group that were further confirmed by IHC studies. CONCLUSIONS:Together these data have important implications for cancer interception in LS, and document for the first-time biological effects of exercise in the immune system of a target organ in patients at-risk for cancer. 10.1158/1078-0432.CCR-23-0088
Neighbor-specific gene expression revealed from physically interacting cells during mouse embryonic development. Proceedings of the National Academy of Sciences of the United States of America Development of multicellular organisms is orchestrated by persistent cell-cell communication between neighboring partners. Direct interaction between different cell types can induce molecular signals that dictate lineage specification and cell fate decisions. Current single-cell RNA-seq technology cannot adequately analyze cell-cell contact-dependent gene expression, mainly due to the loss of spatial information. To overcome this obstacle and resolve cell-cell contact-specific gene expression during embryogenesis, we performed RNA sequencing of physically interacting cells (PIC-seq) and assessed them alongside similar single-cell transcriptomes derived from developing mouse embryos between embryonic day (E) 7.5 and E9.5. Analysis of the PIC-seq data identified gene expression signatures that were dependent on the presence of specific neighboring cell types. Our computational predictions, validated experimentally, demonstrated that neural progenitor (NP) cells upregulate and genes, when exclusively interacting with definitive endoderm (DE) cells. Moreover, there was a reciprocal impact on the transcriptome of DE cells, as they tend to upregulate and when in contact with NP cells. Using individual cell transcriptome data, we formulated a means of computationally predicting the impact of one cell type on the transcriptome of its neighboring cell types. We have further developed a distinctive spatial-t-distributed stochastic neighboring embedding to display the pseudospatial distribution of cells in a 2-dimensional space. In summary, we describe an innovative approach to study contact-specific gene regulation during embryogenesis. 10.1073/pnas.2205371120
Single cell transcriptomics reveals spatial and temporal dynamics of gene expression in the developing mouse spinal cord. Development (Cambridge, England) The coordinated spatial and temporal regulation of gene expression in the vertebrate neural tube determines the identity of neural progenitors and the function and physiology of the neurons they generate. Progress has been made deciphering the gene regulatory programmes that are responsible for this process; however, the complexity of the tissue has hampered the systematic analysis of the network and the underlying mechanisms. To address this, we used single cell mRNA sequencing to profile cervical and thoracic regions of the developing mouse neural tube between embryonic days 9.5-13.5. We confirmed that the data accurately recapitulates neural tube development, allowing us to identify new markers for specific progenitor and neuronal populations. In addition, the analysis highlighted a previously underappreciated temporal component to the mechanisms that generate neuronal diversity, and revealed common features in the sequence of transcriptional events that lead to the differentiation of specific neuronal subtypes. Together, the data offer insight into the mechanisms that are responsible for neuronal specification and provide a compendium of gene expression for classifying spinal cord cell types that will support future studies of neural tube development, function and disease. 10.1242/dev.173807
The proteome and transcriptome of stress granules and P bodies during human T lymphocyte activation. Cell reports Stress granules (SGs) and processing bodies (PBs) are membraneless cytoplasmic assemblies regulating mRNAs under environmental stress such as viral infections, neurological disorders, or cancer. Upon antigen stimulation, T lymphocytes mediate their immune functions under regulatory mechanisms involving SGs and PBs. However, the impact of T cell activation on such complexes in terms of formation, constitution, and relationship remains unknown. Here, by combining proteomic, transcriptomic, and immunofluorescence approaches, we simultaneously characterized the SGs and PBs from primary human T lymphocytes pre and post stimulation. The identification of the proteomes and transcriptomes of SGs and PBs indicate an unanticipated molecular and functional complementarity. Notwithstanding, these granules keep distinct spatial organizations and abilities to interact with mRNAs. This comprehensive characterization of the RNP granule proteomic and transcriptomic landscapes provides a unique resource for future investigations on SGs and PBs in T lymphocytes. 10.1016/j.celrep.2023.112211
A harmonized atlas of mouse spinal cord cell types and their spatial organization. Nature communications Single-cell RNA sequencing data can unveil the molecular diversity of cell types. Cell type atlases of the mouse spinal cord have been published in recent years but have not been integrated together. Here, we generate an atlas of spinal cell types based on single-cell transcriptomic data, unifying the available datasets into a common reference framework. We report a hierarchical structure of postnatal cell type relationships, with location providing the highest level of organization, then neurotransmitter status, family, and finally, dozens of refined populations. We validate a combinatorial marker code for each neuronal cell type and map their spatial distributions in the adult spinal cord. We also show complex lineage relationships among postnatal cell types. Additionally, we develop an open-source cell type classifier, SeqSeek, to facilitate the standardization of cell type identification. This work provides an integrated view of spinal cell types, their gene expression signatures, and their molecular organization. 10.1038/s41467-021-25125-1
Spatial transcriptomic analysis of virtual prostate biopsy reveals confounding effect of tissue heterogeneity on genomic signatures. Molecular cancer Genetic signatures have added a molecular dimension to prognostics and therapeutic decision-making. However, tumour heterogeneity in prostate cancer and current sampling methods could confound accurate assessment. Based on previously published spatial transcriptomic data from multifocal prostate cancer, we created virtual biopsy models that mimic conventional biopsy placement and core size. We then analysed the gene expression of different prognostic signatures (OncotypeDx®, Decipher®, Prostadiag®) using a step-wise approach with increasing resolution from pseudo-bulk analysis of the whole biopsy, to differentiation by tissue subtype (benign, stroma, tumour), followed by distinct tumour grade and finally clonal resolution. The gene expression profile of virtual tumour biopsies revealed clear differences between grade groups and tumour clones, compared to a benign control, which were not reflected in bulk analyses. This suggests that bulk analyses of whole biopsies or tumour-only areas, as used in clinical practice, may provide an inaccurate assessment of gene profiles. The type of tissue, the grade of the tumour and the clonal composition all influence the gene expression in a biopsy. Clinical decision making based on biopsy genomics should be made with caution while we await more precise targeting and cost-effective spatial analyses. 10.1186/s12943-023-01863-2
A spatial vascular transcriptomic, proteomic, and phosphoproteomic atlas unveils an angiocrine Tie-Wnt signaling axis in the liver. Inverso Donato,Shi Jingjing,Lee Ki Hong,Jakab Moritz,Ben-Moshe Shani,Kulkarni Shubhada R,Schneider Martin,Wang Guanxiong,Komeili Marziyeh,Vélez Paula Argos,Riedel Maria,Spegg Carleen,Ruppert Thomas,Schaeffer-Reiss Christine,Helm Dominic,Singh Indrabahadur,Boutros Michael,Chintharlapalli Sudhakar,Heikenwalder Mathias,Itzkovitz Shalev,Augustin Hellmut G Developmental cell Single-cell transcriptomics (scRNA-seq) has revolutionized the understanding of the spatial architecture of tissue structure and function. Advancing the "transcript-centric" view of scRNA-seq analyses is presently restricted by the limited resolution of proteomics and genome-wide techniques to analyze post-translational modifications. Here, by combining spatial cell sorting with transcriptomics and quantitative proteomics/phosphoproteomics, we established the spatially resolved proteome landscape of the liver endothelium, yielding deep mechanistic insight into zonated vascular signaling mechanisms. Phosphorylation of receptor tyrosine kinases was detected preferentially in the central vein area, resulting in an atypical enrichment of tyrosine phosphorylation. Prototypic biological validation identified Tie receptor signaling as a selective and specific regulator of vascular Wnt activity orchestrating angiocrine signaling, thereby controlling hepatocyte function during liver regeneration. Taken together, the study has yielded fundamental insight into the spatial organization of liver endothelial cell signaling. Spatial sorting may be employed as a universally adaptable strategy for multiomic analyses of scRNA-seq-defined cellular (sub)-populations. 10.1016/j.devcel.2021.05.001
Stomach encyclopedia: Combined single-cell and spatial transcriptomics reveal cell diversity and homeostatic regulation of human stomach. Cell reports The stomach is an important digestive organ with various biological functions. However, because of the complexity of its cellular and glandular composition, its precise cellular biology has yet to be elucidated. In this study, we conducted single-cell RNA sequencing (scRNA-seq) and subcellular-level spatial transcriptomics analysis of the human stomach and constructed the largest dataset to date: a stomach encyclopedia. This dataset consists of approximately 380,000 cells from scRNA-seq and the spatial transcriptome, enabling integrated analyses of transcriptional and spatial information of gastric and metaplastic cells. This analysis identified LEFTY1 as an uncharacterized stem cell marker, which was confirmed through lineage tracing analysis. A wide variety of cell-cell interactions between epithelial and stromal cells, including PDGFRABMP4WNT5A fibroblasts, was highlighted in the developmental switch of intestinal metaplasia. Our extensive dataset will function as a fundamental resource in investigations of the stomach, including studies of development, aging, and carcinogenesis. 10.1016/j.celrep.2023.113236
Immune checkpoint blockade induces distinct alterations in the microenvironments of primary and metastatic brain tumors. The Journal of clinical investigation In comparison with responses in recurrent glioblastoma (rGBM), the intracranial response of brain metastases (BrM) to immune checkpoint blockade (ICB) is less well studied. Here, we present an integrated single-cell RNA-Seq (scRNA-Seq) study of 19 ICB-naive and 9 ICB-treated BrM samples from our own and published data sets. We compared them with our previously published scRNA-Seq data from rGBM and found that ICB led to more prominent T cell infiltration into BrM than rGBM. These BrM-infiltrating T cells exhibited a tumor-specific phenotype and displayed greater activated/exhausted features. We also used multiplex immunofluorescence and spatial transcriptomics to reveal that ICB reduced a distinct CD206+ macrophage population in the perivascular space, which may modulate T cell entry into BrM. Furthermore, we identified a subset of progenitor exhausted T cells that correlated with longer overall survival in BrM patients. Our study provides a comprehensive immune cellular landscape of ICB's effect on metastatic brain tumors and offers insights into potential strategies for improving ICB efficacy for brain tumor patients. 10.1172/JCI169314
Single-cell transcriptomes of the human skin reveal age-related loss of fibroblast priming. Solé-Boldo Llorenç,Raddatz Günter,Schütz Sabrina,Mallm Jan-Philipp,Rippe Karsten,Lonsdorf Anke S,Rodríguez-Paredes Manuel,Lyko Frank Communications biology Fibroblasts are an essential cell population for human skin architecture and function. While fibroblast heterogeneity is well established, this phenomenon has not been analyzed systematically yet. We have used single-cell RNA sequencing to analyze the transcriptomes of more than 5,000 fibroblasts from a sun-protected area in healthy human donors. Our results define four main subpopulations that can be spatially localized and show differential secretory, mesenchymal and pro-inflammatory functional annotations. Importantly, we found that this fibroblast 'priming' becomes reduced with age. We also show that aging causes a substantial reduction in the predicted interactions between dermal fibroblasts and other skin cells, including undifferentiated keratinocytes at the dermal-epidermal junction. Our work thus provides evidence for a functional specialization of human dermal fibroblasts and identifies the partial loss of cellular identity as an important age-related change in the human dermis. These findings have important implications for understanding human skin aging and its associated phenotypes. 10.1038/s42003-020-0922-4
Spatial transcriptomics of dorsal root ganglia identifies molecular signatures of human nociceptors. Science translational medicine Nociceptors are specialized sensory neurons that detect damaging or potentially damaging stimuli and are found in the dorsal root ganglia (DRG) and trigeminal ganglia. These neurons are critical for the generation of neuronal signals that ultimately create the perception of pain. Nociceptors are also primary targets for treating acute and chronic pain. Single-cell transcriptomics on mouse nociceptors has transformed our understanding of pain mechanisms. We sought to generate equivalent information for human nociceptors with the goal of identifying transcriptomic signatures of nociceptors, identifying species differences and potential drug targets. We used spatial transcriptomics to molecularly characterize transcriptomes of single DRG neurons from eight organ donors. We identified 12 clusters of human sensory neurons, 5 of which are C nociceptors, as well as 1 C low-threshold mechanoreceptors (LTMRs), 1 Aβ nociceptor, 2 Aδ, 2 Aβ, and 1 proprioceptor subtypes. By focusing on expression profiles for ion channels, G protein-coupled receptors (GPCRs), and other pharmacological targets, we provided a rich map of potential drug targets in the human DRG with direct comparison to mouse sensory neuron transcriptomes. We also compared human DRG neuronal subtypes to nonhuman primates showing conserved patterns of gene expression among many cell types but divergence among specific nociceptor subsets. Last, we identified sex differences in human DRG subpopulation transcriptomes, including a marked increase in calcitonin-related polypeptide alpha () expression in female pruritogen receptor-enriched nociceptors. This comprehensive spatial characterization of human nociceptors might open the door to development of better treatments for acute and chronic pain disorders. 10.1126/scitranslmed.abj8186
Integrating spatial gene expression and breast tumour morphology via deep learning. He Bryan,Bergenstråhle Ludvig,Stenbeck Linnea,Abid Abubakar,Andersson Alma,Borg Åke,Maaskola Jonas,Lundeberg Joakim,Zou James Nature biomedical engineering Spatial transcriptomics allows for the measurement of RNA abundance at a high spatial resolution, making it possible to systematically link the morphology of cellular neighbourhoods and spatially localized gene expression. Here, we report the development of a deep learning algorithm for the prediction of local gene expression from haematoxylin-and-eosin-stained histopathology images using a new dataset of 30,612 spatially resolved gene expression data matched to histopathology images from 23 patients with breast cancer. We identified over 100 genes, including known breast cancer biomarkers of intratumoral heterogeneity and the co-localization of tumour growth and immune activation, the expression of which can be predicted from the histopathology images at a resolution of 100 µm. We also show that the algorithm generalizes well to The Cancer Genome Atlas and to other breast cancer gene expression datasets without the need for re-training. Predicting the spatially resolved transcriptome of a tissue directly from tissue images may enable image-based screening for molecular biomarkers with spatial variation. 10.1038/s41551-020-0578-x
A Spatiotemporal Organ-Wide Gene Expression and Cell Atlas of the Developing Human Heart. Asp Michaela,Giacomello Stefania,Larsson Ludvig,Wu Chenglin,Fürth Daniel,Qian Xiaoyan,Wärdell Eva,Custodio Joaquin,Reimegård Johan,Salmén Fredrik,Österholm Cecilia,Ståhl Patrik L,Sundström Erik,Åkesson Elisabet,Bergmann Olaf,Bienko Magda,Månsson-Broberg Agneta,Nilsson Mats,Sylvén Christer,Lundeberg Joakim Cell The process of cardiac morphogenesis in humans is incompletely understood. Its full characterization requires a deep exploration of the organ-wide orchestration of gene expression with a single-cell spatial resolution. Here, we present a molecular approach that reveals the comprehensive transcriptional landscape of cell types populating the embryonic heart at three developmental stages and that maps cell-type-specific gene expression to specific anatomical domains. Spatial transcriptomics identified unique gene profiles that correspond to distinct anatomical regions in each developmental stage. Human embryonic cardiac cell types identified by single-cell RNA sequencing confirmed and enriched the spatial annotation of embryonic cardiac gene expression. In situ sequencing was then used to refine these results and create a spatial subcellular map for the three developmental phases. Finally, we generated a publicly available web resource of the human developing heart to facilitate future studies on human cardiogenesis. 10.1016/j.cell.2019.11.025
Dissecting mammalian reproduction with spatial transcriptomics. Human reproduction update BACKGROUND:Mammalian reproduction requires the fusion of two specialized cells: an oocyte and a sperm. In addition to producing gametes, the reproductive system also provides the environment for the appropriate development of the embryo. Deciphering the reproductive system requires understanding the functions of each cell type and cell-cell interactions. Recent single-cell omics technologies have provided insights into the gene regulatory network in discrete cellular populations of both the male and female reproductive systems. However, these approaches cannot examine how the cellular states of the gametes or embryos are regulated through their interactions with neighboring somatic cells in the native tissue environment owing to tissue disassociations. Emerging spatial omics technologies address this challenge by preserving the spatial context of the cells to be profiled. These technologies hold the potential to revolutionize our understanding of mammalian reproduction. OBJECTIVE AND RATIONALE:We aim to review the state-of-the-art spatial transcriptomics (ST) technologies with a focus on highlighting the novel biological insights that they have helped to reveal about the mammalian reproductive systems in the context of gametogenesis, embryogenesis, and reproductive pathologies. We also aim to discuss the current challenges of applying ST technologies in reproductive research and provide a sneak peek at what the field of spatial omics can offer for the reproduction community in the years to come. SEARCH METHODS:The PubMed database was used in the search for peer-reviewed research articles and reviews using combinations of the following terms: 'spatial omics', 'fertility', 'reproduction', 'gametogenesis', 'embryogenesis', 'reproductive cancer', 'spatial transcriptomics', 'spermatogenesis', 'ovary', 'uterus', 'cervix', 'testis', and other keywords related to the subject area. All relevant publications until April 2023 were critically evaluated and discussed. OUTCOMES:First, an overview of the ST technologies that have been applied to studying the reproductive systems was provided. The basic design principles and the advantages and limitations of these technologies were discussed and tabulated to serve as a guide for researchers to choose the best-suited technologies for their own research. Second, novel biological insights into mammalian reproduction, especially human reproduction revealed by ST analyses, were comprehensively reviewed. Three major themes were discussed. The first theme focuses on genes with non-random spatial expression patterns with specialized functions in multiple reproductive systems; The second theme centers around functionally interacting cell types which are often found to be spatially clustered in the reproductive tissues; and the thrid theme discusses pathological states in reproductive systems which are often associated with unique cellular microenvironments. Finally, current experimental and computational challenges of applying ST technologies to studying mammalian reproduction were highlighted, and potential solutions to tackle these challenges were provided. Future directions in the development of spatial omics technologies and how they will benefit the field of human reproduction were discussed, including the capture of cellular and tissue dynamics, multi-modal molecular profiling, and spatial characterization of gene perturbations. WIDER IMPLICATIONS:Like single-cell technologies, spatial omics technologies hold tremendous potential for providing significant and novel insights into mammalian reproduction. Our review summarizes these novel biological insights that ST technologies have provided while shedding light on what is yet to come. Our review provides reproductive biologists and clinicians with a much-needed update on the state of art of ST technologies. It may also facilitate the adoption of cutting-edge spatial technologies in both basic and clinical reproductive research. 10.1093/humupd/dmad017
Molecular cartography uncovers evolutionary and microenvironmental dynamics in sporadic colorectal tumors. Cell Colorectal cancer exhibits dynamic cellular and genetic heterogeneity during progression from precursor lesions toward malignancy. Analysis of spatial multi-omic data from 31 human colorectal specimens enabled phylogeographic mapping of tumor evolution that revealed individualized progression trajectories and accompanying microenvironmental and clonal alterations. Phylogeographic mapping ordered genetic events, classified tumors by their evolutionary dynamics, and placed clonal regions along global pseudotemporal progression trajectories encompassing the chromosomal instability (CIN+) and hypermutated (HM) pathways. Integrated single-cell and spatial transcriptomic data revealed recurring epithelial programs and infiltrating immune states along progression pseudotime. We discovered an immune exclusion signature (IEX), consisting of extracellular matrix regulators DDR1, TGFBI, PAK4, and DPEP1, that charts with CIN+ tumor progression, is associated with reduced cytotoxic cell infiltration, and shows prognostic value in independent cohorts. This spatial multi-omic atlas provides insights into colorectal tumor-microenvironment co-evolution, serving as a resource for stratification and targeted treatments. 10.1016/j.cell.2023.11.006
Multiplex digital spatial profiling of proteins and RNA in fixed tissue. Merritt Christopher R,Ong Giang T,Church Sarah E,Barker Kristi,Danaher Patrick,Geiss Gary,Hoang Margaret,Jung Jaemyeong,Liang Yan,McKay-Fleisch Jill,Nguyen Karen,Norgaard Zach,Sorg Kristina,Sprague Isaac,Warren Charles,Warren Sarah,Webster Philippa J,Zhou Zoey,Zollinger Daniel R,Dunaway Dwayne L,Mills Gordon B,Beechem Joseph M Nature biotechnology Digital Spatial Profiling (DSP) is a method for highly multiplex spatial profiling of proteins or RNAs suitable for use on formalin-fixed, paraffin-embedded (FFPE) samples. The approach relies on (1) multiplexed readout of proteins or RNAs using oligonucleotide tags; (2) oligonucleotide tags attached to affinity reagents (antibodies or RNA probes) through a photocleavable (PC) linker; and (3) photocleaving light projected onto the tissue sample to release PC oligonucleotides in any spatial pattern across a region of interest (ROI) covering 1 to ~5,000 cells. DSP is capable of single-cell sensitivity within an ROI using the antibody readout, with RNA detection feasible down to ~600 individual mRNA transcripts. We show spatial profiling of up to 44 proteins and 96 genes (928 RNA probes) in lymphoid, colorectal tumor and autoimmune tissues by using the nCounter system and 1,412 genes (4,998 RNA probes) by using next-generation sequencing (NGS). DSP may be used to profile not only proteins and RNAs in biobanked samples but also immune markers in patient samples, with potential prognostic and predictive potential for clinical decision-making. 10.1038/s41587-020-0472-9
A spatially resolved single-cell genomic atlas of the adult human breast. Nature The adult human breast is comprised of an intricate network of epithelial ducts and lobules that are embedded in connective and adipose tissue. Although most previous studies have focused on the breast epithelial system, many of the non-epithelial cell types remain understudied. Here we constructed the comprehensive Human Breast Cell Atlas (HBCA) at single-cell and spatial resolution. Our single-cell transcriptomics study profiled 714,331 cells from 126 women, and 117,346 nuclei from 20 women, identifying 12 major cell types and 58 biological cell states. These data reveal abundant perivascular, endothelial and immune cell populations, and highly diverse luminal epithelial cell states. Spatial mapping using four different technologies revealed an unexpectedly rich ecosystem of tissue-resident immune cells, as well as distinct molecular differences between ductal and lobular regions. Collectively, these data provide a reference of the adult normal breast tissue for studying mammary biology and diseases such as breast cancer. 10.1038/s41586-023-06252-9
An integrated single cell and spatial transcriptomic map of human white adipose tissue. Nature communications To date, single-cell studies of human white adipose tissue (WAT) have been based on small cohort sizes and no cellular consensus nomenclature exists. Herein, we performed a comprehensive meta-analysis of publicly available and newly generated single-cell, single-nucleus, and spatial transcriptomic results from human subcutaneous, omental, and perivascular WAT. Our high-resolution map is built on data from ten studies and allowed us to robustly identify >60 subpopulations of adipocytes, fibroblast and adipogenic progenitors, vascular, and immune cells. Using these results, we deconvolved spatial and bulk transcriptomic data from nine additional cohorts to provide spatial and clinical dimensions to the map. This identified cell-cell interactions as well as relationships between specific cell subtypes and insulin resistance, dyslipidemia, adipocyte volume, and lipolysis upon long-term weight changes. Altogether, our meta-map provides a rich resource defining the cellular and microarchitectural landscape of human WAT and describes the associations between specific cell types and metabolic states. 10.1038/s41467-023-36983-2
Dissecting mammalian spermatogenesis using spatial transcriptomics. Cell reports Single-cell RNA sequencing has revealed extensive molecular diversity in gene programs governing mammalian spermatogenesis but fails to delineate their dynamics in the native context of seminiferous tubules, the spatially confined functional units of spermatogenesis. Here, we use Slide-seq, a spatial transcriptomics technology, to generate an atlas that captures the spatial gene expression patterns at near-single-cell resolution in the mouse and human testis. Using Slide-seq data, we devise a computational framework that accurately localizes testicular cell types in individual seminiferous tubules. Unbiased analysis systematically identifies spatially patterned genes and gene programs. Combining Slide-seq with targeted in situ RNA sequencing, we demonstrate significant differences in the cellular compositions of spermatogonial microenvironment between mouse and human testes. Finally, a comparison of the spatial atlas generated from the wild-type and diabetic mouse testis reveals a disruption in the spatial cellular organization of seminiferous tubules as a potential mechanism of diabetes-induced male infertility. 10.1016/j.celrep.2021.109915
The power of combined spatial transcriptomics and genome wide association studies (GWAS) approaches to heritable prostate cancer. Journal of cell communication and signaling 10.1007/s12079-023-00724-3
The end of the beginning: application of single-cell sequencing to chronic lymphocytic leukemia. Blood Single-cell analysis has emerged over the past decade as a transformative technology informative for the systematic analysis of complex cell populations such as in cancers and the tumor immune microenvironment. The methodologic and analytical advancements in this realm have evolved rapidly, scaling from but a few cells at its outset to the current capabilities of processing and analyzing hundreds of thousands of individual cells at a time. The types of profiling attainable at individual cell resolution now range from genetic and transcriptomic characterization and extend to epigenomic and spatial analysis. Additionally, the increasing ability to achieve multiomic integration of these data layers now yields ever richer insights into diverse molecular disease subtypes and the patterns of cellular circuitry on a per-cancer basis. Over the years, chronic lymphocytic leukemia (CLL) consistently has been at the forefront of genomic investigation, given the ready accessibility of pure leukemia cells and immune cells from circulating blood of patients with this disease. Herein, we review the recent forays into the application of single-cell analysis to CLL, which are already revealing a new understanding of the natural progression of CLL, the impact of novel therapies, and the interactions with coevolving nonmalignant immune cell populations. As we emerge from the end of the beginning of this technologic revolution, CLL stands poised to reap the benefits of single-cell analysis from the standpoints of uncovering fresh fundamental biological knowledge and of providing a path to devising regimens of personalized diagnosis, treatment, and monitoring. 10.1182/blood.2021014669
The Desmoplastic Stroma of Pancreatic Cancer: Multilayered Levels of Heterogeneity, Clinical Significance, and Therapeutic Opportunities. Cancers Pancreatic cancer remains one of the most lethal malignancies and is becoming a dramatically increasing cause of cancer-related mortality worldwide. Abundant desmoplastic stroma is a histological hallmark of pancreatic ductal adenocarcinoma. Emerging evidence suggests a promising therapeutic effect of several stroma-modifying therapies that target desmoplastic stromal elements in the pancreatic cancer microenvironment. The evidence also unveils multifaceted roles of cancer-associated fibroblasts (CAFs) in manipulating pancreatic cancer progression, immunity, and chemotherapeutic response. Current state-of-the-art technologies, including single-cell transcriptomics and multiplexed tissue imaging techniques, have provided a more profound knowledge of CAF heterogeneity in real-world specimens from pancreatic cancer patients, as well as in genetically engineered mouse models. In this review, we describe recent advances in the understanding of the molecular pathology bases of pancreatic cancer desmoplastic stroma at multilayered levels of heterogeneity, namely, (1) variations in cellular and non-cellular members, including CAF subtypes and extracellular matrix (ECM) proteins; (2) geographical heterogeneity in relation to cell-cell interactions and signaling pathways at niche levels and spatial heterogeneity at locoregional levels or organ levels; and (3) intertumoral stromal heterogeneity at individual levels. This review further discusses the clinicopathological significance of desmoplastic stroma and the potential opportunities for stroma-targeted therapies against this lethal malignancy. 10.3390/cancers14133293
Technological advances in cancer immunity: from immunogenomics to single-cell analysis and artificial intelligence. Signal transduction and targeted therapy Immunotherapies play critical roles in cancer treatment. However, given that only a few patients respond to immune checkpoint blockades and other immunotherapeutic strategies, more novel technologies are needed to decipher the complicated interplay between tumor cells and the components of the tumor immune microenvironment (TIME). Tumor immunomics refers to the integrated study of the TIME using immunogenomics, immunoproteomics, immune-bioinformatics, and other multi-omics data reflecting the immune states of tumors, which has relied on the rapid development of next-generation sequencing. High-throughput genomic and transcriptomic data may be utilized for calculating the abundance of immune cells and predicting tumor antigens, referring to immunogenomics. However, as bulk sequencing represents the average characteristics of a heterogeneous cell population, it fails to distinguish distinct cell subtypes. Single-cell-based technologies enable better dissection of the TIME through precise immune cell subpopulation and spatial architecture investigations. In addition, radiomics and digital pathology-based deep learning models largely contribute to research on cancer immunity. These artificial intelligence technologies have performed well in predicting response to immunotherapy, with profound significance in cancer therapy. In this review, we briefly summarize conventional and state-of-the-art technologies in the field of immunogenomics, single-cell and artificial intelligence, and present prospects for future research. 10.1038/s41392-021-00729-7
Subcellular Transcriptomics and Proteomics: A Comparative Methods Review. Molecular & cellular proteomics : MCP The internal environment of cells is molecularly crowded, which requires spatial organization via subcellular compartmentalization. These compartments harbor specific conditions for molecules to perform their biological functions, such as coordination of the cell cycle, cell survival, and growth. This compartmentalization is also not static, with molecules trafficking between these subcellular neighborhoods to carry out their functions. For example, some biomolecules are multifunctional, requiring an environment with differing conditions or interacting partners, and others traffic to export such molecules. Aberrant localization of proteins or RNA species has been linked to many pathological conditions, such as neurological, cancer, and pulmonary diseases. Differential expression studies in transcriptomics and proteomics are relatively common, but the majority have overlooked the importance of subcellular information. In addition, subcellular transcriptomics and proteomics data do not always colocate because of the biochemical processes that occur during and after translation, highlighting the complementary nature of these fields. In this review, we discuss and directly compare the current methods in spatial proteomics and transcriptomics, which include sequencing- and imaging-based strategies, to give the reader an overview of the current tools available. We also discuss current limitations of these strategies as well as future developments in the field of spatial -omics. 10.1016/j.mcpro.2021.100186
Single-cell genomics and spatial transcriptomics: Discovery of novel cell states and cellular interactions in liver physiology and disease biology. Saviano Antonio,Henderson Neil C,Baumert Thomas F Journal of hepatology Transcriptome analysis enables the study of gene expression in human tissues and is a valuable tool to characterise liver function and gene expression dynamics during liver disease, as well as to identify prognostic markers or signatures, and to facilitate discovery of new therapeutic targets. In contrast to whole tissue RNA sequencing analysis, single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics enables the study of transcriptional activity at the single cell or spatial level. ScRNA-seq has paved the way for the discovery of previously unknown cell types and subtypes in normal and diseased liver, facilitating the study of rare cells (such as liver progenitor cells) and the functional roles of non-parenchymal cells in chronic liver disease and cancer. By adding spatial information to scRNA-seq data, spatial transcriptomics has transformed our understanding of tissue functional organisation and cell-to-cell interactions in situ. These approaches have recently been applied to investigate liver regeneration, organisation and function of hepatocytes and non-parenchymal cells, and to profile the single cell landscape of chronic liver diseases and cancer. Herein, we review the principles and technologies behind scRNA-seq and spatial transcriptomic approaches, highlighting the recent discoveries and novel insights these methodologies have yielded in both liver physiology and disease biology. 10.1016/j.jhep.2020.06.004
Single-cell and spatial analyses of cancer cells: toward elucidating the molecular mechanisms of clonal evolution and drug resistance acquisition. Nagasawa Satoi,Kashima Yukie,Suzuki Ayako,Suzuki Yutaka Inflammation and regeneration Even within a single type of cancer, cells of various types exist and play interrelated roles. Each of the individual cells resides in a distinct microenvironment and behaves differently. Such heterogeneity is the most cumbersome nature of cancers, which is occasionally uncountable when effective prevention or total elimination of cancers is attempted. To understand the heterogeneous nature of each cell, the use of conventional methods for the analysis of "bulk" cells is insufficient. Although some methods are high-throughput and compressive regarding the genes being detected, the obtained data would be from the cell mass, and the average of a large number of the component cells would no longer be measured. Single-cell analysis, which has developed rapidly in recent years, is causing a drastic change. Genome, transcriptome, and epigenome analyses at single-cell resolution currently target cancer cells, cancer-associated fibroblasts, endothelial cells of vessels, and circulating and infiltrating immune cells. In fact, surprisingly diverse features of clonal evolution of cancer cells, during the development of cancer or acquisition of drug resistance, accompanied by corresponding gene expression changes in the circumstantial stromal cells, appeared in recent single-cell analyses. Based on the obtained novel insights, better optimal drug selection and new drug administration sequences were started. Even a remaining concern of the single cell analyses is being addressed. Until very recently, it was impossible to obtain positional information of cells in cancer via single-cell analysis because such information is lost during preparation of single-cell suspensions. A new method, collectively called spatial transcriptome (ST) analysis, has been developed and rapidly applied to various clinical specimens. In this review, we first outline the recent achievements of single-cell cancer analysis in analyzing the molecular basis underlying the acquisition of drug resistance, particularly focusing on the latest anti-epidermal growth factor receptor tyrosine kinase inhibitor, osimertinib. Further, we review the currently available ST analysis methods and introduce our recent attempts regarding the respective topics. 10.1186/s41232-021-00170-x
From bench to bedside: Single-cell analysis for cancer immunotherapy. Cancer cell Single-cell technologies are emerging as powerful tools for cancer research. These technologies characterize the molecular state of each cell within a tumor, enabling new exploration of tumor heterogeneity, microenvironment cell-type composition, and cell state transitions that affect therapeutic response, particularly in the context of immunotherapy. Analyzing clinical samples has great promise for precision medicine but is technically challenging. Successfully identifying predictors of response requires well-coordinated, multi-disciplinary teams to ensure adequate sample processing for high-quality data generation and computational analysis for data interpretation. Here, we review current approaches to sample processing and computational analysis regarding their application to translational cancer immunotherapy research. 10.1016/j.ccell.2021.07.004
Fourth Generation of Next-Generation Sequencing Technologies: Promise and Consequences. Ke Rongqin,Mignardi Marco,Hauling Thomas,Nilsson Mats Human mutation In this review, we discuss the emergence of the fourth-generation sequencing technologies that preserve the spatial coordinates of RNA and DNA sequences with up to subcellular resolution, thus enabling back mapping of sequencing reads to the original histological context. This information is used, for example, in two current large-scale projects that aim to unravel the function of the brain. Also in cancer research, fourth-generation sequencing has the potential to revolutionize the field. Cancer Research UK has named "Mapping the molecular and cellular tumor microenvironment in order to define new targets for therapy and prognosis" one of the grand challenges in tumor biology. We discuss the advantages of sequencing nucleic acids directly in fixed cells over traditional next-generation sequencing (NGS) methods, the limitations and challenges that these new methods have to face to become broadly applicable, and the impact that the information generated by the combination of in situ sequencing and NGS methods will have in research and diagnostics. 10.1002/humu.23051
Biological Misinterpretation of Transcriptional Signatures in Tumor Samples Can Unknowingly Undermine Mechanistic Understanding and Faithful Alignment with Preclinical Data. Clinical cancer research : an official journal of the American Association for Cancer Research PURPOSE:Precise mechanism-based gene expression signatures (GES) have been developed in appropriate in vitro and in vivo model systems, to identify important cancer-related signaling processes. However, some GESs originally developed to represent specific disease processes, primarily with an epithelial cell focus, are being applied to heterogeneous tumor samples where the expression of the genes in the signature may no longer be epithelial-specific. Therefore, unknowingly, even small changes in tumor stroma percentage can directly influence GESs, undermining the intended mechanistic signaling. EXPERIMENTAL DESIGN:Using colorectal cancer as an exemplar, we deployed numerous orthogonal profiling methodologies, including laser capture microdissection, flow cytometry, bulk and multiregional biopsy clinical samples, single-cell RNA sequencing and finally spatial transcriptomics, to perform a comprehensive assessment of the potential for the most widely used GESs to be influenced, or confounded, by stromal content in tumor tissue. To complement this work, we generated a freely-available resource, ConfoundR; https://confoundr.qub.ac.uk/, that enables users to test the extent of stromal influence on an unlimited number of the genes/signatures simultaneously across colorectal, breast, pancreatic, ovarian and prostate cancer datasets. RESULTS:Findings presented here demonstrate the clear potential for misinterpretation of the meaning of GESs, due to widespread stromal influences, which in-turn can undermine faithful alignment between clinical samples and preclinical data/models, particularly cell lines and organoids, or tumor models not fully recapitulating the stromal and immune microenvironment. CONCLUSIONS:Efforts to faithfully align preclinical models of disease using phenotypically-designed GESs must ensure that the signatures themselves remain representative of the same biology when applied to clinical samples. 10.1158/1078-0432.CCR-22-1102
RNA sequencing: the teenage years. Stark Rory,Grzelak Marta,Hadfield James Nature reviews. Genetics Over the past decade, RNA sequencing (RNA-seq) has become an indispensable tool for transcriptome-wide analysis of differential gene expression and differential splicing of mRNAs. However, as next-generation sequencing technologies have developed, so too has RNA-seq. Now, RNA-seq methods are available for studying many different aspects of RNA biology, including single-cell gene expression, translation (the translatome) and RNA structure (the structurome). Exciting new applications are being explored, such as spatial transcriptomics (spatialomics). Together with new long-read and direct RNA-seq technologies and better computational tools for data analysis, innovations in RNA-seq are contributing to a fuller understanding of RNA biology, from questions such as when and where transcription occurs to the folding and intermolecular interactions that govern RNA function. 10.1038/s41576-019-0150-2