Natural RNA Polymerase Aptamers Regulate Transcription in E. coli.

In search for RNA signals that modulate transcription via direct interaction with RNA polymerase (RNAP), we deep sequenced an E. coli genomic library enriched for RNAP-binding RNAs. Many natural RNAP-binding aptamers, termed RAPs, were mapped to the genome. Over 60% of E. coli genes carry RAPs in their mRNA. Combining in vitro and in vivo approaches, we characterized a subset of inhibitory RAPs (iRAPs) that promote Rho-dependent transcription termination. A representative iRAP within the coding region of the essential gene, nadD, greatly reduces its transcriptional output in stationary phase and under oxidative stress, demonstrating that iRAPs control gene expression in response to changing environment. The mechanism of iRAPs involves active uncoupling of transcription and translation, making nascent RNA accessible to Rho. iRAPs encoded in the antisense strand also promote gene expression by reducing transcriptional interference. In essence, our work uncovers a broad class of cis-acting RNA signals that globally control bacterial transcription.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions.

BACKGROUND: Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown. RESULTS: Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels. CONCLUSIONS: Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/ ) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

A cadherin switch marks germ layer formation in the diploblastic sea anemone Nematostella vectensis.

Morphogenesis is a shape-building process during development of multicellular organisms. During this process, the establishment and modulation of cell-cell contacts play an important role. Cadherins, the major cell adhesion molecules, form adherens junctions connecting epithelial cells. Numerous studies of Bilateria have shown that cadherins are associated with the regulation of cell differentiation, cell shape changes, cell migration and tissue morphogenesis. To date, the role of cadherins in non-bilaterians is unknown. Here, we study the expression and function of two paralogous classical cadherins, Cadherin 1 and Cadherin 3, in a diploblastic animal, the sea anemone Nematostella vectensis We show that a cadherin switch accompanies the formation of germ layers. Using specific antibodies, we show that both cadherins are localized to adherens junctions at apical and basal positions in ectoderm and endoderm. During gastrulation, partial epithelial-to-mesenchymal transition of endodermal cells is marked by stepwise downregulation of Cadherin 3 and upregulation of Cadherin 1. Knockdown experiments show that both cadherins are required for maintenance of tissue integrity and tissue morphogenesis. Thus, both sea anemones and bilaterians use independently duplicated cadherins combinatorially for tissue morphogenesis and germ layer differentiation.

Conservation and turnover of miRNAs and their highly complementary targets in early branching animals.

MicroRNAs (miRNAs) are crucial post-transcriptional regulators that have been extensively studied in Bilateria, a group comprising the majority of extant animals, where more than 30 conserved miRNA families have been identified. By contrast, bilaterian miRNA targets are largely not conserved. Cnidaria is the sister group to Bilateria and thus provides a unique opportunity for comparative studies. Strikingly, like their plant counterparts, cnidarian miRNAs have been shown to predominantly have highly complementary targets leading to transcript cleavage by Argonaute proteins. Here, we assess the conservation of miRNAs and their targets by small RNA sequencing followed by miRNA target prediction in eight species of Anthozoa (sea anemones and corals), the earliest-branching cnidarian class. We uncover dozens of novel miRNAs but only a few conserved ones. Further, given their high complementarity, we were able to computationally identify miRNA targets in each species. Besides evidence for conservation of specific miRNA target sites, which are maintained between sea anemones and stony corals across 500 Myr of evolution, we also find indications for convergent evolution of target regulation by different miRNAs. Our data indicate that cnidarians have only few conserved miRNAs and corresponding targets, despite their high complementarity, suggesting a high evolutionary turnover.

Characterization of the piRNA pathway during development of the sea anemone Nematostella vectensis.

