Expression Programs Research Articles

Spatial Transcriptomics of IPMN Reveals Divergent Indolent and Malignant Lineages.

Intraductal papillary mucinous neoplasms (IPMN) occur in 5-10% of the population, but only a small minority progress to pancreatic ductal adenocarcinoma (PDAC). The lack of accurate predictors of high-risk disease leads both to unnecessary operations for indolent neoplasms as well as missed diagnoses of PDAC. Digital spatial RNA profiling (DSP-RNA) provides an opportunity to define and associate transcriptomic states with cancer risk. Whole-transcriptome DSP-RNA profiling was performed on 10 IPMN specimens encompassing the spectrum of dysplastic changes from normal duct to cancer. Ductal epithelial regions within each tissue were annotated as normal duct (NL), low-grade dysplasia (LGD), high-grade dysplasia (HGD), or invasive carcinoma (INV). Gene expression count data was generated by Illumina sequencing and analyzed with R/Bioconductor. Dimension reduction analysis exposed three clusters reflecting IPMN transcriptomic states denoted "normal-like" ( cNL ), "low-risk" ( cLR ) and "high-risk" ( cHR ). In addition to specific marker genes, the three states exhibited significant enrichment for the exocrine, classical, and basal-like programs in PDAC. Specifically, exocrine function diminished in cHR , classical activation distinguished neoplasia from cNL , and basal-like genes were specifically upregulated in cHR . Intriguingly, markers of cHR were detected in NL and LGD regions from specimens with PDAC but not low-grade IPMN. DSP-RNA of IPMN revealed low-risk (indolent) and high-risk (malignant) expression programs that correlated with the activity of exocrine and basal-like PDAC signatures, respectively, and distinguished pathologically low-grade from malignant specimens. These findings contextualize IPMN pathogenesis and have the potential to transform existing risk stratification models. Current consensus guidelines for management of intraductal papillary mucinous neoplasms (IPMN) of the pancreas utilize clinical and radiographic criteria for risk stratification. Unfortunately, the estimated positive predictive value of these criteria for IPMN-associated pancreatic ductal adenocarcinoma (PDAC) is under 50%, indicating that over half of pancreatectomies are performed for benign disease. Moreover, nearly 15% of patients who were deemed "low risk" by the same criteria harbored PDAC. Surgical resection of IPMN has maximal benefit when performed prior to the development of PDAC, as evidence of carcinoma has been associated with a high rate of recurrence and poor overall survival. Thus, the development of molecular diagnostics that improve the accuracy of IPMN risk classification would have immediate relevance for patient care, both in terms of better selecting patients for potentially curative operations, as well as sparing patients with low-risk lesions from invasive procedures.

Single-cell sequencing analysis reveals the dynamic tumour ecosystems of primary and metastatic lymph nodes in nasopharyngeal carcinoma.

Lymph node metastasis contributed to the leading cause and treatment failure in nasopharyngeal carcinoma (NPC). The microenvironment and the cellular communications of lymph node metastasized tumours determine the tumour progression and therapeutic effect, but the ecosystems about the lymph node metastasis (LNM) for NPC patients remain poorly characterized. Here, we integrated the transcriptomes of 47,618 single cells from eight samples related to NPC LNM. The dynamic immune ecosystems and immunosuppressive microenvironment including T cells, myeloid cells and B cells were observed in the lymph node metastatic samples compared with primary tumours. Additionally, the heterogeneity of epithelial cells was also revealed, and several clusters with expression programs that were associated with the progression-free survival of NPC patients were identified. Additionally, our data revealed the complex intercellular communications from primary to lymph node metastasis. The rewiring of CCL signalling which plays an important role in tumour metastasis was further identified. Altogether, we systematically characterized the ecosystem of NPC primary and lymph node metastasized tumours, which may shed light on the development of a therapeutic strategy to improve clinical outcomes of NPC patients with lymph node metastasis.

DsRNA formation leads to preferential nuclear export and gene expression

