Addiction-Associated Genetic Variants Implicate Brain Cell Type- and Region-Specific Cis-Regulatory Elements in Addiction Neurobiology.

Recent large genome-wide association studies have identified multiple confident risk loci linked to addiction-associated behavioral traits. Most genetic variants linked to addiction-associated traits lie in noncoding regions of the genome, likely disrupting cis-regulatory element (CRE) function. CREs tend to be highly cell type-specific and may contribute to the functional development of the neural circuits underlying addiction. Yet, a systematic approach for predicting the impact of risk variants on the CREs of specific cell populations is lacking. To dissect the cell types and brain regions underlying addiction-associated traits, we applied stratified linkage disequilibrium score regression to compare genome-wide association studies to genomic regions collected from human and mouse assays for open chromatin, which is associated with CRE activity. We found enrichment of addiction-associated variants in putative CREs marked by open chromatin in neuronal (NeuN+) nuclei collected from multiple prefrontal cortical areas and striatal regions known to play major roles in reward and addiction. To further dissect the cell type-specific basis of addiction-associated traits, we also identified enrichments in human orthologs of open chromatin regions of female and male mouse neuronal subtypes: cortical excitatory, D1, D2, and PV. Last, we developed machine learning models to predict mouse cell type-specific open chromatin, enabling us to further categorize human NeuN+ open chromatin regions into cortical excitatory or striatal D1 and D2 neurons and predict the functional impact of addiction-associated genetic variants. Our results suggest that different neuronal subtypes within the reward system play distinct roles in the variety of traits that contribute to addiction.SIGNIFICANCE STATEMENT We combine statistical genetic and machine learning techniques to find that the predisposition to for nicotine, alcohol, and cannabis use behaviors can be partially explained by genetic variants in conserved regulatory elements within specific brain regions and neuronal subtypes of the reward system. Our computational framework can flexibly integrate open chromatin data across species to screen for putative causal variants in a cell type- and tissue-specific manner for numerous complex traits.

A cis-regulatory atlas in maize at single-cell resolution.

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Tumor Susceptibility Gene 101 Regulates the Glucocorticoid Receptor through Disorder-Mediated Allostery.

Tumor susceptibility gene 101 (TSG101) is involved in endosomal maturation and has been implicated in the transcriptional regulation of several steroid hormone receptors, although a detailed characterization of such regulation has yet to be conducted. Here we directly measure binding of TSG101 to one steroid hormone receptor, the glucocorticoid receptor (GR). Using biophysical and cellular assays, we show that the coiled-coil domain of TSG101 (1) binds and folds the disordered N-terminal domain of the GR, (2) upon binding improves the DNA binding of the GR in vitro, and (3) enhances the transcriptional activity of the GR in vivo. Our findings suggest that TSG101 is a bona fide transcriptional co-regulator of the GR and reveal how the underlying thermodynamics affect the function of the GR.

The LysR-Type Transcriptional Regulator CrgA Negatively Regulates the Flagellar Master Regulator flhDC in Ralstonia solanacearum GMI1000.

The invasion and colonization of host plants by the destructive pathogen Ralstonia solanacearum rely on its cell motility, which is controlled by multiple factors. Here, we report that the LysR-type transcriptional regulator CrgA (RS_RS16695) represses cell motility in R. solanacearum GMI1000. CrgA possesses common features of a LysR-type transcriptional regulator and contains an N-terminal helix-turn-helix motif as well as a C-terminal LysR substrate-binding domain. Deletion of crgA results in an enhanced swim ring and increased transcription of flhDC In addition, the ΔcrgA mutant possesses more polar flagella than wild-type GMI1000 and exhibits higher expression of the flagellin gene fliC Despite these alterations, the ΔcrgA mutant did not have a detectable growth defect in culture. Yeast one-hybrid and electrophoretic mobility shift assays revealed that CrgA interacts directly with the flhDC promoter. Expressing the ß-glucuronidase (GUS) reporter under the control of the crgA promoter showed that crgA transcription is dependent on cell density. Soil-soaking inoculation with the crgA mutant caused wilt symptoms on tomato (Solanum lycopersicum L. cv. Hong yangli) plants earlier than inoculation with the wild-type GMI1000 but resulted in lower disease severity. We conclude that the R. solanacearum regulator CrgA represses flhDC expression and consequently affects the expression of fliC to modulate cell motility, thereby conditioning disease development in host plants.IMPORTANCERalstonia solanacearum is a widely distributed soilborne plant pathogen that causes bacterial wilt disease on diverse plant species. Motility is a critical virulence attribute of R. solanacearum because it allows this pathogen to efficiently invade and colonize host plants. In R. solanacearum, motility-defective strains are markedly affected in pathogenicity, which is coregulated with multiple virulence factors. In this study, we identified a new LysR-type transcriptional regulator (LTTR), CrgA, that negatively regulates motility. The mutation of the corresponding gene leads to the precocious appearance of wilt symptoms on tomato plants when the pathogen is introduced using soil-soaking inoculation. This study indicates that the regulation of R. solanacearum motility is more complex than previously thought and enhances our understanding of flagellum regulation in R. solanacearum.

