High-Spatial-Resolution Multi-Omics Sequencing via Deterministic Barcoding in Tissue.

We present deterministic barcoding in tissue for spatial omics sequencing (DBiT-seq) for co-mapping of mRNAs and proteins in a formaldehyde-fixed tissue slide via next-generation sequencing (NGS). Parallel microfluidic channels were used to deliver DNA barcodes to the surface of a tissue slide, and crossflow of two sets of barcodes, A1-50 and B1-50, followed by ligation in situ, yielded a 2D mosaic of tissue pixels, each containing a unique full barcode AB. Application to mouse embryos revealed major tissue types in early organogenesis as well as fine features like microvasculature in a brain and pigmented epithelium in an eye field. Gene expression profiles in 10-µm pixels conformed into the clusters of single-cell transcriptomes, allowing for rapid identification of cell types and spatial distributions. DBiT-seq can be adopted by researchers with no experience in microfluidics and may find applications in a range of fields including developmental biology, cancer biology, neuroscience, and clinical pathology.

Precision analysis of mutant U2AF1 activity reveals deployment of stress granules in myeloid malignancies.

Splicing factor mutations are common among cancers, recently emerging as drivers of myeloid malignancies. U2AF1 carries hotspot mutations in its RNA-binding motifs; however, how they affect splicing and promote cancer remain unclear. The U2AF1/U2AF2 heterodimer is critical for 3' splice site (3'SS) definition. To specifically unmask changes in U2AF1 function in vivo, we developed a crosslinking and immunoprecipitation procedure that detects contacts between U2AF1 and the 3'SS AG at single-nucleotide resolution. Our data reveal that the U2AF1 S34F and Q157R mutants establish new 3'SS contacts at -3 and +1 nucleotides, respectively. These effects compromise U2AF2-RNA interactions, resulting predominantly in intron retention and exon exclusion. Integrating RNA binding, splicing, and turnover data, we predicted that U2AF1 mutations directly affect stress granule components, which was corroborated by single-cell RNA-seq. Remarkably, U2AF1-mutant cell lines and patient-derived MDS/AML blasts displayed a heightened stress granule response, pointing to a novel role for biomolecular condensates in adaptive oncogenic strategies.

m6A Modification Prevents Formation of Endogenous Double-Stranded RNAs and Deleterious Innate Immune Responses during Hematopoietic Development.

N6-methyladenosine (m6A) is the most abundant RNA modification, but little is known about its role in mammalian hematopoietic development. Here, we show that conditional deletion of the m6A writer METTL3 in murine fetal liver resulted in hematopoietic failure and perinatal lethality. Loss of METTL3 and m6A activated an aberrant innate immune response, mediated by the formation of endogenous double-stranded RNAs (dsRNAs). The aberrantly formed dsRNAs were long, highly m6A modified in their native state, characterized by low folding energies, and predominantly protein coding. We identified coinciding activation of pattern recognition receptor pathways normally tasked with the detection of foreign dsRNAs. Disruption of the aberrant immune response via abrogation of downstream Mavs or Rnasel signaling partially rescued the observed hematopoietic defects in METTL3-deficient cells in vitro and in vivo. Our results suggest that m6A modification protects against endogenous dsRNA formation and a deleterious innate immune response during mammalian hematopoietic development.

HuD Is a Neural Translation Enhancer Acting on mTORC1-Responsive Genes and Counteracted by the Y3 Small Non-coding RNA.

The RNA-binding protein HuD promotes neurogenesis and favors recovery from peripheral axon injury. HuD interacts with many mRNAs, altering both stability and translation efficiency. We generated a nucleotide resolution map of the HuD RNA interactome in motor neuron-like cells, identifying HuD target sites in 1,304 mRNAs, almost exclusively in the 3' UTR. HuD binds many mRNAs encoding mTORC1-responsive ribosomal proteins and translation factors. Altered HuD expression correlates with the translation efficiency of these mRNAs and overall protein synthesis, in a mTORC1-independent fashion. The predominant HuD target is the abundant, small non-coding RNA Y3, amounting to 70% of the HuD interaction signal. Y3 functions as a molecular sponge for HuD, dynamically limiting its recruitment to polysomes and its activity as a translation and neuron differentiation enhancer. These findings uncover an alternative route to the mTORC1 pathway for translational control in motor neurons that is tunable by a small non-coding RNA.