When mRNAs have been transcribed and processed in the nucleus, they are exported to the cytoplasm for translation. This export is mediated by the export receptor heterodimer Mex67–Mtr2 in the yeast Saccharomyces cerevisiae (TAP–p15 in humans)1,2. Interestingly, many long non-coding RNAs (lncRNAs) also leave the nucleus but it is currently unclear why they move to the cytoplasm3. Here we show that antisense RNAs (asRNAs) accelerate mRNA export by annealing with their sense counterparts through the helicase Dbp2. These double-stranded RNAs (dsRNAs) dominate export compared with single-stranded RNAs (ssRNAs) because they have a higher capacity and affinity for the export receptor Mex67. In this way, asRNAs boost gene expression, which is beneficial for cells. This is particularly important when the expression program changes. Consequently, the degradation of dsRNA, or the prevention of its formation, is toxic for cells. This mechanism illuminates the general cellular occurrence of asRNAs and explains their nuclear export.

The Effect of an Emotion Recognition and Expression Program on the Alexithymia, Emotion Expression Skills and Positive and Negative Symptoms of Patients with Schizophrenia in a Community Mental Health Center

This study aimed to examine the effect of an emotion recognition and expression program (EREP) on the alexithymia, emotion expression skills and positive and negative symptoms of patients with schizophrenia. The study had a non-randomized, quasi-experimental design including a pretest, post-test, and follow-up test. It was conducted with 36 patients with schizophrenia (n = 18 intervention group, n = 18 control group) who regularly visited a Community Mental Health Center (CMHC) in Türkiye and participated voluntarily. The EREP was applied to the intervention group for eight weeks. "Personal Information Form", "Emotion Expression Scale (EES)", "Toronto Alexithymia Scale (TAS)", and "Positive Negative Syndrome Scale (PANSS)" were applied to all participants in the pretest, post-test, and follow-up test. The follow-up test was applied one month after the end of the sessions. Number, percentage, chi-square test, and repeated measures analysis of variance were used for data evaluation. In the total alexithymia score, there was a significant difference in the group interaction by time in the intervention group compared to the control group. In terms of total alexithymia score, the post-test and follow-up test mean scores of the intervention group were lower than the control group (p < 0.05; η2 = 0.122). There was a significant time*group interaction in the positive emotion subscale of the EES (p < 0.05; η2 = 0.121). The findings of our study indicated that the EREP had a positive effect on the alexithymia scores of patients with schizophrenia. We found that the EREP used in our study contributed to the reduction of alexithymia levels in patients with schizophrenia.

Dynamic immune recovery process after liver transplantation revealed by single-cell multi-omics analysis

Elucidating the temporal process of immune remodeling under immunosuppressive treatment after liver transplantation (LT) is critical for precise clinical management strategies. Here, we performed a single-cell multi-omics analysis of peripheral blood mononuclear cells (PBMCs) collected from LT patients (with and without acute cellular rejection [ACR]) at 13 time points. Validation was performed in two independent cohorts with additional LT patients and healthy controls. Our study revealed a four-phase recovery process after LT and delineated changes in immune cell composition, expression programs, and interactions along this process. The intensity of the immune response differs between the ACR and non-ACR patients. Notably, the newly identified inflamed NK cells, CD14+RNASE2+ monocytes, and FOS-expressing monocytes emerged as predictive indicators of ACR. This study illuminates the longitudinal evolution of the immune cell landscape under tacrolimus-based immunosuppressive treatment during LT recovery, providing a four-phase framework that aids the clinical management of LT patients.

Transcriptomic Maps of Colorectal Liver Metastasis: Machine Learning of Gene Activation Patterns and Epigenetic Trajectories in Support of Precision Medicine.

