Search

Type III CRISPR immune systems bind viral or plasmid RNA transcripts and activate Csm3/Cmr4 and Cas10 nucleases to uniquely cleave both invader RNA and DNA, respectively. Additionally, type III effector complexes generate cyclic oligoadenylate (cOA) signaling molecules to activate trans-acting, auxiliary Csm6/Csx1 ribonucleases, previously proposed to be non-specific in their in vivo RNA cleavage preference. Despite extensive in vitro studies, the nuclease requirements of type III systems in their native contexts remain poorly understood. Here we systematically investigated the in vivo roles for immunity of each of the three Streptococcus thermophilus (Sth) type III-A Cas nucleases and cOA signaling by challenging nuclease defective mutant strains with plasmid and phage infections. Our results reveal that RNA cleavage by Csm6 is both sufficient and essential for maintaining wild-type levels of immunity. Importantly, Csm6 RNase activity leads to immunity against even high levels of phage challenge without causing host cell dormancy or death. Transcriptomic analyses during phage infection indicated Csm6-mediated and crRNA-directed preferential cleavage of phage transcripts. Our findings highlight the critical role of Csm6 RNase activity in type III immunity and demonstrate specificity for invader RNA transcripts by Csm6 to ensure host cell survival upon phage infection., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

RNA interactome studies have revealed that hundreds of zinc-finger proteins (ZFPs) are candidate RNA-binding proteins (RBPs), yet their RNA substrates and functional significance remain largely uncharacterized. Here, we present a systematic multi-omics analysis of the DNA- and RNA-binding targets and regulatory roles of more than 100 ZFPs representing 37 zinc-finger families. We show that multiple ZFPs are previously unknown regulators of RNA splicing, alternative polyadenylation, stability, or translation. The examined ZFPs show widespread sequence-specific RNA binding and preferentially bind proximal to transcription start sites. Additionally, several ZFPs associate with their targets at both the DNA and RNA levels. We highlight ZNF277, a C2H2 ZFP that binds thousands of RNA targets and acts as a multi-functional RBP. We also show that ZNF473 is a DNA/RNA-associated protein that regulates the expression and splicing of cell cycle genes. Our results reveal diverse roles for ZFPs in transcriptional and post-transcriptional gene regulation., Competing Interests: Declaration of interests G.W.Y. is a scientific advisory board (SAB) member of Jumpcode Genomics and a co-founder, member of the board of directors, on the SAB, equity holder, and paid consultant for Eclipse BioInnovations. G.W.Y. is a Distinguished Visiting Professor at the National University of Singapore. G.W.Y.’s interests have been reviewed and approved by UCSD in accordance with its conflict-of-interest policies. B.R.G. is a co-founder and SAB member for RNAConnect., (Copyright © 2024 University of California San Diego. Published by Elsevier Inc. All rights reserved.)

The prevalence of obesity is increasing in older adults and contributes to age-related decline. Caloric restriction (CR) alleviates obesity phenotypes and delays the onset of age-related changes. However, how late in life organisms benefit from switching from a high-(H) to a low-calorie (L) diet is unclear. We transferred male flies from a H to a L (HL) diet or vice versa (LH) at different times during life. Both shifts immediately change fly rate of aging even when applied late in life. HL shift rapidly reduces fly mortality rate to briefly lower rate than in flies on a constant L diet, and extends lifespan. Transcriptomic analysis uncovers that flies aged on H diet have acquired increased stress response, which may have temporal advantage over flies aged on L diet and leads to rapid decrease in mortality rate after HL switch. Conversely, a LH shift increases mortality rate, which is temporarily higher than in flies aged on a H diet, and shortens lifespan. Unexpectedly, more abundant transcriptomic changes accompanied LH shift, including increase in ribosome biogenesis, stress response and growth. These changes reflect protection from sudden release of ROS, energy storage, and use of energy to growth, which all likely contribute to higher mortality rate. As the beneficial effects of CR on physiology and lifespan are conserved across many organisms, our study provides framework to study underlying mechanisms of CR interventions that counteract the detrimental effects of H diets and reduce rate of aging even when initiated later in life., Competing Interests: Competing interests statement:The authors declare no competing interest.

CRISPR-Cas systems provide heritable immunity against viruses and other mobile genetic elements by incorporating fragments of invader DNA into the host CRISPR array as spacers. Integration of new spacers is localized to the 5' end of the array, and in certain Gram-negative Bacteria this polarized localization is accomplished by the integration host factor. For most other Bacteria and Archaea, the mechanism for 5' end localization is unknown. Here we show that archaeal histones play a key role in directing integration of CRISPR spacers. In Pyrococcus furiosus, deletion of either histone A or B impairs integration. In vitro, purified histones are sufficient to direct integration to the 5' end of the CRISPR array. Archaeal histone tetramers and bacterial integration host factor induce similar U-turn bends in bound DNA. These findings indicate a co-evolution of CRISPR arrays with chromosomal DNA binding proteins and a widespread role for binding and bending of DNA to facilitate accurate spacer integration., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

