Selected Model Organisms Research Articles

INTREPPPID-an orthologue-informed quintuplet network for cross-species prediction of protein-protein interaction.

An overwhelming majority of protein-protein interaction (PPI) studies are conducted in a select few model organisms largely due to constraints in time and cost of the associated 'wet lab' experiments. In silico PPI inference methods are ideal tools to overcome these limitations, but often struggle with cross-species predictions. We present INTREPPPID, a method that incorporates orthology data using a new 'quintuplet' neural network, which is constructed with five parallel encoders with shared parameters. INTREPPPID incorporates both a PPI classification task and an orthologous locality task. The latter learns embeddings of orthologues that have small Euclidean distances between them and large distances between embeddings of all other proteins. INTREPPPID outperforms all other leading PPI inference methods tested on both the intraspecies and cross-species tasks using strict evaluation datasets. We show that INTREPPPID's orthologous locality loss increases performance because of the biological relevance of the orthologue data and not due to some other specious aspect of the architecture. Finally, we introduce PPI.bio and PPI Origami, a web server interface for INTREPPPID and a software tool for creating strict evaluation datasets, respectively. Together, these two initiatives aim to make both the use and development of PPI inference tools more accessible to the community.

Decoding functional proteome information in model organisms using protein language models.

Protein language models have been tested and proved to be reliable when used on curated datasets but have not yet been applied to full proteomes. Accordingly, we tested how two different machine learning-based methods performed when decoding functional information from the proteomes of selected model organisms. We found that protein language models are more precise and informative than deep learning methods for all the species tested and across the three gene ontologies studied, and that they better recover functional information from transcriptomicexperiments. The results obtained indicate that these language models are likely to be suitable for large-scale annotation and downstream analyses, and we recommend a guide for their use.

The Glutathione System: A Journey from Cyanobacteria to Higher Eukaryotes.

From bacteria to plants and humans, the glutathione system plays a pleiotropic role in cell defense against metabolic, oxidative and metal stresses. Glutathione (GSH), the γ-L-glutamyl-L-cysteinyl-glycine nucleophile tri-peptide, is the central player of this system that acts in redox homeostasis, detoxification and iron metabolism in most living organisms. GSH directly scavenges diverse reactive oxygen species (ROS), such as singlet oxygen, superoxide anion, hydrogen peroxide, hydroxyl radical, nitric oxide and carbon radicals. It also serves as a cofactor for various enzymes, such as glutaredoxins (Grxs), glutathione peroxidases (Gpxs), glutathione reductase (GR) and glutathione-S-transferases (GSTs), which play crucial roles in cell detoxication. This review summarizes what is known concerning the GSH-system (GSH, GSH-derived metabolites and GSH-dependent enzymes) in selected model organisms (Escherichia coli, Saccharomyces cerevisiae, Arabidopsis thaliana and human), emphasizing cyanobacteria for the following reasons. Cyanobacteria are environmentally crucial and biotechnologically important organisms that are regarded as having evolved photosynthesis and the GSH system to protect themselves against the ROS produced by their active photoautotrophic metabolism. Furthermore, cyanobacteria synthesize the GSH-derived metabolites, ergothioneine and phytochelatin, that play crucial roles in cell detoxication in humans and plants, respectively. Cyanobacteria also synthesize the thiol-less GSH homologs ophthalmate and norophthalmate that serve as biomarkers of various diseases in humans. Hence, cyanobacteria are well-suited to thoroughly analyze the role/specificity/redundancy of the players of the GSH-system using a genetic approach (deletion/overproduction) that is hardly feasible with other model organisms (E. coli and S. cerevisiae do not synthesize ergothioneine, while plants and humans acquire it from their soil and their diet, respectively).

Evolution of proteins involved in the final steps of juvenile hormone synthesis

Juvenile hormone (JH), a sesquiterpenoid produced by the insect corpus allatum gland (CA), is a key regulator of insect metamorphosis, reproduction, caste differentiation, and polyphenism. The first part of JH biosynthesis occurs via the universal eukaryotic mevalonate pathway. The final steps involve epoxidation and methylation. However, the sequence of these steps might not be conserved among all insects and Crustacea. Therefore, we used available genomic and transcriptomic data and identified JH acid methyltransferase (JHAMT), analyzed their genomic duplications in selected model organisms, and reconstructed their phylogeny. We have further reconstructed phylogeny of FAMeT proteins and show that evolution of this protein group is more complicated than originally appreciated. The analysis delineates important milestones in the evolution of several JH biosynthetic enzymes in arthropods, reviews major literature data on the last steps of JH synthesis, and defines questions and some hypotheses worth pursuing experimentally.