Candida Biofilms: An Update on Developmental Mechanisms and Therapeutic Challenges.

Fungi of the genus Candida are important etiological agents of superficial and life-threatening infections in individuals with a compromised immune system. One of the main characteristics of Candida is its ability to form highly drug tolerance biofilms in the human host. Biofilms are a dynamic community of multiple cell types whose formation over time is orchestrated by a network of transcription regulators. In this brief review, we provide an update of the processes involved in biofilm formation by Candida spp. (formation, treatment, and control), as well as the transcriptional circuitry that regulates its development and interactions with other microorganisms. Candida albicans is known to build mixed species biofilms with other Candida species and with various other bacterial species in different host niches. Taken together, these properties play a key role in Candida pathogenesis. In addition, this review gathers recent studies with new insights and perspectives for the treatment and control of Candida biofilms.

Low affinity binding sites in an activating CRM mediate negative autoregulation of the Drosophila Hox gene Ultrabithorax.

Specification of cell identity and the proper functioning of a mature cell depend on precise regulation of gene expression. Both binary ON/OFF regulation of transcription, as well as more fine-tuned control of transcription levels in the ON state, are required to define cell types. The Drosophila melanogaster Hox gene, Ultrabithorax (Ubx), exhibits both of these modes of control during development. While ON/OFF regulation is needed to specify the fate of the developing wing (Ubx OFF) and haltere (Ubx ON), the levels of Ubx within the haltere differ between compartments along the proximal-distal axis. Here, we identify and molecularly dissect the novel contribution of a previously identified Ubx cis-regulatory module (CRM), anterobithorax (abx), to a negative auto-regulatory loop that decreases Ubx expression in the proximal compartment of the haltere as compared to the distal compartment. We find that Ubx, in complex with the known Hox cofactors, Homothorax (Hth) and Extradenticle (Exd), acts through low-affinity Ubx-Exd binding sites to reduce the levels of Ubx transcription in the proximal compartment. Importantly, we also reveal that Ubx-Exd-binding site mutations sufficient to result in de-repression of abx activity in a transgenic context are not sufficient to de-repress Ubx expression when mutated at the endogenous locus, suggesting the presence of multiple mechanisms through which Ubx-mediated repression occurs. Our results underscore the complementary nature of CRM analysis through transgenic reporter assays and genome modification of the endogenous locus; but, they also highlight the increasing need to understand gene regulation within the native context to capture the potential input of multiple genomic elements on gene control.

Regional Gene Repression by DNA Double-Strand Breaks in G1 Phase Cells.

DNA damage responses (DDR) to double-strand breaks (DSBs) alter cellular transcription programs at the genome-wide level. Through processes that are less well understood, DSBs also alter transcriptional responses locally, which may be important for efficient DSB repair. Here, we developed an approach to elucidate the cis-acting responses to DSBs in G1 phase cells. We found that DSBs within a gene body silence its expression, as well as the transcription of local undamaged genes at a distance defined by the spread of γ-H2AX from the DSB. Importantly, DSBs not only repress ongoing transcription but also block the inducible expression of regional genes. DSB-mediated transcriptional repression depends on DDR signaling but does not require the generation of inaccessible chromatin. Our findings demonstrate that in G1 phase cells, DDR signaling establishes a robust and extensive region of transcriptional repression spreading from DSB sites and introduce an approach to study the mechanistic impact of targeted DNA breaks in nearly any chromatin environment.