RNA binding proteins (RBPs) are central regulators of gene expression implicated in all facets of RNA metabolism. As such, they play key roles in cellular physiology and disease etiology. Since different steps of post-transcriptional gene expression tend to occur in specific regions of the cell, including nuclear or cytoplasmic locations, defining the subcellular distribution properties of RBPs is an important step in assessing their potential functions. Here, we present the RBP Image Database, a resource that details the subcellular localization features of 301 RBPs in the human HepG2 and HeLa cell lines, based on the results of systematic immuno-fluorescence studies conducted using a highly validated collection of RBP antibodies and a panel of 12 markers for specific organelles and subcellular structures. The unique features of the RBP Image Database include: (i) hosting of comprehensive representative images for each RBP-marker pair, with ∼250,000 microscopy images; (ii) a manually curated controlled vocabulary of annotation terms detailing the localization features of each factor; and (iii) a user-friendly interface allowing the rapid querying of the data by target or annotation. The RBP Image Database is freely available at https://rnabiology.ircm.qc.ca/RBPImage/., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

CRISPR-Cas systems are functionally diverse prokaryotic antiviral defense systems, which encompass six distinct types (I-VI) that each encode different effector Cas nucleases with distinct nucleic acid cleavage specificities. By harnessing the unique attributes of the various CRISPR-Cas systems, a range of innovative CRISPR-based DNA and RNA targeting tools and technologies have been developed. Here, we exploit the ability of type III-A CRISPR-Cas systems to carry out RNA-guided and sequence-specific target RNA cleavage for establishment of research tools for post-transcriptional control of gene expression. Type III-A systems from three bacterial species ( L. lactis , S. epidermidis , and S. thermophilus ) were each expressed on a single plasmid in E. coli , and the efficiency and specificity of gene knockdown was assessed by northern blot and transcriptomic analysis. We show that engineered type III-A modules can be programmed using tailored CRISPR RNAs to efficiently knock down gene expression of both coding and noncoding RNAs in vivo. Moreover, simultaneous degradation of multiple cellular mRNA transcripts can be directed by utilizing a CRISPR array expressing corresponding gene-targeting crRNAs. Our results demonstrate the utility of distinct type III-A modules to serve as specific and effective gene knockdown platforms in heterologous cells. This transcriptome engineering technology has the potential to be further refined and exploited for key applications including gene discovery and gene pathway analyses in additional prokaryotic and perhaps eukaryotic cells and organisms., (© 2022 Woodside et al.; Published by Cold Spring Harbor Laboratory Press for the RNA Society.)