PIWI-interacting RNAs (piRNAs) and associated proteins comprise a conserved pathway for silencing transposons in metazoan germlines. piRNA pathway components are also expressed in multipotent somatic stem cells in various organisms. piRNA functions have been extensively explored in bilaterian model systems, however, comprehensive studies in non-bilaterian phyla remain limited. Here we investigate the piRNA pathway during the development of Nematostella vectensis, a well-established model system belonging to Cnidaria, the sister group to Bilateria. To date, no population of somatic stem cells has been identified in this organism, despite its long life-span and regenerative capacities that require a constant cell-renewal. We show that Nematostella piRNA pathway components are broadly expressed in early developmental stages, while piRNAs themselves show differential expression, suggesting specific developmental roles of distinct piRNA families. In adults, piRNA associated proteins are enriched in the germline but also expressed in somatic cells, indicating putative stem cell properties. Furthermore, we provide experimental evidence that Nematostella piRNAs cleave transposable elements as well as protein-coding genes. Our results demonstrate that somatic expression of piRNA associated proteins as well as the roles of piRNAs in transposon repression and gene regulation are likely ancestral features that evolved before the split between Cnidaria and Bilateria.

The double-stranded transcriptome of Escherichia coli.

Advances in high-throughput transcriptome analyses have revealed hundreds of antisense RNAs (asRNAs) for many bacteria, although few have been characterized, and the number of functional asRNAs remains unknown. We have developed a genome-wide high-throughput method to identify functional asRNAs in vivo. Most mechanisms of gene regulation via asRNAs require an RNA-RNA interaction with its target RNA, and we hypothesized that a functional asRNA would be found in a double strand (dsRNA), duplexed with its cognate RNA in a single cell. We developed a method of isolating dsRNAs from total RNA by immunoprecipitation with a ds-RNA specific antibody. Total RNA and immunoprecipitated dsRNA from Escherichia coli RNase III WT and mutant strains were deep-sequenced. A statistical model was applied to filter for biologically relevant dsRNA regions, which were subsequently categorized by location relative to annotated genes. A total of 316 potentially functional asRNAs were identified in the RNase III mutant strain and are encoded primarily opposite to the 5' ends of transcripts, but are also found opposite ncRNAs, gene junctions, and the 3' ends. A total of 21 sense/antisense RNA pairs identified in dsRNAs were confirmed by Northern blot analyses. Most of the RNA steady-state levels were higher or detectable only in the RNase III mutant strain. Taken together, our data indicate that a significant amount of dsRNA is formed in the cell, that RNase III degrades or processes these dsRNAs, and that dsRNA plays a major role in gene regulation in E. coli.

Analysis of SMAD1/5 target genes in a sea anemone reveals ZSWIM4-6 as a novel BMP signaling modulator.

BMP signaling has a conserved function in patterning the dorsal-ventral body axis in Bilateria and the directive axis in anthozoan cnidarians. So far, cnidarian studies have focused on the role of different BMP signaling network components in regulating pSMAD1/5 gradient formation. Much less is known about the target genes downstream of BMP signaling. To address this, we generated a genome-wide list of direct pSMAD1/5 target genes in the anthozoan Nematostella vectensis, several of which were conserved in Drosophila and Xenopus. Our ChIP-seq analysis revealed that many of the regulatory molecules with documented bilaterally symmetric expression in Nematostella are directly controlled by BMP signaling. We identified several so far uncharacterized BMP-dependent transcription factors and signaling molecules, whose bilaterally symmetric expression may be indicative of their involvement in secondary axis patterning. One of these molecules is zswim4-6, which encodes a novel nuclear protein that can modulate the pSMAD1/5 gradient and potentially promote BMP-dependent gene repression.

Phase I Study of the Liposomal Formulation of Eribulin (E7389-LF): Results from the Advanced Gastric Cancer Expansion Cohort.

