Marcela Dávila López

Computational screen for spliceosomal RNA genes aids in defining the phylogenetic distribution of major and minor spliceosomal components

The RNA molecules of the spliceosome are critical for specificity and catalysis during splicing of eukaryotic pre-mRNA. In order to examine the evolution and phylogenetic distribution of these RNAs, we analyzed 149 eukaryotic genomes representing a broad range of phylogenetic groups. RNAs were predicted using high-sensitivity local alignment methods and profile HMMs in combination with covariance models. The results provide the most comprehensive view so far of the phylogenetic distribution of spliceosomal RNAs. RNAs were predicted in many phylogenetic groups where these RNA were not previously reported. Examples are RNAs of the major (U2-type) spliceosome in all fungal lineages, in lower metazoa and many protozoa. We also identified the minor (U12-type) spliceosomal U11 and U6atac RNAs in Acanthamoeba castellanii, where U12 spliceosomal RNA as well as minor introns were reported recently. In addition, minor-spliceosome-specific RNAs were identified in a number of phylogenetic groups where previously such RNAs were not observed, including the nematode Trichinella spiralis, the slime mold Physarum polycephalum and the fungal lineages Zygomycota and Chytridiomycota. The detailed map of the distribution of the U12-type RNA genes supports an early origin of the minor spliceosome and points to a number of occasions during evolution where it was lost.

Polyglucosan body myopathy caused by defective ubiquitin ligase RBCK1.

Glycogen storage diseases are important causes of myopathy and cardiomyopathy. We describe 10 patients from 8 families with childhood or juvenile onset of myopathy, 8 of whom also had rapidly progressive cardiomyopathy, requiring heart transplant in 4. The patients were homozygous or compound heterozygous for missense or truncating mutations in RBCK1, which encodes for a ubiquitin ligase, and had extensive polyglucosan accumulation in skeletal muscle and in the heart in cases of cardiomyopathy. We conclude that RBCK1 deficiency is a frequent cause of polyglucosan storage myopathy associated with progressive muscle weakness and cardiomyopathy. Ann Neurol 2013;74:914–919

Inventory and analysis of the protein subunits of the ribonucleases P and MRP provides further evidence of homology between the yeast and human enzymes

The RNases P and MRP are involved in tRNA and rRNA processing, respectively. Both enzymes in eukaryotes are composed of an RNA molecule and 9–12 protein subunits. Most of the protein subunits are shared between RNases P and MRP. We have here performed a computational analysis of the protein subunits in a broad range of eukaryotic organisms using profile-based searches and phylogenetic methods. A number of novel homologues were identified, giving rise to a more complete inventory of RNase P/MRP proteins. We present evidence of a relationship between fungal Pop8 and the protein subunit families Rpp14/Pop5 as well as between fungal Pop6 and metazoan Rpp25. These relationships further emphasize a structural and functional similarity between the yeast and human P/MRP complexes. We have also identified novel P and MRP RNAs and analysis of all available sequences revealed a K-turn motif in a large number of these RNAs. We suggest that this motif is a binding site for the Pop3/Rpp38 proteins and we discuss other structural features of the RNA subunit and possible relationships to the protein subunit repertoire.

Small RNA deep sequencing discriminates subsets of extracellular vesicles released by melanoma cells – Evidence of unique microRNA cargos