In vivo Hox binding specificity revealed by systematic changes to a single cis regulatory module.

Hox proteins belong to a family of transcription factors with similar DNA binding specificities that control animal differentiation along the antero-posterior body axis. Hox proteins are expressed in partially overlapping regions where each one is responsible for the formation of particular organs and structures through the regulation of specific direct downstream targets. Thus, explaining how each Hox protein can selectively control its direct targets from those of another Hox protein is fundamental to understand animal development. Here we analyse a cis regulatory module directly regulated by seven different Drosophila Hox proteins and uncover how different Hox class proteins differentially control its expression. We find that regulation by one or another Hox protein depends on the combination of three modes: Hox-cofactor dependent DNA-binding specificity; Hox-monomer binding sites; and interaction with positive and negative Hox-collaborator proteins. Additionally, we find that similar regulation can be achieved by Amphioxus orthologs, suggesting these three mechanisms are conserved from insects to chordates.

Dissection of acute stimulus-inducible nucleosome remodeling in mammalian cells.

Accessibility of the genomic regulatory information is largely controlled by the nucleosome-organizing activity of transcription factors (TFs). While stimulus-induced TFs bind to genomic regions that are maintained accessible by lineage-determining TFs, they also increase accessibility of thousands of cis-regulatory elements. Nucleosome remodeling events underlying such changes and their interplay with basal positioning are unknown. Here, we devised a novel quantitative framework discriminating different types of nucleosome remodeling events in micrococcal nuclease ChIP-seq (chromatin immunoprecipitation [ChIP] combined with high-throughput sequencing) data sets and used it to analyze nucleosome dynamics at stimulus-regulated cis-regulatory elements. At enhancers, remodeling preferentially affected poorly positioned nucleosomes while sparing well-positioned nucleosomes flanking the enhancer core, indicating that inducible TFs do not suffice to overrule basal nucleosomal organization maintained by lineage-determining TFs. Remodeling events appeared to be combinatorially driven by multiple TFs, with distinct TFs showing, however, different remodeling efficiencies. Overall, these data provide a systematic view of the impact of stimulation on nucleosome organization and genome accessibility in mammalian cells.

Distinct Genome Replication and Transcription Strategies within the Growing Filovirus Family.

Research on filoviruses has historically focused on the highly pathogenic ebola- and marburgviruses. Indeed, until recently, these were the only two genera in the filovirus family. Recent advances in sequencing technologies have facilitated the discovery of not only a new ebolavirus, but also three new filovirus genera and a sixth proposed genus. While two of these new genera are similar to the ebola- and marburgviruses, the other two, discovered in saltwater fishes, are considerably more diverse. Nonetheless, these viruses retain a number of key features of the other filoviruses. Here, we review the key characteristics of filovirus replication and transcription, highlighting similarities and differences between the viruses. In particular, we focus on key regulatory elements in the genomes, replication and transcription strategies, and the conservation of protein domains and functions among the viruses. In addition, using computational analyses, we were able to identify potential homology and functions for some of the genes of the novel filoviruses with previously unknown functions. Although none of the newly discovered filoviruses have yet been isolated, initial studies of some of these viruses using minigenome systems have yielded insights into their mechanisms of replication and transcription. In general, the Cuevavirus and proposed Dianlovirus genera appear to follow the transcription and replication strategies employed by the ebola- and marburgviruses, respectively. While our knowledge of the fish filoviruses is currently limited to sequence analysis, the lack of certain conserved motifs and even entire genes necessitates that they have evolved distinct mechanisms of replication and transcription.

Epigenetic memory in gene regulation and immune response.