The molecular mechanisms of the liver metastasis of colorectal cancer (CRLM) remain poorly understood. Here, we applied machine learning and bioinformatics trajectory inference to analyze a gene expression dataset of CRLM. We studied the co-regulation patterns at the gene level, the potential paths of tumor development, their functional context, and their prognostic relevance. Our analysis confirmed the subtyping of five liver metastasis subtypes (LMS). We provide gene-marker signatures for each LMS, and a comprehensive functional characterization that considers both the hallmarks of cancer and the tumor microenvironment. The ordering of CRLMs along a pseudotime-tree revealed a continuous shift in expression programs, suggesting a developmental relationship between the subtypes. Notably, trajectory inference and personalized analysis discovered a range of epigenetic states that shape and guide metastasis progression. By constructing prognostic maps that divided the expression landscape into regions associated with favorable and unfavorable prognoses, we derived a prognostic expression score. This was associated with critical processes such as epithelial-mesenchymal transition, treatment resistance, and immune evasion. These factors were associated with responses to neoadjuvant treatment and the formation of an immuno-suppressive, mesenchymal state. Our machine learning-based molecular profiling provides an in-depth characterization of CRLM heterogeneity with possible implications for treatment and personalized diagnostics.

Decoding and Recoding of mRNA Sequences by the Ribosome.

Faithful translation of messenger RNA (mRNA) into protein is essential to maintain protein homeostasis in the cell. Spontaneous translation errors are very rare due to stringent selection of cognate aminoacyl transfer RNAs (tRNAs) and the tight control of the mRNA reading frame by the ribosome. Recoding events, such as stop codon readthrough, frameshifting, and translational bypassing, reprogram the ribosome to make intentional mistakes and produce alternative proteins from the same mRNA. The hallmark of recoding is the change of ribosome dynamics. The signals for recoding are built into the mRNA, but their reading depends on the genetic makeup of the cell, resulting in cell-specific changes in expression programs. In this review, I discuss the mechanisms of canonical decoding and tRNA-mRNA translocation; describe alternative pathways leading to recoding; and identify the links among mRNA signals, ribosome dynamics, and recoding.

Temporal single-cell atlas of non-neuronal retinal cells reveals dynamic, coordinated multicellular responses to central nervous system injury.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2+ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Application of single-cell RNA sequencing in head and neck squamous cell carcinoma.

Single-cell RNA sequencing has been broadly applied to head and neck squamous cell carcinoma (HNSCC) for characterizing the heterogeneity and genomic mutations of HNSCC benefiting from the advantage of single-cell resolution. We summarized most of the current studies and aimed to explore their research methods and ideas, as well as how to transform them into clinical applications. Through single-cell RNA sequencing, we found the differences in tumor cells' expression programs and differentiation tracks. The studies of immune microenvironment allowed us to distinguish immune cell subpopulations, the extensive expression of immune checkpoints, and the complex crosstalk network between immune cells and non-immune cells. For cancer-associated fibroblasts (CAFs), single-cell RNA sequencing had made an irreplaceable contribution to the exploration of their differentiation status, specific CAFs markers, and the interaction with tumor cells and immune cells. In addition, we demonstrated in detail how single-cell RNA sequencing explored the HNSCC epithelial-to-mesenchymal transition (EMT) model and the mechanism of drug resistance, as well as its clinical value.

Microenvironmental sensing by fibroblasts controls macrophage population size

Animal tissues comprise diverse cell types. However, the mechanisms controlling the number of each cell type within tissue compartments remain poorly understood. Here, we report that different cell types utilize distinct strategies to control population numbers. Proliferation of fibroblasts, stromal cells important for tissue integrity, is limited by space availability. In contrast, proliferation of macrophages, innate immune cells involved in defense, repair, and homeostasis, is constrained by growth factor availability. Examination of density-dependent gene expression in fibroblasts revealed that Hippo and TGF-β target genes are both regulated by cell density. We found YAP1, the transcriptional coactivator of the Hippo signaling pathway, directly regulates expression of Csf1, the lineage-specific growth factor for macrophages, through an enhancer of Csf1 that is specifically active in fibroblasts. Activation of YAP1 in fibroblasts elevates Csf1 expression and is sufficient to increase the number of macrophages at steady state. Our data also suggest that expression programs in fibroblasts that change with density may result from sensing of mechanical force through actin-dependent mechanisms. Altogether, we demonstrate that two different modes of population control are connected and coordinated to regulate cell numbers of distinct cell types. Sensing of the tissue environment may serve as a general strategy to control tissue composition.

Comparative chromatin accessibility upon BDNF stimulation delineates neuronal regulatory elements