Melanoma cells release different types of extracellular vesicles (EVs) into the extracellular milieu that are involved with communication and signaling in the tumor microenvironment. Subsets of EVs include exosomes, microvesicles, and apoptotic bodies that carry protein and genetic (RNA) cargos. To define the contribution of the RNA cargo of melanoma cell derived EVs we performed small RNA sequencing to identify different small RNAs in the EV subsets. Using validated centrifugation protocols, we separated these EV subsets released by the melanoma cell line MML-1, and performed RNA sequencing with the Ion Torrent platform. Various, but different, non-coding RNAs were detected in the EV subsets, including microRNA, mitochondrial associated tRNA, small nucleolar RNA, small nuclear RNA, Ro associated Y-RNA, vault RNA and Y-RNA. We identified in total 1041 miRNAs in cells and EV subsets. Hierarchical clustering showed enrichment of specific miRNAs in exosomes, including hsa-miR-214-3p, hsa-miR-199a-3p and hsa-miR-155-5p, all being associated with melanoma progression. Comparison of exosomal miRNAs with miRNAs in clinical melanoma samples indicate that multiple miRNAs in exosomes also are expressed specifically in melanoma tissues, but not in benign naevi. This study shows for the first time the presence of distinct small RNAs in subsets of EVs released by melanoma cells, with significant similarities to clinical melanoma tissue, and provides unique insights into the contribution of EV associated extracellular RNA in cancer.

Conserved and variable domains of RNase MRP RNA

Ribonuclease MRP is a eukaryotic ribonucleoprotein complex consisting of one RNA molecule and 7-10 protein subunits. One important function of MRP is to catalyze an endonucleolytic cleavage during processing of rRNA precursors. RNase MRP is evolutionary related to RNase P which is critical for tRNA processing. A large number of MRP RNA sequences that now are available have been used to identify conserved primary and secondary structure features of the molecule. MRP RNA has structural features in common with P RNA such as a conserved catalytic core, but it also has unique features and is characterized by a domain highly variable between species. Information regarding primary and secondary structure features is of interest not only in basic studies of the function of MRP RNA, but also because mutations in the RNA give rise to human genetic diseases such as cartilage-hair hypoplasia.

Histone modifications influence mediator interactions with chromatin

The Mediator complex transmits activation signals from DNA bound transcription factors to the core transcription machinery. Genome wide localization studies have demonstrated that Mediator occupancy not only correlates with high levels of transcription, but that the complex also is present at transcriptionally silenced locations. We provide evidence that Mediator localization is guided by an interaction with histone tails, and that this interaction is regulated by their post-translational modifications. A quantitative, high-density genetic interaction map revealed links between Mediator components and factors affecting chromatin structure, especially histone deacetylases. Peptide binding assays demonstrated that pure wild-type Mediator forms stable complexes with the tails of Histone H3 and H4. These binding assays also showed Mediator—histone H4 peptide interactions are specifically inhibited by acetylation of the histone H4 lysine 16, a residue critical in transcriptional silencing. Finally, these findings were validated by tiling array analysis that revealed a broad correlation between Mediator and nucleosome occupancy in vivo, but a negative correlation between Mediator and nucleosomes acetylated at histone H4 lysine 16. Our studies show that chromatin structure and the acetylation state of histones are intimately connected to Mediator localization.

Analysis of gene order conservation in eukaryotes identifies transcriptionally and functionally linked genes.

The order of genes in eukaryotes is not entirely random. Studies of gene order conservation are important to understand genome evolution and to reveal mechanisms why certain neighboring genes are more difficult to separate during evolution. Here, genome-wide gene order information was compiled for 64 species, representing a wide variety of eukaryotic phyla. This information is presented in a browser where gene order may be displayed and compared between species. Factors related to non-random gene order in eukaryotes were examined by considering pairs of neighboring genes. The evolutionary conservation of gene pairs was studied with respect to relative transcriptional direction, intergenic distance and functional relationship as inferred by gene ontology. The results show that among gene pairs that are conserved the divergently and co-directionally transcribed genes are much more common than those that are convergently transcribed. Furthermore, highly conserved pairs, in particular those of fungi, are characterized by a short intergenic distance. Finally, gene pairs of metazoa and fungi that are evolutionary conserved and that are divergently transcribed are much more likely to be related by function as compared to poorly conserved gene pairs. One example is the ribosomal protein gene pair L13/S16, which is unusual as it occurs both in fungi and alveolates. A specific functional relationship between these two proteins is also suggested by the fact that they are part of the same operon in both eubacteria and archaea. In conclusion, factors associated with non-random gene order in eukaryotes include relative gene orientation, intergenic distance and functional relationships. It seems likely that certain pairs of genes are conserved because the genes involved have a transcriptional and/or functional relationship. The results also indicate that studies of gene order conservation aid in identifying genes that are related in terms of transcriptional control.