Estimating Total Quantitative Protein Content in Escherichia coli, Saccharomyces cerevisiae, and HeLa Cells.

The continuous improvement of proteomic techniques, most notably mass spectrometry, has generated quantified proteomes of many organisms with unprecedented depth and accuracy. However, there is still a significant discrepancy in the reported numbers of total protein molecules per specific cell type. In this article, we explore the results of proteomic studies of Escherichia coli, Saccharomyces cerevisiae, and HeLa cells in terms of total protein copy numbers per cell. We observe up to a ten-fold difference between reported values. Investigating possible reasons for this discrepancy, we conclude that neither an unmeasured fraction of the proteome nor biases in the quantification of individual proteins can explain the observed discrepancy. We normalize protein copy numbers in each study using a total protein amount per cell as reported in the literature and create integrated proteome maps of the selected model organisms. Our results indicate that cells contain from one to three million protein molecules per µm3 and that protein copy density decreases with increasing organism complexity.

Heart rate and behavioral responses in three phylogenetically distant aquatic model organisms exposed to environmental concentrations of carbaryl and fenitrothion

Carbaryl and fenitrothion are two insecticides sharing a common mode of action, the inhibition of the acetylcholinesterase (AChE) activity. Their use is now regulated or banned in different countries, and the environmental levels of both compounds in aquatic ecosystems have decreased to the range of pg/L to ng/L. As these concentrations are below the non-observed-adverse-effect-concentrations (NOAEC) for AChE inhibition reported for both compounds in aquatic organisms, there is a general agreement that the current levels of these two chemicals are safe for aquatic organisms. In this study we have exposed zebrafish, Japanese medaka and Daphnia magna to concentrations of carbaryl and fenitrothion under their NOAECs for 24-h, and the effects on heart rate (HR), basal locomotor activity (BLA), visual motor response (VMR), startle response (SR) and its habituation have been evaluated. Both pesticides increased the HR in the three selected model organisms, although the intensity of this effect was chemical-, concentration- and organism-dependent. The exposure to both pesticides also led to a decrease in BLA and an increase in VMR in all three species, although this effect was only significant in zebrafish larvae. For SR and its habituation, the response profile was more species- and concentration-specific. The results presented in this manuscript demonstrate that concentrations of carbaryl and fenitrothion well below their respective NOAECs induce tachycardia and the impairment of ecologically relevant behaviors in phylogenetically distinct aquatic model organisms, both vertebrates and invertebrates, emphasizing the need to include this range of concentrations in the environmental risk assessment.

Zebrafish and Medaka as model organisms for climate change research: Global literature scientometric analysis

Determining how climate change affects aquatic organisms, such as fish is vital, since this could directly or indirectly impact food and protein sources that are important for human nutrition. Thus, identifying suitable organisms for studying the impacts of climate change on aquatic species is essential. It is most effective to select model organisms for climate change study and determine how each organism might adapt within the diversity of organisms present. This study aimed to review the current development and frontiers of climate change’s model organism based on the literature. We conducted a scientometric analysis by differentiating between publications on different model species, the number and origin of authors and affiliations involved, the citation analysis, and the most common keywords used. Increased publication numbers for Zebrafish and Medaka were detected during the analysis of the networks. Our results showed that both species are among the most important aquatic model organisms for climate change related research. Furthermore, we found that these model organisms, especially the Zebrafish are becoming increasingly important towards climate change related studies, because of their simple anatomy and established biological studies. Our analysis could be on the forefront for disseminating and communicating scientific knowledge and impactful discoveries to researchers, academics, policymakers, and to the public worldwide for future contribution to the community resources preservation.

Government funding of research beyond biomedicine: challenges and opportunities for neuroethology

Curiosity-driven research is fundamental for neuroethology and depends crucially on governmental funding. Here, we highlight similarities and differences in funding of curiosity-driven research across countries by comparing two major funding agencies—the National Science Foundation (NSF) in the United States and the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG). We interviewed representatives from each of the two agencies, focusing on general funding trends, levels of young investigator support, career-life balance, and international collaborations. While our analysis revealed a negative trend in NSF funding of biological research, including curiosity-driven research, German researchers in these areas have benefited from a robust positive trend in DFG funding. The main reason for the decrease in curiosity-driven research in the US is that the NSF has only partially been able to compensate for the funding gap resulting from the National Institutes of Health restricting their support to biomedical research using select model organisms. Notwithstanding some differences in funding programs, particularly those relevant for scientists in the postdoctoral phase, both the NSF and DFG clearly support curiosity-driven research.

