Open-ST: High-resolution spatial transcriptomics in 3D.

Spatial transcriptomics (ST) methods unlock molecular mechanisms underlying tissue development, homeostasis, or disease. However, there is a need for easy-to-use, high-resolution, cost-efficient, and 3D-scalable methods. Here, we report Open-ST, a sequencing-based, open-source experimental and computational resource to address these challenges and to study the molecular organization of tissues in 2D and 3D. In mouse brain, Open-ST captured transcripts at subcellular resolution and reconstructed cell types. In primary head-and-neck tumors and patient-matched healthy/metastatic lymph nodes, Open-ST captured the diversity of immune, stromal, and tumor populations in space, validated by imaging-based ST. Distinct cell states were organized around cell-cell communication hotspots in the tumor but not the metastasis. Strikingly, the 3D reconstruction and multimodal analysis of the metastatic lymph node revealed spatially contiguous structures not visible in 2D and potential biomarkers precisely at the 3D tumor/lymph node boundary. All protocols and software are available at https://rajewsky-lab.github.io/openst.

TP53 regulates miRNA association with AGO2 to remodel the miRNA-mRNA interaction network.

DNA damage activates TP53-regulated surveillance mechanisms that are crucial in suppressing tumorigenesis. TP53 orchestrates these responses directly by transcriptionally modulating genes, including microRNAs (miRNAs), and by regulating miRNA biogenesis through interacting with the DROSHA complex. However, whether the association between miRNAs and AGO2 is regulated following DNA damage is not yet known. Here, we show that, following DNA damage, TP53 interacts with AGO2 to induce or reduce AGO2's association of a subset of miRNAs, including multiple let-7 family members. Furthermore, we show that specific mutations in TP53 decrease rather than increase the association of let-7 family miRNAs, reducing their activity without preventing TP53 from interacting with AGO2. This is consistent with the oncogenic properties of these mutants. Using AGO2 RIP-seq and PAR-CLIP-seq, we show that the DNA damage-induced increase in binding of let-7 family members to the RISC complex is functional. We unambiguously determine the global miRNA-mRNA interaction networks involved in the DNA damage response, validating them through the identification of miRNA-target chimeras formed by endogenous ligation reactions. We find that the target complementary region of the let-7 seed tends to have highly fixed positions and more variable ones. Additionally, we observe that miRNAs, whose cellular abundance or differential association with AGO2 is regulated by TP53, are involved in an intricate network of regulatory feedback and feedforward circuits. TP53-mediated regulation of AGO2-miRNA interaction represents a new mechanism of miRNA regulation in carcinogenesis.

ARGONAUTE2 cooperates with SWI/SNF complex to determine nucleosome occupancy at human Transcription Start Sites.

Argonaute (AGO) proteins have a well-established role in post-transcriptional regulation of gene expression as key component of the RNA silencing pathways. Recent evidence involves AGO proteins in mammalian nuclear processes such as transcription and splicing, though the mechanistic aspects of AGO nuclear functions remain largely elusive. Here, by SILAC-based interaction proteomics, we identify the chromatin-remodelling complex SWI/SNF as a novel AGO2 interactor in human cells. Moreover, we show that nuclear AGO2 is loaded with a novel class of Dicer-dependent short RNAs (sRNAs), that we called swiRNAs, which map nearby the Transcription Start Sites (TSSs) bound by SWI/SNF. The knock-down of AGO2 decreases nucleosome occupancy at the first nucleosome located downstream of TSSs in a swiRNA-dependent manner. Our findings indicate that in human cells AGO2 binds SWI/SNF and a novel class of sRNAs to establish nucleosome occupancy on target TSSs.

The role of TP53 in miRNA loading onto AGO2 and in remodelling the miRNA-mRNA interaction network.