PURPOSE: In the dose-expansion part of this open-label, phase I study, we explored the efficacy and safety of E7389-LF (liposomal formulation of eribulin) in Japanese patients with advanced gastric cancer. PATIENTS AND METHODS: Patients with advanced gastric cancer who had been previously treated with ≥2 lines of chemotherapy received E7389-LF 2.0 mg/m2 every 3 weeks (the previously determined maximum tolerated dose, the primary objective of Study 114). Secondary objectives included objective response rate (ORR), progression-free survival (PFS), and safety; exploratory objectives included disease control rate (DCR) and clinical benefit rate (CBR), as well as pharmacodynamic measurements of serum biomarkers. RESULTS: As of June 24, 2021, 34 patients were enrolled and treated (10 from the original dose-expansion cohort, expanded to include 24 additional patients). Six patients had partial responses, for an ORR of 17.6% [95% confidence interval (CI), 6.8-34.5], and the median PFS was 3.7 months (95% CI, 2.7-4.8). The DCR was 79.4% (95% CI, 62.1-91.3), and the CBR was 32.4% (95% CI, 17.4-50.5). Overall, 32 patients (94.1%) experienced treatment-related adverse events, and 26 patients (76.5%) experienced grade ≥3 events, most commonly neutropenia (41.2%) and leukopenia (29.4%). Of the 8 endothelial cell/vasculature markers tested in this study, 7 were significantly increased among patients treated with E7389-LF; these changes were generally consistent regardless of best overall response. CONCLUSIONS: E7389-LF 2.0 mg/m2 every 3 weeks was tolerable and showed preliminary activity for the treatment of patients with gastric cancer.

Muscle cell-type diversification is driven by bHLH transcription factor expansion and extensive effector gene duplications.

Animals are typically composed of hundreds of different cell types, yet mechanisms underlying the emergence of new cell types remain unclear. Here we address the origin and diversification of muscle cells in the non-bilaterian, diploblastic sea anemone Nematostella vectensis. We discern two fast and two slow-contracting muscle cell populations, which differ by extensive sets of paralogous structural protein genes. We find that the regulatory gene set of the slow cnidarian muscles is remarkably similar to the bilaterian cardiac muscle, while the two fast muscles differ substantially from each other in terms of transcription factor profiles, though driving the same set of structural protein genes and having similar physiological characteristics. We show that anthozoan-specific paralogs of Paraxis/Twist/Hand-related bHLH transcription factors are involved in the formation of fast and slow muscles. Our data suggest that the subsequent recruitment of an entire effector gene set from the inner cell layer into the neural ectoderm contributes to the evolution of a novel muscle cell type. Thus, we conclude that extensive transcription factor gene duplications and co-option of effector modules act as an evolutionary mechanism underlying cell type diversification during metazoan evolution.

Topological structures and syntenic conservation in sea anemone genomes.

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

An ancestral Wnt-Brachyury feedback loop in axial patterning and recruitment of mesoderm-determining target genes.

Transcription factors are crucial drivers of cellular differentiation during animal development and often share ancient evolutionary origins. The T-box transcription factor Brachyury plays a pivotal role as an early mesoderm determinant and neural repressor in vertebrates; yet, the ancestral function and key evolutionary transitions of the role of this transcription factor remain obscure. Here, we present a genome-wide target-gene screen using chromatin immunoprecipitation sequencing in the sea anemone Nematostella vectensis, an early branching non-bilaterian, and the sea urchin Strongylocentrotus purpuratus, a representative of the sister lineage of chordates. Our analysis reveals an ancestral gene regulatory feedback loop connecting Brachyury, FoxA and canonical Wnt signalling involved in axial patterning that predates the cnidarian-bilaterian split about 700 million years ago. Surprisingly, we also found that part of the gene regulatory network controlling the fate of neuromesodermal progenitors in vertebrates was already present in the common ancestor of cnidarians and bilaterians. However, while several endodermal and neuronal Brachyury target genes are ancestrally shared, hardly any of the key mesodermal downstream targets in vertebrates are found in the sea anemone or the sea urchin. Our study suggests that a limited number of target genes involved in mesoderm formation were newly acquired in the vertebrate lineage, leading to a dramatic shift in the function of this ancestral developmental regulator.

Emergence of novel cephalopod gene regulation and expression through large-scale genome reorganization.