Surviving Mars: new insights into the persistence of facultative anaerobic microbes from analogue sites

AbstractMars analogue environments are some of the most extreme locations on Earth. Their unique combination of multiples extremes (e.g. high salinity, anoxia and low nutrient availability) make them valuable sources for finding new polyextremophilic microbes, and for exploring the limits of life. Mars, especially at its surface, is still considered to be very hostile to life but it probably possesses geological subsurface niches where the occurrence of (polyextremophilic) life is conceivable. Despite their well-recognized relevance, current knowledge on the capability of (facultative) anaerobic microbes to withstand extraterrestrial/Martian conditions, either as single strains or in communities, is still very sparse. Therefore, space experiments simulating the Martian environmental conditions by using space as a tool for astrobiological research are needed to substantiate the hypotheses of habitability of Mars. Addressing this knowledge gap is one of the main goals of the project MEXEM (Mars EXposed Extremophiles Mixture), where selected model organisms will be subjected to space for a period of 3 months. These experiments will take place on the Exobiology facility (currently under development and implementation), located outside the International Space Station. Such space experiments require a series of preliminary tests and ground data collection for the selected microbial strains. Here, we report on the survivability of Salinisphaera shabanensis and Buttiauxella sp. MASE-IM-9 after exposure to Mars-relevant stress factors (such as desiccation and ultraviolet (UV) radiation under anoxia). Both organisms showed survival after anoxic desiccation for up to 3 months but this could be further extended (nearly doubled) by adding artificial Mars regolith (MGS-1S; 0.5% wt/v) and sucrose (0.1 M). Survival after desiccation was also observed when both organisms were mixed before treatment. Mixing also positively influenced survival after exposure to polychromatic Mars-like UV radiation (200–400 nm) up to 12 kJ m−2, both in suspension and in a desiccated form.

Polarized Dishevelled dissolution and reassembly drives embryonic axis specification in sea star oocytes

The organismal body axes that are formed during embryogenesis are intimately linked to intrinsic asymmetries established at the cellular scale in oocytes.1 However, the mechanisms that generate cellular asymmetries within the oocyte and then transduce that polarity to organismal scale body axes are poorly understood outside of select model organisms. Here, we report an axis-defining event in meiotic oocytes of the sea star Patiria miniata. Dishevelled (Dvl) is a cytoplasmic Wnt pathway effector required for axis development in diverse species,2-4 but the mechanisms governing its function and distribution remain poorly defined. Using time-lapse imaging, we find that Dvl localizes uniformly to puncta throughout the cell cortex in Prophase I-arrested oocytes but becomes enriched at the vegetal pole following meiotic resumption through a dissolution-reassembly mechanism. This process is driven by an initial disassembly phase of Dvl puncta, followed by selective reformation of Dvl assemblies at the vegetal pole. Rather than being driven by Wnt signaling, this localization behavior is coupled to meiotic cell cycle progression and influenced by Lamp1+ endosome association and Frizzled receptors pre-localized within the oocyte cortex. Our results reveal a cell cycle-linked mechanism by which maternal cellular polarity is transduced to the embryo through spatially regulated Dvl dynamics.

Explaining Regeneration: Cells and Limbs as Complex Living Systems, Learning From History.

Regeneration has been investigated since Aristotle, giving rise to many ways of explaining what this process is and how it works. Current research focuses on gene expression and cell signaling of regeneration within individual model organisms. We tend to look to model organisms on the reasoning that because of evolution, information gained from other species must in some respect be generalizable. However, for all that we have uncovered about how regeneration works within individual organisms, we have yet to translate what we have gleaned into achieving the goal of regenerative medicine: to harness and enhance our own regenerative abilities. Turning to history may provide a crucial perspective in advancing us toward this goal. History gives perspective, allowing us to reflect on how our predecessors did their work and what assumptions they made, thus also revealing limitations. History, then, may show us how we can move from our current reductionist thinking focused on particular selected model organisms toward generalizations about this crucial process that operates across complex living systems and move closer to repairing our own damaged bodies.

Developing Functional Genomics Platforms for Fungi.