BACKGROUND: DNA damage transactivates tumour protein p53 (TP53)-regulated surveillance, crucial in suppressing tumorigenesis. TP53 mediates this process directly by transcriptionally modulating gene and microRNA (miRNA) expression and indirectly by regulating miRNA biogenesis. However, the role of TP53 in regulating miRNA-AGO2 loading and global changes in AGO2 binding to its gene targets in response to DNA damage are unknown. These processes might be novel mechanisms by which TP53 regulates miRNAs in response to DNA damage. METHODS: To show the network of miRNA-mRNA interactions that occur in response to DNA damage, we stimulated TP53 wild-type and null cell-lines with doxorubicin and performed RNA sequencing from total RNA (RNA-Seq) and AGO2-immunoprecipitated RNA (AGO2-RIP-Seq). We used a combined AGO2 RIP-seq and AGO2 PAR-CLIP-seq (photoactivatable-ribonucleoside-enhanced cross-linking and immunoprecipitation) approach to determine the exact sites of interaction between the AGO2-bound miRNAs and their mRNA targets. FINDINGS: TP53 directly associated with AGO2, and induced and reduced loading of a subset of miRNAs, including the lethal 7 (let-7) miRNA family members, onto AGO2 in response to DNA damage. Although mutated TP53 maintained its capacity to interact with AGO2, it mediated unloading instead of loading of let-7 family miRNAs, thereby reducing their activity. We determined the miRNA-mRNA interaction networks involved in the response to DNA damage both in the presence and absence of TP53. Furthermore, we showed that miRNAs whose cellular abundance or differential loading onto AGO2 was regulated by TP53 were involved in an intricate network of regulatory feedback and feedforward circuits that fine-tuned gene expression levels in response to DNA damage to permit the repair of DNA damage or initiation of programmed cell death. INTERPRETATION: Control of AGO2 loading by TP53 is a new mechanism of miRNA regulation in carcinogenesis. FUNDING: UK Medical Research Council, Action Against Cancer.

Increased chronic lymphocytic leukemia proliferation upon IgM stimulation is sustained by the upregulation of miR-132 and miR-212.

To assess the involvement of microRNAs (miRNAs) in B-cell receptor (BCR) stimulation, we first evaluated miRNA profiling following IgM cross-linking in chronic lymphocytic leukemia (CLL) cells and in normal B lymphocytes. Second, we combined miRNA and gene expression data to identify putative miRNA functional networks. miRNA profiling showed distinctive patterns of regulation after stimulation in leukemic versus normal B lymphocytes and identified a differential responsiveness to BCR engagement in CLL subgroups according to the immunoglobulin heavy chain variable region mutational status and clinical outcome. The most significantly modulated miRNAs in stimulated CLL are miR-132 and miR-212. Notably, these miRNAs appeared regulated in progressive but not in stable CLL. Accordingly, gene profiling showed a significant transcriptional response to stimulation exclusively in progressive CLL. Based on these findings, we combined miRNA and gene expression data to investigate miR-132 and miR-212 candidate interactions in this CLL subgroup. Correlation analysis pointed to a link between these miRNAs and RB/E2F and TP53 cascades with proproliferative effects, as corroborated by functional analyses. Finally, basal levels of miR-132 and miR-212 were measured in an independent cohort of 20 unstimulated CLL cases and both showed lower expression in progressive compared to stable patients, suggesting an association between the expression of these molecules and disease prognosis. Overall, our results support a model involving miR-132 and miR-212 upregulation in sustaining disease progression in CLL. These miRNAs may therefore provide new valuable strategies for therapeutic intervention.

An emerging player in the adaptive immune response: microRNA-146a is a modulator of IL-2 expression and activation-induced cell death in T lymphocytes.