Neuronal stimulation induced by the brain‐derived neurotrophic factor (BDNF) triggers gene expression, which is crucial for neuronal survival, differentiation, synaptic plasticity, memory formation, and neurocognitive health. However, its role in chromatin regulation is unclear. Here, using temporal profiling of chromatin accessibility and transcription in mouse primary cortical neurons upon either BDNF stimulation or depolarization (KCl), we identify features that define BDNF‐specific chromatin‐to‐gene expression programs. Enhancer activation is an early event in the regulatory control of BDNF‐treated neurons, where the bZIP motif‐binding Fos protein pioneered chromatin opening and cooperated with co‐regulatory transcription factors (Homeobox, EGRs, and CTCF) to induce transcription. Deleting cis‐regulatory sequences affect BDNF‐mediated Arc expression, a regulator of synaptic plasticity. BDNF‐induced accessible regions are linked to preferential exon usage by neurodevelopmental disorder‐related genes and the heritability of neuronal complex traits, which were validated in human iPSC‐derived neurons. Thus, we provide a comprehensive view of BDNF‐mediated genome regulatory features using comparative genomic approaches to dissect mammalian neuronal stimulation.

Discovery of Unannotated Small Open Reading Frames in Streptococcus pneumoniae D39 Involved in Quorum Sensing and Virulence Using Ribosome Profiling.

ABSTRACTStreptococcus pneumoniae, an opportunistic human pathogen, is the leading cause of community-acquired pneumonia and an agent of otitis media, septicemia, and meningitis. Although genomic and transcriptomic studies of S. pneumoniae have provided detailed perspectives on gene content and expression programs, they have lacked information pertaining to the translational landscape, particularly at a resolution that identifies commonly overlooked small open reading frames (sORFs), whose importance is increasingly realized in metabolism, regulation, and virulence. To identify protein-coding sORFs in S. pneumoniae, antibiotic-enhanced ribosome profiling was conducted. Using translation inhibitors, 114 novel sORFs were detected, and the expression of a subset of them was experimentally validated. Two loci associated with virulence and quorum sensing were examined in deeper detail. One such sORF, rio3, overlaps with the noncoding RNA srf-02 that was previously implicated in pathogenesis. Targeted mutagenesis parsing rio3 from srf-02 revealed that rio3 is responsible for the fitness defect seen in a murine nasopharyngeal colonization model. Additionally, two novel sORFs located adjacent to the quorum sensing receptor rgg1518 were found to impact regulatory activity. Our findings emphasize the importance of sORFs present in the genomes of pathogenic bacteria and underscore the utility of ribosome profiling for identifying the bacterial translatome.

Abstract 1297: Chemotherapy-induced senescence activates robust, parallel programs of immune checkpoint expression that can be targeted with immunotherapy

Abstract We and others have previously shown that the breast cancers most difficult to eradicate with chemotherapy are TP53 wild-type, and these are among the most lethal breast cancers. For instance, chemotherapy-treated TNBC patients with TP53 wild-type tumors have a median overall survival of 45 months, vs 263 months for TP53 mutant tumors. Treatments fail to eradicate these cancers for two reasons: the tumor cells 1) avoid intrinsic cell death, and 2) escape immune clearance. We have shown that p53 drives a program of senescence that arrests the cell cycle to prevent intrinsic modes of cell death such as mitotic catastrophe, apoptosis and nutrient deprivation. However, it is not currently understood how breast cancer cells that enter senescence to survive chemotherapy escape immune clearance. Here, we show tumor cells from mice and human patients that survive chemotherapy activate complex programs of immune modulation. Surprisingly, the surviving, senescent tumor cells are highly enriched for many antigen presentation genes and positive regulators of T cell activation, suggesting they have switched from immunologically “cold” to “hot”. Unfortunately, these senescent tumors concurrently upregulate redundant expression of many T cell inhibitory checkpoint genes, including CD80 and PD-L1. scRNA-seq and cell imaging revealed that CD80 and PD-L1 each mark unique populations of cells in the treated tumor, typified by p53 signaling or interferon signaling, respectively. In p53 wild-type, syngeneic, orthotopic mouse mammary tumor models that recapitulate human breast cancer response to chemotherapy, treatment of tumors with chemotherapy followed by targeting of the PD-L1 and/or CD80 axes improved response, including complete eradication in some instances. Unfortunately, however, even combination strategies failed to elicit a cure in the majority of cases. Our findings reveal the formidable challenge of eliminating residual disease populated by senescent cells that express multiple redundant immune inhibitory pathways and suggest rational strategies are needed based on the specific checkpoint pathways expressed in residual disease. Citation Format: Ashkan Shahbandi, Fang-Yen Chiu, Nathan A. Ungerleider, Ashlyn Y. Anderson, Heather L. Machado, Zachary F. Pursell, Raegan Kvadas, Sonia G. Rao, James G. Jackson. Chemotherapy-induced senescence activates robust, parallel programs of immune checkpoint expression that can be targeted with immunotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1297.