ABSTRACTFungi are responsible for diseases that result in the deaths of over a million individuals each year and devastating crop infestations that threaten global food supplies. However, outside of a select few model organisms, the majority of fungal genes are uncharacterized. The roles of these genes in the biology of the organism, pathogenesis, and mediating interactions with the environment and other microbes are unknown. Historically, fungal gene characterization has primarily relied on classical genetic screens. However, advances in sequencing technology have enabled more rapid methods of gene functional characterization. Large-scale transcriptional profiling projects are one solution to generating hypotheses about fungal gene function. Together with other ‘omics techniques and newer tools that enable massively parallel mutant screens, knowledge of fungal gene function will be substantially improved. Understanding the function of fungal genes will be instrumental in increasing global food security, protecting ecosystems, and improving health outcomes.

A Cre-Dependent CRISPR/dCas9 System for Gene Expression Regulation in Neurons.

Site-specific genetic and epigenetic targeting of distinct cell populations is a central goal in molecular neuroscience and is crucial to understand the gene regulatory mechanisms that underlie complex phenotypes and behaviors. While recent technological advances have enabled unprecedented control over gene expression, many of these approaches are focused on selected model organisms and/or require labor-intensive customization for different applications. The simplicity and modularity of clustered regularly interspaced short palindromic repeats (CRISPR)-based systems have transformed genome editing and expanded the gene regulatory toolbox. However, there are few available tools for cell-selective CRISPR regulation in neurons. We designed, validated, and optimized CRISPR activation (CRISPRa) and CRISPR interference (CRISPRi) systems for Cre recombinase-dependent gene regulation. Unexpectedly, CRISPRa systems based on a traditional double-floxed inverted open reading frame (DIO) strategy exhibited leaky target gene induction even without Cre. Therefore, we developed an intron-containing Cre-dependent CRISPRa system (SVI-DIO-dCas9-VPR) that alleviated leaky gene induction and outperformed the traditional DIO system at endogenous genes in HEK293T cells and rat primary neuron cultures. Using gene-specific CRISPR sgRNAs, we demonstrate that SVI-DIO-dCas9-VPR can activate numerous rat or human genes (GRM2, Tent5b, Fos, Sstr2, and Gadd45b) in a Cre-specific manner. To illustrate the versatility of this tool, we created a parallel CRISPRi construct that successfully inhibited expression from a luciferase reporter in HEK293T cells only in the presence of Cre. These results provide a robust framework for Cre-dependent CRISPR-dCas9 approaches across different model systems, and enable cell-specific targeting when combined with common Cre driver lines or Cre delivery via viral vectors.

Twelve quick steps for genome assembly and annotation in the classroom.

Eukaryotic genome sequencing and de novo assembly, once the exclusive domain of well-funded international consortia, have become increasingly affordable, thus fitting the budgets of individual research groups. Third-generation long-read DNA sequencing technologies are increasingly used, providing extensive genomic toolkits that were once reserved for a few select model organisms. Generating high-quality genome assemblies and annotations for many aquatic species still presents significant challenges due to their large genome sizes, complexity, and high chromosome numbers. Indeed, selecting the most appropriate sequencing and software platforms and annotation pipelines for a new genome project can be daunting because tools often only work in limited contexts. In genomics, generating a high-quality genome assembly/annotation has become an indispensable tool for better understanding the biology of any species. Herein, we state 12 steps to help researchers get started in genome projects by presenting guidelines that are broadly applicable (to any species), sustainable over time, and cover all aspects of genome assembly and annotation projects from start to finish. We review some commonly used approaches, including practical methods to extract high-quality DNA and choices for the best sequencing platforms and library preparations. In addition, we discuss the range of potential bioinformatics pipelines, including structural and functional annotations (e.g., transposable elements and repetitive sequences). This paper also includes information on how to build a wide community for a genome project, the importance of data management, and how to make the data and results Findable, Accessible, Interoperable, and Reusable (FAIR) by submitting them to a public repository and sharing them with the research community.

CAPt'n of Actin Dynamics: Recent Advances in the Molecular, Developmental and Physiological Functions of Cyclase-Associated Protein (CAP).