Cells must fine-tune their gene expression programs for optimal cellular activities in their natural growth conditions. Transcriptional memory, a unique transcriptional response, plays a pivotal role in faster reactivation of genes upon environmental changes, and is facilitated if genes were previously in an active state. Hyper-activation of gene expression by transcriptional memory is critical for cellular differentiation, development, and adaptation. TREM (Transcriptional REpression Memory), a distinct type of transcriptional memory, promoting hyper-repression of unnecessary genes, upon environmental changes has been recently reported. These two transcriptional responses may optimize specific gene expression patterns, in rapidly changing environments. Emerging evidence suggests that they are also critical for immune responses. In addition to memory B and T cells, innate immune cells are transcriptionally hyperactivated by restimulation, with the same or different pathogens known as trained immunity. In this review, we briefly summarize recent progress in chromatin-based regulation of transcriptional memory, and its potential role in immune responses. [BMB Reports 2019; 52(2): 127-132].

High-throughput sequencing identifies aetiology-dependent differences in ductular reaction in human chronic liver disease.

Ductular reaction (DR) represents the activation of hepatic progenitor cells (HPCs) and has been associated with features of advanced chronic liver disease; yet it is not clear whether these cells contribute to disease progression and how the composition of their micro-environment differs depending on the aetiology. This study aimed to identify HPC-associated signalling pathways relevant in different chronic liver diseases using a high-throughput sequencing approach. DR/HPCs were isolated using laser microdissection from patient samples diagnosed with HCV or primary sclerosing cholangitis (PSC), as models for hepatocellular or biliary regeneration. Key signals were validated at the protein level for a cohort of 56 patients (20 early and 36 advanced stage). In total, 330 genes were significantly differentially expressed between the HPCs in HCV and PSC. Recruitment and homing of inflammatory cells were distinctly different depending on the aetiology. HPCs in PSC were characterised by a response to oxidative stress (e.g. JUN, VNN1) and neutrophil-attractant chemokines (CXCL5, CXCL6, IL-8), whereas HPCs in HCV were identified by T- and B-lymphocyte infiltration. Moreover, we found that communication between HPCs and macrophages was aetiology driven. In PSC, a high frequency of CCL28-positive macrophages was observed in the portal infiltrate, already in early disease in the absence of advanced fibrosis, while in HCV, HPCs showed a strong expression of the macrophage scavenger receptor MARCO. Interestingly, DR/HPCs in PSC showed more deposition of ECM (e.g. FN1, LAMC2, collagens) compared to HCV, where an increase of pro-invasive genes (e.g. PDGFRA, IGF2) was observed. Additionally, endothelial cells in the vicinity of DR/HPCs showed differential immunopositivity (e.g. IGF2 and INHBA expression). In conclusion, our data shine light on the role of DR/HPCs in immune signalling, fibrogenesis and angiogenesis in chronic liver disease. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

A Developmental Program Truncates Long Transcripts to Temporally Regulate Cell Signaling.

Rapid mitotic divisions and a fixed transcription rate limit the maximal length of transcripts in early Drosophila embryos. Previous studies suggested that transcription of long genes is initiated but aborted, as early nuclear divisions have short interphases. Here, we identify long genes that are expressed during short nuclear cycles as truncated transcripts. The RNA binding protein Sex-lethal physically associates with transcripts for these genes and is required to support early termination to specify shorter transcript isoforms in early embryos of both sexes. In addition, one truncated transcript for the gene short-gastrulation encodes a product in embryos that functionally relates to a previously characterized dominant-negative form, which maintains TGF-ß signaling in the off-state. In summary, our results reveal a developmental program of short transcripts functioning to help temporally regulate Drosophila embryonic development, keeping cell signaling at early stages to a minimum in order to support its proper initiation at cellularization.

CosR is an oxidative stress sensing a MarR-type transcriptional repressor in Corynebacterium glutamicum.