Coleoid cephalopods (squid, cuttlefish, octopus) have the largest nervous system among invertebrates that together with many lineage-specific morphological traits enables complex behaviors. The genomic basis underlying these innovations remains unknown. Using comparative and functional genomics in the model squid Euprymna scolopes, we reveal the unique genomic, topological, and regulatory organization of cephalopod genomes. We show that coleoid cephalopod genomes have been extensively restructured compared to other animals, leading to the emergence of hundreds of tightly linked and evolutionary unique gene clusters (microsyntenies). Such novel microsyntenies correspond to topological compartments with a distinct regulatory structure and contribute to complex expression patterns. In particular, we identify a set of microsyntenies associated with cephalopod innovations (MACIs) broadly enriched in cephalopod nervous system expression. We posit that the emergence of MACIs was instrumental to cephalopod nervous system evolution and propose that microsyntenic profiling will be central to understanding cephalopod innovations.

Genomic SELEX: a discovery tool for genomic aptamers.

Genomic SELEX is a discovery tool for genomic aptamers, which are genomically encoded functional domains in nucleic acid molecules that recognize and bind specific ligands. When combined with genomic libraries and using RNA-binding proteins as baits, Genomic SELEX used with high-throughput sequencing enables the discovery of genomic RNA aptamers and the identification of RNA-protein interaction networks. Here we describe how to construct and analyze genomic libraries, how to choose baits for selections, how to perform the selection procedure and finally how to analyze the enriched sequences derived from deep sequencing. As a control procedure, we recommend performing a "Neutral" SELEX experiment in parallel to the selection, omitting the selection step. This control experiment provides a background signal for comparison with the positively selected pool. We also recommend deep sequencing the initial library in order to facilitate the final in silico analysis of enrichment with respect to the initial levels. Counter selection procedures, using modified or inactive baits, allow strengthening the binding specificity of the winning selected sequences.

Cnidarian-bilaterian comparison reveals the ancestral regulatory logic of the ß-catenin dependent axial patterning.

In animals, body axis patterning is based on the concentration-dependent interpretation of graded morphogen signals, which enables correct positioning of the anatomical structures. The most ancient axis patterning system acting across animal phyla relies on ß-catenin signaling, which directs gastrulation, and patterns the main body axis. However, within Bilateria, the patterning logic varies significantly between protostomes and deuterostomes. To deduce the ancestral principles of ß-catenin-dependent axial patterning, we investigate the oral-aboral axis patterning in the sea anemone Nematostella-a member of the bilaterian sister group Cnidaria. Here we elucidate the regulatory logic by which more orally expressed ß-catenin targets repress more aborally expressed ß-catenin targets, and progressively restrict the initially global, maternally provided aboral identity. Similar regulatory logic of ß-catenin-dependent patterning in Nematostella and deuterostomes suggests a common evolutionary origin of these processes and the equivalence of the cnidarian oral-aboral and the bilaterian posterior-anterior body axes.

RNA polymerase II-binding aptamers in human ACRO1 satellites disrupt transcription in cis.

Transcription elongation is a highly regulated process affected by many proteins, RNAs and the underlying DNA. Here we show that the nascent RNA can interfere with transcription in human cells, extending our previous findings from bacteria and yeast. We identified a variety of Pol II-binding aptamers (RAPs), prominent in repeat elements such as ACRO1 satellites, LINE1 retrotransposons and CA simple repeats, and also in several protein-coding genes. ACRO1 repeat, when translated in silico, exhibits ~50% identity with the Pol II CTD sequence. Taken together with a recent proposal that proteins in general tend to interact with RNAs similar to their cognate mRNAs, this suggests a mechanism for RAP binding. Using a reporter construct, we show that ACRO1 potently inhibits Pol II elongation in cis. We propose a novel mode of transcriptional regulation in humans, in which the nascent RNA binds Pol II to silence its own expression.

Ancient animal genome architecture reflects cell type identities.

The level of conservation of ancient metazoan gene order (synteny) is remarkable. Despite this, the functionality of the vast majority of such regions in metazoan genomes remains elusive. Utilizing recently published single-cell expression data from several anciently diverging metazoan species, we reveal the level of correspondence between cell types and genomic synteny, identifying genomic regions conferring ancient cell type identity.

Nascent RNA signaling to yeast RNA Pol II during transcription elongation.