Widespread association of the Argonaute protein AGO2 with meiotic chromatin suggests a distinct nuclear function in mammalian male reproduction.

Argonaute 2 (AGO2) is a ubiquitously expressed protein critical for regulation of mRNA translation and vital to animal development. AGO2 protein is found in both cytoplasmic and nuclear compartments, and although its cytoplasmic role is well studied, the biological relevance of nuclear AGO2 is unclear. Here, we address this problem in vivo using spermatogenic cells as a model. We find that AGO2 transiently binds both chromatin and nucleus-specific mRNA transcripts of hundreds of genes required for sperm production during male meiosis in mice, and that germline conditional knockout (cKO) of Ago2 causes depletion of the encoded proteins. Correspondingly, Ago2 cKO males show abnormal sperm head morphology and reduced sperm count, along with reduced postnatal viability of offspring. Together, our data reveal an unexpected nuclear role for AGO2 in enhancing expression of developmentally important genes during mammalian male reproduction.

Pre-mRNA splicing factor U2AF2 recognizes distinct conformations of nucleotide variants at the center of the pre-mRNA splice site signal.

The essential pre-mRNA splicing factor U2AF2 (also called U2AF65) identifies polypyrimidine (Py) tract signals of nascent transcripts, despite length and sequence variations. Previous studies have shown that the U2AF2 RNA recognition motifs (RRM1 and RRM2) preferentially bind uridine-rich RNAs. Nonetheless, the specificity of the RRM1/RRM2 interface for the central Py tract nucleotide has yet to be investigated. We addressed this question by determining crystal structures of U2AF2 bound to a cytidine, guanosine, or adenosine at the central position of the Py tract, and compared U2AF2-bound uridine structures. Local movements of the RNA site accommodated the different nucleotides, whereas the polypeptide backbone remained similar among the structures. Accordingly, molecular dynamics simulations revealed flexible conformations of the central, U2AF2-bound nucleotide. The RNA binding affinities and splicing efficiencies of structure-guided mutants demonstrated that U2AF2 tolerates nucleotide substitutions at the central position of the Py tract. Moreover, enhanced UV-crosslinking and immunoprecipitation of endogenous U2AF2 in human erythroleukemia cells showed uridine-sensitive binding sites, with lower sequence conservation at the central nucleotide positions of otherwise uridine-rich, U2AF2-bound splice sites. Altogether, these results highlight the importance of RNA flexibility for protein recognition and take a step towards relating splice site motifs to pre-mRNA splicing efficiencies.

Hypomorphic mutation of the mouse Huntington's disease gene orthologue.

Rare individuals with inactivating mutations in the Huntington's disease gene (HTT) exhibit variable abnormalities that imply essential HTT roles during organ development. Here we report phenotypes produced when increasingly severe hypomorphic mutations in the murine HTT orthologue Htt, (HdhneoQ20, HdhneoQ50, HdhneoQ111), were placed over a null allele (Hdhex4/5). The most severe hypomorphic allele failed to rescue null lethality at gastrulation, while the intermediate, though still severe, alleles yielded recessive perinatal lethality and a variety of fetal abnormalities affecting body size, skin, skeletal and ear formation, and transient defects in hematopoiesis. Comparative molecular analysis of wild-type and Htt-null retinoic acid-differentiated cells revealed gene network dysregulation associated with organ development that nominate polycomb repressive complexes and miRNAs as molecular mediators. Together these findings demonstrate that Htt is required both pre- and post-gastrulation to support normal development.

Low iron promotes megakaryocytic commitment of megakaryocytic-erythroid progenitors in humans and mice.