Abstract 1380: Setdb1-loss reactivates ERV expression and interferon signaling to induce immune-mediated melanoma clearance

Abstract While immunotherapies have been revolutionary for the treatment of many cancers including melanoma, most patients do not respond, or develop resistance to these therapies. Several epigenetic modulators have been identified as potential targets for augmenting anti-tumor immunity, however the full breadth of their role is not well-understood. To identify modulators of anti-tumor immunity in melanoma, we performed a whole-genome CRISPR screen using the YUMMER-G model and identified members of the HUSH complex and its associated methyltransferase, Setdb1, as drop-out hits. To validate these findings, we used CRISPR-Cas9 to knockout Setdb1 individually and observed that Setdb1-/- tumors are cleared by their host mice in a CD8+ T-cell dependent manner. To determine the source of immunogenicity in Setdb1-/- cells, we performed RNA-sequencing, cytokine profiling, and evaluated antigen presentation by flow cytometry. These assays revealed increased MHC-I levels, elevated secretion of Ifn-, and induction of an interferon-stimulated gene (ISG) expression program, consistent with the observed immunogenicity. We further determined that treating Setdb1-/- cells with IFNAR-blocking antibodies restores MHC-I levels to wild-type levels, suggestive of a tumor-intrinsic interferon loop. We hypothesized that type-I interferon signaling is generated through sensing of expressed endogenous retroviral elements (ERVs), and consistently observe upregulation of ERVs in Setdb1-/- cells. Previous studies have described the potential for ERV-derived antigens to act as substrates for anti-tumor activity by CD8+ T cells. We identified T cells in the tumor microenvironment of Setdb1-/- tumors that bind tetramer carrying the P15E peptide derived from ERV envelope. Together, these findings suggest that tumors lacking expression of Setdb1-/- upregulate ERV transcripts, driving expression of type-I interferons which generate an inflamed tumor microenvironment and potentiate efficient tumor clearance of Setdb1-/- tumors. Overall, this study emphasizes a critical role for regulators of ERV expression and type-I interferon expression, and suggests their potential as therapeutic targets for augmenting anti-cancer immune responses. Citation Format: Meaghan K. McGeary, William Damsky, Drew Daniels, Goran Micevic, Eric Song, Hua Jane Lou, Clotilde Calderwood, Sateja Paradkar, Akiko Iwasaki, David Calderwood, Benjamin Turk, Marcus Bosenberg. Setdb1-loss reactivates ERV expression and interferon signaling to induce immune-mediated melanoma clearance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1380.

Abstract 6304: Targeting HIF2a with siRNA: From preclinical models to the clinic