Activation of the T cell-mediated immune response has been associated with changes in the expression of specific microRNAs (miRNAs). However, the role of miRNAs in the development of an effective immune response is just beginning to be explored. This study focuses on the functional role of miR-146a in T lymphocyte-mediated immune response and provides interesting clues on the transcriptional regulation of miR-146a during T-cell activation. We show that miR-146a is low in human naive T cells and is abundantly expressed in human memory T cells; consistently, miR-146a is induced in human primary T lymphocytes upon T-cell receptor (TCR) stimulation. Moreover, we identified NF-kB and c-ETS binding sites as required for the induction of miR-146a transcription upon TCR engagement. Our results demonstrate that several signaling pathways, other than inflammation, are influenced by miR-146a. In particular, we provide experimental evidence that miR-146a modulates activation-induced cell death (AICD), acting as an antiapoptotic factor, and that Fas-associated death domain (FADD) is a target of miR-146a. Furthermore, miR-146a enforced expression impairs both activator protein 1 (AP-1) activity and interleukin-2 (IL-2) production induced by TCR engagement, thus suggesting a role of this miRNA in the modulation of adaptive immunity.

Epigenomic landscape of human colorectal cancer unveils an aberrant core of pan-cancer enhancers orchestrated by YAP/TAZ.

Cancer is characterized by pervasive epigenetic alterations with enhancer dysfunction orchestrating the aberrant cancer transcriptional programs and transcriptional dependencies. Here, we epigenetically characterize human colorectal cancer (CRC) using de novo chromatin state discovery on a library of different patient-derived organoids. By exploring this resource, we unveil a tumor-specific deregulated enhancerome that is cancer cell-intrinsic and independent of interpatient heterogeneity. We show that the transcriptional coactivators YAP/TAZ act as key regulators of the conserved CRC gained enhancers. The same YAP/TAZ-bound enhancers display active chromatin profiles across diverse human tumors, highlighting a pan-cancer epigenetic rewiring which at single-cell level distinguishes malignant from normal cell populations. YAP/TAZ inhibition in established tumor organoids causes extensive cell death unveiling their essential role in tumor maintenance. This work indicates a common layer of YAP/TAZ-fueled enhancer reprogramming that is key for the cancer cell state and can be exploited for the development of improved therapeutic avenues.

Characterization of B- and T-lineage acute lymphoblastic leukemia by integrated analysis of MicroRNA and mRNA expression profiles.

Acute lymphoblastic leukemia (ALL) is an heterogeneous disease comprising several subentities that differ for both immunophenotypic and molecular characteristics. Over the years, the biological understanding of this neoplasm has largely increased. Gene expression profiling has allowed to identify specific signatures for the different ALL subsets and permitted the identification of pathways deregulated by a given lesion. MicroRNAs (miRNAs) are small noncoding RNAs, which play a pivotal role in several cellular functions. In this study, we investigated miRNAs expression profiles in a series of adult ALL cases by microarray analysis. Unsupervised hierarchical clustering largely recapitulated ALL subgroups. Furthermore, we identified miR-148, miR-151, and miR-424 as discriminative of T-lineage versus B-lineage ALL; ANOVA highlighted a set of six miRNAs-namely miR-425-5p, miR-191, miR-146b, miR-128, miR-629, and miR-126-that can discriminate B-lineage ALL subgroups harboring specific molecular lesions. These results were confirmed and extended by quantitative-PCR on a further cohort of cases. Finally, we used Pearson correlation analysis to combine miRNA and gene expression profiles. The distribution of correlation coefficients generated by comparing the expression of every miRNA/gene pair in our data set shows enrichment of both positively and negatively correlated pairs over background distributions obtained using randomized data. Moreover, a clear enrichment for predicted miRNA:target pairs is observed at negative correlation coefficient intervals. Signal-to-noise ratio highlighted several miRNA/gene pairs with a possible role in the disease. In fact, gene set enrichment analysis of genes composing the selected miRNA/gene pairs displays over-representation of functional categories related to cancer and cell-cycle regulation.

Quelling: post-transcriptional gene silencing guided by small RNAs in Neurospora crassa.