Regulation of DNA replication at the end of the mitochondrial D-loop involves the helicase TWINKLE and a conserved sequence element

The majority of mitochondrial DNA replication events are terminated prematurely. The nascent DNA remains stably associated with the template, forming a triple-stranded displacement loop (D-loop) structure. However, the function of the D-loop region of the mitochondrial genome remains poorly understood. Using a comparative genomics approach we here identify two closely related 15 nt sequence motifs of the D-loop, strongly conserved among vertebrates. One motif is at the D-loop 5′-end and is part of the conserved sequence block 1 (CSB1). The other motif, here denoted coreTAS, is at the D-loop 3′-end. Both these sequences may prevent transcription across the D-loop region, since light and heavy strand transcription is terminated at CSB1 and coreTAS, respectively. Interestingly, the replication of the nascent D-loop strand, occurring in a direction opposite to that of heavy strand transcription, is also terminated at coreTAS, suggesting that coreTAS is involved in termination of both transcription and replication. Finally, we demonstrate that the loading of the helicase TWINKLE at coreTAS is reversible, implying that this site is a crucial component of a switch between D-loop formation and full-length mitochondrial DNA replication.

Mediator Promotes CENP-A Incorporation at Fission Yeast Centromeres.

ABSTRACT At Schizosaccharomyces pombe centromeres, heterochromatin formation is required for de novo incorporation of the histone H3 variant CENP-ACnp1, which in turn directs kinetochore assembly and ultimately chromosome segregation during mitosis. Noncoding RNAs (ncRNAs) transcribed by RNA polymerase II (Pol II) directs heterochromatin formation through not only the RNA interference (RNAi) machinery but also RNAi-independent RNA processing factors. Control of centromeric ncRNA transcription is therefore a key factor for proper centromere function. We here demonstrate that Mediator directs ncRNA transcription and regulates centromeric heterochromatin formation in fission yeast. Mediator colocalizes with Pol II at centromeres, and loss of the Mediator subunit Med20 causes a dramatic increase in pericentromeric transcription and desilencing of the core centromere. As a consequence, heterochromatin formation is impaired via both the RNAi-dependent and -independent pathways, resulting in loss of CENP-ACnp1 from the core centromere, a defect in kinetochore function, and a severe chromosome segregation defect. Interestingly, the increased centromeric transcription observed in med20Δ cells appears to directly block CENP-ACnp1 incorporation since inhibition of Pol II transcription can suppress the observed phenotypes. Our data thus identify Mediator as a crucial regulator of ncRNA transcription at fission yeast centromeres and add another crucial layer of regulation to centromere function.

Cyclin-Dependent Kinase 8 Regulates Mitotic Commitment in Fission Yeast

ABSTRACT Temporal changes in transcription programs are coupled to control of cell growth and division. We here report that Mediator, a conserved coregulator of eukaryotic transcription, is part of a regulatory pathway that controls mitotic entry in fission yeast. The Mediator subunit cyclin-dependent kinase 8 (Cdk8) phosphorylates the forkhead 2 (Fkh2) protein in a periodic manner that coincides with gene activation during mitosis. Phosphorylation prevents degradation of the Fkh2 transcription factor by the proteasome, thus ensuring cell cycle-dependent variations in Fkh2 levels. Interestingly, Cdk8-dependent phosphorylation of Fkh2 controls mitotic entry, and mitotic entry is delayed by inactivation of the Cdk8 kinase activity or mutations replacing the phosphorylated serine residues of Fkh2. In addition, mutations in Fkh2, which mimic protein phosphorylation, lead to premature mitotic entry. Therefore, Fkh2 regulates not only the onset of mitotic transcription but also the correct timing of mitotic entry via effects on the Wee1 kinase. Our findings thus establish a new pathway linking the Mediator complex to control of mitotic transcription and regulation of mitotic entry in fission yeast.