Abstract Hypoxia-inducible factor 2 alpha (HIF2a) is arguably the most important driver of kidney cancer. HIF2a is constitutively activated following von Hippel-Lindau (VHL) gene inactivation, the signature event of the most common type of kidney cancer, clear cell renal cell carcinoma (ccRCC). HIF2a functions as a heterodimeric transcription factor and regulates a program of gene expression that promotes cell proliferation, stemness, and angiogenesis. Using a highly specific inhibitor designed to target a structural vulnerability in HIF2a (PT2399), we previously showed that approximately 50% of ccRCCs are dependent on HIF2a. However, prolonged drug exposure results in resistance and the acquisition of gatekeeper mutations, which we reported first in patient-derived xenografts (PDXs) and subsequently in humans. Using the same PDX platform that previously validated PT2399, we show that HIF2a can be effectively inhibited using a tumor-directed siRNA (siHIF2). Referring herein to both first- and second-generation (ARO-HIF2) siRNA drugs, siHIF2 is specifically taken up by human ccRCC tumors transplanted in mice, where it depletes HIF2a inhibiting target gene expression and tumor growth. Through orthogonal RNA-seq studies integrating both PT2399 and siHIF2 in PDXs, we provide unprecedented detail on the HIF2a effector transcriptome, which we further dissect by incorporating ChIP-seq. A PDX line was generated from a ccRCC patient who had paraneoplastic polycythemia (a HIF2a dependent syndrome due to erythropoietin [Epo] secretion by the tumor) and participated in the phase I trial of ARO-HIF2 (NCT04169711). We show that siHIF2 effectively depleted HIF2a in both the PDX as well as in the patient, that it normalized Epo (and hemoglobin), and that it inhibited tumor growth. siHIF2 has activity against both wild-type and drug-resistant mutant HIF2a and is expected to be active in patients progressing on PT2977 (belzutifan), a PT2399-related drug recently approved by the FDA. To our knowledge, this is the first example of functional inactivation of an oncoprotein with a tumor-directed siRNA in humans. In summary, these studies provide unique insight into HIF2a (the only known core dependency in ccRCC), illustrate how it can be effectively inhibited by an siRNA drug, and establish a paradigm for the development of tumor directed siRNA-based therapeutics. Citation Format: Yuanqing Ma, Christina Stevens, Olivia Brandenburg, Vanina Toffessi Tcheuyap, Quyen N. Do, Faeze Saatchi, Tanner Hardy, Oluwatomilade Fatunde, Alyssa Macchiaroli, Jeffrey Miyata, Deyssy Carrillo, Thomas Schluep, So Wong, Alana Christie, Payal Kapur, Ivan Pedrosa, James Hamilton, James Brugarolas. Targeting HIF2a with siRNA: From preclinical models to the clinic [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 6304.

The biological role of primary and secondary plants metabolites

Metabolism consists of closely coordinated series of enzyme-mediated chemical reactions that take place in the plant organism, resulting in the synthesis and use of a wide variety of molecules in the category of carbohydrates, amino acids, fatty acids, nucleotides and polymers derived from them (polysaccharides, proteins, lipids, DNA, RNA, etc.). All these processes are defined as primary metabolism and the respective compounds, which are essential for the survival of the plant, are described as primary metabolites. In addition to the primary metabolites, which play a role in maintaining the viability of the plant (proteins, carbohydrates and lipids), a number of compounds such as terpenes, steroids, anthocyanins, anthraquinones, phenols and polyphenols, which belong to the "secondary metabolism", are also synthesized. Secondary metabolites (SMs) are present only in certain species, often manifesting specificity of organ or tissue, can be identified only at a certain stage of growth and development within a species, or can be activated only during periods of stress caused by the attack. microorganisms or nutrient depletion. Their synthesis seems to have no direct significance for the synthesizing cell, but may be decisive for the development and functioning of the body as a whole. Their synthesis is not a vital part of the gene expression and developmental program, these metabolites are not simple catabolic products, have a diversified structure and can be frequently re-included in metabolic processes. The boundary between primary and secondary metabolism is uncertain, as many primary metabolism intermediates play similar roles in secondary metabolism. Some obscure amino acids are infallibly SMs, while sterols are essential structural compounds of many organisms and should therefore be considered primary metabolites.

Abstract PD6-03: Spatio-molecular dissection of the breast cancer metastatic microenvironment

Abstract Metastatic breast cancer (MBC) remains incurable due to inevitable development of therapeutic resistance. Although tumor cell intrinsic mechanisms of resistance in MBC are beginning to be elucidated by bulk sequencing studies, the roles of the tumor microenvironment and intratumor heterogeneity in therapeutic resistance remain underexplored due to both technological barriers and limited availability of samples. To comprehensively capture these characteristics we have adapted a research biopsy protocol to collect tissue for an array of single-cell and spatio-molecular assays whose performance we have optimized for MBC, including single-cell and single-nucleus RNA sequencing, Slide-Seq, Multiplexed Error-Robust FISH (MERFISH), Expansion Sequencing (ExSEQ), Co-detection by Indexing (CODEX) and Multiplexed Ion Beam Imaging (MIBI). To date, we have successfully performed single-cell or single-nucleus RNAseq in 67 MBC biopsies and generated detailed accompanying clinical annotations for each. These samples provide a representation of the clinicopathological diversity of MBC including different breast cancer subtypes (44 HR+/HER2-, 3 HR-/HER2+, 3 HR+/HER2+, 16 TNBC, 1 unknown), common anatomic sites of metastasis (37 liver, 9 axilla, 7 breast, 5 bone, 3 chest wall, 3 neck, 1 brain, 1 lung, 1 skin), metastatic presentations (53 recurrent, 14 de novo) and histologic subtypes in the breast (45 IDC, 7 ILC, 6 mixed, 3 DCIS, 1 mucinous, 5 unknown/NA). Following optimization, both single-cell and single-nucleus RNA seq perform well in these MBC biopsies recovering all expected cell types including the malignant, stromal (e.g. fibroblasts, endothelial cells), myeloid (e.g. monocytes, macrophages) and lymphoid compartments (e.g. T cells, B cells, NK cells) as well as relevant oncogenic programs (e.g. cell cycle programs in all compartments; EMT-like and ER signaling programs in the malignant compartment, immune checkpoint programs in the lymphoid compartment; and fibroblast activation and vascular homeostasis programs in the stromal compartment). In addition to differences between the two techniques, these data demonstrate substantial intratumor heterogeneity in cell type composition. For example in liver biopsies the average number of cells per sample compartment by single nucleus RNA-seq was 6745 malignant (56%, SD 4216), 4637 stromal (41%, SD 3727), 1196 lymphoid (8%, SD 1617) and 874 myeloid (6%, SD 852); in breast biopsies the average number of cells per compartment by single nucleus RNA-seq was 6421 malignant (70%, SD 3497), 1628 stromal (24%, SD 117), 333 lymphoid (4%, SD 170) and 213 myeloid (3%, SD 117). Additionally, we find both inter- and intra-tumor heterogeneity in expression patterns and programs including, for example, expression of ER, PR and HER2 within clinical receptor subtypes (log normalized counts for ER expression in tumor cells by single cell RNA-seq: HR+/HER2- 0.921 (SD 0.714); HR+/HER2+ 0.768 (SD 0.624); HR-/HER2+ 0.018 (SD 0.122); and HR-/HER2- 0.005 (SD 0.066). For a subset of 13 biopsies we are also completing the spatiomolecular characterization methods on serial sections of a single adjacent biopsy. This unique experimental setup was designed to enable efficient comparison and integration of these assays. In spite of differences between experimental techniques and readouts, cell typing can be approached by annotation transfer from matching single cell or single nucleus RNAseq data, enabling exploratory analyses including evaluation of spatial phenotypes and cell type colocalization. Overall, these single cell and spatial data afford a comprehensive atlas including cell types, cell states/programs, cell interactions and spatial organization in MBC lesions. Future analyses will include serial biopsies over time and integration of clinicopathologic data including therapeutic response and resistance. Citation Format: Daniel L Abravanel, Johanna Klughammer, Timothy Blosser, Yury Goltsev, Sizun Jiang, Yunjao Bai, Evan Murray, Shahar Alon, Yi Cui, Daniel R Goodwin, Anubhav Sinha, Ofir Cohen, Michal Slyper, Orr Ashenberg, Danielle Dionne, Judit Jané-Valbuena, Caroline BM Porter, Asa Segerstolpe, Julia Waldman, Sébastien Vigneau, Karla Helvie, Allison Frangieh, Laura DelloStritto, Miraj Patel, Jingyi We, Kathleen Pfaff, Nicole Cullen, Ana Lako, Madison Turner, Isaac Wakiro, Sara Napolitano, Abhay Kanodia, Rebecca Ortiz, Colin MacKichan, Stephanie Inga, Judy Chen, Aaron R Thorner, Asaf Rotem, Scott Rodig, Fei Chen, Edward S Boyden, Garry P Nolan, Xiaowei Zhuang, Orit Rozenblatt-Rosen, Bruce E Johnson, Aviv Regev, Nikhil Wagle. Spatio-molecular dissection of the breast cancer metastatic microenvironment [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr PD6-03.

Posttranscriptional regulation of cellulose synthase genes by small RNAs derived from cellulose synthase antisense transcripts.

Transcriptional regulatory mechanisms governing plant cell wall biosynthesis are incomplete. Expression programs that activate wall biosynthesis are well understood, but mechanisms that control the attenuation of gene expression networks remain elusive. Previous work has shown that small RNAs (sRNAs) derived from the HvCESA6 (Hordeum vulgare, Hv) antisense transcripts are naturally produced and are capable of regulating aspects of wall biosynthesis. Here, we further test the hypothesis that CESA‐derived sRNAs generated from CESA antisense transcripts are involved in the regulation of cellulose and broader cell wall biosynthesis. Antisense transcripts were detected for some but not all members of the CESA gene family in both barley and Brachypodium distachyon. Phylogenetic analysis indicates that antisense transcripts are detected for most primary cell wall CESA genes, suggesting a possible role in the transition from primary to secondary cell wall biosynthesis. Focusing on one antisense transcript, HvCESA1 shows dynamic expression throughout development, is correlated with corresponding sRNAs over the same period and is anticorrelated with HvCESA1 mRNA expression. To assess the broader impacts of CESA‐derived sRNAs on the regulation of cell wall biosynthesis, transcript profiling was performed on barley tissues overexpressing CESA‐derived sRNAs. Together, the data support the hypothesis that CESA antisense transcripts function through an RNA‐induced silencing mechanism, to degrade cis transcripts, and may also trigger trans‐acting silencing on related genes to alter the expression of cell wall gene networks.

Comparative transcriptome analysis reveals distinct gene expression profiles in Brachypodium distachyon infected by two fungal pathogens

BackgroundThe production of cereal crops is frequently affected by diseases caused by Fusarium graminearum and Magnaporthe oryzae, two devastating fungal pathogens. To improve crop resistance, many studies have focused on understanding the mechanisms of host defense against these two fungi individually. However, our knowledge of the common and different host defenses against these pathogens is very limited.ResultsIn this study, we employed Brachypodium distachyon as a model for cereal crops and performed comparative transcriptomics to study the dynamics of host gene expression at different infection stages. We found that infection with either F. graminearum or M. oryzae triggered massive transcriptomic reprogramming in the diseased tissues. Numerous defense-related genes were induced with dynamic changes during the time course of infection, including genes that function in pattern detection, MAPK cascade, phytohormone signaling, transcription, protein degradation, and secondary metabolism. In particular, the expression of jasmonic acid signaling genes and proteasome component genes were likely specifically inhibited or manipulated upon infection by F. graminearum.ConclusionsOur analysis showed that, although the affected host pathways are similar, their expression programs and regulations are distinct during infection by F. graminearum and M. oryzae. The results provide valuable insight into the interactions between B. distachyon and two important cereal pathogens.

A cellular and spatial map of the choroid plexus across brain ventricles and ages

The choroid plexus (ChP) in each brain ventricle produces cerebrospinal fluid (CSF) and forms the blood-CSF barrier. Here, we construct a single-cell and spatial atlas of each ChP in the developing, adult, and aged mouse brain. We delineate diverse cell types, subtypes, cell states, and expression programs in epithelial and mesenchymal cells across ages and ventricles. In the developing ChP, we predict a common progenitor pool for epithelial and neuronal cells, validated by lineage tracing. Epithelial and fibroblast cells show regionalized expression by ventricle, starting at embryonic stages and persisting with age, with a dramatic transcriptional shift with maturation, and a smaller shift in each aged cell type. With aging, epithelial cells upregulate host-defense programs, and resident macrophages upregulate interleukin-1β (IL-1β) signaling genes. Our atlas reveals cellular diversity, architecture and signaling across ventricles during development, maturation, and aging of the ChP-brain barrier.