The filamentous fungus Neurospora crassa is a model organism for the study of gene silencing. The most characterized gene silencing mechanism in this ascomycete is quelling, which occurs at the post-transcriptional level. Quelling is triggered by the introduction of transgenes and results in silencing of both transgenes and cognate endogenous mRNAs. Quelling is related to co-suppression, observed in plants, and RNA interference in animals; it requires an Argonaute protein and acts by generating small RNA molecules (about 25 nt long), which in turn target mRNAs to be silenced. It has been recently shown that quelling is needed for the taming of transposons but, unlike other model organisms, does not seem to play any role in heterochromatin assembly and maintenance.

Human Cortical Organoids Expose a Differential Function of GSK3 on Cortical Neurogenesis.

The regulation of the proliferation and polarity of neural progenitors is crucial for the development of the brain cortex. Animal studies have implicated glycogen synthase kinase 3 (GSK3) as a pivotal regulator of both proliferation and polarity, yet the functional relevance of its signaling for the unique features of human corticogenesis remains to be elucidated. We harnessed human cortical brain organoids to probe the longitudinal impact of GSK3 inhibition through multiple developmental stages. Chronic GSK3 inhibition increased the proliferation of neural progenitors and caused massive derangement of cortical tissue architecture. Single-cell transcriptome profiling revealed a direct impact on early neurogenesis and uncovered a selective role of GSK3 in the regulation of glutamatergic lineages and outer radial glia output. Our dissection of the GSK3-dependent transcriptional network in human corticogenesis underscores the robustness of the programs determining neuronal identity independent of tissue architecture.

Small interfering RNAs that trigger posttranscriptional gene silencing are not required for the histone H3 Lys9 methylation necessary for transgenic tandem repeat stabilization in Neurospora crassa.

In Neurospora crassa, the introduction of a transgene can lead to small interfering RNA (siRNA)-mediated posttranscriptional gene silencing (PTGS) of homologous genes. siRNAs can also guide locus-specific methylation of Lys9 of histone H3 (Lys9H3) in Schizosaccharomyces pombe. Here we tested the hypothesis that transgenically derived siRNAs may contemporaneously both activate the PTGS mechanism and induce chromatin modifications at the transgene and the homologous endogenous gene. We carried out chromatin immunoprecipitation using a previously characterized albino-1 (al-1) silenced strain but detected no alterations in the pattern of histone modifications at the endogenous al-1 locus, suggesting that siRNAs produced from the transgenic locus do not trigger modifications in trans of those histones tested. Instead, we found that the transgenic locus was hypermethylated at Lys9H3 in our silenced strain and remained hypermethylated in the quelling defective mutants (qde), further demonstrating that the PTGS machinery is dispensable for Lys9H3 methylation. However, we found that a mutant in the histone Lys9H3 methyltransferase dim-5 was unable to maintain PTGS, with transgenic copies being rapidly lost, resulting in reversion of the silenced phenotype. These results indicate that the defect in PTGS of the Deltadim-5 strain is due to the inability to maintain the transgene in tandem, suggesting a role for DIM-5 in stabilizing such repeated sequences. We conclude that in Neurospora, siRNAs produced from the transgenic locus are used in the RNA-induced silencing complex-mediated PTGS pathway and do not communicate with an RNAi-induced initiation of transcriptional gene silencing complex to effect chromatin-based silencing.

Characterization of Transcription Termination-Associated RNAs: New Insights into their Biogenesis, Tailing, and Expression in Primary Tumors.