Cyclase-associated protein (CAP) has been discovered three decades ago in budding yeast as a protein that associates with the cyclic adenosine monophosphate (cAMP)-producing adenylyl cyclase and that suppresses a hyperactive RAS2 variant. Since that time, CAP has been identified in all eukaryotic species examined and it became evident that the activity in RAS-cAMP signaling is restricted to a limited number of species. Instead, its actin binding activity is conserved among eukaryotes and actin cytoskeleton regulation emerged as its primary function. However, for many years, the molecular functions as well as the developmental and physiological relevance of CAP remained unknown. In the present article, we will compile important recent progress on its molecular functions that identified CAP as a novel key regulator of actin dynamics, i.e., the spatiotemporally controlled assembly and disassembly of actin filaments (F-actin). These studies unraveled a cooperation with ADF/Cofilin and Twinfilin in F-actin disassembly, a nucleotide exchange activity on globular actin monomers (G-actin) that is required for F-actin assembly and an inhibitory function towards the F-actin assembly factor INF2. Moreover, by focusing on selected model organisms, we will review current literature on its developmental and physiological functions, and we will present studies implicating CAP in human pathologies. Together, this review article summarizes and discusses recent achievements in understanding the molecular, developmental and physiological functions of CAP, which led this protein emerge as a novel CAPt’n of actin dynamics.

DeepAdd: Protein function prediction from k-mer embedding and additional features

With the application of new high throughput sequencing technology, a large number of protein sequences is becoming available. Determination of the functional characteristics of these proteins by experiments is an expensive endeavor that requires a lot of time. Furthermore, at the organismal level, such kind of experimental functional analyses can be conducted only for a very few selected model organisms. Computational function prediction methods can be used to fill this gap. The functions of proteins are classified by Gene Ontology (GO), which contains more than 40,000 classifications in three domains, Molecular Function (MF), Biological Process (BP), and Cellular Component (CC). Additionally, since proteins have many functions, function prediction represents a multi-label and multi-class problem. We developed a new method to predict protein function from sequence. To this end, natural language model was used to generate word embedding of sequence and learn features from it by deep learning, and additional features to locate every protein. Our method uses the dependencies between GO classes as background information to construct a deep learning model. We evaluate our method using the standards established by the Computational Assessment of Function Annotation (CAFA) and have noticeable improvement over several algorithms, such as FFPred, DeepGO, GoFDR and other methods compared on the CAFA3 datasets.

The Genome3D Consortium for Structural Annotations of Selected Model Organisms.

Genome3D consortium is a collaborative project involving protein structure prediction and annotation resources developed by six world-leading structural bioinformatics groups, based in the United Kingdom (namely Blundell, Murzin, Gough, Sternberg, Orengo, and Jones). The main objective of Genome3D serves as a common portal to provide both predicted models and annotations of proteins in model organisms, using several resources developed by these labs such as CATH-Gene3D, DOMSERF, pDomTHREADER, PHYRE, SUPERFAMILY, FUGUE/TOCATTA, and VIVACE. These resources primarily use SCOP- and/or CATH-based protein domain assignments. Another objective of Genome3D is to compare structural classifications of protein domains in CATH and SCOP databases and to provide a consensus mapping of CATH and SCOP protein superfamilies. CATH/SCOP mapping analyses led to the identification of total of 1429 consensus superfamilies.Currently, Genome3D provides structural annotations for ten model organisms, including Homo sapiens, Arabidopsis thaliana, Mus musculus, Escherichia coli, Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Plasmodium falciparum, Staphylococcus aureus, and Schizosaccharomyces pombe. Thus, Genome3D serves as a common gateway to each structure prediction/annotation resource and allows users to perform comparative assessment of the predictions. It, thus, assists researchers to broaden their perspective on structure/function predictions of their query protein of interest in selected model organisms.

Depletion of Hemoglobin Transcripts and Long-Read Sequencing Improves the Transcriptome Annotation of the Polar Bear (Ursus maritimus).