Although reverse-transcriptase (RT) enzymes are critical reagents for research and biotechnology, their mechanical properties are not well understood. In particular, we know little about their relative speed and response to structural obstacles in the template. Commercial retroviral RTs stop at many positions along mixed sequence templates, resulting in truncated cDNA products that complicate downstream analysis. By contrast, group II intron-encoded RTs appear to copy long RNAs with high processivity and minimal stops. However, their speed, consistency and pausing behavior have not been explored. Here, we analyze RT velocity as the enzyme moves through heterogeneous sequences and structures that are embedded within a long noncoding RNA transcript. We observe that heterogeneities in the template are highly disruptive to primer extension by retroviral RTs. However, sequence composition and template structure have negligible effects on behavior of group II intron RTs, such as MarathonRT (MRT). Indeed, MRT copies long RNAs in a single pass, and displays synchronized primer extension at a constant speed of 25 nt/sec. In addition, it passes through stable RNA structural motifs without perturbation of velocity. Taken together, the results demonstrate that consistent, robust translocative behavior is a hallmark of group II intron-encoded RTs, some of which operate at high velocity., (© The Author(s) 2022. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Type III CRISPR-Cas systems have a unique mode of interference, involving crRNA-guided recognition of nascent RNA and leading to DNA and RNA degradation. How type III systems acquire new CRISPR spacers is currently not well understood. Here, we characterize CRISPR spacer uptake by a type III-A system within its native host, Streptococcus thermophilus. Adaptation by the type II-A system in the same host provided a basis for comparison. Cas1 and Cas2 proteins were critical for type III adaptation but deletion of genes responsible for crRNA biogenesis or interference did not detectably change spacer uptake patterns, except those related to host counter-selection. Unlike the type II-A system, type III spacers are acquired in a PAM- and orientation-independent manner. Interestingly, certain regions of plasmids and the host genome were particularly well-sampled during type III-A, but not type II-A, spacer uptake. These regions included the single-stranded origins of rolling-circle replicating plasmids, rRNA and tRNA encoding gene clusters, promoter regions of expressed genes and 5' UTR regions involved in transcription attenuation. These features share the potential to form DNA secondary structures, suggesting a preferred substrate for type III adaptation. Lastly, the type III-A system adapted to and protected host cells from lytic phage infection., (© The Author(s) 2021. Published by Oxford University Press on behalf of Nucleic Acids Research.)

The human and mouse genomes contain instructions that specify RNAs and proteins and govern the timing, magnitude, and cellular context of their production. To better delineate these elements, phase III of the Encyclopedia of DNA Elements (ENCODE) Project has expanded analysis of the cell and tissue repertoires of RNA transcription, chromatin structure and modification, DNA methylation, chromatin looping, and occupancy by transcription factors and RNA-binding proteins. Here we summarize these efforts, which have produced 5,992 new experimental datasets, including systematic determinations across mouse fetal development. All data are available through the ENCODE data portal (https://www.encodeproject.org), including phase II ENCODE 1 and Roadmap Epigenomics 2 data. We have developed a registry of 926,535 human and 339,815 mouse candidate cis-regulatory elements, covering 7.9 and 3.4% of their respective genomes, by integrating selected datatypes associated with gene regulation, and constructed a web-based server (SCREEN; http://screen.encodeproject.org) to provide flexible, user-defined access to this resource. Collectively, the ENCODE data and registry provide an expansive resource for the scientific community to build a better understanding of the organization and function of the human and mouse genomes.

The Encylopedia of DNA Elements (ENCODE) Project launched in 2003 with the long-term goal of developing a comprehensive map of functional elements in the human genome. These included genes, biochemical regions associated with gene regulation (for example, transcription factor binding sites, open chromatin, and histone marks) and transcript isoforms. The marks serve as sites for candidate cis-regulatory elements (cCREs) that may serve functional roles in regulating gene expression 1 . The project has been extended to model organisms, particularly the mouse. In the third phase of ENCODE, nearly a million and more than 300,000 cCRE annotations have been generated for human and mouse, respectively, and these have provided a valuable resource for the scientific community.

Many proteins regulate the expression of genes by binding to specific regions encoded in the genome 1 . Here we introduce a new data set of RNA elements in the human genome that are recognized by RNA-binding proteins (RBPs), generated as part of the Encyclopedia of DNA Elements (ENCODE) project phase III. This class of regulatory elements functions only when transcribed into RNA, as they serve as the binding sites for RBPs that control post-transcriptional processes such as splicing, cleavage and polyadenylation, and the editing, localization, stability and translation of mRNAs. We describe the mapping and characterization of RNA elements recognized by a large collection of human RBPs in K562 and HepG2 cells. Integrative analyses using five assays identify RBP binding sites on RNA and chromatin in vivo, the in vitro binding preferences of RBPs, the function of RBP binding sites and the subcellular localization of RBPs, producing 1,223 replicated data sets for 356 RBPs. We describe the spectrum of RBP binding throughout the transcriptome and the connections between these interactions and various aspects of RNA biology, including RNA stability, splicing regulation and RNA localization. These data expand the catalogue of functional elements encoded in the human genome by the addition of a large set of elements that function at the RNA level by interacting with RBPs.

CRISPR-Cas adaptive immune systems are used by prokaryotes to defend against invaders like viruses and other mobile genetic elements. Immune memories are stored in the form of 'spacers' which are short DNA sequences that are captured from invaders and added to the CRISPR array during a process called 'adaptation'. Spacers are transcribed and the resulting CRISPR (cr)RNAs assemble with different Cas proteins to form effector complexes that recognize matching nucleic acid and destroy it ('interference'). Adaptation can be 'naïve', i.e. independent of any existing spacer matches, or it can be 'primed', i.e. spurred by the crRNA-mediated detection of a complete or partial match to an invader sequence. Here we show that primed adaptation occurs in Pyrococcus furiosus. Although P. furiosus has three distinct CRISPR-Cas interference systems (I-B, I-A and III-B), only the I-B system and Cas3 were necessary for priming. Cas4, which is important for selection and processing of new spacers in naïve adaptation, was also essential for priming. Loss of either the I-B effector proteins or Cas3 reduced naïve adaptation. However, when Cas3 and all crRNP genes were deleted, uptake of correctly processed spacers was observed, indicating that none of these interference proteins are necessary for naïve adaptation., (© The Author(s) 2020. Published by Oxford University Press on behalf of Nucleic Acids Research.)

We report here the complete genome sequence of Streptococcus ratti strain JH145. Streptococcus ratti is a cariogenic species of mutans streptococcus that has been isolated from rat and human teeth. The strain JH145, derived from strain BHT-2, is interesting for oral health because it does not produce cariogenic lactic acid but shows robust biofilm production., (Copyright © 2020 Garrett et al.)

CRISPR-Cas systems provide heritable immunity against viruses by capturing short invader DNA sequences, termed spacers, and incorporating them into the CRISPR loci of the prokaryotic host genome. Here, we investigate DNA elements that control accurate spacer uptake in the type II-A CRISPR locus of Streptococcus thermophilus. We determined that purified Cas1 and Cas2 proteins catalyze spacer integration with high specificity for CRISPR repeat junctions. We show that 10 bp of the CRISPR leader sequence is critical for stimulating polarized integration preferentially at the repeat proximal to the leader. Spacer integration proceeds through a two-step transesterification reaction where the 3' hydroxyl groups of the spacer target both repeat borders on opposite strands. The leader-proximal end of the repeat is preferentially targeted for the first site of integration through recognition of sequences spanning the leader-repeat junction. Subsequently, second-site integration at the leader-distal end of the repeat is specified by multiple determinants including a length-defining mechanism relying on a repeat element proximal to the second site of integration. Our results highlight the intrinsic ability of type II Cas1/Cas2 proteins to coordinate directional and site-specific spacer integration into the CRISPR locus to ensure precise duplication of the repeat required for CRISPR immunity., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Acquiring foreign spacer DNA into the CRISPR locus is an essential primary step of the CRISPR-Cas pathway in prokaryotes for developing host immunity to mobile genetic elements. Here, we investigate spacer integration in vitro using proteins from Pyrococcus furiosus and demonstrate that Cas1 and Cas2 are sufficient to accurately integrate spacers into a minimal CRISPR locus. Using high-throughput sequencing, we identified high frequency spacer integration occurring at the same CRISPR repeat border sites utilized in vivo, as well as at several non-CRISPR plasmid sequences which share features with repeats. Analysis of non-CRISPR integration sites revealed that Cas1 and Cas2 are directed to catalyze full-site spacer integration at specific DNA stretches where guanines and/or cytosines are 30 base pairs apart and the intervening sequence harbors several positionally conserved bases. Moreover, assaying a series of CRISPR repeat mutations, followed by sequencing of the integration products, revealed that the specificity of integration is primarily directed by sequences at the leader-repeat junction as well as an adenine-rich sequence block in the mid-repeat. Together, our results indicate that P. furiosus Cas1 and Cas2 recognize multiple sequence features distributed over a 30 base pair DNA region for accurate spacer integration at the CRISPR repeat., (© The Author(s) 2019. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Increasing evidence suggests that transcriptional control and chromatin activities at large involve regulatory RNAs, which likely enlist specific RNA-binding proteins (RBPs). Although multiple RBPs have been implicated in transcription control, it has remained unclear how extensively RBPs directly act on chromatin. We embarked on a large-scale RBP ChIP-seq analysis, revealing widespread RBP presence in active chromatin regions in the human genome. Like transcription factors (TFs), RBPs also show strong preference for hotspots in the genome, particularly gene promoters, where their association is frequently linked to transcriptional output. Unsupervised clustering reveals extensive co-association between TFs and RBPs, as exemplified by YY1, a known RNA-dependent TF, and RBM25, an RBP involved in splicing regulation. Remarkably, RBM25 depletion attenuates all YY1-dependent activities, including chromatin binding, DNA looping, and transcription. We propose that various RBPs may enhance network interaction through harnessing regulatory RNAs to control transcription., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Translation and decay of eukaryotic mRNAs is controlled by shortening of the poly(A) tail and release of the poly(A)-binding protein Pab1/PABP. The Ccr4-Not complex contains two exonucleases-Ccr4 and Caf1/Pop2-that mediate mRNA deadenylation. Here, using a fully reconstituted biochemical system with proteins from the fission yeast Schizosaccharomyces pombe, we show that Pab1 interacts with Ccr4-Not, stimulates deadenylation, and differentiates the roles of the nuclease enzymes. Surprisingly, Pab1 release relies on Ccr4 activity. In agreement with this, in vivo experiments in budding yeast show that Ccr4 is a general deadenylase that acts on all mRNAs. In contrast, Caf1 only trims poly(A) not bound by Pab1. As a consequence, Caf1 is a specialized deadenylase required for the selective deadenylation of transcripts with lower rates of translation elongation and reduced Pab1 occupancy. These findings reveal a coupling between the rates of translation and deadenylation that is dependent on Pab1 and Ccr4-Not., (Copyright © 2018 MRC Laboratory of Molecular Biology. Published by Elsevier Inc. All rights reserved.)

RNA binding proteins (RBPs) orchestrate the production, processing, and function of mRNAs. Here, we present the affinity landscapes of 78 human RBPs using an unbiased assay that determines the sequence, structure, and context preferences of these proteins in vitro by deep sequencing of bound RNAs. These data enable construction of "RNA maps" of RBP activity without requiring crosslinking-based assays. We found an unexpectedly low diversity of RNA motifs, implying frequent convergence of binding specificity toward a relatively small set of RNA motifs, many with low compositional complexity. Offsetting this trend, however, we observed extensive preferences for contextual features distinct from short linear RNA motifs, including spaced "bipartite" motifs, biased flanking nucleotide composition, and bias away from or toward RNA structure. Our results emphasize the importance of contextual features in RNA recognition, which likely enable targeting of distinct subsets of transcripts by different RBPs that recognize the same linear motif., (Copyright © 2018 The Authors. Published by Elsevier Inc. All rights reserved.)

To achieve adaptive and heritable immunity against viruses and other mobile genetic elements, CRISPR-Cas systems must capture and store short DNA fragments (spacers) from these foreign elements into host genomic CRISPR arrays. This process is catalyzed by conserved Cas1/Cas2 integration complexes, but the specific roles of another highly conserved protein linked to spacer acquisition, the Cas4 nuclease, are just now emerging. Here, we show that two Cas4 nucleases (Cas4-1 and Cas4-2) play critical roles in CRISPR spacer acquisition in Pyrococcus furiosus. The nuclease activities of both Cas4 proteins are required to process protospacers to the correct size. Cas4-1 specifies the upstream PAM (protospacer adjacent motif), while Cas4-2 specifies the conserved downstream motif. Both Cas4 proteins ensure CRISPR spacer integration in a defined orientation leading to CRISPR immunity. Collectively, these findings provide in vivo evidence for critical roles of Cas4 nucleases in protospacer generation and functional spacer integration at CRISPR arrays., (Copyright © 2018 Elsevier Inc. All rights reserved.)

We report here the complete genome sequence of Streptococcus thermophilus DGCC 7710. S. thermophilus is widely used in industrial dairy production., (Copyright © 2018 Hatmaker et al.)

The generation of inhibitory interneuron progenitors from human embryonic stem cells (ESCs) is of great interest due to their potential use in transplantation therapies designed to treat central nervous system disorders. The medial ganglionic eminence (MGE) is a transient embryonic structure in the ventral telencephalon that is a major source of cortical GABAergic inhibitory interneuron progenitors. These progenitors migrate tangentially to sites in the cortex and differentiate into a variety of interneuron subtypes, forming local synaptic connections with excitatory projection neurons to modulate activity of the cortical circuitry. The homeobox domain-containing transcription factor NKX2.1 is highly expressed in the MGE and pre-optic area of the ventral subpallium and is essential for specifying cortical interneuron fate. Using a combination of growth factor agonists and antagonists to specify ventral telencephalic fates, we previously optimized a protocol for the efficient generation of NKX2.1-positive MGE-like neural progenitors from human ESCs. To establish their identity, we now characterize the transcriptome of these MGE-like neural progenitors using RNA sequencing and demonstrate the capacity of these cells to differentiate into inhibitory interneurons in vitro using a neuron-astrocyte co-culture system. These data provide information on the potential origin of interneurons in the human brain.

In eukaryotic cells, RNAs exist as ribonucleoprotein particles (RNPs). Despite the importance of these complexes in many biological processes, including splicing, polyadenylation, stability, transportation, localization, and translation, their compositions are largely unknown. We affinity-purified 20 distinct RNA-binding proteins (RBPs) from cultured Drosophila melanogaster cells under native conditions and identified both the RNA and protein compositions of these RNP complexes. We identified "high occupancy target" (HOT) RNAs that interact with the majority of the RBPs we surveyed. HOT RNAs encode components of the nonsense-mediated decay and splicing machinery, as well as RNA-binding and translation initiation proteins. The RNP complexes contain proteins and mRNAs involved in RNA binding and post-transcriptional regulation. Genes with the capacity to produce hundreds of mRNA isoforms, ultracomplex genes, interact extensively with heterogeneous nuclear ribonuclear proteins (hnRNPs). Our data are consistent with a model in which subsets of RNPs include mRNA and protein products from the same gene, indicating the widespread existence of auto-regulatory RNPs. From the simultaneous acquisition and integrative analysis of protein and RNA constituents of RNPs, we identify extensive cross-regulatory and hierarchical interactions in post-transcriptional control., (© 2015 Stoiber et al.; Published by Cold Spring Harbor Laboratory Press.)

Alternative splicing is regulated by RNA binding proteins (RBPs) that recognize pre-mRNA sequence elements and activate or repress adjacent exons. Here, we used RNA interference and RNA-seq to identify splicing events regulated by 56 Drosophila proteins, some previously unknown to regulate splicing. Nearly all proteins affected alternative first exons, suggesting that RBPs play important roles in first exon choice. Half of the splicing events were regulated by multiple proteins, demonstrating extensive combinatorial regulation. We observed that SR and hnRNP proteins tend to act coordinately with each other, not antagonistically. We also identified a cross-regulatory network where splicing regulators affected the splicing of pre-mRNAs encoding other splicing regulators. This large-scale study substantially enhances our understanding of recent models of splicing regulation and provides a resource of thousands of exons that are regulated by 56 diverse RBPs., (© 2015 Brooks et al.; Published by Cold Spring Harbor Laboratory Press.)

In species with a heterogametic sex, population genetics theory predicts that DNA sequences on the X chromosome can evolve faster than comparable sequences on autosomes. Both neutral and nonneutral evolutionary processes can generate this pattern. Complex traits like gene expression are not predicted to have accelerated evolution by these theories, yet a "faster-X" pattern of gene expression divergence has recently been reported for both Drosophila and mammals. Here, we test the hypothesis that accelerated adaptive evolution of cis-regulatory sequences on the X chromosome is responsible for this pattern by comparing the relative contributions of cis- and trans-regulatory changes to patterns of faster-X expression divergence observed between strains and species of Drosophila with a range of divergence times. We find support for this hypothesis, especially among male-biased genes, when comparing different species. However, we also find evidence that trans-regulatory differences contribute to a faster-X pattern of expression divergence both within and between species. This contribution is surprising because trans-acting regulators of X-linked genes are generally assumed to be randomly distributed throughout the genome. We found, however, that X-linked transcription factors appear to preferentially regulate expression of X-linked genes, providing a potential mechanistic explanation for this result. The contribution of trans-regulatory variation to faster-X expression divergence was larger within than between species, suggesting that it is more likely to result from neutral processes than positive selection. These data show how accelerated evolution of both coding and noncoding sequences on the X chromosome can lead to accelerated expression divergence on the X chromosome relative to autosomes., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Animal transcriptomes are dynamic, with each cell type, tissue and organ system expressing an ensemble of transcript isoforms that give rise to substantial diversity. Here we have identified new genes, transcripts and proteins using poly(A)+ RNA sequencing from Drosophila melanogaster in cultured cell lines, dissected organ systems and under environmental perturbations. We found that a small set of mostly neural-specific genes has the potential to encode thousands of transcripts each through extensive alternative promoter usage and RNA splicing. The magnitudes of splicing changes are larger between tissues than between developmental stages, and most sex-specific splicing is gonad-specific. Gonads express hundreds of previously unknown coding and long non-coding RNAs (lncRNAs), some of which are antisense to protein-coding genes and produce short regulatory RNAs. Furthermore, previously identified pervasive intergenic transcription occurs primarily within newly identified introns. The fly transcriptome is substantially more complex than previously recognized, with this complexity arising from combinatorial usage of promoters, splice sites and polyadenylation sites.

Accurate gene model annotation of reference genomes is critical for making them useful. The modENCODE project has improved the D. melanogaster genome annotation by using deep and diverse high-throughput data. Since transcriptional activity that has been evolutionarily conserved is likely to have an advantageous function, we have performed large-scale interspecific comparisons to increase confidence in predicted annotations. To support comparative genomics, we filled in divergence gaps in the Drosophila phylogeny by generating draft genomes for eight new species. For comparative transcriptome analysis, we generated mRNA expression profiles on 81 samples from multiple tissues and developmental stages of 15 Drosophila species, and we performed cap analysis of gene expression in D. melanogaster and D. pseudoobscura. We also describe conservation of four distinct core promoter structures composed of combinations of elements at three positions. Overall, each type of genomic feature shows a characteristic divergence rate relative to neutral models, highlighting the value of multispecies alignment in annotating a target genome that should prove useful in the annotation of other high priority genomes, especially human and other mammalian genomes that are rich in noncoding sequences. We report that the vast majority of elements in the annotation are evolutionarily conserved, indicating that the annotation will be an important springboard for functional genetic testing by the Drosophila community., (© 2014 Chen et al.; Published by Cold Spring Harbor Laboratory Press.)

Anopheles darlingi is the principal neotropical malaria vector, responsible for more than a million cases of malaria per year on the American continent. Anopheles darlingi diverged from the African and Asian malaria vectors ∼100 million years ago (mya) and successfully adapted to the New World environment. Here we present an annotated reference A. darlingi genome, sequenced from a wild population of males and females collected in the Brazilian Amazon. A total of 10 481 predicted protein-coding genes were annotated, 72% of which have their closest counterpart in Anopheles gambiae and 21% have highest similarity with other mosquito species. In spite of a long period of divergent evolution, conserved gene synteny was observed between A. darlingi and A. gambiae. More than 10 million single nucleotide polymorphisms and short indels with potential use as genetic markers were identified. Transposable elements correspond to 2.3% of the A. darlingi genome. Genes associated with hematophagy, immunity and insecticide resistance, directly involved in vector-human and vector-parasite interactions, were identified and discussed. This study represents the first effort to sequence the genome of a neotropical malaria vector, and opens a new window through which we can contemplate the evolutionary history of anopheline mosquitoes. It also provides valuable information that may lead to novel strategies to reduce malaria transmission on the South American continent. The A. darlingi genome is accessible at www.labinfo.lncc.br/index.php/anopheles-darlingi.

Many investigators who do high-throughput sequencing of mRNA (mRNA-Seq) use kits for library preparation purchased from Illumina. Recognizing that these kits are continually being updated and improved, we provide here some background information, tips, and troubleshooting advice for the kits available at the time of this writing.

Drosophila melanogaster is one of the most well studied genetic model organisms; nonetheless, its genome still contains unannotated coding and non-coding genes, transcripts, exons and RNA editing sites. Full discovery and annotation are pre-requisites for understanding how the regulation of transcription, splicing and RNA editing directs the development of this complex organism. Here we used RNA-Seq, tiling microarrays and cDNA sequencing to explore the transcriptome in 30 distinct developmental stages. We identified 111,195 new elements, including thousands of genes, coding and non-coding transcripts, exons, splicing and editing events, and inferred protein isoforms that previously eluded discovery using established experimental, prediction and conservation-based approaches. These data substantially expand the number of known transcribed elements in the Drosophila genome and provide a high-resolution view of transcriptome dynamics throughout development.

Drosophila melanogaster cell lines are important resources for cell biologists. Here, we catalog the expression of exons, genes, and unannotated transcriptional signals for 25 lines. Unannotated transcription is substantial (typically 19% of euchromatic signal). Conservatively, we identify 1405 novel transcribed regions; 684 of these appear to be new exons of neighboring, often distant, genes. Sixty-four percent of genes are expressed detectably in at least one line, but only 21% are detected in all lines. Each cell line expresses, on average, 5885 genes, including a common set of 3109. Expression levels vary over several orders of magnitude. Major signaling pathways are well represented: most differentiation pathways are "off" and survival/growth pathways "on." Roughly 50% of the genes expressed by each line are not part of the common set, and these show considerable individuality. Thirty-one percent are expressed at a higher level in at least one cell line than in any single developmental stage, suggesting that each line is enriched for genes characteristic of small sets of cells. Most remarkable is that imaginal disc-derived lines can generally be assigned, on the basis of expression, to small territories within developing discs. These mappings reveal unexpected stability of even fine-grained spatial determination. No two cell lines show identical transcription factor expression. We conclude that each line has retained features of an individual founder cell superimposed on a common "cell line" gene expression pattern.

To gain insight into how genomic information is translated into cellular and developmental programs, the Drosophila model organism Encyclopedia of DNA Elements (modENCODE) project is comprehensively mapping transcripts, histone modifications, chromosomal proteins, transcription factors, replication proteins and intermediates, and nucleosome properties across a developmental time course and in multiple cell lines. We have generated more than 700 data sets and discovered protein-coding, noncoding, RNA regulatory, replication, and chromatin elements, more than tripling the annotated portion of the Drosophila genome. Correlated activity patterns of these elements reveal a functional regulatory network, which predicts putative new functions for genes, reveals stage- and tissue-specific regulators, and enables gene-expression prediction. Our results provide a foundation for directed experimental and computational studies in Drosophila and related species and also a model for systematic data integration toward comprehensive genomic and functional annotation.

The collection of components required to carry out the intricate processes involved in generating and maintaining a living, breathing and, sometimes, thinking organism is staggeringly complex. Where do all of the parts come from? Early estimates stated that about 100,000 genes would be required to make up a mammal; however, the actual number is less than one-quarter of that, barely four times the number of genes in budding yeast. It is now clear that the 'missing' information is in large part provided by alternative splicing, the process by which multiple different functional messenger RNAs, and therefore proteins, can be synthesized from a single gene.

Alternative splicing is a powerful means of controlling gene expression and increasing protein diversity. Most genes express a limited number of mRNA isoforms, but there are several examples of genes that use alternative splicing to generate hundreds, thousands and even tens of thousands of isoforms. Collectively such genes are considered to undergo complex alternative splicing. The best example is the Drosophila Down syndrome cell adhesion molecule (Dscam) gene, which can generate 38,016 isoforms by the alternative splicing of 95 variable exons. In this review, we will describe several genes that use complex alternative splicing to generate large repertoires of mRNAs and what is known about the mechanisms by which they do so.

Drosophila Dscam encodes 38,016 distinct axon guidance receptors through the mutually exclusive alternative splicing of 95 variable exons. Importantly, known mechanisms that ensure the mutually exclusive splicing of pairs of exons cannot explain this phenomenon in Dscam. I have identified two classes of conserved elements in the Dscam exon 6 cluster, which contains 48 alternative exons--the docking site, located in the intron downstream of constitutive exon 5, and the selector sequences, which are located upstream of each exon 6 variant. Strikingly, each selector sequence is complementary to a portion of the docking site, and this pairing juxtaposes one, and only one, alternative exon to the upstream constitutive exon. The mutually exclusive nature of the docking site:selector sequence interactions suggests that the formation of these competing RNA structures is a central component of the mechanism guaranteeing that only one exon 6 variant is included in each Dscam mRNA.

SR proteins are essential metazoan splicing factors that contain an RNA-binding domain and an arginine/serine-rich domain that functions to promote assembly of the spliceosome. The prevailing model over the past several years suggests that the RS domains function as protein-interaction domains. However, two new papers from Green et al. demonstrate that these RS domains directly contact the pre-mRNA within the functional spliceosome. The sequential character of these contacts suggests that RS domain interactions with RNA promote spliceosome assembly.

Single-strand conformational polymorphism analysis has been used successfully to identify single nucleotide changes within sequences based on the fact that multidetection enhancement gels will separate molecules based on their conformation rather than their size. We have expanded the utility of this technique to analyze easily the alternative splicing of pre-mRNAs containing multiple mutually exclusive exons of the same size. We have used this technique to study the Caenorhabditis elegans let-2 gene containing two alternative exons and the Drosophilia melanogaster Dscam gene, which contains 12 mutually exclusive exons. The ease and the quantitative nature of this technique should be very useful.

RNA interference (RNAi) is a useful tool for degrading targeted messenger RNAs (mRNAs) and thus "knocking down" the abundance of the encoded protein. We have been using RNAi in cultured Drosophila cells to evaluate the effect of "knocking down" numerous mRNA processing factors on the alternative splicing of specific pre-mRNAs. This relatively simple technique has allowed us to identify a number of splicing factors that impact the alternative splicing of particular alternatively spliced exons. This approach can be extended to examine the splicing of nearly any gene.

SR proteins are essential pre-mRNA splicing factors that have been shown to bind a number of exonic splicing enhancers where they function to stimulate the splicing of adjacent introns. Members of the SR protein family contain one or two N-terminal RNA binding domains, as well as a C-terminal arginine-serine (RS) rich domain. The RS domains mediate protein-protein interactions with other RS domain containing proteins and are essential for many, but not all, SR protein functions. Hybrid proteins containing an RS domain fused to the bacteriophage MS2 coat protein are sufficient to activate enhancer-dependent splicing in HeLa cell nuclear extract when bound to the pre-mRNA. Here we report progress towards determining the protein sequence requirements for RS domain function. We show that the RS domains from non-SR proteins can also function as splicing activation domains when tethered to the pre-mRNA. Truncation experiments with the RS domain of the human SR protein 9G8 identified a 29 amino acid segment, containing 26 arginine or serine residues, that is sufficient to activate splicing when fused to MS2. We also show that synthetic domains composed solely of RS dipeptides are capable of activating splicing, although their potency is proportional to their size.

Alternative splicing is an important means of regulating the expression of eukaryotic genes and enhancing protein diversity. A detailed examination of the Drosophila Sex-lethal gene has led to two significant discoveries-the role of the splicing factor SPF45 in defining the site of exon ligation, and that alternative splicing can be regulated at the second step.

The essential splicing factor U2AF (U2 auxiliary factor) is a heterodimer composed of 65-kDa (U2AF(65)) and 35-kDa (U2AF(35)) subunits. U2AF(35) has multiple functions in pre-mRNA splicing. First, U2AF(35) has been shown to function by directly interacting with the AG at the 3' splice site. Second, U2AF(35) is thought to play a role in the recruitment of U2AF(65) by serine-arginine-rich (SR) proteins in enhancer-dependent splicing. It has been proposed that the physical interaction between the arginine-serine-rich (RS) domain of U2AF(35) and SR proteins is important for this activity. However, other data suggest that this may not be the case. Here, we report the identification of a mammalian gene that encodes a 26-kDa protein bearing strong sequence similarity to U2AF(35), designated U2AF(26). The N-terminal 187 amino acids of U2AF(35) and U2AF(26) are nearly identical. However, the C-terminal domain of U2AF(26) lacks many characteristics of the U2AF(35) RS domain and, therefore, might be incapable of interacting with SR proteins. We show that U2AF(26) can associate with U2AF(65) and can functionally substitute for U2AF(35) in both constitutive and enhancer-dependent splicing, demonstrating that the RS domain of the small U2AF subunit is not required for splicing enhancer function. Finally, we show that U2AF(26) functions by enhancing the binding of U2AF(65) to weak 3' splice sites. These studies identify U2AF(26) as a mammalian splicing factor and demonstrate that distinct U2AF complexes can participate in pre-mRNA splicing. Based on its sequence and functional similarity to U2AF(35), U2AF(26) may play a role in regulating alternative splicing.