Search

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.

Supplementary Figure S5 shows pharmacodynamic changes in biomarkers over time (A) and changes by best overall response (B) in all study 114 cohorts.

Supplementary Table S3 provides a summary of infusion-related reactions during cycle 1, per premedication.

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.

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. Among the so far uncharacterized BMP-dependent transcription factors and signaling molecules we identified several, whose bilaterally symmetric expression may be indicative of their involvement in secondary axis patterning. One of these molecules, zswim4-6, encodes a novel nuclear modulator of the pSMAD1/5 gradient potentially promoting BMP-dependent gene repression. Strikingly, overexpression of the zebrafish homologue zswim5 suggests that its effect on the pSMAD1/5 gradient is conserved between anthozoan Cnidaria and Bilateria.

Draft genome sequences of non-bilaterian species have provided important insights into the evolution of the metazoan gene repertoire. However, there is little information about the evolution of gene clusters, genome architectures and karyotypes during animal evolution. Here we report chromosome-level genome assemblies of two related anthozoan cnidarians, the sea anemones, Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence of the chromosomes between the two species. We show that, in contrast to Bilateria, Hox and NK clusters of investigated cnidarians are disintegrated, indicating that microsynteny conservation is largely lost. In line with that, we find no evidence for topologically associated domains, suggesting fundamental difference in long-range gene regulation compared to vertebrates. However, both sea anemone genomes show remarkable chromosomal conservation with other cnidarians, several bilaterians and the sponge Ephydatia muelleri, allowing us to reconstruct the putative cnidarian and metazoan chromosomes, consisting of 19 and 16 ancestral linkage groups, respectively. These data suggest that large parts of the ancestral metazoan genome have been retained in chromosomes of some extant lineages, yet, higher order gene regulation may have evolved only after the cnidarian-bilaterian split.

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 ChIP-seq 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 signaling 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 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.

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.

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.