The mechanisms underlying thrombocytosis in patients with iron deficiency anemia remain unknown. Here, we present findings that support the hypothesis that low iron biases the commitment of megakaryocytic (Mk)-erythroid progenitors (MEPs) toward the Mk lineage in both human and mouse. In MEPs of transmembrane serine protease 6 knockout (Tmprss6-/-) mice, which exhibit iron deficiency anemia and thrombocytosis, we observed a Mk bias, decreased labile iron, and decreased proliferation relative to wild-type (WT) MEPs. Bone marrow transplantation assays suggest that systemic iron deficiency, rather than a local role for Tmprss6-/- in hematopoietic cells, contributes to the MEP lineage commitment bias observed in Tmprss6-/- mice. Nontransgenic mice with acquired iron deficiency anemia also show thrombocytosis and Mk-biased MEPs. Gene expression analysis reveals that messenger RNAs encoding genes involved in metabolic, vascular endothelial growth factor, and extracellular signal-regulated kinase (ERK) pathways are enriched in Tmprss6-/- vs WT MEPs. Corroborating our findings from the murine models of iron deficiency anemia, primary human MEPs exhibit decreased proliferation and Mk-biased commitment after knockdown of transferrin receptor 2, a putative iron sensor. Signal transduction analyses reveal that both human and murine MEP have lower levels of phospho-ERK1/2 in iron-deficient conditions compared with controls. These data are consistent with a model in which low iron in the marrow environment affects MEP metabolism, attenuates ERK signaling, slows proliferation, and biases MEPs toward Mk lineage commitment.

Hydrogen peroxide is a neuronal alarmin that triggers specific RNAs, local translation of Annexin A2, and cytoskeletal remodeling in Schwann cells.

Schwann cells are key players in neuro-regeneration: They sense "alarm" signals released by degenerating nerve terminals and differentiate toward a proregenerative phenotype, with phagocytosis of nerve debris and nerve guidance. At the murine neuromuscular junction, hydrogen peroxide (H2O2) is a key signal of Schwann cells' activation in response to a variety of nerve injuries. Here we report that Schwann cells exposed to low doses of H2O2 rewire the expression of several RNAs at both transcriptional and translational levels. Among the genes positively regulated at both levels, we identified an enriched cluster involved in cytoskeleton remodeling and cell migration, with the Annexin (Anxa) proteins being the most represented family. We show that both Annexin A2 (Anxa2) transcript and protein accumulate at the tips of long pseudopods that Schwann cells extend upon H2O2 exposure. Interestingly, Schwann cells reply to this signal and to nerve injury by locally translating Anxa2 in pseudopods, and undergo an extensive cytoskeleton remodeling. Our results show that, similarly to neurons, Schwann cells take advantage of local protein synthesis to change shape and move toward damaged axonal terminals to facilitate axonal regeneration.

riboWaltz: Optimization of ribosome P-site positioning in ribosome profiling data.

Ribosome profiling is a powerful technique used to study translation at the genome-wide level, generating unique information concerning ribosome positions along RNAs. Optimal localization of ribosomes requires the proper identification of the ribosome P-site in each ribosome protected fragment, a crucial step to determine the trinucleotide periodicity of translating ribosomes, and draw correct conclusions concerning where ribosomes are located. To determine the P-site within ribosome footprints at nucleotide resolution, the precise estimation of its offset with respect to the protected fragment is necessary. Here we present riboWaltz, an R package for calculation of optimal P-site offsets, diagnostic analysis and visual inspection of ribosome profiling data. Compared to existing tools, riboWaltz shows improved accuracies for P-site estimation and neat ribosome positioning in multiple case studies. riboWaltz was implemented in R and is available as an R package at https://github.com/LabTranslationalArchitectomics/RiboWaltz.

Identification and dynamic changes of RNAs isolated from RALY-containing ribonucleoprotein complexes.

RALY is a member of the heterogeneous nuclear ribonucleoprotein family (hnRNP), a large family of RNA-binding proteins involved in many aspects of RNA metabolism. Although RALY interactome has been recently characterized, a comprehensive global analysis of RALY-associated RNAs is lacking and the biological function of RALY remains elusive. Here, we performed RIP-seq analysis to identify RALY interacting RNAs and assessed the role of RALY in gene expression. We demonstrate that RALY binds specific coding and non-coding RNAs and associates with translating mRNAs of mammalian cells. Among the identified transcripts, we focused on ANXA1 and H1FX mRNAs, encoding for Annexin A1 and for the linker variant of the histone H1X, respectively. Both proteins are differentially expressed by proliferating cells and are considered as markers for tumorigenesis. We demonstrate that cells lacking RALY expression exhibit changes in the levels of H1FX and ANXA1 mRNAs and proteins in an opposite manner. We also provide evidence for a direct binding of RALY to the U-rich elements present within the 3΄UTR of both transcripts. Thus, our results identify RALY as a poly-U binding protein and as a regulator of H1FX and ANXA1 in mammalian cells.