Next-generation sequencing has uncovered novel classes of small RNAs (sRNAs) in eukaryotes, in addition to the well-known miRNAs, siRNAs, and piRNAs. In particular, sRNA species arise from transcription start sites (TSSs) and the transcription termination sites (TTSs) of genes. However, a detailed characterization of these new classes of sRNAs is still lacking. Here, we present a comprehensive study of sRNAs derived from TTSs of expressed genes (TTSa-RNAs) in human cell lines and primary tissues. Taking advantage of sRNA-sequencing, we show that TTSa-RNAs are present in the nuclei of human cells, are loaded onto both AGO1 and AGO2, and their biogenesis does not require DICER and AGO2 endonucleolytic activity. TTSa-RNAs display a strong bias against a G residue in the first position at 5' end, a known feature of AGO-bound sRNAs, and a peculiar oligoA tail at 3' end. AGO-bound TTSa-RNAs derive from genes involved in cell cycle progression regulation and DNA integrity checkpoints. Finally, we provide evidence that TTSa-RNAs can be detected by sRNA-Seq in primary human tissue, and their expression increases in tumor samples as compared to nontumor tissues, suggesting that in the future, TTSa-RNAs might be explored as biomarker for diagnosis or prognosis of human malignancies.

Redundancy of the two dicer genes in transgene-induced posttranscriptional gene silencing in Neurospora crassa.

RNA interference (RNAi) in animals, cosuppression in plants, and quelling in fungi are homology-dependent gene silencing mechanisms in which the introduction of either double-stranded RNA (dsRNA) or transgenes induces sequence-specific mRNA degradation. These phenomena share a common genetic and mechanistic basis. The accumulation of short interfering RNA (siRNA) molecules that guide sequence-specific mRNA degradation is a common feature in both silencing mechanisms, as is the component of the RNase complex involved in mRNA cleavage. During RNAi in animal cells, dsRNA is processed into siRNA by an RNase III enzyme called Dicer. Here we show that elimination of the activity of two Dicer-like genes by mutation in the fungus Neurospora crassa eliminates transgene-induced gene silencing (quelling) and the processing of dsRNA to an siRNA form. The two Dicer-like genes appear redundant because single mutants are quelling proficient. This first demonstration of the involvement of Dicer in gene silencing induced by transgenes supports a model by which a dsRNA produced by the activity of cellular RNA-dependent RNA polymerases on transgenic transcripts is an essential intermediate of silencing.

In vitro phosphorylation and kinase assays in Neurospora crassa.

Phosphorylation assay is a widespread technique usually necessary for the identification of a specific kinase substrate and/or for the measurement of kinase activity. As an example of the technique, here we describe an assay aimed to test the phosphorylation of the myelin basic protein (MBP) by protein kinase C (PKC), which is overexpressed and purified from Neurospora. The kinase is immunopurified from Neurospora using the expression vector pMYX2 and the FLAG epitope. The purified PKC and the MBP are then incubated in the presence of radioactive ATP, and the phosphorylated product is separated using the polyacrylamide gel electrophoresis technique.

The post-transcriptional gene silencing machinery functions independently of DNA methylation to repress a LINE1-like retrotransposon in Neurospora crassa.

Post-transcriptional gene silencing (PTGS) involving small interfering RNA (siRNA)-directed degradation of RNA transcripts and transcriptional silencing via DNA methylation have each been proposed as mechanisms of genome defence against invading nucleic acids, such as transposons and viruses. Furthermore, recent data from plants indicates that many transposons are silenced via a combination of the two mechanisms, and siRNAs can direct methylation of transposon sequences. We investigated the contribution of DNA methylation and the PTGS pathway to transposon control in the filamentous fungus Neurospora crassa. We found that repression of the LINE1-like transposon, Tad, requires the Argonaute protein QDE2 and Dicer, each of which are required for transgene-induced PTGS (quelling) in N.crassa. Interestingly, unlike quelling, the RNA-dependent RNA polymerase QDE1 and the RecQ DNA helicase QDE3 were not required for Tad control, suggesting the existence of specialized silencing pathways for diverse kinds of repetitive elements. In contrast, Tad elements were not significantly methylated and the DIM2 DNA methyltransferase, responsible for all known DNA methylation in Neurospora, had no effect on Tad control. Thus, an RNAi-related transposon silencing mechanism operates during the vegetative phase of N.crassa that is independent of DNA methylation, highlighting a major difference between this organism and other methylation-proficient species.