The MarR family is unique to both bacteria and archaea. The members of this family, one of the most prevalent families of transcriptional regulators in bacteria, enable bacteria to adapt to changing environmental conditions, such as the presence of antibiotics, toxic chemicals, or reactive oxygen species (ROS), mainly by thiol-disulfide switches. Although the genome of Corynebacterium glutamicum encodes a large number of the putative MarR-type transcriptional regulators, their physiological and biochemical functions have so far been limited to only two proteins, regulator of oxidative stress response RosR and quinone oxidoreductase regulator QosR. Here, we report that the ncgl2617 gene (cosR) of C. glutamicum encoding an MarR-type transcriptional regulator plays an important role in oxidative stress resistance. The cosR null mutant is found to be more resistant to various oxidants and antibiotics, accompanied by a decrease in ROS production and protein carbonylation levels under various stresses. Protein biochemical function analysis shows that two Cys residues presenting at 49 and 62 sites in CosR are redox-active. They form intermolecular disulfide bonds in CosR under oxidative stress. This CosR oxidation leads to its dissociation from promoter DNA, depression of the target DNA, and increased oxidative stress resistance of C. glutamicum. Together, the results reveal that CosR is a redox-sensitive regulator that senses peroxide stress to mediate oxidative stress resistance in C. glutamicum.

Genomic analysis of transcriptional networks directing progression of cell states during MGE development.

BACKGROUND: Homeodomain (HD) transcription factor (TF) NKX2-1 critical for the regional specification of the medial ganglionic eminence (MGE) as well as promoting the GABAergic and cholinergic neuron fates via the induction of TFs such as LHX6 and LHX8. NKX2-1 defines MGE regional identity in large part through transcriptional repression, while specification and maturation of GABAergic and cholinergic fates is mediated in part by transcriptional activation via TFs such as LHX6 and LHX8. Here we analyze the signaling and TF pathways, downstream of NKX2-1, required for GABAergic and cholinergic neuron fate maturation. METHODS: Differential ChIP-seq analysis was used to identify regulatory elements (REs) where chromatin state was sensitive to change in the Nkx2-1cKO MGE at embryonic day (E) 13.5. TF motifs in the REs were identified using RSAT. CRISPR-mediated genome editing was used to generate enhancer knockouts. Differential gene expression in these knockouts was analyzed through RT-qPCR and in situ hybridization. Functional analysis of motifs within hs623 was analyzed via site directed mutagenesis and reporter assays in primary MGE cultures. RESULTS: We identified 4782 activating REs (aREs) and 6391 repressing REs (rREs) in the Nkx2-1 conditional knockout (Nkx2-1cKO) MGE. aREs are associated with basic-Helix-Loop-Helix (bHLH) TFs. Deletion of hs623, an intragenic Tcf12 aRE, caused a reduction of Tcf12 expression in the sub-ventricular zone (SVZ) and mantle zone (MZ) of the MGE. Mutation of LHX, SOX and octamers, within hs623, caused a reduction of hs623 activity in MGE primary cultures. CONCLUSIONS: Tcf12 expression in the SVZ of the MGE is mediated through aRE hs623. The activity of hs623 is dependent on LHX6, SOX and octamers. Thus, maintaining the expression of Tcf12 in the SVZ involves on TF pathways parallel and genetically downstream of NKX2-1.

Lateral inhibition: Two modes of non-autonomous negative autoregulation by neuralized.

Developmental patterning involves the progressive subdivision of tissue into different cell types by invoking different genetic programs. In particular, cell-cell signaling is a universally deployed means of specifying distinct cell fates in adjacent cells. For this mechanism to be effective, it is essential that an asymmetry be established in the signaling and responding capacities of the participating cells. Here we focus on the regulatory mechanisms underlying the role of the neuralized gene and its protein product in establishing and maintaining asymmetry of signaling through the Notch pathway. The context is the classical process of "lateral inhibition" within Drosophila proneural clusters, which is responsible for distinguishing the sensory organ precursor (SOP) and non-SOP fates among adjacent cells. We find that neur is directly regulated in proneural clusters by both proneural transcriptional activators and Enhancer of split basic helix-loop-helix repressors (bHLH-Rs), via two separate cis-regulatory modules within the neur locus. We show that this bHLH-R regulation is required to prevent the early, pre-SOP expression of neur from being maintained in a subset of non-SOPs following SOP specification. Lastly, we demonstrate that Neur activity in the SOP is required to inhibit, in a cell non-autonomous manner, both neur expression and Neur function in non-SOPs, thus helping to secure the robust establishment of distinct cell identities within the developing proneural cluster.