The hnRNP RALY regulates transcription and cell proliferation by modulating the expression of specific factors including the proliferation marker E2F1.

The heterogeneous nuclear ribonucleoproteins (hnRNP) form a large family of RNA-binding proteins that exert numerous functions in RNA metabolism. RALY is a member of the hnRNP family that binds poly-U-rich elements within several RNAs and regulates the expression of specific transcripts. RALY is up-regulated in different types of cancer, and its down-regulation impairs cell cycle progression. However, the RALY's role in regulating RNA levels remains elusive. Here, we show that numerous genes coding for factors involved in transcription and cell cycle regulation exhibit an altered expression in RALY-down-regulated HeLa cells, consequently causing impairments in transcription, cell proliferation, and cell cycle progression. Interestingly, by comparing the list of RALY targets with the list of genes affected by RALY down-regulation, we found an enrichment of RALY mRNA targets in the down-regulated genes upon RALY silencing. The affected genes include the E2F transcription factor family. Given its role as proliferation-promoting transcription factor, we focused on E2F1. We demonstrate that E2F1 mRNA stability and E2F1 protein levels are reduced in cells lacking RALY expression. Finally, we also show that RALY interacts with transcriptionally active chromatin in both an RNA-dependent and -independent manner and that this association is abolished in the absence of active transcription. Taken together, our results highlight the importance of RALY as an indirect regulator of transcription and cell cycle progression through the regulation of specific mRNA targets, thus strengthening the possibility of a direct gene expression regulation exerted by RALY.

The Typhoid Toxin Promotes Host Survival and the Establishment of a Persistent Asymptomatic Infection.

Bacterial genotoxins, produced by several Gram-negative bacteria, induce DNA damage in the target cells. While the responses induced in the host cells have been extensively studied in vitro, the role of these effectors during the course of infection remains poorly characterized. To address this issue, we assessed the effects of the Salmonella enterica genotoxin, known as typhoid toxin, in in vivo models of murine infection. Immunocompetent mice were infected with isogenic S. enterica, serovar Typhimurium (S. Typhimurium) strains, encoding either a functional or an inactive typhoid toxin. The presence of the genotoxic subunit was detected 10 days post-infection in the liver of infected mice. Unexpectedly, its expression promoted the survival of the host, and was associated with a significant reduction of severe enteritis in the early phases of infection. Immunohistochemical and transcriptomic analysis confirmed the toxin-mediated suppression of the intestinal inflammatory response. The presence of a functional typhoid toxin further induced an increased frequency of asymptomatic carriers. Our data indicate that the typhoid toxin DNA damaging activity increases host survival and favours long-term colonization, highlighting a complex cross-talk between infection, DNA damage response and host immune response. These findings may contribute to understand why such effectors have been evolutionary conserved and horizontally transferred among Gram-negative bacteria.

RNAcommender: genome-wide recommendation of RNA-protein interactions.

MOTIVATION: Information about RNA-protein interactions is a vital pre-requisite to tackle the dissection of RNA regulatory processes. Despite the recent advances of the experimental techniques, the currently available RNA interactome involves a small portion of the known RNA binding proteins. The importance of determining RNA-protein interactions, coupled with the scarcity of the available information, calls for in silico prediction of such interactions. RESULTS: We present RNAcommender, a recommender system capable of suggesting RNA targets to unexplored RNA binding proteins, by propagating the available interaction information taking into account the protein domain composition and the RNA predicted secondary structure. Our results show that RNAcommender is able to successfully suggest RNA interactors for RNA binding proteins using little or no interaction evidence. RNAcommender was tested on a large dataset of human RBP-RNA interactions, showing a good ranking performance (average AUC ROC of 0.75) and significant enrichment of correct recommendations for 75% of the tested RBPs. RNAcommender can be a valid tool to assist researchers in identifying potential interacting candidates for the majority of RBPs with uncharacterized binding preferences. AVAILABILITY AND IMPLEMENTATION: The software is freely available at http://rnacommender.disi.unitn.it CONTACT: gianluca.corrado@unitn.it or andrea.passerini@unitn.itSupplementary information: Supplementary data are available at Bioinformatics online.