Transcription as the key step in gene expression is a highly regulated process. The speed of transcription elongation depends on the underlying gene sequence and varies on a gene by gene basis. The reason for this sequence dependence is not known in detail. Recently, our group studied the cross talk between the nascent RNA and the transcribing RNA polymerase by screening the Escherichia coli genome for RNA sequences with high affinity to RNA Pol by performing genomic SELEX. This approach led to the identification of RNA polymerase-binding APtamers termed "RAPs". RAPs can have positive and negative effects on gene expression. A subgroup is able to downregulate transcription via the activity of the termination factor Rho. In this study, we used a similar SELEX setup using yeast genomic DNA as source of RNA sequences and highly purified yeast RNA Pol II as bait and obtained almost 1300 yeast-derived RAPs. Yeast RAPs are found throughout the genome within genes and antisense to genes, they are overrepresented in the non-transcribed strand of yeast telomeres and underrepresented in intergenic regions. Genes harbouring a RAP are more likely to show lower mRNA levels. By determining the endogenous expression levels as well as using a reporter system, we show that RAPs located within coding regions can reduce the transcript level downstream of the RAP. Here we demonstrate that RAPs represent a novel type of regulatory RNA signal in Saccharomyces cerevisiae that act in cis and interfere with the elongating transcription machinery to reduce the transcriptional output.

Finding aptamers and small ribozymes in unexpected places.

The discovery of the catalytic properties of RNAs was a milestone for our view of how life emerged and forced us to reformulate many of our dogmas. The urge to grasp the whole spectrum of potential activities of RNA molecules stimulated two decades of fervent research resulting in a deep understanding of RNA-based phenomena. Most ribozymes were discovered by serendipity during the analysis of chemical processes, whereas RNA aptamers were identified through meticulous design and selection even before their discovery in nature. The desire to obtain aptamers led to the development of sophisticated technology and the design of efficient strategies. With the new notion that transcriptomes cover a major part of genomes and determine the identity of cells, it is reasonable to speculate that many more aptamers and ribozymes are awaiting their discovery in unexpected places. Now, in the genomic era with the development of powerful bioinformatics and sequencing methods, we are overwhelmed with tools for studying the genomes of all living and possibly even extinct organisms. Genomic SELEX (systematic evolution of ligands by exponential enrichment) coupled with deep sequencing and sophisticated computational analysis not only gives access to unexplored parts of sequenced genomes but also allows screening metagenomes in an unbiased manner.

Monitoring genomic sequences during SELEX using high-throughput sequencing: neutral SELEX.

BACKGROUND: SELEX is a well established in vitro selection tool to analyze the structure of ligand-binding nucleic acid sequences called aptamers. Genomic SELEX transforms SELEX into a tool to discover novel, genomically encoded RNA or DNA sequences binding a ligand of interest, called genomic aptamers. Concerns have been raised regarding requirements imposed on RNA sequences undergoing SELEX selection. METHODOLOGY/PRINCIPAL FINDINGS: To evaluate SELEX and assess the extent of these effects, we designed and performed a Neutral SELEX experiment omitting the selection step, such that the sequences are under the sole selective pressure of SELEX's amplification steps. Using high-throughput sequencing, we obtained thousands of full-length sequences from the initial genomic library and the pools after each of the 10 rounds of Neutral SELEX. We compared these to sequences obtained from a Genomic SELEX experiment deriving from the same initial library, but screening for RNAs binding with high affinity to the E. coli regulator protein Hfq. With each round of Neutral SELEX, sequences became less stable and changed in nucleotide content, but no sequences were enriched. In contrast, we detected substantial enrichment in the Hfq-selected set with enriched sequences having structural stability similar to the neutral sequences but with significantly different nucleotide selection. CONCLUSIONS/SIGNIFICANCE: Our data indicate that positive selection in SELEX acts independently of the neutral selective requirements imposed on the sequences. We conclude that Genomic SELEX, when combined with high-throughput sequencing of positively and neutrally selected pools, as well as the gnomic library, is a powerful method to identify genomic aptamers.