Transcriptome studies evaluating whole blood and tissues are often confounded by overrepresentation of highly abundant transcripts. These abundant transcripts are problematic, as they compete with and prevent the detection of rare RNA transcripts, obscuring their biological importance. This issue is more pronounced when using long-read sequencing technologies for isoform-level transcriptome analysis, as they have relatively lower throughput compared to short-read sequencers. As a result, long-read based transcriptome analysis is prohibitively expensive for non-model organisms. While there are off-the-shelf kits available for select model organisms capable of depleting highly abundant transcripts for alpha (HBA) and beta (HBB) hemoglobin, they are unsuitable for non-model organisms. To address this, we have adapted the recent CRISPR/Cas9-based depletion method (depletion of abundant sequences by hybridization) for long-read full-length cDNA sequencing approaches that we call Long-DASH. Using a recombinant Cas9 protein with appropriate guide RNAs, full-length hemoglobin transcripts can be depleted in vitro prior to performing any short- and long-read sequencing library preparations. Using this method, we sequenced depleted full-length cDNA in parallel using both our Oxford Nanopore Technology (ONT) based R2C2 long-read approach, as well as the Illumina short-read based Smart-seq2 approach. To showcase this, we have applied our methods to create an isoform-level transcriptome from whole blood samples derived from three polar bears (Ursus maritimus). Using Long-DASH, we succeeded in depleting hemoglobin transcripts and generated deep Smart-seq2 Illumina datasets and 3.8 million R2C2 full-length cDNA consensus reads. Applying Long-DASH with our isoform identification pipeline, Mandalorion, we discovered ∼6,000 high-confidence isoforms and a number of novel genes. This indicates that there is a high diversity of gene isoforms within U. maritimus not yet reported. This reproducible and straightforward approach has not only improved the polar bear transcriptome annotations but will serve as the foundation for future efforts to investigate transcriptional dynamics within the 19 polar bear subpopulations around the Arctic.

Exploratory locomotion, a predictor of addiction vulnerability, is oligogenic in rats selected for this phenotype

Artificially selected model organisms can reveal hidden features of the genetic architecture of the complex disorders that they model. Addictions are disease phenotypes caused by different intermediate phenotypes and pathways and thereby are potentially highly polygenic. High responder (bHR) and low responder (bLR) rat lines have been selectively bred (b) for exploratory locomotion (EL), a behavioral phenotype correlated with novelty-seeking, impulsive response to reward, and vulnerability to addiction, and is inversely correlated with spontaneous anxiety and depression-like behaviors. The rapid response to selection indicates loci of large effect for EL. Using exome sequencing of HR and LR rats, we identified alleles in gene-coding regions that segregate between the two lines. Quantitative trait locus (QTL) analysis in F2 rats derived from a bHR × bLR intercross confirmed that these regions harbored genes affecting EL. The combined effects of the seven genome-wide significant QTLs accounted for approximately one-third of the total variance in EL, and two-thirds of the variance attributable to genetic factors, consistent with an oligogenic architecture of EL estimated both from the phenotypic distribution of F2 animals and rapid response to selection. Genetic association in humans linked APBA2, the ortholog of the gene at the center of the strongest QTL, with substance use disorders and related behavioral phenotypes. Our finding is also convergent with molecular and animal behavioral studies implicating Apba2 in locomotion. These results provide multilevel evidence for genes/loci influencing EL. They shed light on the genetic architecture of oligogenicity in animals artificially selected for a phenotype modeling a more complex disorder in humans.

Nucleomorph Small RNAs in Cryptophyte and Chlorarachniophyte Algae.

The regulation of gene expression and RNA maturation underlies fundamental processes such as cell homeostasis, development, and stress acclimation. The biogenesis and modification of RNA is tightly controlled by an array of regulatory RNAs and nucleic acid-binding proteins. While the role of small RNAs (sRNAs) in gene expression has been studied in-depth in select model organisms, little is known about sRNA biology across the eukaryotic tree of life. We used deep sequencing to explore the repertoires of sRNAs encoded by the miniaturized, endosymbiotically derived “nucleomorph” genomes of two single-celled algae, the cryptophyte Guillardia theta and the chlorarachniophyte Bigelowiella natans. A total of 32.3 and 35.3 million reads were generated from G. theta and B. natans, respectively. In G. theta, we identified nucleomorph U1, U2, and U4 spliceosomal small nuclear RNAs (snRNAs) as well as 11 C/D box small nucleolar RNAs (snoRNAs), five of which have potential plant and animal homologs. The snoRNAs are predicted to perform 2′-O methylation of rRNA (but not snRNA). In B. natans, we found the previously undetected 5S rRNA as well as six orphan sRNAs. Analysis of chlorarachniophyte snRNAs shed light on the removal of the miniature 18–21 nt introns found in B. natans nucleomorph genes. Neither of the nucleomorph genomes appears to encode RNA pseudouridylation machinery, and U5 snRNA cannot be found in the cryptophyte G. theta. Considering the central roles of U5 snRNA and RNA modifications in other organisms, cytoplasm-to-nucleomorph RNA shuttling in cryptophyte algae is a distinct possibility.