RiboAbacus: a model trained on polyribosome images predicts ribosome density and translational efficiency from mammalian transcriptomes.

Fluctuations in mRNA levels only partially contribute to determine variations in mRNA availability for translation, producing the well-known poor correlation between transcriptome and proteome data. Recent advances in microscopy now enable researchers to obtain high resolution images of ribosomes on transcripts, providing precious snapshots of translation in vivo. Here we propose RiboAbacus, a mathematical model that for the first time incorporates imaging data in a predictive model of transcript-specific ribosome densities and translational efficiencies. RiboAbacus uses a mechanistic model of ribosome dynamics, enabling the quantification of the relative importance of different features (such as codon usage and the 5' ramp effect) in determining the accuracy of predictions. The model has been optimized in the human Hek-293 cell line to fit thousands of images of human polysomes obtained by atomic force microscopy, from which we could get a reference distribution of the number of ribosomes per mRNA with unmatched resolution. After validation, we applied RiboAbacus to three case studies of known transcriptome-proteome datasets for estimating the translational efficiencies, resulting in an increased correlation with corresponding proteomes. RiboAbacus is an intuitive tool that allows an immediate estimation of crucial translation properties for entire transcriptomes, based on easily obtainable transcript expression levels.

Global translation variations in host cells upon attack of lytic and sublytic Staphylococcus aureus α-haemolysin.

Genome-wide analyses of translation can provide major contributions in our understanding of the complex interplay between virulent factors and host cells. So far, the activation of host translational control mechanisms by bacterial toxins, owing to specific recruitment of mRNAs, RNA-binding proteins (RBPs) and ncRNAs (non-coding RNAs), are far from being understood. In the present study, we characterize for the first time the changes experienced by the translational control system of host cells in response to the well-known Staphylococcus aureus α-haemolysin (AHL) under both sublytic and lytic conditions. By comparing variations occurring in the cellular transcriptome and translatome, we give evidence that global gene expression is primarily rewired at the translational level, with the contribution of the RBP ELAVL1 (HuR) in the sublytic response. These results reveal the importance of translational control during host-pathogen interaction, opening new approaches for AHL-induced diseases.

Translational downregulation of HSP90 expression by iron chelators in neuroblastoma cells.

Iron is an essential cellular nutrient, being a critical cofactor of several proteins involved in cell growth and replication. Compared with normal cells, neoplastic cells have been shown to require a greater amount of iron, thus laying the basis for the promising anticancer activity of iron chelators. In this work, we evaluated the effects of molecules with iron chelation activity on neuroblastoma (NB) cell lines. Of the 17 iron chelators tested, six reduced cell viability of two NB cell lines with an inhibition of growth of 50% below 10 µM; four of the six molecules-ciclopirox olamine (CPX), piroctone, 8-hydroxyquinoline, and deferasirox-were also shown to efficiently chelate intracellular iron within minutes after addition. Effects on cell viability of one of the compounds, CPX, were indeed dependent on chelation of intracellular iron and mediated by both G0/G1 cell cycle block and induction of apoptosis. By combined transcriptome and translatome profiling we identified early translational downregulation of several members of the heat shock protein group as a specific effect of CPX treatment. We functionally confirmed iron-dependent depletion of HSP90 and its client proteins at pharmacologically achievable concentrations of CPX, and we extended this effect to piroctone, 8-hydroxyquinoline, and deferasirox. Given the documented sensitivity of NB cells to HSP90 inhibition, we propose CPX and other iron chelators as investigational antitumor agents in NB therapy.

Whole-genome cartography of p53 response elements ranked on transactivation potential.

BACKGROUND: Many recent studies using ChIP-seq approaches cross-referenced to trascriptome data and also to potentially unbiased in vitro DNA binding selection experiments are detailing with increasing precision the p53-directed gene regulatory network that, nevertheless, is still expanding. However, most experiments have been conducted in established cell lines subjected to specific p53-inducing stimuli, both factors potentially biasing the results. RESULTS: We developed p53retriever, a pattern search algorithm that maps p53 response elements (REs) and ranks them according to predicted transactivation potentials in five classes. Besides canonical, full site REs, we developed specific pattern searches for non-canonical half sites and 3/4 sites and show that they can mediate p53-dependent responsiveness of associated coding sequences. Using ENCODE data, we also mapped p53 REs in about 44,000 distant enhancers and identified a 16-fold enrichment for high activity REs within those sites in the comparison with genomic regions near transcriptional start sites (TSS). Predictions from our pattern search were cross-referenced to ChIP-seq, ChIP-exo, expression, and various literature data sources. Based on the mapping of predicted functional REs near TSS, we examined expression changes of thirteen genes as a function of different p53-inducing conditions, providing further evidence for PDE2A, GAS6, E2F7, APOBEC3H, KCTD1, TRIM32, DICER, HRAS, KITLG and TGFA p53-dependent regulation, while MAP2K3, DNAJA1 and potentially YAP1 were identified as new direct p53 target genes. CONCLUSIONS: We provide a comprehensive annotation of canonical and non-canonical p53 REs in the human genome, ranked on predicted transactivation potential. We also establish or corroborate direct p53 transcriptional control of thirteen genes. The entire list of identified and functionally classified p53 REs near all UCSC-annotated genes and within ENCODE mapped enhancer elements is provided. Our approach is distinct from, and complementary to, existing methods designed to identify p53 response elements. p53retriever is available as an R package at: http://tomateba.github.io/p53retriever .

tRanslatome: an R/Bioconductor package to portray translational control.

UNLABELLED: High-throughput technologies have led to an explosion of genomic data available for automated analysis. The consequent possibility to simultaneously sample multiple layers of variation along the gene expression flow requires computational methods integrating raw information from different '-omics'. It has been recently demonstrated that translational control is a widespread phenomenon, with profound and still underestimated regulation capabilities. Although detecting changes in the levels of total messenger RNAs (mRNAs; the transcriptome), of polysomally loaded mRNAs (the translatome) and of proteins (the proteome) is experimentally feasible in a high-throughput way, the integration of these levels is still far from being robustly approached. Here we introduce tRanslatome, a new R/Bioconductor package, which is a complete platform for the simultaneous pairwise analysis of transcriptome, translatome and proteome data. The package includes most of the available statistical methods developed for the analysis of high-throughput data, allowing the parallel comparison of differentially expressed genes and the corresponding differentially enriched biological themes. Notably, it also enables the prediction of translational regulatory elements on mRNA sequences. The utility of this tool is demonstrated with two case studies. AVAILABILITY AND IMPLEMENTATION: tRanslatome is available in Bioconductor.

TrkA is amplified in malignant melanoma patients and induces an anti-proliferative response in cell lines.

BACKGROUND: The nerve growth factor (NGF) receptor tyrosine-kinase TrkA is a well-known determinant of the melanocytic lineage, through modulation of the MAPK and AKT cascades. While TrkA gene is frequently rearranged in cancers, its involvement in malignant melanoma (MM) development is still unclear. METHODS: We analyzed a dataset of primary cutaneous MM (n = 31) by array comparative genomic hybridization (aCGH), to identify genomic amplifications associated with tumor progression. The analysis was validated by genomic quantitative PCR (qPCR) on an extended set of cases (n = 64) and the results were correlated with the clinical outcome. To investigate TrkA molecular pathways and cellular function, we generated inducible activation of the NGF-TrkA signaling in human MM cell lines. RESULTS: We identified amplification of 1q23.1, where the TrkA locus resides, as a candidate hotspot implicated in the progression of MM. Across 40 amplicons detected, segmental amplification of 1q23.1 showed the strongest association with tumor thickness. By validation of the analysis, TrkA gene amplification emerged as a frequent event in primary melanomas (50 % of patients), and correlated with worse clinical outcome. However, experiments in cell lines revealed that induction of the NGF-TrkA signaling produced a phenotype of dramatic suppression of cell proliferation through inhibition of cell division and pronounced intracellular vacuolization, in a way straightly dependent on NGF activation of TrkA. These events were triggered via MAPK activity but not via AKT, and involved p21(cip1) protein increase, compatibly with a mechanism of oncogene-induced growth arrest. CONCLUSIONS: Taken together, our findings point to TrkA as a candidate oncogene in MM and support a model in which the NGF-TrkA-MAPK pathway may mediate a trade-off between neoplastic transformation and adaptive anti-proliferative response.