The oncogenic role of circPVT1 in head and neck squamous cell carcinoma is mediated through the mutant p53/YAP/TEAD transcription-competent complex.

BACKGROUND: Circular RNAs are a class of endogenous RNAs with various functions in eukaryotic cells. Worthy of note, circular RNAs play a critical role in cancer. Currently, nothing is known about their role in head and neck squamous cell carcinoma (HNSCC). The identification of circular RNAs in HNSCC might become useful for diagnostic and therapeutic strategies in HNSCC. RESULTS: Using samples from 115 HNSCC patients, we find that circPVT1 is over-expressed in tumors compared to matched non-tumoral tissues, with particular enrichment in patients with TP53 mutations. circPVT1 up- and down-regulation determine, respectively, an increase and a reduction of the malignant phenotype in HNSCC cell lines. We show that circPVT1 expression is transcriptionally enhanced by the mut-p53/YAP/TEAD complex. circPVT1 acts as an oncogene modulating the expression of miR-497-5p and genes involved in the control of cell proliferation. CONCLUSIONS: This study shows the oncogenic role of circPVT1 in HNSCC, extending current knowledge about the role of circular RNAs in cancer.

RNAi-dependent and RNAi-independent mechanisms contribute to the silencing of RIPed sequences in Neurospora crassa.

RNA interference (RNAi) can silence genes at the transcriptional level by targeting locus-specific Lys9H3 methylation or at the post-transcriptional level by targeting mRNA degradation. Here we have cloned and sequenced genomic regions methylated in Lys9H3 in Neurospora crassa to test the requirements for components of the RNAi pathway in this modification. We find that 90% of clones map to repeated sequences and relics of transposons that have undergone repeat-induced point mutations (RIP). We find siRNAs derived from transposon relics indicating that the RNAi machinery targets these regions. This is confirmed by the fact that the presence of these siRNAs depends on components of the RNAi pathway such as the RdRP (QDE-1), the putative RecQ helicase (QDE-3) and the two Dicer enzymes. We show that Lys9H3 methylation of RIP sequences is not affected in mutants of the RNAi pathway indicating that the RNAi machinery is not involved in transcriptional gene silencing in Neurospora. We find that RIP regions are transcribed and that the transcript level increases in the mutants of the RNAi pathway. These data suggest that the biological function of the Neurospora RNAi machinery is to control transposon relics and repeated sequences by targeting degradation of transcripts derived from these regions.

The RNA-dependent RNA polymerase, QDE-1, is a rate-limiting factor in post-transcriptional gene silencing in Neurospora crassa.

The RNA-dependent RNA polymerase (RdRP) qde-1 is an essential component of post-transcriptional gene silencing (PTGS), termed 'quelling' in the fungus Neurospora crassa. Here we show that the overexpression of QDE-1 results in a dramatic increase in the efficiency of quelling, with a concomitant net increase in the quantity of al-1 siRNAs. Moreover, in overexpressed strains there is a significant reduction in the number of transgenes required to induce quelling, and an increase in the phenotypic stability despite progressive loss of tandemly repeated transgenes, which normally determines reversion of a silenced phenotype to wild type. These data suggest that the activation and maintenance of silencing in Neurospora appear to rely both on the cellular amount of QDE-1 and the amount of transgenic copies producing RNA molecules that act as a substrate for the RdRP, implicating QDE-1 as a rate-limiting factor in PTGS.

Microarray dataset of Jurkat cells following miR-93 over-expression.

The dataset presented here represents a microarray experiment of Jurkat cell line over-expressing miR-93 after lentiviral transgenic construct transduction. Three biological replicates have been performed. We further provide normalized and processed data, log2 Fold Change based ranked list and GOterms resulting table. The raw microarray data are available in the ArrayExpress database (www.ebi.ac.uk/arrayexpress) under accession number ArrayExpress: E-MTAB-4588.