Metagenomic mining of regulatory elements enables programmable species-selective gene expression.

Robust and predictably performing synthetic circuits rely on the use of well-characterized regulatory parts across different genetic backgrounds and environmental contexts. Here we report the large-scale metagenomic mining of thousands of natural 5' regulatory sequences from diverse bacteria, and their multiplexed gene expression characterization in industrially relevant microbes. We identified sequences with broad and host-specific expression properties that are robust in various growth conditions. We also observed substantial differences between species in terms of their capacity to utilize exogenous regulatory sequences. Finally, we demonstrate programmable species-selective gene expression that produces distinct and diverse output patterns in different microbes. Together, these findings provide a rich resource of characterized natural regulatory sequences and a framework that can be used to engineer synthetic gene circuits with unique and tunable cross-species functionality and properties, and also suggest the prospect of ultimately engineering complex behaviors at the community level.

Crystal structure of an inactive variant of the quorum-sensing master regulator HapR from the protease-deficient non-O1, non-O139 Vibrio cholerae strain V2.

HapR is a TetR-family transcriptional regulator that controls quorum sensing in Vibrio cholerae, the causative agent of cholera. HapR regulates the expression of hemagglutinin protease, virulence and biofilm genes. The crystal structure of wild-type HapR from V. cholerae strain O1 El Tor C6706 has previously been solved. In this study, the structure of a DNA-binding-deficient variant of HapR (HapRV2) derived from the protease-deficient V. cholerae serotype O37 strain V2 is reported. The structure reveals no structural differences compared with wild-type HapR. However, structural alignment of HapRV2 with the TetR-family member QacR in complex with its operator DNA suggests that the aspartate residue located between the regulatory and DNA-binding domains may clash with and electrostatically repel the phosphate backbone of DNA to prevent binding.

Concomitant Regulation by a LacI-Type Transcriptional Repressor XylR on Genes Involved in Xylan and Xylose Metabolism and the Type III Secretion System in Rice Pathogen Xanthomonas oryzae pv. oryzae.

The hypersensitive response and pathogenicity (hrp) genes of Xanthomonas oryzae pv. oryzae, the causal agent of bacterial leaf blight of rice, encode components of the type III secretion system and are essential for virulence. Expression of hrp genes is regulated by two key hrp regulators, HrpG and HrpX; HrpG regulates hrpX and hrpA, and HrpX regulates the other hrp genes on hrpB-hrpF operons. We previously reported the sugar-dependent quantitative regulation of HrpX; the regulator highly accumulates in the presence of xylose, followed by high hrp gene expression. Here, we found that, in a mutant lacking the LacI-type transcriptional regulator XylR, HrpX accumulation and hrp gene expression were high even in the medium without xylose, reaching the similar levels present in the wild type incubated in the xylose-containing medium. XylR also negatively regulated one of two xylose isomerase genes (xylA2 but not xylA1) by binding to the motif sequence in the upstream region of the gene. Xylose isomerase is an essential enzyme in xylose metabolism and interconverts between xylose and xylulose. Our results suggest that, in the presence of xylose, inactivation of XylR leads to greater xylan and xylose utilization and, simultaneously, to higher accumulation of HrpX, followed by higher hrp gene expression in the bacterium.

The molecular dynamics of long noncoding RNA control of transcription in PTEN and its pseudogene.

RNA has been found to interact with chromatin and modulate gene transcription. In human cells, little is known about how long noncoding RNAs (lncRNAs) interact with target loci in the context of chromatin. We find here, using the phosphatase and tensin homolog (PTEN) pseudogene as a model system, that antisense lncRNAs interact first with a 5' UTR-containing promoter-spanning transcript, which is then followed by the recruitment of DNA methyltransferase 3a (DNMT3a), ultimately resulting in the transcriptional and epigenetic control of gene expression. Moreover, we find that the lncRNA and promoter-spanning transcript interaction are based on a combination of structural and sequence components of the antisense lncRNA. These observations suggest, on the basis of this one example, that evolutionary pressures may be placed on RNA structure more so than sequence conservation. Collectively, the observations presented here suggest a much more complex and vibrant RNA regulatory world may be operative in the regulation of gene expression.