Androgen receptor signaling inhibitors (ARSi) and taxanes are mainstays for patients with metastatic castration-resistant prostate cancer (mCRPC). However, patient response is heterogeneous, and the molecular underpinnings of treatment resistance are not well elucidated. To identify clinically meaningful mechanisms of treatment resistance, we performed transcriptome analysis of circulating tumor cells (CTCs) isolated from mCRPC patients enrolled in two independent prospective clinical trials: PROPHECY, a clinical study of patients (n=118) treated with abiraterone or enzalutamide followed by docetaxel; and TAXYNERGY where patients were randomized to docetaxel or cabazitaxel treatment. CTCs were obtained at baseline (before treatment), on treatment and at progression and their comprehensive transcriptomic analysis was correlated with clinical outcomes. To uncover potential involvement of the circulating immune macroenvironment (CIME) in treatment resistance, we performed transcriptomic analysis of matching peripheral blood mononuclear cells (PBMCs) using an established, rigorous, blood-derived transcriptional modular framework. In PROPHECY, CTC RNA-seq identified that RB loss concurrently with enhanced E2F signaling networks were associated with intrinsic ARSi resistance. Using single sample GSEA (ssGSEA) score, we identified that the RB/E2F common signature at baseline was associated with short PFS (median PFS=6.5 months) and OS (median OS=24.5 months) (hazard ratio (HR) = 3.5; 95% CI 1.5-8.2) in men with mCRPC. We further developed a BRCA-loss transcriptional signature, and validated it in the SU2C mCRPC patient cohort, expanding the identification of patients with BRCA-loss phenotypes beyond genomic loss. Applying this signature to PROPHECY baseline samples, we showed that men with high BRCA-loss scores experienced shorter OS (HR=2.42; 95% CI=1-5.9). Through the comparison of CTC transcriptomic profiles at progression with baseline, we identified an inflammatory response signature in CTCs which was significantly associated with acquired ARSi resistance. Transcriptomic PBMC analysis further identified enrichment of inflammasome gene signatures at progression, with concurrent downregulation of CD8+ T and NK cells. Furthermore, preliminary data from both clinical trials, showed a significant upregulation of TGF-β1 and corresponding TGFβ-Receptor signaling pathway in CTCs from patients at progression following taxane treatment, suggesting a role for TGFβ pathway in clinical response to taxane chemotherapy. Taken together, these data demonstrate that liquid biopsy transcriptomics of both tumor cells and immune cells can identify molecular pathways associated with treatment resistance paving the way for treatment optimization and the development of novel precision therapies in patients with mCRPC. Citation Format: Jiaren Zhang, Bob Zimmermann, Giuseppe Galletti, Susan Halabi, Ada Gjyrezi, Qian Yang, Santosh Gupta, Akanksha Verma, Andrea Sboner, Monika Anand, Daniel J. George, Simon G. Gregory, Prerna Mahtani, Seunghee Hong, Virginia Pascual, Clio P. Mavragani, Emmanuel S. Antonarakis, David M. Nanus, Scott T. Tagawa, Olivier Elemento, Andrew J. Armstrong, Paraskevi Giannakakou. Liquid biopsy transcriptomics identify pathways associated with poor outcomes and immune phenotypes in men with mCRPC [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 646.

139 Background: Androgen receptor signaling inhibitors (ARSi) are a mainstay for patients with metastatic castration-resistant prostate cancer (mCRPC). However, patient response is heterogeneous and the molecular underpinnings of ARSi resistance are not well elucidated. Methods: We performed transcriptome analysis of circulating tumor cells (CTCs) and peripheral blood mononuclear cells (PBMC) in the context of PROPHECY, a prospective clinical trial of men (n = 118) with mCRPC treated with abiraterone (Abi) or enzalutamide (Enza). We obtained CTCs at baseline (before treatment) and at the time of progression on Abi/Enza, performed a comprehensive transcriptomic analysis of CTC patient samples (n = 40) and correlated with clinical outcomes to identify mechanisms of ARSi resistance. In addition, we also performed a transcriptomic analysis of matching peripheral blood mononuclear cells (PBMCs) in order to uncover potential involvement of the circulating immune macroenvironment (CIME) in ARSi resistance. The proportional hazard model was used to determine the prognostic significance of these signatures in predicting overall survival (OS) and progression-free survival (PFS). Results: CTC RNA-sequencing identified that RB loss concurrently with enhanced E2F signaling transcriptional networks were associated with intrinsic ARSi resistance. Using single sample GSEA (ssGSEA) score, we identified that the RB/E2F common signature at baseline was associated with short PFS (median PFS = 6.5 months) and OS (median OS = 24.5 months) (hazard ratio (HR) = 3.5; 95% CI 1.5-8.2) in men with mCRPC. We further developed a BRCA loss transcriptional signature which we validated in the SU2C mCRPC patient cohort, by showing that BRCA loss transcriptional network reflected BRCA genomic alterations as it was significantly enriched in the SU2C BRCA-altered patients vs unaltered patients. Generating BRCA loss ssGSEA scores in the PROPHECY cohort we observed that patients with high BRCA loss scores at baseline experienced shorter OS (HR = 2.42; 95% CI = 1-5.9). Through the comparison of CTC transcriptomic profiles at progression with baseline, we identified an inflammatory response signature in CTCs which was significantly associated with acquired ARSi resistance. Transcriptomic analysis of matching PBMCs identified enrichment of inflammasome gene signatures indicative of activated innate immunity at progression, with concurrent downregulation of CD8 T and NK cells. Importantly, CTC gene signatures of RB loss/E2F signaling had a significant positive association with this CIME signatures. Conclusions: Taken together, these data demonstrate that liquid biopsy transcriptomics of both tumor cells and immune cells can identify molecular pathways associated with clinical ARSi resistance paving the way for treatment optimization and the development of novel precision therapies in patients with mCRPC.

Cnidaria (sea anemones, jellyfish, corals and hydra) form a close sister group to Bilateria. Within this clade, the sea anemone Nematostella vectensis has emerged as a slow evolving model for investigating characteristics of the cnidarian-bilaterian common ancestor, which diverged near the Cambrian explosion. Here, using long read sequencing and high throughput chromosome conformation capture, we generate high quality chromosome-level genome assemblies for N. vectensis and the closely related edwardsiid sea anemone, Scolanthus callimorphus. In both cases we find a robust set of 15 chromosomes comprising a stable linkage group detectable within all major clades of sequenced cnidarian genomes. Further, both genomes show remarkable chromosomal conservation with chordates. In contrast with Bilateria, we report that extended Hox and NK gene clusters are chromosomally linked but do not retain a tight spatial conservation. Accordingly, there is a lack of evidence for topologically associated domains, which have been implicated in the evolutionary pressure to retain tight microsyntenic gene clusters. We also uncover ultra-conserved noncoding elements at levels previously undetected in non-chordate lineages. Both genomes are accessible through an actively updated genome browser and database at https://simrbase.stowers.org

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 investigated the oral-aboral axis patterning in the sea anemoneNematostella- 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 inNematostellaand deuterostomes suggests a common evolutionary origin of these processes.

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.

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.

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

How dispersal strategies impact the distribution of species and subsequent speciation events is a fundamental question in evolutionary biology. Sedentary benthic marine organisms, such as corals or sea anemones usually rely on motile larval stages for dispersal and therefore have a relatively restricted distribution along coasts. Edwardsiella lineata and Edwardsiella carnea are virtually indistinguishable edwardsiid sea anemones native to the east American and the Northern European coast, respectively. E. lineata is a facultative parasite to the ctenophore Mnemiopsis leidyi, while the life cycle of E. carnea is unknown. Recently M. leidyi was found in the Skagerrak carrying Edwardsiella sp. parasites, which raised the intriguing possibility that the invasive comb jellies acted as cargo for the facultative E. lineata parasites to establish a new population in Northern Europe. Here, we assessed the genetic differences between these two cryptic Edwardsiella species and isolated parasites from the invasive comb jelly M. leidyi in Sweden by comparing rRNA, whole transcriptomes, SNPs, ITS2 sequences and the gene complements of key developmental regulators, the Wnt gene family. We show that E. carnea and the parasite transcriptomes are more than 99% identical, hence demonstrating that E. carnea has a previously unknown parasitic life stage. ITS2 sequence analysis of E. carnea and E. lineata suggest that they may not be reproductively isolated. The transcriptomes of E. lineata and E. carnea are ∼97% identical. We also estimate that the species diverged between 18.7 and 21.6 million years ago.

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 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. WIREs RNA 2012, 3:73–91. doi: 10.1002/wrna.105 For further resources related to this article, please visit the WIREs website.

An unexpectedly high number of regulatory RNAs have been recently discovered that fine-tune the function of genes at all levels of expression. We employed Genomic SELEX, a method to identify protein-binding RNAs encoded in the genome, to search for further regulatory RNAs in Escherichia coli. We used the global regulator protein Hfq as bait, because it can interact with a large number of RNAs, promoting their interaction. The enriched SELEX pool was subjected to deep sequencing, and 8865 sequences were mapped to the E. coli genome. These short sequences represent genomic Hfq-aptamers and are part of potential regulatory elements within RNA molecules. The motif 5'-AAYAAYAA-3' was enriched in the selected RNAs and confers low-nanomolar affinity to Hfq. The motif was confirmed to bind Hfq by DMS footprinting. The Hfq aptamers are 4-fold more frequent on the antisense strand of protein coding genes than on the sense strand. They were enriched opposite to translation start sites or opposite to intervening sequences between ORFs in operons. These results expand the repertoire of Hfq targets and also suggest that Hfq might regulate the expression of a large number of genes via interaction with cis-antisense RNAs.

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.

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.

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.